# -*- coding: utf-8 -*- require "numru/ggraph" include NumRu # added IntegPeriod = 4320000.0e0 tn = 50 dt = IntegPeriod / tn ##for t in [3600,4200] #for t in timearray ##for t in timearray[0..2] ##for t in timearray[6*24*5..timearray.size-1] #for timestep in 1..tn #timestep = 20 # p t # p timestep ## system('ruby draw_core.rb ' + file0 + ' ' + var0 + ' ' + file + ' ' + var + ' ' + t.to_s ) ## system('ruby draw_core3.rb ' + t.to_s ) ## system('ruby draw_core4.rb ' + t.to_s ) # system('ruby draw_core4_mod.rb ' + t.to_s ) # system('ruby draw_xt.rb ' + timestep.to_s + ' ' + tn.to_s ) # system('ruby draw_tz.rb ' + timestep.to_s ) ### # system('ruby draw_xz_thermalexpansionterm.rb' ) # system('ruby draw_z_density-thermalexpansiontermmeanx.rb' ) # system('ruby draw_z_exner-thermalexpansiontermmeanx.rb' ) # system('ruby draw_t_watervapormixingratio_atSfc.rb' ) # system('ruby draw_t_exner_atSfc.rb' ) # system('ruby draw_tz_velxmeanx.rb ' ) # system('ruby draw_tz_velzmeanx.rb ' ) # system('ruby draw_tz_exnermeanx.rb ' ) # system('ruby draw_tz_tempmeanx.rb ' ) system('ruby draw_tz_watervapormixingratiomeanx.rb ' ) # system('ruby draw_tz_condensationmeanx.rb ' ) # system('ruby draw_tz_condensationmeanx.rb ' + timestep.to_s ) # system('ruby draw_tz_cloudwatermixingratiomeanx.rb ' ) # system('ruby draw_tz_rainwatermixingratiomeanx.rb ' ) # system('ruby draw_tz_rerativehumiditymeanx.rb ' ) # system('ruby draw_tz_2mai.rb' ) # system('ruby draw_xt_kousui.rb ' + timestep.to_s ) # system('ruby draw_xt_jouhatsu.rb ' + timestep.to_s ) # system('ruby draw_xt_velx_at500m.rb ' + timestep.to_s ) # system('ruby draw_xt_temp_at500m.rb ' + timestep.to_s ) # system('ruby draw_xt_watervapormixingratio_at500m.rb ' + timestep.to_s ) # system('ruby draw_xt_watervapormixingratio_at1000m.rb ' + timestep.to_s ) # system('ruby draw_xt_watervapormixingratio_at2000m.rb ' + timestep.to_s ) # system('ruby draw_xt_velz_at500m.rb ' + timestep.to_s ) # system('ruby draw_xt_rerativehumidity_at500m.rb ' + timestep.to_s ) ### ##system("pstoimg -flip r270 -type gif dcl.ps") # system("convert -rotate 90 dcl.ps dcl_H2O-l-Rain_FallFluxAtLB_xt" + (timestep-1).to_s + ".gif" ) system("convert -rotate 90 dcl.ps dcl.gif" ) ### ## 熱膨張項の密度への寄与 (密度) # com = sprintf( "mv dcl.gif zzz-figs/dcl_DensThermExp_z_1-2_all_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_DensThermExp_z_1-2_all_z=0-12km_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_DensThermExp_z_1-2_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_DensThermExp_z_1-2_x256km_z24km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_DensThermExp_z_1-2_x256km_z36km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_DensThermExp_z_1-2_x256km_z48km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_DensThermExp_z_1-3_all_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_DensThermExp_z_1-3_all_z=0-12km_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_DensThermExp_z_1-3_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_DensThermExp_z_1-3_x256km_z24km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_DensThermExp_z_1-3_x256km_z36km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_DensThermExp_z_1-3_x256km_z48km.gif" ) ## 熱膨張項の密度への寄与 (エクスナー関数) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExp_z_1-2_all_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExp_z_1-2_all_z=0-12km_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExp_z_1-2_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExp_z_1-2_x256km_z24km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExp_z_1-2_x256km_z36km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExp_z_1-2_x256km_z48km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExp_z_1-3_all_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExp_z_1-3_all_z=0-12km_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExp_z_1-3_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExp_z_1-3_x256km_z24km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExp_z_1-3_x256km_z36km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExp_z_1-3_x256km_z48km.gif" ) ## 地表面での水蒸気混合比 # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_t_atSfc.gif" ) ## 地表面でのエクスナー関数 # com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner_t_atSfc.gif" ) ## 水平平均した水平速度 # 各計算設定において同じトーンを使う # com = sprintf( "mv dcl.gif zzz-figs/dcl_VelX_tz.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_VelZ_tz_dx200dz100.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_VelZ_tz_dx200dz200.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_VelZ_tz_dx400dz100.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_VelZ_tz_dx400dz200.gif" ) # 計算設定ごとにトーンを貼り替え # com = sprintf( "mv dcl.gif zzz-figs/dcl_VelX_tz_omakase.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_VelX_tz%0#{3}d.gif", timestep ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_VelX_tz-contour%0#{3}d.gif", timestep ) ## 水平平均したエクスナー関数 # com = sprintf( "mv dcl.gif zzz-figs/tz-exner_meanx/dcl_Exner_tz.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner_tz.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner_tz2.gif" ) ## 水平平均した温度 # com = sprintf( "mv dcl.gif zzz-figs/dcl_Temp_tz.gif") ## 水平平均した水蒸気混合比 # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_tz.gif" ) com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_tz_z=0-10km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_tz_t=0-3day.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_tz_t=0-3day_z=0-15000.gif" ) ## 水平平均した雲水混合比 # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-l-CloudAll_tz.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-l-CloudAll_tz_t=0-3day.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-l-CloudAll_tz_t=0-3day_z=0-15000.gif" ) ## 水平平均した雨水混合比 # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-l-RainAll_tz.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-l-RainAll_tz_t=0-3day.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-l-RainAll_tz_t=0-3day_z=0-15000.gif" ) ## 水平平均した正味の凝結量 # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-g_Cond_tz.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-g_Cond_tz%0#{3}d.gif", timestep ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-g_Cond_tz001.gif" ) ## 水平平均した相対湿度 # com = sprintf( "mv dcl.gif zzz-figs/dcl_RevHumid_tz.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_RevHumid_t=0-3day_tz.gif" ) ## 降水量 # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-l-Rain_FallFluxAtLB_xt%0#{3}d.gif", timestep ) ## 海面からの蒸発量 # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-g_SfcMassFlux_xt%0#{3}d.gif", timestep ) ## 海面からの熱フラックス # com = sprintf( "mv dcl.gif zzz-figs/dcl_SfcHeatFlux_xt%0#{3}d.gif", timestep ) ## 高度 500 m あたり (450 m) での水平風 # com = sprintf( "mv dcl.gif zzz-figs/dcl_VelX_xt-at500m%0#{3}d.gif", timestep ) ## 高度 500 m あたり (450 m) での温度 # com = sprintf( "mv dcl.gif zzz-figs/dcl_Temp_xt-at500m%0#{3}d.gif", timestep ) ## 高度 500 m あたり (450 m) での水蒸気混合比 # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_xt-at500m%0#{3}d.gif", timestep ) ## 高度 1000 m あたり (1050 m) での水蒸気混合比 # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_xt-at1000m%0#{3}d.gif", timestep ) ## 高度 2000 m あたり (1950 m) での水蒸気混合比 # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_xt-at2000m%0#{3}d.gif", timestep ) ## 高度 5000 m あたり (4950 m) での鉛直流 # com = sprintf( "mv dcl.gif zzz-figs/dcl_VelZ_xt-at5000m%0#{3}d.gif", timestep ) ## 高度 500 m あたり (450 m) での鉛直流 # com = sprintf( "mv dcl.gif zzz-figs/dcl_VelZ_xt-at500m%0#{3}d.gif", timestep ) ## 高度 500 m あたり (450 m) での相対湿度 # com = sprintf( "mv dcl.gif zzz-figs/dcl_RevHumid_xt-at500m%0#{3}d.gif", timestep ) system(com) #end