# 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



#########################################################################
## wet days frequency > 1 mm
# cdo
./cdo eca_rr1,1 tpr_daily_yes.nc tpr_daily_yes_freq_1.nc
./cdo eca_rr1,1 tpr_daily_no.nc tpr_daily_no_freq_1.nc
./cdo eca_rr1,1 tpr_daily_yes_MAM.nc tpr_daily_yes_MAM_freq_1.nc
./cdo eca_rr1,1 tpr_daily_no_MAM.nc tpr_daily_no_MAM_freq_1.nc
./cdo eca_rr1,1 tpr_daily_yes_JJA.nc tpr_daily_yes_JJA_freq_1.nc
./cdo eca_rr1,1 tpr_daily_no_JJA.nc tpr_daily_no_JJA_freq_1.nc
./cdo eca_rr1,1 tpr_daily_yes_SON.nc tpr_daily_yes_SON_freq_1.nc
./cdo eca_rr1,1 tpr_daily_no_SON.nc tpr_daily_no_SON_freq_1.nc
./cdo eca_rr1,1 tpr_daily_yes_DJF.nc tpr_daily_yes_DJF_freq_1.nc
./cdo eca_rr1,1 tpr_daily_no_DJF.nc tpr_daily_no_DJF_freq_1.nc

#GrADS  
clear
reinit
sdfopen tpr_daily_yes_freq_1.nc
sdfopen tpr_daily_no_freq_1.nc
set gxout shaded

page q1
d  wet_days_index_.1/20
cbarn
draw title YES
page q2
d  wet_days_index_.2/20
cbarn
draw title NO
page q4
set clevs -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80
d (wet_days_index_.1- wet_days_index_.2)/20
cbarn
draw title YES-NO
print tpr_freq_1.eps


# GrADS seasonal
clear
reinit

sdfopen tpr_daily_yes_MAM_freq_1.nc
sdfopen tpr_daily_no_MAM_freq_1.nc
set gxout shaded
set grid off
page q1
set clevs -30 -25 -20 -15 -10 -5 5 10 15 20 25 30
d (wet_days_index_.1-wet_days_index_.2)/20
cbarn
draw title MAM
close 2
close 1

sdfopen tpr_daily_yes_JJA_freq_1.nc
sdfopen tpr_daily_no_JJA_freq_1.nc
set gxout shaded
set grid off
page q2
set clevs -30 -25 -20 -15 -10 -5 5 10 15 20 25 30
d (wet_days_index_.1-wet_days_index_.2)/20
cbarn
draw title JJA
close 2
close 1

sdfopen tpr_daily_yes_SON_freq_1.nc
sdfopen tpr_daily_no_SON_freq_1.nc
set gxout shaded
set grid off
page q3
set clevs -30 -25 -20 -15 -10 -5 5 10 15 20 25 30
d (wet_days_index_.1-wet_days_index_.2)/20
cbarn
draw title SON
close 2
close 1

sdfopen tpr_daily_yes_DJF_freq_1.nc
sdfopen tpr_daily_no_DJF_freq_1.nc
set gxout shaded
set grid off
page q4
set clevs -30 -25 -20 -15 -10 -5 5 10 15 20 25 30
d (wet_days_index_.1-wet_days_index_.2)/20
cbarn
draw title DJF

print tpr_freq_1_seasonal.eps


#########################################################################
## wet days frequency > 5 mm
# cdo
./cdo eca_rr1,5 tpr_daily_yes.nc tpr_daily_yes_freq_5.nc
./cdo eca_rr1,5 tpr_daily_no.nc tpr_daily_no_freq_5.nc
./cdo eca_rr1,5 tpr_daily_yes_MAM.nc tpr_daily_yes_MAM_freq_5.nc
./cdo eca_rr1,5 tpr_daily_no_MAM.nc tpr_daily_no_MAM_freq_5.nc
./cdo eca_rr1,5 tpr_daily_yes_JJA.nc tpr_daily_yes_JJA_freq_5.nc
./cdo eca_rr1,5 tpr_daily_no_JJA.nc tpr_daily_no_JJA_freq_5.nc
./cdo eca_rr1,5 tpr_daily_yes_SON.nc tpr_daily_yes_SON_freq_5.nc
./cdo eca_rr1,5 tpr_daily_no_SON.nc tpr_daily_no_SON_freq_5.nc
./cdo eca_rr1,5 tpr_daily_yes_DJF.nc tpr_daily_yes_DJF_freq_5.nc
./cdo eca_rr1,5 tpr_daily_no_DJF.nc tpr_daily_no_DJF_freq_5.nc

#GrADS  
clear
reinit
sdfopen tpr_daily_yes_freq_5.nc
sdfopen tpr_daily_no_freq_5.nc
set gxout shaded

page q1
set clevs 0 20 40 60 80 100 120 140 160 180
d  wet_days_index_.1/20
cbarn
draw title YES
page q2
d  wet_days_index_.2/20
cbarn
draw title NO
page q4
set clevs -50 -40 -30 -20 -10 10 20 30 40 50
d (wet_days_index_.1- wet_days_index_.2)/20
cbarn
draw title YES-NO
print tpr_freq_5.eps

# GrADS seasonal
clear
reinit

sdfopen tpr_daily_yes_MAM_freq_5.nc
sdfopen tpr_daily_no_MAM_freq_5.nc
set gxout shaded
set grid off
page q1
set clevs -30 -25 -20 -15 -10 -5 5 10 15 20 25 30
d (wet_days_index_.1-wet_days_index_.2)/20
cbarn
draw title MAM
close 2
close 1

sdfopen tpr_daily_yes_JJA_freq_5.nc
sdfopen tpr_daily_no_JJA_freq_5.nc
set gxout shaded
set grid off
page q2
set clevs -30 -25 -20 -15 -10 -5 5 10 15 20 25 30
d (wet_days_index_.1-wet_days_index_.2)/20
cbarn
draw title JJA
close 2
close 1

sdfopen tpr_daily_yes_SON_freq_5.nc
sdfopen tpr_daily_no_SON_freq_5.nc
set gxout shaded
set grid off
page q3
set clevs -30 -25 -20 -15 -10 -5 5 10 15 20 25 30
d (wet_days_index_.1-wet_days_index_.2)/20
cbarn
draw title SON
close 2
close 1

sdfopen tpr_daily_yes_DJF_freq_5.nc
sdfopen tpr_daily_no_DJF_freq_5.nc
set gxout shaded
set grid off
page q4
set clevs -30 -25 -20 -15 -10 -5 5 10 15 20 25 30
d (wet_days_index_.1-wet_days_index_.2)/20
cbarn
draw title DJF

print tpr_freq_5_seasonal.eps

#########################################################################
## wet days frequency > 10 mm
# cdo
./cdo eca_rr1,10 tpr_daily_yes.nc tpr_daily_yes_freq_10.nc
./cdo eca_rr1,10 tpr_daily_no.nc tpr_daily_no_freq_10.nc
./cdo eca_rr1,10 tpr_daily_yes_MAM.nc tpr_daily_yes_MAM_freq_10.nc
./cdo eca_rr1,10 tpr_daily_no_MAM.nc tpr_daily_no_MAM_freq_10.nc
./cdo eca_rr1,10 tpr_daily_yes_JJA.nc tpr_daily_yes_JJA_freq_10.nc
./cdo eca_rr1,10 tpr_daily_no_JJA.nc tpr_daily_no_JJA_freq_10.nc
./cdo eca_rr1,10 tpr_daily_yes_SON.nc tpr_daily_yes_SON_freq_10.nc
./cdo eca_rr1,10 tpr_daily_no_SON.nc tpr_daily_no_SON_freq_10.nc
./cdo eca_rr1,10 tpr_daily_yes_DJF.nc tpr_daily_yes_DJF_freq_10.nc
./cdo eca_rr1,10 tpr_daily_no_DJF.nc tpr_daily_no_DJF_freq_10.nc

#GrADS  
clear
reinit
sdfopen tpr_daily_yes_freq_10.nc
sdfopen tpr_daily_no_freq_10.nc
set gxout shaded

page q1
set clevs 0 20 40 60 80 100 120
d  wet_days_index_.1/20
cbarn
draw title YES
page q2
set clevs 0 20 40 60 80 100 120
d  wet_days_index_.2/20
cbarn
draw title NO
page q4
set clevs -30 -25 -20 -15 -10 -5 5 10 15 20 25 30
d (wet_days_index_.1- wet_days_index_.2)/20
cbarn
draw title YES-NO
print tpr_freq_10.eps

# GrADS seasonal
clear
reinit

sdfopen tpr_daily_yes_MAM_freq_10.nc
sdfopen tpr_daily_no_MAM_freq_10.nc
set gxout shaded
set grid off
page q1
set clevs -12 -10 -8 -6 -4 -2 2 4 6 8 10 12
d (wet_days_index_.1-wet_days_index_.2)/20
cbarn
draw title MAM
close 2
close 1

sdfopen tpr_daily_yes_JJA_freq_10.nc
sdfopen tpr_daily_no_JJA_freq_10.nc
set gxout shaded
set grid off
page q2
set clevs -12 -10 -8 -6 -4 -2 2 4 6 8 10 12
d (wet_days_index_.1-wet_days_index_.2)/20
cbarn
draw title JJA
close 2
close 1

sdfopen tpr_daily_yes_SON_freq_10.nc
sdfopen tpr_daily_no_SON_freq_10.nc
set gxout shaded
set grid off
page q3
set clevs -12 -10 -8 -6 -4 -2 2 4 6 8 10 12
d (wet_days_index_.1-wet_days_index_.2)/20
cbarn
draw title SON
close 2
close 1

sdfopen tpr_daily_yes_DJF_freq_10.nc
sdfopen tpr_daily_no_DJF_freq_10.nc
set gxout shaded
set grid off
page q4
set clevs -12 -10 -8 -6 -4 -2 2 4 6 8 10 12
d (wet_days_index_.1-wet_days_index_.2)/20
cbarn
draw title DJF

print tpr_freq_10_seasonal.eps


#########################################################################
## wet days frequency > 15 mm
# cdo
./cdo eca_rr1,15 tpr_daily_yes.nc tpr_daily_yes_freq_15.nc
./cdo eca_rr1,15 tpr_daily_no.nc tpr_daily_no_freq_15.nc
./cdo eca_rr1,15 tpr_daily_yes_MAM.nc tpr_daily_yes_MAM_freq_15.nc
./cdo eca_rr1,15 tpr_daily_no_MAM.nc tpr_daily_no_MAM_freq_15.nc
./cdo eca_rr1,15 tpr_daily_yes_JJA.nc tpr_daily_yes_JJA_freq_15.nc
./cdo eca_rr1,15 tpr_daily_no_JJA.nc tpr_daily_no_JJA_freq_15.nc
./cdo eca_rr1,15 tpr_daily_yes_SON.nc tpr_daily_yes_SON_freq_15.nc
./cdo eca_rr1,15 tpr_daily_no_SON.nc tpr_daily_no_SON_freq_15.nc
./cdo eca_rr1,15 tpr_daily_yes_DJF.nc tpr_daily_yes_DJF_freq_15.nc
./cdo eca_rr1,15 tpr_daily_no_DJF.nc tpr_daily_no_DJF_freq_15.nc

#GrADS  
clear
reinit
sdfopen tpr_daily_yes_freq_15.nc
sdfopen tpr_daily_no_freq_15.nc
set gxout shaded

page q1
set clevs 0 10 20 30 40 50 60 70 80
d  wet_days_index_.1/20
cbarn
draw title YES
page q2
set clevs 0 10 20 30 40 50 60 70 80
d  wet_days_index_.2/20
cbarn
draw title NO
page q4
set clevs -30 -25 -20 -15 -10 -5 5 10 15 20 25 30
d (wet_days_index_.1- wet_days_index_.2)/20
cbarn
draw title YES-NO
print tpr_freq_15.eps

# GrADS seasonal
clear
reinit

sdfopen tpr_daily_yes_MAM_freq_15.nc
sdfopen tpr_daily_no_MAM_freq_15.nc
set gxout shaded
set grid off
page q1
set clevs -12 -10 -8 -6 -4 -2 2 4 6 8 10 12
d (wet_days_index_.1-wet_days_index_.2)/20
cbarn
draw title MAM
close 2
close 1

sdfopen tpr_daily_yes_JJA_freq_15.nc
sdfopen tpr_daily_no_JJA_freq_15.nc
set gxout shaded
set grid off
page q2
set clevs -12 -10 -8 -6 -4 -2 2 4 6 8 10 12
d (wet_days_index_.1-wet_days_index_.2)/20
cbarn
draw title JJA
close 2
close 1

sdfopen tpr_daily_yes_SON_freq_15.nc
sdfopen tpr_daily_no_SON_freq_15.nc
set gxout shaded
set grid off
page q3
set clevs -12 -10 -8 -6 -4 -2 2 4 6 8 10 12
d (wet_days_index_.1-wet_days_index_.2)/20
cbarn
draw title SON
close 2
close 1

sdfopen tpr_daily_yes_DJF_freq_15.nc
sdfopen tpr_daily_no_DJF_freq_15.nc
set gxout shaded
set grid off
page q4
set clevs -12 -10 -8 -6 -4 -2 2 4 6 8 10 12
d (wet_days_index_.1-wet_days_index_.2)/20
cbarn
draw title DJF

print tpr_freq_15_seasonal.eps




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

上一篇：estimating rainfall PDFs and densities in Africa simulations
下一篇：SVD analysis of two climate fields