# cdo
./cdo selname,tpr daily_no_land.nc tpr_daily_no.nc
./cdo selname,tpr daily_yes_land.nc tpr_daily_yes.nc
./cdo splitseas tpr_daily_no.nc tpr_daily_no_
./cdo splitseas tpr_daily_yes.nc tpr_daily_yes_


# cdo  WET 5 days

./cdo eca_cwd,1 tpr_daily_yes.nc tpr_daily_yes_wet5d.nc
./cdo eca_cwd,1 tpr_daily_no.nc tpr_daily_no_wet5d.nc
./cdo eca_cwd,1 tpr_daily_yes_MAM.nc tpr_daily_yes_MAM_wet5d.nc
./cdo eca_cwd,1 tpr_daily_no_MAM.nc tpr_daily_no_MAM_wet5d.nc
./cdo eca_cwd,1 tpr_daily_yes_JJA.nc tpr_daily_yes_JJA_wet5d.nc
./cdo eca_cwd,1 tpr_daily_no_JJA.nc tpr_daily_no_JJA_wet5d.nc
./cdo eca_cwd,1 tpr_daily_yes_SON.nc tpr_daily_yes_SON_wet5d.nc
./cdo eca_cwd,1 tpr_daily_no_SON.nc tpr_daily_no_SON_wet5d.nc
./cdo eca_cwd,1 tpr_daily_yes_DJF.nc tpr_daily_yes_DJF_wet5d.nc
./cdo eca_cwd,1 tpr_daily_no_DJF.nc tpr_daily_no_DJF_wet5d.nc

# GrADS annual mean
clear
reinit
sdfopen tpr_daily_yes_wet5d.nc
sdfopen tpr_daily_no_wet5d.nc
set gxout shaded

page q1
d number_of_cwd_p.1/20
cbarn
draw title YES
page q2
d number_of_cwd_p.2/20
cbarn
draw title NO
page q3
set clevs -5 -4 -3 -2 -1 0 1 2 3 4 5
d (number_of_cwd_p.1-number_of_cwd_p.2)/20
cbarn
draw title YES-NO

print trp5d.eps

# GrADS WET days
clear
reinit

sdfopen tpr_daily_yes_MAM_wet5d.nc
sdfopen tpr_daily_no_MAM_wet5d.nc
set gxout shaded
set grid off
page q1
set clevs -50 -40 -30 -20 -10 0 10 20 30 40 50
d number_of_cwd_p.1-number_of_cwd_p.2
cbarn
draw title MAM
close 2
close 1

sdfopen tpr_daily_yes_JJA_wet5d.nc
sdfopen tpr_daily_no_JJA_wet5d.nc
set gxout shaded
set grid off
page q2
set clevs -50 -40 -30 -20 -10 0 10 20 30 40 50
d number_of_cwd_p.1-number_of_cwd_p.2
cbarn
draw title JJA
close 2
close 1

sdfopen tpr_daily_yes_SON_wet5d.nc
sdfopen tpr_daily_no_SON_wet5d.nc
set gxout shaded
set grid off
page q3
set clevs -50 -40 -30 -20 -10 0 10 20 30 40 50
d number_of_cwd_p.1-number_of_cwd_p.2
cbarn
draw title SON
close 2
close 1

sdfopen tpr_daily_yes_DJF_wet5d.nc
sdfopen tpr_daily_no_DJF_wet5d.nc
set gxout shaded
set grid off
page q4
set clevs -50 -40 -30 -20 -10 0 10 20 30 40 50
d number_of_cwd_p.1-number_of_cwd_p.2
cbarn
draw title DJF

print trp5d_seasonal.eps


转载本文请联系原作者获取授权，同时请注明本文来自王国杰科学网博客。
链接地址：https://wap.sciencenet.cn/blog-569118-626170.html

上一篇：estimating PDF with kernels
下一篇：estimating rainfall PDFs and densities in Africa simulations