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/fig7/air_validation_scatter_2000_2009_basin.pro

    r94 r97  
    9090reinitplt, /z,/invert 
    9191key_portrait = 1 
    92  
     92; 
    9393openps, FILENAME = project_od_env+'air_validation_scatter_2000_2009_'+ocean+'.ps' 
    94  
    95 ;; Give the location of mooring for validation of basic meteorological variables 
    96  
    97  
     94; 
     95; Give the location of mooring for validation of basic meteorological variables 
     96; 
     97; 
    9898sitelist=['8s67e','12s55e', '8s55e', '8s80.5e', '1.5s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $ 
    9999           '0n90e', '1.5n90e', '4n90e','8n90e','12n90e', '15n90e', '5s95e', $ 
    … …  
    106106          '0n0e', '0n10w', '0n23w', '0n35w', '10s10w', '12n23w', '12n38w', '14s32w', '15n38w', '19s34w', '20n38w', $ 
    107107          '21n23w', '4n23w', '4n38w', '6s10w', '8n38w', '8s30w'] 
    108  
     108; 
    109109ocean='global' 
    110  
     110; 
    111111sitelist=['8s67e','12s55e', '8s55e', '8s80.5e', '1.5s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $ 
    112112           '0n90e', '1.5n90e', '4n90e','8n90e','12n90e', '15n90e', '5s95e'] 
    113113ocean='IO' 
    114  
     114; 
    115115;sitelist=['8s165e', '8s180w',  '8s155w', '8s125w', '8s110w', '8s95w',  '5s156e', '5s165e', '5s180w', '5s170w', $ 
    116116;          '5s155w', '5s140w', '5s125w', '5s110w', '5s95w', '2s156e', '2s165e', '2s180w', '2s170w', '2s155w', '2s140w', $ 
    … …  
    120120;          '8n156e', '8n165e', '8n180w', '8n170w', '9n140w', '8n125w', '8n110w', '8n95w'] 
    121121;ocean='Pac' 
    122  
     122; 
    123123;sitelist=['0n0e', '0n10w', '0n23w', '0n35w', '10s10w', '12n23w', '12n38w', '14s32w', '15n38w', '19s34w', '20n38w', $ 
    124124;          '21n23w', '4n23w', '4n38w', '6s10w', '8n38w', '8s30w'] 
    125125;ocean='AO' 
    126  
     126; 
    127127;sitelist=['8s165e', '8s180w',  '8s155w','5s156e', '5s165e', '5s180w', '5s170w','2s156e', '2s165e', '2s180w', $ 
    128128;          '2s170w','0n147e', '0n156e', '0n165e', '0n180w', '0n170w', '2n147e', '2n156e', '2n165e', '2n180w', $ 
    … …  
    130130;          '8n125w', '8n110w', '8n95w'] 
    131131;ocean='Pac_warmpool' 
    132  
     132; 
    133133;sitelist=['8s125w', '8s110w', '8s95w','5s155w', '5s140w', '5s125w', '5s110w', '5s95w','2s155w', '2s140w', $ 
    134134;          '2s125w', '2s110w', '2s95w', '0n155w', '0n140w', '0n125w', '0n110w', '0n95w', '2n155w', '2n140w', $ 
    135135;         '2n125w', '2n110w', '2n95w',  '5n155w', '5n140w', '5n125w', '5n110w', '5n95w'] 
    136136;ocean='Pac_coldtongue' 
    137  
     137; 
    138138;sitelist=['0n0e','0n35w', '4n23w', '4n38w', '8n38w'] 
    139139;ocean='AO_warmpool' 
    140  
     140; 
    141141;sitelist=['0n10w','0n23w','10s10w','12n23w','12n38w','14s32w','15n38w', '19s34w','20n38w', $ 
    142142;          '21n23w', '6s10w', '8s30w'] 
    143143;ocean='AO_coldtongue' 
    144  
    145 nsmooth=1.    ;; statistics are with 7 day smoothed 
    146 ;;   This program will create the following text files with statistics of respective variables 
     144; 
     145nsmooth=1.    ; statistics are with 7 day smoothed 
     146;   This program will create the following text files with statistics of respective variables 
    147147close,/all 
    148  
     148; 
    149149fi_air_erai=project_id_env+'air_2000_2009_erai_'+ocean+'.txt' 
    150150openw,1,fi_air_erai 
    … …  
    157157fi_air_ncep1=project_id_env+'air_2000_2009_ncep1_'+ocean+'.txt' 
    158158openw,5,fi_air_ncep1 
    159  
    160  
     159; 
     160; 
    161161printf,1, 'x     y      cor    bias     std     rmsd    mean_tao' 
    162162printf,2, 'x     y      cor    bias     std     rmsd    mean_tao' 
    … …  
    164164printf,4, 'x     y      cor    bias     std     rmsd    mean_tao' 
    165165printf,5, 'x     y      cor    bias     std     rmsd    mean_tao' 
    166  
    167 ;; first reading the whole ERAI uncorrected and corrected data 
    168  
     166; 
     167; first reading the whole ERAI uncorrected and corrected data 
     168; 
    169169file=project_id_env+'erai_t2m_19890101_20091231_oafluxgrid.nc' 
    170170initncdf, file 
    171 unc=read_ncdf('t2m',date1,date2,file=file,/nostr)  
     171unc=read_ncdf('t2m',date1,date2,file=file,/nostr) 
    172172unc=unc-273.15 
    173173help, unc 
    174  
     174; 
    175175;save, unc, cor, oaf, nce, nce1, filename="airt_unc_cor_oaf_nce_nce1_2000_2009_global.idl" 
    176176restore, file="/Users/pkb/work/MY_SAXO/airt_unc_cor_oaf_nce_nce1_2000_2009_global.idl" 
    177177help, unc, cor, oaf, nce, nce1 
    178  
     178; 
    179179nn=n_elements(sitelist) 
    180180for n=0, nn-1 do begin 
    181  
    182 ;; reading data from mooring 
    183  
    184     site=sitelist(n)  
     181; 
     182; reading data from mooring 
     183; 
     184    site=sitelist(n) 
    185185csite=site 
    186186    print, csite 
    … …  
    188188    y=y_site_location(site) 
    189189    if (y ge 0. and y le 30.) then y=y+360. 
    190     dx=0.5  
    191     dy=0.5  
     190    dx=0.5 
     191    dy=0.5 
    192192    box=[y-dy, y+dy, x-dx, x+dx] 
    193193    read_variables_v2, csite,date1,date2,nsmooth, $ 
    194194         at, sw,rh,sst,wu,wv,ws,lh 
    195  
    196  
    197 ;    air=rh_to_sphum(rh,at,1008)               ;;        at  ->  air temperature 
    198     air=at                                               ;;        air ->  sea surface temperature 
    199                                                 ;;        lh  ->  latent heat flux 
    200                                                 ;;        rh  ->  relative humidity 
    201                                                 ;;   wu,wv,ws ->  wind speed 
    202 ;; 
    203  
    204  
    205 ;;  extracting the corrected and uncorrected ERAI data at the locations 
     195; 
     196; 
     197;    air=rh_to_sphum(rh,at,1008)               ;        at  ->  air temperature 
     198    air=at                                               ;        air ->  sea surface temperature 
     199                                                ;        lh  ->  latent heat flux 
     200                                                ;        rh  ->  relative humidity 
     201                                                ;   wu,wv,ws ->  wind speed 
     202; 
     203; 
     204; 
     205;  extracting the corrected and uncorrected ERAI data at the locations 
    206206    nsmooth=1. 
    207  
     207; 
    208208    extract_flux_tropflux,unc,box, $ 
    209209        tropflux 
    210210    uncr=tropflux 
    211  
     211; 
    212212     extract_flux_tropflux,cor,box, $ 
    213213        tropflux 
    214214     corr=tropflux 
    215  
     215; 
    216216     extract_flux_tropflux,oaf,box, $ 
    217217        tropflux 
    218218     oafl=tropflux 
    219  
     219; 
    220220     extract_flux_tropflux,nce,box, $ 
    221221        tropflux 
    222222     ncep=tropflux 
    223  
     223; 
    224224     extract_flux_tropflux,nce1,box, $ 
    225225        tropflux 
    226226     ncep1=tropflux 
    227  
    228  
    229     ind=where(finite(air))  
    230     air=air(ind)  
    231     uncr_air=uncr(ind)   
     227; 
     228; 
     229    ind=where(finite(air)) 
     230    air=air(ind) 
     231    uncr_air=uncr(ind) 
    232232    corr_air=corr(ind) 
    233     oafl=oafl(ind)  
    234     ncep=ncep(ind)  
     233    oafl=oafl(ind) 
     234    ncep=ncep(ind) 
    235235    ncep1=ncep1(ind) 
    236236    mean_tao=total(air,/nan)/n_elements(ind) 
    237  
     237; 
    238238    statistics_3var_v1, air, uncr_air, corr_air, $ 
    239239         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    240  
     240; 
    241241    printf, 1, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f6.2)' 
    242242    printf, 2, x, y, cor2, bias2, std2, rmsd2, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f6.2)' 
    243  
     243; 
    244244    statistics_3var_v1, air, oafl, ncep, $ 
    245245         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    246246    printf, 3, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f6.2)' 
    247247    printf, 4, x, y, cor2, bias2, std2, rmsd2, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f6.2,3x,f4.2,3x,f4.2,3x,f6.2)' 
    248  
     248; 
    249249    statistics_3var_v1, air, ncep1, ncep, $ 
    250250         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    251251    printf, 5, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f6.2)' 
    252  
     252; 
    253253endfor 
    254254close,/all 
    255  
     255; 
    256256fi_air_erai=project_id_env+'air_2000_2009_erai_'+ocean+'.txt' 
    257257res=read_ascii(fi_air_erai,data_start=1) 
    … …  
    259259lat=reform(ff(0,*)) 
    260260lon=reform(ff(1,*)) 
    261 cor_era=reform(ff(2,*))    
     261cor_era=reform(ff(2,*)) 
    262262cor_erai=total(cor_era)/n_elements(cor_era) 
    263 bias_era=reform(ff(3,*))   
     263bias_era=reform(ff(3,*)) 
    264264bias_erai=total(bias_era)/n_elements(bias_era) 
    265 std_era=reform(ff(4,*))    
     265std_era=reform(ff(4,*)) 
    266266std_erai=total(std_era)/n_elements(std_era) 
    267 rmsd_era=reform(ff(5,*))   
     267rmsd_era=reform(ff(5,*)) 
    268268rmsd_erai=total(rmsd_era)/n_elements(rmsd_era) 
    269 mean_tao=reform(ff(6,*))   
     269mean_tao=reform(ff(6,*)) 
    270270mean_erai=bias_era+mean_tao 
    271  
     271; 
    272272print, '' 
    273273print, 'ERAI' 
    … …  
    279279xyouts, 22.3,30.2, cstat, charsize=0.9 
    280280xyouts, 22.3,29.5, 'cor   bias   std   rmsd', charsize=0.9 
    281  
     281; 
    282282oplot, [22,31], [22,31] 
    283283ab=linfit(mean_tao, mean_erai,yfit=yfit) 
    284 a=float(ab(0))  
     284a=float(ab(0)) 
    285285b=float(ab(1)) 
    286286oplot, mean_tao, yfit, color=250, thick=2 
    287  
     287; 
    288288fi_air_trop=project_id_env+'air_2000_2009_trop_'+ocean+'.txt' 
    289289res=read_ascii(fi_air_trop,data_start=1) 
    … …  
    291291lat=reform(ff(0,*)) 
    292292lon=reform(ff(1,*)) 
    293 cor_tro=reform(ff(2,*))    
     293cor_tro=reform(ff(2,*)) 
    294294cor_trop=total(cor_tro)/n_elements(cor_tro) 
    295 bias_tro=reform(ff(3,*))   
     295bias_tro=reform(ff(3,*)) 
    296296bias_trop=total(bias_tro)/n_elements(bias_tro) 
    297 std_tro=reform(ff(4,*))    
     297std_tro=reform(ff(4,*)) 
    298298std_trop=total(std_tro)/n_elements(std_tro) 
    299 rmsd_tro=reform(ff(5,*))   
     299rmsd_tro=reform(ff(5,*)) 
    300300rmsd_trop=total(rmsd_tro)/n_elements(rmsd_tro) 
    301 mean_tao=reform(ff(6,*))   
     301mean_tao=reform(ff(6,*)) 
    302302mean_trop=bias_tro+mean_tao 
    303  
     303; 
    304304print, '' 
    305305print, 'TropFlux' 
    306306print, cor_trop, bias_trop, std_trop, rmsd_trop 
    307307cstat=string(cor_trop, bias_trop, std_trop, rmsd_trop, format='(f4.2,1x,f6.2,1x,f4.2,1x,f4.2)') 
    308  
     308; 
    309309splot, mean_tao, mean_trop, title='T2m - TAO Vs TropFlux', subtitle='', $ 
    310310     charsize=1.1, xtitle='TAO T2m', ytitle='TropFlux T2m', small=[2,3,2],/noer, psym=2, $ 
    … …  
    313313xyouts, 22.3,30.2, cstat, charsize=0.9 
    314314xyouts, 22.3,29.5, 'cor   bias   std   rmsd', charsize=0.9 
    315  
     315; 
    316316ab=linfit(mean_tao, mean_trop,yfit=yfit) 
    317 a=float(ab(0))  
     317a=float(ab(0)) 
    318318b=float(ab(1)) 
    319319oplot, mean_tao, yfit, color=250, thick=2 
    320  
     320; 
    321321fi_air_oaflx=project_id_env+'air_2000_2009_oaflx_'+ocean+'.txt' 
    322322res=read_ascii(fi_air_oaflx,data_start=1) 
    … …  
    324324lat=reform(ff(0,*)) 
    325325lon=reform(ff(1,*)) 
    326 cor_oaf=reform(ff(2,*))    
     326cor_oaf=reform(ff(2,*)) 
    327327cor_oafl=total(cor_oaf)/n_elements(cor_oaf) 
    328 bias_oaf=reform(ff(3,*))   
     328bias_oaf=reform(ff(3,*)) 
    329329bias_oafl=total(bias_oaf)/n_elements(bias_oaf) 
    330 std_oaf=reform(ff(4,*))    
     330std_oaf=reform(ff(4,*)) 
    331331std_oafl=total(std_oaf)/n_elements(std_oaf) 
    332 rmsd_oaf=reform(ff(5,*))   
     332rmsd_oaf=reform(ff(5,*)) 
    333333rmsd_oafl=total(rmsd_oaf)/n_elements(rmsd_oaf) 
    334 mean_tao=reform(ff(6,*))   
     334mean_tao=reform(ff(6,*)) 
    335335mean_oafl=bias_oaf+mean_tao 
    336  
     336; 
    337337print, '' 
    338338print, 'OAFlux' 
    339339print, cor_oafl, bias_oafl, std_oafl, rmsd_oafl 
    340340cstat=string(cor_oafl, bias_oafl, std_oafl, rmsd_oafl, format='(f4.2,1x,f6.2,1x,f4.2,1x,f4.2)') 
    341  
     341; 
    342342splot, mean_tao, mean_oafl, title='T2m - TAO Vs OAFlux', subtitle='', $ 
    343343     charsize=1.1, xtitle='TAO T2m', ytitle='OAFlux T2m', small=[2,3,3],/noer, psym=2, $ 
    … …  
    346346xyouts, 22.3,30.2, cstat, charsize=0.9 
    347347xyouts, 22.3,29.5, 'cor   bias   std   rmsd', charsize=0.9 
    348  
     348; 
    349349ab=linfit(mean_tao, mean_oafl,yfit=yfit) 
    350 a=float(ab(0))  
     350a=float(ab(0)) 
    351351b=float(ab(1)) 
    352352oplot, mean_tao, yfit, color=250, thick=2 
    353  
    354  
    355  
     353; 
     354; 
     355; 
    356356fi_air_ncep=project_id_env+'air_2000_2009_ncep2_'+ocean+'.txt' 
    357357res=read_ascii(fi_air_ncep,data_start=1) 
    … …  
    359359lat=reform(ff(0,*)) 
    360360lon=reform(ff(1,*)) 
    361 cor_nce=reform(ff(2,*))    
     361cor_nce=reform(ff(2,*)) 
    362362cor_ncep=total(cor_nce)/n_elements(cor_nce) 
    363 bias_nce=reform(ff(3,*))   
     363bias_nce=reform(ff(3,*)) 
    364364bias_ncep=total(bias_nce)/n_elements(bias_nce) 
    365 std_nce=reform(ff(4,*))    
     365std_nce=reform(ff(4,*)) 
    366366std_ncep=total(std_nce)/n_elements(std_nce) 
    367 rmsd_nce=reform(ff(5,*))   
     367rmsd_nce=reform(ff(5,*)) 
    368368rmsd_ncep=total(rmsd_nce)/n_elements(rmsd_nce) 
    369 mean_tao=reform(ff(6,*))   
     369mean_tao=reform(ff(6,*)) 
    370370mean_ncep=bias_nce+mean_tao 
    371  
     371; 
    372372print, '' 
    373373print, 'NCEP2' 
    374374print, cor_ncep, bias_ncep, std_ncep, rmsd_ncep 
    375375cstat=string(cor_ncep, bias_ncep, std_ncep, rmsd_ncep, format='(f4.2,1x,f6.2,1x,f4.2,1x,f4.2)') 
    376  
     376; 
    377377splot, mean_tao, mean_ncep, title='T2m - TAO Vs NCEP2', subtitle='', $ 
    378378     charsize=1.1, xtitle='TAO T2m', ytitle='NCEP2 T2m', small=[2,3,4],/noer, psym=2, $ 
    … …  
    381381xyouts, 22.3,30.2, cstat, charsize=0.9 
    382382xyouts, 22.3,29.5, 'cor   bias   std   rmsd', charsize=0.9 
    383  
     383; 
    384384ab=linfit(mean_tao, mean_ncep,yfit=yfit) 
    385 a=float(ab(0))  
     385a=float(ab(0)) 
    386386b=float(ab(1)) 
    387387oplot, mean_tao, yfit, color=250, thick=2 
    388  
     388; 
    389389fi_air_ncep1=project_id_env+'air_2000_2009_ncep1_'+ocean+'.txt' 
    390390res=read_ascii(fi_air_ncep1,data_start=1) 
    … …  
    392392lat=reform(ff(0,*)) 
    393393lon=reform(ff(1,*)) 
    394 cor_nce=reform(ff(2,*))    
     394cor_nce=reform(ff(2,*)) 
    395395cor_ncep=total(cor_nce)/n_elements(cor_nce) 
    396 bias_nce=reform(ff(3,*))   
     396bias_nce=reform(ff(3,*)) 
    397397bias_ncep=total(bias_nce)/n_elements(bias_nce) 
    398 std_nce=reform(ff(4,*))    
     398std_nce=reform(ff(4,*)) 
    399399std_ncep=total(std_nce)/n_elements(std_nce) 
    400 rmsd_nce=reform(ff(5,*))   
     400rmsd_nce=reform(ff(5,*)) 
    401401rmsd_ncep=total(rmsd_nce)/n_elements(rmsd_nce) 
    402 mean_tao=reform(ff(6,*))   
     402mean_tao=reform(ff(6,*)) 
    403403mean_ncep=bias_nce+mean_tao 
    404  
     404; 
    405405print, '' 
    406406print, 'NCEP1' 
    407407print, cor_ncep, bias_ncep, std_ncep, rmsd_ncep 
    408408cstat=string(cor_ncep, bias_ncep, std_ncep, rmsd_ncep, format='(f4.2,1x,f6.2,1x,f4.2,1x,f4.2)') 
    409  
     409; 
    410410splot, mean_tao, mean_ncep, title='T2m - TAO Vs NCEP_NCAR', subtitle='', $ 
    411411     charsize=1.1, xtitle='TAO T2m', ytitle='NCEP', small=[2,3,5],/noer, psym=2, $ 
    … …  
    414414xyouts, 22.3,30.2, cstat, charsize=0.9 
    415415xyouts, 22.3,29.5, 'cor   bias   std   rmsd', charsize=0.9 
    416  
     416; 
    417417ab=linfit(mean_tao, mean_ncep,yfit=yfit) 
    418 a=float(ab(0))  
     418a=float(ab(0)) 
    419419b=float(ab(1)) 
    420420oplot, mean_tao, yfit, color=250, thick=2 
    421  
     421; 
    422422closeps 
    423  
     423; 
    424424end 
Note: See TracChangeset for help on using the changeset viewer.