Ignore:
Timestamp:
08/30/11 14:38:29 (13 years ago)
Author:
pinsard
Message:

suppress blank lines trailing blank

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/paper01/fig14/whots_validation_net_2004_2007_all_v10.pro

    r94 r97  
    9292key_portrait = 1 
    9393coefpalit=.9 
    94  
     94; 
    9595openps, FILENAME = project_od_env+"WHOTS_validation_net_2004_2007_all_v10.ps" 
    9696; partie a changer 
    97  
    98 ;; WHOTS in the NC pacific  2004-09 
     97; 
     98; WHOTS in the NC pacific  2004-09 
    9999;22.45N, 157.54W (202.46E) 
    100100box=[202, 203, 22., 23.] 
    101  
    102 ;; first leg 
    103  
     101; 
     102; first leg 
     103; 
    104104fi=project_id_env+"met_data_WHOTS_2004_2005.txt" 
    105105res=read_ascii(fi,data_start=1) 
    … …  
    110110lwr_04=reform(ff(8,*)) 
    111111net_04=[lhf_04+shf_04+swr_04+lwr_04] 
    112  
     112; 
    113113fi=project_id_env+"met_data_WHOTS_2005_2006.txt" 
    114114res=read_ascii(fi,data_start=1) 
    … …  
    119119lwr_05=reform(ff(8,*)) 
    120120net_05=[lhf_05+shf_05+swr_05+lwr_05] 
    121  
     121; 
    122122fi=project_id_env+"met_data_WHOTS_2006_2007.txt" 
    123123res=read_ascii(fi,data_start=1) 
    … …  
    128128lwr_06=reform(ff(8,*)) 
    129129net_06=[lhf_06+shf_06+swr_06+lwr_06] 
    130  
     130; 
    131131net=[net_04, net_05, net_06] 
    132  
     132; 
    133133help, net 
    134  
    135 ;; reading other data sets 
    136  
    137 st1=20040813  
     134; 
     135; reading other data sets 
     136; 
     137st1=20040813 
    138138en1=20050725 
    139 st2=20050728  
     139st2=20050728 
    140140en2=20060624 
    141 st3=20060627  
     141st3=20060627 
    142142en3=20070628 
    143 ;st4=20070626  
     143;st4=20070626 
    144144;en4=20080606 
    145 ;st5=20080605  
     145;st5=20080605 
    146146; en5=20090715 
    147  
     147; 
    148148file=project_id_env+'erai_lhf_19890101_20091231_oafluxgrid.nc' 
    149149initncdf, file 
    150 lhf1=-1*read_ncdf('lhf',st1,en1,file=file,/nostr, box=box)  
    151 lhf1=grossemoyenne(lhf1, "xy", box=box,/nan)  
     150lhf1=-1*read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) 
     151lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) 
    152152lhf1=reform(lhf1) 
    153 lhf2=-1*read_ncdf('lhf',st2,en2,file=file,/nostr, box=box)  
    154 lhf2=grossemoyenne(lhf2, "xy", box=box,/nan)  
     153lhf2=-1*read_ncdf('lhf',st2,en2,file=file,/nostr, box=box) 
     154lhf2=grossemoyenne(lhf2, "xy", box=box,/nan) 
    155155lhf2=reform(lhf2) 
    156 lhf3=-1*read_ncdf('lhf',st3,en3,file=file,/nostr, box=box)  
    157 lhf3=grossemoyenne(lhf3, "xy", box=box,/nan)  
     156lhf3=-1*read_ncdf('lhf',st3,en3,file=file,/nostr, box=box) 
     157lhf3=grossemoyenne(lhf3, "xy", box=box,/nan) 
    158158lhf3=reform(lhf3) 
    159  
     159; 
    160160file=project_id_env+'erai_shf_19890101_20091231_oafluxgrid.nc' 
    161161initncdf, file 
    162 shf1=-1*read_ncdf('shf',st1,en1,file=file,/nostr, box=box)  
    163 shf1=grossemoyenne(shf1, "xy", box=box,/nan)  
     162shf1=-1*read_ncdf('shf',st1,en1,file=file,/nostr, box=box) 
     163shf1=grossemoyenne(shf1, "xy", box=box,/nan) 
    164164shf1=reform(shf1) 
    165 shf2=-1*read_ncdf('shf',st2,en2,file=file,/nostr, box=box)  
    166 shf2=grossemoyenne(shf2, "xy", box=box,/nan)  
     165shf2=-1*read_ncdf('shf',st2,en2,file=file,/nostr, box=box) 
     166shf2=grossemoyenne(shf2, "xy", box=box,/nan) 
    167167shf2=reform(shf2) 
    168 shf3=-1*read_ncdf('shf',st3,en3,file=file,/nostr, box=box)  
    169 shf3=grossemoyenne(shf3, "xy", box=box,/nan)  
     168shf3=-1*read_ncdf('shf',st3,en3,file=file,/nostr, box=box) 
     169shf3=grossemoyenne(shf3, "xy", box=box,/nan) 
    170170shf3=reform(shf3) 
    171  
     171; 
    172172file=project_id_env+'erai_lwr_19890101_20091231_oafluxgrid.nc' 
    173173initncdf, file 
    174 lwr1=read_ncdf('lwr',st1-1,en1,file=file,/nostr, box=box)  
    175 lwr1=grossemoyenne(lwr1, "xy", box=box,/nan)  
     174lwr1=read_ncdf('lwr',st1-1,en1,file=file,/nostr, box=box) 
     175lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) 
    176176lwr1=reform(lwr1) 
    177 lwr2=read_ncdf('lwr',st2-1,en2,file=file,/nostr, box=box)  
    178 lwr2=grossemoyenne(lwr2, "xy", box=box,/nan)  
     177lwr2=read_ncdf('lwr',st2-1,en2,file=file,/nostr, box=box) 
     178lwr2=grossemoyenne(lwr2, "xy", box=box,/nan) 
    179179lwr2=reform(lwr2) 
    180 lwr3=read_ncdf('lwr',st3-1,en3,file=file,/nostr, box=box)  
    181 lwr3=grossemoyenne(lwr3, "xy", box=box,/nan)  
     180lwr3=read_ncdf('lwr',st3-1,en3,file=file,/nostr, box=box) 
     181lwr3=grossemoyenne(lwr3, "xy", box=box,/nan) 
    182182lwr3=reform(lwr3) 
    183  
     183; 
    184184file=project_id_env+"erai_swr_19910101_20091231_oafluxgrid.nc" 
    185185initncdf, file 
    186 swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box)  
    187 swr1=grossemoyenne(swr1, "xy", box=box,/nan)  
     186swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) 
     187swr1=grossemoyenne(swr1, "xy", box=box,/nan) 
    188188swr1=reform(swr1) 
    189 swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box)  
    190 swr2=grossemoyenne(swr2, "xy", box=box,/nan)  
     189swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box) 
     190swr2=grossemoyenne(swr2, "xy", box=box,/nan) 
    191191swr2=reform(swr2) 
    192 swr3=read_ncdf('swr',st3,en3,file=file,/nostr, box=box)  
    193 swr3=grossemoyenne(swr3, "xy", box=box,/nan)  
     192swr3=read_ncdf('swr',st3,en3,file=file,/nostr, box=box) 
     193swr3=grossemoyenne(swr3, "xy", box=box,/nan) 
    194194swr3=reform(swr3) 
    195  
    196 net1=swr1+(lwr1+shf1+lhf1)  
    197 net2=swr2+(lwr2+shf2+lhf2)  
     195; 
     196net1=swr1+(lwr1+shf1+lhf1) 
     197net2=swr2+(lwr2+shf2+lhf2) 
    198198net3=swr3+(lwr3+shf3+lhf3) 
    199199era=[net1,net2,net3] 
    200 lhf1=0.  
    201 lhf2=0.  
     200lhf1=0. 
     201lhf2=0. 
    202202lhf3=0. 
    203 shf1=0.  
    204 shf2=0.  
     203shf1=0. 
     204shf2=0. 
    205205shf3=0. 
    206 lwr1=0.  
    207 lwr2=0.  
     206lwr1=0. 
     207lwr2=0. 
    208208lwr3=0. 
    209 swr1=0.  
    210 swr2=0.  
     209swr1=0. 
     210swr2=0. 
    211211swr3=0. 
    212 net1=0.  
    213 net2=0.  
     212net1=0. 
     213net2=0. 
    214214net3=0. 
    215215help, era 
    216  
     216; 
    217217file=project_id_env+"TropFlux_19890101_20091231.nc" 
    218218initncdf, file 
    219  
    220 lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box)  
    221 lhf1=grossemoyenne(lhf1, "xy", box=box,/nan)  
     219; 
     220lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) 
     221lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) 
    222222lhf1=reform(lhf1) 
    223 lhf2=read_ncdf('lhf',st2,en2,file=file,/nostr, box=box)  
    224 lhf2=grossemoyenne(lhf2, "xy", box=box,/nan)  
     223lhf2=read_ncdf('lhf',st2,en2,file=file,/nostr, box=box) 
     224lhf2=grossemoyenne(lhf2, "xy", box=box,/nan) 
    225225lhf2=reform(lhf2) 
    226 lhf3=read_ncdf('lhf',st3,en3,file=file,/nostr, box=box)  
    227 lhf3=grossemoyenne(lhf3, "xy", box=box,/nan)  
     226lhf3=read_ncdf('lhf',st3,en3,file=file,/nostr, box=box) 
     227lhf3=grossemoyenne(lhf3, "xy", box=box,/nan) 
    228228lhf3=reform(lhf3) 
    229  
    230 shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box)  
    231 shf1=grossemoyenne(shf1, "xy", box=box,/nan)  
     229; 
     230shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box) 
     231shf1=grossemoyenne(shf1, "xy", box=box,/nan) 
    232232shf1=reform(shf1) 
    233 shf2=read_ncdf('shf',st2,en2,file=file,/nostr, box=box)  
    234 shf2=grossemoyenne(shf2, "xy", box=box,/nan)  
     233shf2=read_ncdf('shf',st2,en2,file=file,/nostr, box=box) 
     234shf2=grossemoyenne(shf2, "xy", box=box,/nan) 
    235235shf2=reform(shf2) 
    236 shf3=read_ncdf('shf',st3,en3,file=file,/nostr, box=box)  
    237 shf3=grossemoyenne(shf3, "xy", box=box,/nan)  
     236shf3=read_ncdf('shf',st3,en3,file=file,/nostr, box=box) 
     237shf3=grossemoyenne(shf3, "xy", box=box,/nan) 
    238238shf3=reform(shf3) 
    239  
    240 lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box)  
    241 lwr1=grossemoyenne(lwr1, "xy", box=box,/nan)  
     239; 
     240lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box) 
     241lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) 
    242242lwr1=reform(lwr1) 
    243 lwr2=read_ncdf('lwr',st2,en2,file=file,/nostr, box=box)  
    244 lwr2=grossemoyenne(lwr2, "xy", box=box,/nan)  
     243lwr2=read_ncdf('lwr',st2,en2,file=file,/nostr, box=box) 
     244lwr2=grossemoyenne(lwr2, "xy", box=box,/nan) 
    245245lwr2=reform(lwr2) 
    246 lwr3=read_ncdf('lwr',st3,en3,file=file,/nostr, box=box)  
    247 lwr3=grossemoyenne(lwr3, "xy", box=box,/nan)  
     246lwr3=read_ncdf('lwr',st3,en3,file=file,/nostr, box=box) 
     247lwr3=grossemoyenne(lwr3, "xy", box=box,/nan) 
    248248lwr3=reform(lwr3) 
    249  
    250 swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box)  
    251 swr1=grossemoyenne(swr1, "xy", box=box,/nan)  
     249; 
     250swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) 
     251swr1=grossemoyenne(swr1, "xy", box=box,/nan) 
    252252swr1=reform(swr1) 
    253 swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box)  
    254 swr2=grossemoyenne(swr2, "xy", box=box,/nan)  
     253swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box) 
     254swr2=grossemoyenne(swr2, "xy", box=box,/nan) 
    255255swr2=reform(swr2) 
    256 swr3=read_ncdf('swr',st3,en3,file=file,/nostr, box=box)  
    257 swr3=grossemoyenne(swr3, "xy", box=box,/nan)  
     256swr3=read_ncdf('swr',st3,en3,file=file,/nostr, box=box) 
     257swr3=grossemoyenne(swr3, "xy", box=box,/nan) 
    258258swr3=reform(swr3) 
    259  
     259; 
    260260file=project_id_env+"TropFlux_swr_19890101_20091231_NRT_v50.nc" 
    261261initncdf, file 
    262 olr1=read_ncdf("sw", st1-1, en1, file=file,/nostr, box=box)  
    263 olr1=grossemoyenne(olr1, "xy",/nan, box=box)  
     262olr1=read_ncdf("sw", st1-1, en1, file=file,/nostr, box=box) 
     263olr1=grossemoyenne(olr1, "xy",/nan, box=box) 
    264264olr1=reform(olr1) 
    265 olr2=read_ncdf("sw", st2-1, en2, file=file,/nostr, box=box)  
    266 olr2=grossemoyenne(olr2, "xy",/nan, box=box)  
     265olr2=read_ncdf("sw", st2-1, en2, file=file,/nostr, box=box) 
     266olr2=grossemoyenne(olr2, "xy",/nan, box=box) 
    267267olr2=reform(olr2) 
    268 olr3=read_ncdf("sw", st3-1, en3, file=file,/nostr, box=box)  
    269 olr3=grossemoyenne(olr3, "xy",/nan, box=box)  
     268olr3=read_ncdf("sw", st3-1, en3, file=file,/nostr, box=box) 
     269olr3=grossemoyenne(olr3, "xy",/nan, box=box) 
    270270olr3=reform(olr3) 
    271  
    272 net1=swr1+lwr1+(shf1+lhf1)  
    273 net2=swr2+lwr2+(shf2+lhf2)  
     271; 
     272net1=swr1+lwr1+(shf1+lhf1) 
     273net2=swr2+lwr2+(shf2+lhf2) 
    274274net3=swr3+lwr3+(shf3+lhf3) 
    275 neto_1=olr1+lwr1+lhf1+shf1  
    276 neto_2=olr2+lwr2+lhf2+shf2  
     275neto_1=olr1+lwr1+lhf1+shf1 
     276neto_2=olr2+lwr2+lhf2+shf2 
    277277neto_3=olr3+lwr3+lhf3+shf3 
    278  
    279 trop=[net1,net2,net3]  
     278; 
     279trop=[net1,net2,net3] 
    280280olr=[neto_1,neto_2,neto_3] 
    281  
    282 lhf1=0.  
    283 lhf2=0.  
     281; 
     282lhf1=0. 
     283lhf2=0. 
    284284lhf3=0. 
    285 shf1=0.  
    286 shf2=0.  
     285shf1=0. 
     286shf2=0. 
    287287shf3=0. 
    288 lwr1=0.  
    289 lwr2=0.  
     288lwr1=0. 
     289lwr2=0. 
    290290lwr3=0. 
    291 swr1=0.  
    292 swr2=0.  
     291swr1=0. 
     292swr2=0. 
    293293swr3=0. 
    294 net1=0.  
    295 net2=0.  
     294net1=0. 
     295net2=0. 
    296296net3=0. 
    297 neto_1=0  
    298 neto_2=0  
     297neto_1=0 
     298neto_2=0 
    299299neto_3=0 
    300  
     300; 
    301301help, trop, olr 
    302  
     302; 
    303303file=project_id_env+'swr_oafluxgrid_1985_2007.nc' 
    304304initncdf, file 
    305 swr1=read_ncdf("swr", st1, en1, file=file,/nostr, box=box)  
    306 swr1=grossemoyenne(swr1, "xy",/nan, box=box)  
     305swr1=read_ncdf("swr", st1, en1, file=file,/nostr, box=box) 
     306swr1=grossemoyenne(swr1, "xy",/nan, box=box) 
    307307swr1=reform(swr1) 
    308 swr2=read_ncdf("swr", st2, en2, file=file,/nostr, box=box)  
    309 swr2=grossemoyenne(swr2, "xy",/nan, box=box)  
     308swr2=read_ncdf("swr", st2, en2, file=file,/nostr, box=box) 
     309swr2=grossemoyenne(swr2, "xy",/nan, box=box) 
    310310swr2=reform(swr2) 
    311 swr3=read_ncdf("swr", st3, en3, file=file,/nostr, box=box)  
    312 swr3=grossemoyenne(swr3, "xy",/nan, box=box)  
     311swr3=read_ncdf("swr", st3, en3, file=file,/nostr, box=box) 
     312swr3=grossemoyenne(swr3, "xy",/nan, box=box) 
    313313swr3=reform(swr3) 
    314  
     314; 
    315315file=project_id_env+'lwr_oafluxgrid_1985_2007.nc' 
    316316initncdf, file 
    317 lwr1=-1*read_ncdf("lwr", st1, en1, file=file,/nostr, box=box)  
    318 lwr1=grossemoyenne(lwr1, "xy",/nan, box=box)  
     317lwr1=-1*read_ncdf("lwr", st1, en1, file=file,/nostr, box=box) 
     318lwr1=grossemoyenne(lwr1, "xy",/nan, box=box) 
    319319lwr1=reform(lwr1) 
    320 lwr2=-1*read_ncdf("lwr", st2, en2, file=file,/nostr, box=box)  
    321 lwr2=grossemoyenne(lwr2, "xy",/nan, box=box)  
     320lwr2=-1*read_ncdf("lwr", st2, en2, file=file,/nostr, box=box) 
     321lwr2=grossemoyenne(lwr2, "xy",/nan, box=box) 
    322322lwr2=reform(lwr2) 
    323 lwr3=-1*read_ncdf("lwr", st3, en3, file=file,/nostr, box=box)  
    324 lwr3=grossemoyenne(lwr3, "xy",/nan, box=box)  
     323lwr3=-1*read_ncdf("lwr", st3, en3, file=file,/nostr, box=box) 
     324lwr3=grossemoyenne(lwr3, "xy",/nan, box=box) 
    325325lwr3=reform(lwr3) 
    326  
     326; 
    327327file=project_id_env+'lhf_oafluxgrid_1985_2009.nc' 
    328328initncdf, file 
    329 lhf1=-1*read_ncdf("lhf", st1, en1, file=file,/nostr, box=box)  
    330 lhf1=grossemoyenne(lhf1, "xy",/nan, box=box)  
     329lhf1=-1*read_ncdf("lhf", st1, en1, file=file,/nostr, box=box) 
     330lhf1=grossemoyenne(lhf1, "xy",/nan, box=box) 
    331331lhf1=reform(lhf1) 
    332 lhf2=-1*read_ncdf("lhf", st2, en2, file=file,/nostr, box=box)  
    333 lhf2=grossemoyenne(lhf2, "xy",/nan, box=box)  
     332lhf2=-1*read_ncdf("lhf", st2, en2, file=file,/nostr, box=box) 
     333lhf2=grossemoyenne(lhf2, "xy",/nan, box=box) 
    334334lhf2=reform(lhf2) 
    335 lhf3=-1*read_ncdf("lhf", st3, en3, file=file,/nostr, box=box)  
    336 lhf3=grossemoyenne(lhf3, "xy",/nan, box=box)  
     335lhf3=-1*read_ncdf("lhf", st3, en3, file=file,/nostr, box=box) 
     336lhf3=grossemoyenne(lhf3, "xy",/nan, box=box) 
    337337lhf3=reform(lhf3) 
    338  
     338; 
    339339file=project_id_env+'shf_oafluxgrid_1985_2009.nc' 
    340340initncdf, file 
    341 shf1=-1*read_ncdf("shf", st1, en1, file=file,/nostr, box=box)  
    342 shf1=grossemoyenne(shf1, "xy",/nan, box=box)  
     341shf1=-1*read_ncdf("shf", st1, en1, file=file,/nostr, box=box) 
     342shf1=grossemoyenne(shf1, "xy",/nan, box=box) 
    343343shf1=reform(shf1) 
    344 shf2=-1*read_ncdf("shf", st2, en2, file=file,/nostr, box=box)  
    345 shf2=grossemoyenne(shf2, "xy",/nan, box=box)  
     344shf2=-1*read_ncdf("shf", st2, en2, file=file,/nostr, box=box) 
     345shf2=grossemoyenne(shf2, "xy",/nan, box=box) 
    346346shf2=reform(shf2) 
    347 shf3=-1*read_ncdf("shf", st3, en3, file=file,/nostr, box=box)  
    348 shf3=grossemoyenne(shf3, "xy",/nan, box=box)  
     347shf3=-1*read_ncdf("shf", st3, en3, file=file,/nostr, box=box) 
     348shf3=grossemoyenne(shf3, "xy",/nan, box=box) 
    349349shf3=reform(shf3) 
    350  
    351 net1=swr1+lwr1+lhf1+shf1  
    352 net2=swr2+lwr2+lhf2+shf2  
     350; 
     351net1=swr1+lwr1+lhf1+shf1 
     352net2=swr2+lwr2+lhf2+shf2 
    353353net3=swr3+lwr3+lhf3+shf3 
    354  
     354; 
    355355oaf=[net1,net2,net3] 
    356 lhf1=0.  
    357 lhf2=0.  
     356lhf1=0. 
     357lhf2=0. 
    358358lhf3=0. 
    359 shf1=0.  
    360 shf2=0.  
     359shf1=0. 
     360shf2=0. 
    361361shf3=0. 
    362 lwr1=0.  
    363 lwr2=0.  
     362lwr1=0. 
     363lwr2=0. 
    364364lwr3=0. 
    365 swr1=0.  
    366 swr2=0.  
     365swr1=0. 
     366swr2=0. 
    367367swr3=0. 
    368 net1=0.  
     368net1=0. 
    369369net2=0. 
    370370net3=0. 
    371371help, oaf 
    372  
     372; 
    373373fi=project_id_env+'swr_ncep2_oaflxgrid_19890101_20091231.nc' 
    374374initncdf, fi 
    375 swr1=read_ncdf("swr", st1, en1, file=fi,/nostr, box=box)  
    376 swr1=grossemoyenne(swr1, "xy",/nan, box=box)  
     375swr1=read_ncdf("swr", st1, en1, file=fi,/nostr, box=box) 
     376swr1=grossemoyenne(swr1, "xy",/nan, box=box) 
    377377swr1=reform(swr1) 
    378 swr2=read_ncdf("swr", st2, en2, file=fi,/nostr, box=box)  
    379 swr2=grossemoyenne(swr2, "xy",/nan, box=box)  
     378swr2=read_ncdf("swr", st2, en2, file=fi,/nostr, box=box) 
     379swr2=grossemoyenne(swr2, "xy",/nan, box=box) 
    380380swr2=reform(swr2) 
    381 swr3=read_ncdf("swr", st3, en3, file=fi,/nostr, box=box)  
    382 swr3=grossemoyenne(swr3, "xy",/nan, box=box)  
     381swr3=read_ncdf("swr", st3, en3, file=fi,/nostr, box=box) 
     382swr3=grossemoyenne(swr3, "xy",/nan, box=box) 
    383383swr3=reform(swr3) 
    384  
     384; 
    385385fi=project_id_env+'NCEP2_flux_19890101_20090729.nc" 
    386386initncdf, fi 
    387 lwr1=-1*read_ncdf("lwr", st1, en1, file=fi,/nostr, box=box)  
    388 lwr1=grossemoyenne(lwr1, "xy",/nan, box=box)  
     387lwr1=-1*read_ncdf("lwr", st1, en1, file=fi,/nostr, box=box) 
     388lwr1=grossemoyenne(lwr1, "xy",/nan, box=box) 
    389389lwr1=reform(lwr1) 
    390 lwr2=-1*read_ncdf("lwr", st2, en2, file=fi,/nostr, box=box)  
    391 lwr2=grossemoyenne(lwr2, "xy",/nan, box=box)  
     390lwr2=-1*read_ncdf("lwr", st2, en2, file=fi,/nostr, box=box) 
     391lwr2=grossemoyenne(lwr2, "xy",/nan, box=box) 
    392392lwr2=reform(lwr2) 
    393 lwr3=-1*read_ncdf("lwr", st3, en3, file=fi,/nostr, box=box)  
    394 lwr3=grossemoyenne(lwr3, "xy",/nan, box=box)  
     393lwr3=-1*read_ncdf("lwr", st3, en3, file=fi,/nostr, box=box) 
     394lwr3=grossemoyenne(lwr3, "xy",/nan, box=box) 
    395395lwr3=reform(lwr3) 
    396  
     396; 
    397397fi=project_id_env+'lhf_ncep2_oafluxgrid_19890101_20091231.nc' 
    398398initncdf, fi 
    399 lhf1=-1*read_ncdf("lhf", st1, en1, file=fi,/nostr, box=box)  
    400 lhf1=grossemoyenne(lhf1, "xy",/nan, box=box)  
     399lhf1=-1*read_ncdf("lhf", st1, en1, file=fi,/nostr, box=box) 
     400lhf1=grossemoyenne(lhf1, "xy",/nan, box=box) 
    401401lhf1=reform(lhf1) 
    402 lhf2=-1*read_ncdf("lhf", st2, en2, file=fi,/nostr, box=box)  
    403 lhf2=grossemoyenne(lhf2, "xy",/nan, box=box)  
     402lhf2=-1*read_ncdf("lhf", st2, en2, file=fi,/nostr, box=box) 
     403lhf2=grossemoyenne(lhf2, "xy",/nan, box=box) 
    404404lhf2=reform(lhf2) 
    405 lhf3=-1*read_ncdf("lhf", st3, en3, file=fi,/nostr, box=box)  
    406 lhf3=grossemoyenne(lhf3, "xy",/nan, box=box)  
     405lhf3=-1*read_ncdf("lhf", st3, en3, file=fi,/nostr, box=box) 
     406lhf3=grossemoyenne(lhf3, "xy",/nan, box=box) 
    407407lhf3=reform(lhf3) 
    408  
     408; 
    409409fi=project_id_env+'shf_ncep2_oafluxgrid_19890101_20091231.nc' 
    410410initncdf, fi 
    411 shf1=-1*read_ncdf("shf", st1, en1, file=fi,/nostr, box=box)  
    412 shf1=grossemoyenne(shf1, "xy",/nan, box=box)  
     411shf1=-1*read_ncdf("shf", st1, en1, file=fi,/nostr, box=box) 
     412shf1=grossemoyenne(shf1, "xy",/nan, box=box) 
    413413shf1=reform(shf1) 
    414 shf2=-1*read_ncdf("shf", st2, en2, file=fi,/nostr, box=box)  
    415 shf2=grossemoyenne(shf2, "xy",/nan, box=box)  
     414shf2=-1*read_ncdf("shf", st2, en2, file=fi,/nostr, box=box) 
     415shf2=grossemoyenne(shf2, "xy",/nan, box=box) 
    416416shf2=reform(shf2) 
    417 shf3=-1*read_ncdf("shf", st3, en3, file=fi,/nostr, box=box)  
    418 shf3=grossemoyenne(shf3, "xy",/nan, box=box)  
     417shf3=-1*read_ncdf("shf", st3, en3, file=fi,/nostr, box=box) 
     418shf3=grossemoyenne(shf3, "xy",/nan, box=box) 
    419419shf3=reform(shf3) 
    420  
    421 net1=swr1+lwr1+lhf1+shf1  
    422 net2=swr2+lwr2+lhf2+shf2  
     420; 
     421net1=swr1+lwr1+lhf1+shf1 
     422net2=swr2+lwr2+lhf2+shf2 
    423423net3=swr3+lwr3+lhf3+shf3 
    424  
     424; 
    425425nce2=[net1,net2,net3] 
    426426help, nce2 
    427  
    428 lhf1=0.  
    429 lhf2=0.  
     427; 
     428lhf1=0. 
     429lhf2=0. 
    430430lhf3=0. 
    431 shf1=0.  
    432 shf2=0.  
     431shf1=0. 
     432shf2=0. 
    433433shf3=0. 
    434 lwr1=0.  
    435 lwr2=0.  
     434lwr1=0. 
     435lwr2=0. 
    436436lwr3=0. 
    437 swr1=0.  
    438 swr2=0.  
     437swr1=0. 
     438swr2=0. 
    439439swr3=0. 
    440 net1=0.  
    441 net2=0.  
     440net1=0. 
     441net2=0. 
    442442net3=0. 
    443  
     443; 
    444444file=project_id_env+"ncep1_flux_19890101_20091231.nc" 
    445445initncdf, file 
    446  
    447 lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box)  
    448 lhf1=grossemoyenne(lhf1, "xy", box=box,/nan)  
     446; 
     447lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) 
     448lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) 
    449449lhf1=reform(lhf1) 
    450 lhf2=read_ncdf('lhf',st2,en2,file=file,/nostr, box=box)  
    451 lhf2=grossemoyenne(lhf2, "xy", box=box,/nan)  
     450lhf2=read_ncdf('lhf',st2,en2,file=file,/nostr, box=box) 
     451lhf2=grossemoyenne(lhf2, "xy", box=box,/nan) 
    452452lhf2=reform(lhf2) 
    453 lhf3=read_ncdf('lhf',st3,en3,file=file,/nostr, box=box)  
    454 lhf3=grossemoyenne(lhf3, "xy", box=box,/nan)  
     453lhf3=read_ncdf('lhf',st3,en3,file=file,/nostr, box=box) 
     454lhf3=grossemoyenne(lhf3, "xy", box=box,/nan) 
    455455lhf3=reform(lhf3) 
    456  
    457 shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box)  
    458 shf1=grossemoyenne(shf1, "xy", box=box,/nan)  
     456; 
     457shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box) 
     458shf1=grossemoyenne(shf1, "xy", box=box,/nan) 
    459459shf1=reform(shf1) 
    460 shf2=read_ncdf('shf',st2,en2,file=file,/nostr, box=box)  
    461 shf2=grossemoyenne(shf2, "xy", box=box,/nan)  
     460shf2=read_ncdf('shf',st2,en2,file=file,/nostr, box=box) 
     461shf2=grossemoyenne(shf2, "xy", box=box,/nan) 
    462462shf2=reform(shf2) 
    463 shf3=read_ncdf('shf',st3,en3,file=file,/nostr, box=box)  
    464 shf3=grossemoyenne(shf3, "xy", box=box,/nan)  
     463shf3=read_ncdf('shf',st3,en3,file=file,/nostr, box=box) 
     464shf3=grossemoyenne(shf3, "xy", box=box,/nan) 
    465465shf3=reform(shf3) 
    466  
    467 lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box)  
    468 lwr1=grossemoyenne(lwr1, "xy", box=box,/nan)  
     466; 
     467lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box) 
     468lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) 
    469469lwr1=reform(lwr1) 
    470 lwr2=read_ncdf('lwr',st2,en2,file=file,/nostr, box=box)  
    471 lwr2=grossemoyenne(lwr2, "xy", box=box,/nan)  
     470lwr2=read_ncdf('lwr',st2,en2,file=file,/nostr, box=box) 
     471lwr2=grossemoyenne(lwr2, "xy", box=box,/nan) 
    472472lwr2=reform(lwr2) 
    473 lwr3=read_ncdf('lwr',st3,en3,file=file,/nostr, box=box)  
    474 lwr3=grossemoyenne(lwr3, "xy", box=box,/nan)  
     473lwr3=read_ncdf('lwr',st3,en3,file=file,/nostr, box=box) 
     474lwr3=grossemoyenne(lwr3, "xy", box=box,/nan) 
    475475lwr3=reform(lwr3) 
    476  
    477 swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box)  
    478 swr1=grossemoyenne(swr1, "xy", box=box,/nan)  
     476; 
     477swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) 
     478swr1=grossemoyenne(swr1, "xy", box=box,/nan) 
    479479swr1=reform(swr1) 
    480 swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box)  
    481 swr2=grossemoyenne(swr2, "xy", box=box,/nan)  
     480swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box) 
     481swr2=grossemoyenne(swr2, "xy", box=box,/nan) 
    482482swr2=reform(swr2) 
    483 swr3=read_ncdf('swr',st3,en3,file=file,/nostr, box=box)  
    484 swr3=grossemoyenne(swr3, "xy", box=box,/nan)  
     483swr3=read_ncdf('swr',st3,en3,file=file,/nostr, box=box) 
     484swr3=grossemoyenne(swr3, "xy", box=box,/nan) 
    485485swr3=reform(swr3) 
    486  
    487 net1=swr1+(lwr1+shf1+lhf1)  
    488 net2=swr2+(lwr2+shf2+lhf2)  
     486; 
     487net1=swr1+(lwr1+shf1+lhf1) 
     488net2=swr2+(lwr2+shf2+lhf2) 
    489489net3=swr3+(lwr3+shf3+lhf3) 
    490490nce1=[net1,net2,net3] 
    491491help, nce1 
    492  
    493 lhf1=0.  
    494 lhf2=0.  
     492; 
     493lhf1=0. 
     494lhf2=0. 
    495495lhf3=0. 
    496 shf1=0.  
    497 shf2=0.  
     496shf1=0. 
     497shf2=0. 
    498498shf3=0. 
    499 lwr1=0.  
    500 lwr2=0.  
     499lwr1=0. 
     500lwr2=0. 
    501501lwr3=0. 
    502 swr1=0.  
    503 swr2=0.  
     502swr1=0. 
     503swr2=0. 
    504504swr3=0. 
    505 net1=0.  
    506 net2=0.  
     505net1=0. 
     506net2=0. 
    507507net3=0. 
    508  
     508; 
    509509ind=where(finite(net)) 
    510 net=net(ind)  
    511 era=era(ind)  
    512 trop=trop(ind)  
     510net=net(ind) 
     511era=era(ind) 
     512trop=trop(ind) 
    513513oaf=oaf(ind) 
    514 nce2=nce2(ind)  
    515 nce1=nce1(ind)  
     514nce2=nce2(ind) 
     515nce1=nce1(ind) 
    516516olr=olr(ind) 
    517  
    518  
     517; 
     518; 
    519519statistics, net, era, $ 
    520520                 cor, bias, std, rmsd 
    521521print, cor, bias, std,  rmsd 
    522522cstat_era=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    523  
     523; 
    524524statistics, net, trop, $ 
    525525                 cor, bias, std, rmsd 
    526526print, cor, bias, std,  rmsd 
    527527cstat_trop=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    528  
     528; 
    529529statistics, net, oaf, $ 
    530530                 cor, bias, std, rmsd 
    531531print, cor, bias, std,  rmsd 
    532532cstat_oaf=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    533  
     533; 
    534534statistics, net, nce1, $ 
    535535                 cor, bias, std, rmsd 
    536536print, cor, bias, std,  rmsd 
    537537cstat_nce1=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    538  
     538; 
    539539statistics, net, nce2, $ 
    540540                 cor, bias, std, rmsd 
    541541print, cor, bias, std,  rmsd 
    542542cstat_nce2=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    543  
     543; 
    544544statistics, net, olr, $ 
    545545                 cor, bias, std, rmsd 
    546546print, cor, bias, std,  rmsd 
    547547cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    548  
    549  
     548; 
     549; 
    550550splot, net, era, charsize=1.1, title='Net Flux - WHOTS Vs ERAI', $ 
    551551     xrange=[-250,250], yrange=[-250,250], small=[2,3,1], psym=2, xmin=1, ymin=1 
    552552oplot, [-250,250], [-250,250] 
    553553xyouts, 0,-200, cstat_era, charsize=1. 
    554  
     554; 
    555555splot, net, trop, charsize=1.1, title='Net Flux - WHOTS Vs TropFlux',/noer, $ 
    556556     xrange=[-250,250], yrange=[-250,250], small=[2,3,2], psym=2, xmin=1, ymin=1 
    557557oplot, [-250,250], [-250,250] 
    558558xyouts, 0,-200, cstat_trop, charsize=1. 
    559  
     559; 
    560560splot, net, oaf, charsize=1.1, title='Net Flux - WHOTS Vs OAFlux',/noer, $ 
    561561     xrange=[-250,250], yrange=[-250,250], small=[2,3,3], psym=2, xmin=1, ymin=1 
    562562oplot, [-250,250], [-250,250] 
    563563xyouts, 0,-200, cstat_oaf, charsize=1. 
    564  
     564; 
    565565splot, net, nce1, charsize=1.1, title='Net Flux - WHOTS Vs NCEP1',/noer, $ 
    566566     xrange=[-250,250], yrange=[-250,250], small=[2,3,4], psym=2, xmin=1, ymin=1 
    567567oplot, [-250,250], [-250,250] 
    568568xyouts, 0,-200, cstat_nce1, charsize=1. 
    569  
     569; 
    570570splot, net, nce2, charsize=1.1, title='Net Flux - WHOTS Vs NCEP2',/noer, $ 
    571571     xrange=[-250,250], yrange=[-250,250], small=[2,3,5], psym=2, xmin=1, ymin=1 
    572572oplot, [-250,250], [-250,250] 
    573573xyouts, 0,-200, cstat_nce2, charsize=1. 
    574  
     574; 
    575575splot, net, olr, charsize=1.1, title='Net Flux - WHOTS Vs TropFlux_NRT',/noer, $ 
    576576     xrange=[-250,250], yrange=[-250,250], small=[2,3,6], psym=2, xmin=1, ymin=1 
    577577oplot, [-250,250], [-250,250] 
    578578xyouts, 0,-200, cstat_olr, charsize=1. 
    579  
     579; 
    580580closeps 
     581; 
    581582end 
Note: See TracChangeset for help on using the changeset viewer.