# -*- coding: utf-8 -*- require "numru/ggraph" include NumRu # added IntegPeriod = 4320000.0e0 dt = 432000 tn = IntegPeriod / dt #for t in timearray #for timestep in 1..tn # p timestep ##### 基本場 (z 図) ##### ## 温度場 # system('ruby draw_z_basicstate.rb ' + "Temp" ) # system("convert -rotate 90 dcl.ps dcl.gif" ) ## com = sprintf( "mv dcl.gif zzz-figs/dcl_TempBZ_z.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_TempBZ_z_z=0-10km.gif" ) # system(com) ## 圧力場 # system('ruby draw_z_basicstate.rb ' + "Press" ) # system("convert -rotate 90 dcl.ps dcl.gif" ) ## com = sprintf( "mv dcl.gif zzz-figs/dcl_PressBZ_z.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_PressBZ_z_z=0-10km.gif" ) # system(com) ## 湿度場 # system('ruby draw_z_basicstate.rb ' + "Hum" ) # system("convert -rotate 90 dcl.ps dcl.gif" ) ## com = sprintf( "mv dcl.gif zzz-figs/dcl_HumBZ_z.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_HumBZ_z_z=0-10km.gif" ) # system(com) ## 水蒸気混合比場 # system('ruby draw_z_basicstate.rb ' + "H2O-g" ) # system("convert -rotate 90 dcl.ps dcl.gif" ) ## com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gBZ_z.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gBZ_z_z=0-10km.gif" ) # system(com) #exit ##### 擾乱場 (xz 図) ##### ## 温度場 # system('ruby draw_xz_disturb.rb ' + "Temp" ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_TempDist_xz.gif" ) # system(com) ## 温位場 # system('ruby draw_xz_disturb.rb ' + "PTemp" ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_PTempDist_xz.gif" ) # system(com) ## 水蒸気混合比場 # system('ruby draw_xz_disturb.rb ' + "H2O-g" ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gDist_xz.gif" ) # system(com) #exit ##### t 図 ##### ## 地表面でのエクスナー関数 # system('ruby draw_t_exner_atSfc.rb' ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner_t_atSfc_hikaku.gif" ) ## com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner_t_atSfc_1-2_all.gif" ) ## com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner_t_atSfc_1-3_all.gif" ) # system(com) ## 地表面での水蒸気混合比 system('ruby draw_t_watervapormixingratio_atSfc.rb' ) system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_t_atSfc.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_t_atSfc_1-2_all.gif" ) com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_t_atSfc_1-3_all.gif" ) system(com) ## 水平平均水蒸気混合比の総量の時間変化 # system('ruby draw_t_budgetofwatervapormeanx.rb' ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Budget-H2O-gAll_t.gif" ) # system(com) ## 6 時間ごとに積算した水平平均降水量 # system('ruby draw_t_int-rainfallfluxatLB.rb' ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Int-H2O-l-Rain_FallFluxAtLB_t.gif" ) # system(com) ## 水平平均降水量の積算値の時間変化 # system('ruby draw_t_sekisan-rainfallfluxatLB.rb' ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Sekisan-H2O-l-Rain_FallFluxAtLB_t.gif" ) # system(com) ## 水平平均した降水量 # system('ruby draw_t_rainfallfluxatLB.rb' ) # system('ruby draw_t_rainfallfluxatLB.rb ' + timestep.to_s ) # system("convert -rotate 90 dcl.ps dcl.gif" ) ## com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-l-Rain_FallFluxAtLB_t.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-l-Rain_FallFluxAtLB_t_1-3_x256km_z48km_%0#{8}d.gif", timestep*dt ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-l-Rain_FallFluxAtLB_t_1-5_x256km_z48km_%0#{8}d.gif", timestep*dt ) # system(com) ## 水平平均蒸発量の積算値の時間変化 # system('ruby draw_t_sekisan-watervaporsfcmassflux.rb' ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Sekisan-H2O-g_SfcMassFlux_t.gif" ) # system(com) ## 水平平均した蒸発量 # system('ruby draw_t_watervapormixingratio_sfcmassflux.rb' ) # system('ruby draw_t_watervapormixingratio_sfcmassflux.rb ' + timestep.to_s ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-g_SfcMassFlux_t.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-g_SfcMassFlux_t_1-3_x256km_z48km_%0#{8}d.gif", timestep*dt ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-g_SfcMassFlux_t_1-5_x256km_z48km_%0#{8}d.gif", timestep*dt ) # system(com) ##### z 図 ##### ## 時間平均した水蒸気混合比の鉛直分布 # system('ruby draw_z_watervapormixingratiomeanx.rb' ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_z.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_z_t=29-30day.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_z_z=0-10km_t=49-50day.gif" ) # system(com) ## 時間平均した温度の鉛直分布 # system('ruby draw_z_tempmeanx.rb' ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Temp_z.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Temp_z_t=49-50day.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Temp_z_z=0-10km_t=49-50day.gif" ) # system(com) ## 時間平均したエクスナー関数の鉛直分布 # system('ruby draw_z_exnermeanx.rb' ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner_z.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner_z_z=0-10km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner_z_t=49-50day.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner_z_z=0-10km_t=49-50day.gif" ) # system(com) ## 時間積分した熱膨張項への寄与 (エクスナー関数) # system('ruby draw_z_exner-thermalexpansiontermmeanx.rb' ) # system("convert -rotate 90 dcl.ps dcl.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.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" ) # system(com) ##### tz 図 ##### ## 水平平均したエクスナー関数 # system('ruby draw_tz_exnermeanx.rb ' ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner_tz.gif" ) # system(com) ## 水平平均した水平速度 # system('ruby draw_tz_velxmeanx.rb ' ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_VelX_tz.gif" ) # system(com) ## 水平平均した鉛直速度 # system('ruby draw_tz_velzmeanx.rb ' ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_VelZ_tz.gif" ) # system(com) ## 水平平均した温度 # system('ruby draw_tz_tempmeanx.rb ' ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Temp_tz.gif") # system(com) ## 水平平均した水蒸気混合比 # system('ruby draw_tz_watervapormixingratiomeanx.rb ' ) # system("convert -rotate 90 dcl.ps dcl.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" ) # system(com) ## 水平平均した正味の凝結量 # system('ruby draw_tz_condensationmeanx.rb ' ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-g_Cond_tz.gif" ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # system(com) ## 水平平均した雲水混合比 # system('ruby draw_tz_cloudwatermixingratiomeanx.rb ' ) # system("convert -rotate 90 dcl.ps dcl.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_z=0-20km.gif" ) # system(com) ## 水平平均した雨水混合比 # system('ruby draw_tz_rainwatermixingratiomeanx.rb ' ) # system("convert -rotate 90 dcl.ps dcl.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_z=0-20km.gif" ) # system(com) ## 水平平均した相対湿度 # system('ruby draw_tz_rerativehumiditymeanx.rb ' ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_RevHumid_tz.gif" ) # system(com) # system('ruby draw_tz_2mai.rb' ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # system(com) ## 熱膨張項の密度への寄与 (エクスナー関数) # system('ruby draw_tz_thermalexpansionterm-dthetadt.rb' ) # system('ruby draw_tz_thermalexpansionterm-dthetadt.rb ' + timestep.to_s ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner-ThermExpTerm_tz_1-3_x256km_z48km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner-ThermExpTerm_tz_1-5_x256km_z48km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner-ThermExpTerm_tz_1-3_x256km_z48km_%0#{8}d.gif", timestep*dt ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner-ThermExpTerm_tz_1-5_x256km_z48km_%0#{8}d.gif", timestep*dt ) # system(com) ##### t 図 ##### ## 熱膨張項の密度への寄与 (密度) # system('ruby draw_z_density-thermalexpansiontermmeanx.rb' ) # 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_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_z48km.gif" ) # system(com) ## 熱膨張項の密度への寄与 (エクスナー関数) system('ruby draw_z_exner-thermalexpansiontermdthetadtmeanx.rb' ) system("convert -rotate 90 dcl.ps dcl.gif" ) com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDthetaDt_z_1-2_all_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDthetaDt_z_1-2_all_z=0-12km_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDthetaDt_z_1-2_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDthetaDt_z_1-2_x256km_z24km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDthetaDt_z_1-2_x256km_z36km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDthetaDt_z_1-2_x256km_z48km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDthetaDt_z_1-2_z=0-12km_x256km_z48km.gif" ) com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDthetaDt_z_1-3_all_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDthetaDt_z_1-3_all_z=0-12km_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDthetaDt_z_1-3_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDthetaDt_z_1-3_x256km_z24km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDthetaDt_z_1-3_x256km_z36km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDthetaDt_z_1-3_x256km_z48km.gif" ) system(com) system('ruby draw_z_exner-thermalexpansiontermdqdtmeanx.rb' ) system("convert -rotate 90 dcl.ps dcl.gif" ) com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDqDt_z_1-2_all_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDqDt_z_1-2_all_z=0-12km_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDqDt_z_1-2_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDqDt_z_1-2_x256km_z24km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDqDt_z_1-2_x256km_z36km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDqDt_z_1-2_x256km_z48km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDqDt_z_1-2_z=0-12km_x256km_z48km.gif" ) com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDqDt_z_1-3_all_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDqDt_z_1-3_all_z=0-12km_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDqDt_z_1-3_x256km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDqDt_z_1-3_x256km_z24km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDqDt_z_1-3_x256km_z36km.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_ExnerThermExpDqDt_z_1-3_x256km_z48km.gif" ) system(com) system('ruby draw_z_exner-thermalexpansiontermmeanx.rb' ) system("convert -rotate 90 dcl.ps dcl.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-2_z=0-12km_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" ) system(com) #end ##### xt 図 ##### #file = 'thermal-moist_PTempAll.nc' #var = 'PTempAll' # #gphys = GPhys::IO.open( file , var ) #timearray = gphys.coord(3).val # #for t in timearray # p t ## 降水量 # system('ruby draw_xt_kousui.rb ' + timestep.to_s ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-l-Rain_FallFluxAtLB_xt%0#{3}d.gif", timestep ) # system(com) ## 海面からの蒸発量 # system('ruby draw_xt_jouhatsu.rb ' + timestep.to_s ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-g_SfcMassFlux_xt%0#{3}d.gif", timestep ) # system(com) ## 高度 500 m あたり (450 m) での水平風 # system('ruby draw_xt_velx_at500m.rb ' + timestep.to_s ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_VelX_xt-at500m%0#{3}d.gif", timestep ) # system(com) ## 高度 500 m あたり (450 m) での温度 # system('ruby draw_xt_temp_at500m.rb ' + timestep.to_s ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_Temp_xt-at500m%0#{3}d.gif", timestep ) # system(com) ## 高度 500 m あたり (450 m) での水蒸気混合比 # system('ruby draw_xt_watervapormixingratio_at500m.rb ' + timestep.to_s ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_xt-at500m%0#{3}d.gif", timestep ) # system(com) ## 高度 1000 m あたり (1050 m) での水蒸気混合比 # system('ruby draw_xt_watervapormixingratio_at1000m.rb ' + timestep.to_s ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_xt-at1000m%0#{3}d.gif", timestep ) # system(com) ## 高度 2000 m あたり (1950 m) での水蒸気混合比 # system('ruby draw_xt_watervapormixingratio_at2000m.rb ' + timestep.to_s ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_xt-at2000m%0#{3}d.gif", timestep ) # system(com) ## 高度 500 m あたり (450 m) での鉛直流 # system('ruby draw_xt_velz_at500m.rb ' + timestep.to_s ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_VelZ_xt-at500m%0#{3}d.gif", timestep ) # system(com) ## 高度 500 m あたり (450 m) での相対湿度 # system('ruby draw_xt_rerativehumidity_at500m.rb ' + timestep.to_s ) # system("convert -rotate 90 dcl.ps dcl.gif" ) # com = sprintf( "mv dcl.gif zzz-figs/dcl_RevHumid_xt-at500m%0#{3}d.gif", timestep ) # system(com) #end