Changeset 43 for trunk/src/paper01/fig3


Ignore:
Timestamp:
04/11/11 18:05:51 (13 years ago)
Author:
pinsard
Message:

add minimal header in paper01 IDL files

Location:
trunk/src/paper01/fig3
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/paper01/fig3/air_validation_scatter_2000_2009_v50.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
    2 pro air_validation_scatter_2000_2009_v50,date1,date2 
     1;+ 
     2; _air_validation_scatter_2000_2009_v50.pro: 
     3; 
     4; ======================================== 
     5; air_validation_scatter_2000_2009_v50.pro 
     6; ======================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> date1=++ 
     23;  IDL> date2=++ 
     24;  IDL> air_validation_scatter_2000_2009_v50, date1, date2 
     25; 
     26; 
     27; EVOLUTIONS 
     28; ========== 
     29; 
     30; $Id$ 
     31; 
     32; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     33; 
     34;   * minimal header 
     35; 
     36;- 
     37 
     38pro air_validation_scatter_2000_2009_v50, date1, date2 
    339@common 
    4 ;------------------------------------------------------------ 
    540 
    641reinitplt, /z,/invert 
    … …  
    843 
    944openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1145;; Before running this program, you have to compile the following subrutines 
    1246;;   1.  .r read_era_total 
    … …  
    3165 
    3266nsmooth=1.    ;; statistics are with 7 day smoothed 
    33 ;------------------------------------------------------------------------------------------------------------------------ 
    3467;;   This program will create the following text files with statistics of respective variables 
    35 ;------------------------------------------------------------------------------------------------------------------------ 
    3668close,/all 
    3769 
    … …  
    5486printf,5, 'x     y      cor    bias     std     rmsd    mean_tao' 
    5587 
    56 ;------------------------------------------------------------------------------------------------------------------------ 
    5788;; first reading the whole ERAI uncorrected and corrected data 
    5889 
    5990file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_t2m_19890101_20091231_oafluxgrid.nc' 
    6091initncdf, file 
    61 unc=read_ncdf('t2m',date1,date2,file=file,/nostr) & unc=unc-273.15  
     92unc=read_ncdf('t2m',date1,date2,file=file,/nostr) & unc=unc-273.15 
    6293help, unc 
    6394 
    … …  
    75106fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/air_2m_ncep2_oafluxgrid_19890101_20091231.nc' 
    76107initncdf, fi 
    77 nce=read_ncdf("air", date1-1, date2, file=fi,/nostr)  
     108nce=read_ncdf("air", date1-1, date2, file=fi,/nostr) 
    78109nce=nce-273.15 
    79110help, nce 
    … …  
    81112file='/Volumes/Iomega_HDD/flux_reconstruction/ncep/t2m_ncep1_19890101_20091231.nc' 
    82113initncdf, file 
    83 nce1=read_ncdf("t2m", date1, date2, file=file,/nostr)  
     114nce1=read_ncdf("t2m", date1, date2, file=file,/nostr) 
    84115nce1=nce1-273.15 
    85116help, nce1 
    86117 
    87 ;------------------------------------------------------------------------------------------------------------------------ 
    88118nn=n_elements(sitelist) 
    89119date1=date1 
    … …  
    107137                                                ;;        lh  ->  latent heat flux 
    108138                                                ;;        rh  ->  relative humidity 
    109                                                 ;;   wu,wv,ws ->  wind speed     
    110 ;;     
     139                                                ;;   wu,wv,ws ->  wind speed 
     140;; 
    111141 
    112142 
    … …  
    118148    uncr=tropflux 
    119149 
    120      extract_flux_tropflux,cor,box, $   
     150     extract_flux_tropflux,cor,box, $ 
    121151        tropflux 
    122152     corr=tropflux 
    123      
     153 
    124154     extract_flux_tropflux,oaf,box, $ 
    125155        tropflux 
    … …  
    135165 
    136166 
    137     ind=where(finite(air)) & air=air(ind) & uncr_air=uncr(ind)  & corr_air=corr(ind)  
     167    ind=where(finite(air)) & air=air(ind) & uncr_air=uncr(ind)  & corr_air=corr(ind) 
    138168    oafl=oafl(ind) & ncep=ncep(ind) & ncep1=ncep1(ind) 
    139169    mean_tao=total(air,/nan)/n_elements(ind) 
    140       
     170 
    141171    statistics_3var_v1, air, uncr_air, corr_air, $ 
    142172         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    143   
     173 
    144174    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)' 
    145175    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)' 
    … …  
    155185 
    156186endfor 
    157 ;----------------------------------------------------------------------------------------- 
    158 close,/all   
     187close,/all 
    159188 
    160189fi_air_erai='/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_erai_v50.txt' 
    … …  
    294323oplot, mean_tao, yfit, color=250, thick=2 
    295324 
    296 ;----------------------------------------------------------- 
    297325closeps 
    298326 
    … …  
    302330return 
    303331end 
    304 ;-------------------------------------------------------------------------- 
    305332function x_site_location, site 
    306333    n1=strpos(site, 's') 
    … …  
    317344return, float(x) 
    318345end 
    319 ;-------------------------------------------------------------------------- 
    320346function y_site_location, site 
    321347    n1=strpos(site, 'e') 
    … …  
    343369return,float(y) 
    344370end 
    345  
    346 ;-------------------------------------------------------------------------- 
    347  

  • trunk/src/paper01/fig3/q2m_validation_scatter_2000_2009_v50.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
    2 pro q2m_validation_scatter_2000_2009_v50,date1,date2 
     1;+ 
     2; .. _q2m_validation_scatter_2000_2009_v50.pro 
     3; 
     4; ======================================== 
     5; q2m_validation_scatter_2000_2009_v50.pro 
     6; ======================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> date1=++ 
     23;  IDL> date2=++ 
     24;  IDL> q2m_validation_scatter_2000_2009_v50, date1, date2 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
     36pro q2m_validation_scatter_2000_2009_v50, date1, date2 
    337@common 
    4 ;------------------------------------------------------------ 
    538 
    639reinitplt, /z,/invert 
    … …  
    841 
    942openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1143;; Before running this program, you have to compile the following subrutines 
    1244;;   1.  .r read_era_total 
    … …  
    3163 
    3264nsmooth=1.    ;; statistics are with 7 day smoothed 
    33 ;------------------------------------------------------------------------------------------------------------------------ 
    3465;;   This program will create the following text files with statistics of respective variables 
    35 ;------------------------------------------------------------------------------------------------------------------------ 
    3666close,/all 
    3767 
    … …  
    5484printf,5, 'x     y      cor    bias     std     rmsd    mean_tao' 
    5585 
    56 ;------------------------------------------------------------------------------------------------------------------------ 
    5786;; first reading the whole ERAI uncorrected and corrected data 
    5887 
    … …  
    6493file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_q2m_19890101_20091231_v20.nc' 
    6594initncdf, file 
    66 cor=read_ncdf('q2m',date1,date2,file=file,/nostr)  
     95cor=read_ncdf('q2m',date1,date2,file=file,/nostr) 
    6796help, cor 
    6897 
    … …  
    74103fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/sphum_ncep2_oafluxgrid_19890101_20091231.nc' 
    75104initncdf, fi 
    76 nce=read_ncdf("shum", date1-1, date2, file=fi,/nostr)  
     105nce=read_ncdf("shum", date1-1, date2, file=fi,/nostr) 
    77106nce=nce*1000 
    78107help, nce 
    … …  
    84113help, nce1 
    85114 
    86 ;------------------------------------------------------------------------------------------------------------------------ 
    87115nn=n_elements(sitelist) 
    88116date1=date1 
    … …  
    106134                                                ;;        lh  ->  latent heat flux 
    107135                                                ;;        rh  ->  relative humidity 
    108                                                 ;;   wu,wv,ws ->  wind speed     
    109 ;;     
     136                                                ;;   wu,wv,ws ->  wind speed 
     137;; 
    110138 
    111139 
    … …  
    117145    uncr=tropflux 
    118146 
    119      extract_flux_tropflux,cor,box, $   
     147     extract_flux_tropflux,cor,box, $ 
    120148        tropflux 
    121149     corr=tropflux 
    122      
     150 
    123151     extract_flux_tropflux,oaf,box, $ 
    124152        tropflux 
    … …  
    133161     ncep1=tropflux 
    134162 
    135     ind=where(finite(q2m)) & q2m=q2m(ind) & uncr_q2m=uncr(ind)  & corr_q2m=corr(ind)  
    136     oafl=oafl(ind) & ncep=ncep(ind) & ncep1=ncep1(ind)  
     163    ind=where(finite(q2m)) & q2m=q2m(ind) & uncr_q2m=uncr(ind)  & corr_q2m=corr(ind) 
     164    oafl=oafl(ind) & ncep=ncep(ind) & ncep1=ncep1(ind) 
    137165    mean_tao=total(q2m,/nan)/n_elements(ind) 
    138       
     166 
    139167    statistics_3var_v1, q2m, uncr_q2m, corr_q2m, $ 
    140168         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    141   
     169 
    142170    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)' 
    143171    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)' 
    … …  
    153181 
    154182endfor 
    155 ;----------------------------------------------------------------------------------------- 
    156 close,/all   
     183close,/all 
    157184 
    158185fi_q2m_erai='/Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_erai_v50.txt' 
    … …  
    291318oplot, mean_tao, yfit, color=250, thick=2 
    292319 
    293 ;----------------------------------------------------------- 
    294320closeps 
    295321 
    … …  
    299325return 
    300326end 
    301 ;-------------------------------------------------------------------------- 
    302327function x_site_location, site 
    303328    n1=strpos(site, 's') 
    … …  
    314339return, float(x) 
    315340end 
    316 ;-------------------------------------------------------------------------- 
    317341function y_site_location, site 
    318342    n1=strpos(site, 'e') 
    … …  
    340364return,float(y) 
    341365end 
    342  
    343 ;-------------------------------------------------------------------------- 
    344  

  • trunk/src/paper01/fig3/read_variables_v2.pro

    r41 r43  
    1 ;------------------------------------------------------------------------------------------------- 
     1;+ 
     2; .. _read_variables_v2.pro 
     3; 
     4; ===================== 
     5; read_variables_v2.pro 
     6; ===================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> csite=++ 
     23;  IDL> ++ 
     24;  IDL> read_variables_v2, csite,date1,date2,nsmooth, at, sw,rh,sst,wu,wv,ws, lh 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
    236pro read_variables_v2, csite,date1,date2,nsmooth, $ 
    337                     at, sw,rh,sst,wu,wv,ws, lh 
    438 
    539; 
    6 ;------------------------------------------------------------------------------------------------- 
    740@common 
    841 
    … …  
    75108    lon(tt0(ind))=lon0(ind) 
    76109  endif 
    77 endif  
     110endif 
    78111 
    79112 
    … …  
    122155ind_at=where(at_q ne 1 and at_q ne 2)   &  ind_ws=where(ws_q ne 1 and ws_q ne 2) 
    123156ind_rh=where(rh_q ne 1 and rh_q ne 2)   &  ind_sst=where(sst_q ne 1 and sst_q ne 2) 
    124 ind_lh=where(lh_q ne 1 and lh_q ne 2)  
    125   
     157ind_lh=where(lh_q ne 1 and lh_q ne 2) 
     158 
    126159ind=union(ind_at, union(ind_rh, union(ind_ws, union(ind_lh, ind_sst)))) 
    127160 
    … …  
    135168;; Replace missing values by "NaN" 
    136169; 
    137 tsvars=['at','sw','rh','sst','wu','wv','ws','lat','lon']  
     170tsvars=['at','sw','rh','sst','wu','wv','ws','lat','lon'] 
    138171vars=[tsvars] 
    139172nn=n_elements(vars) 
    … …  
    150183ws=smooth(ws,nsmooth,/nan) & sw=smooth(sw,nsmooth,/nan) 
    151184 
    152 ;------------------------------------------------------------------------------------------------- 
    153185end 
    154 ;------------------------------------------------------------------------------------------------- 
    155186 
    156187 
    157 ;------------------------------------------------------------------------------------------------- 
    158188function time_lec, fi 
    159189tt=ncdf_lec(fi,var='time') 
    … …  
    167197return, tt 
    168198end 
    169  
    170 ;------------------------------------------------------------------------------------------------- 

  • trunk/src/paper01/fig3/sst_validation_scatter_2000_2009_v50.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
    2 pro sst_validation_scatter_2000_2009_v50,date1,date2 
     1;+ 
     2; .. _sst_validation_scatter_2000_2009_v50.pro 
     3; 
     4; ======================================== 
     5; sst_validation_scatter_2000_2009_v50.pro 
     6; ======================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> date1=++ 
     23;  IDL> date2=++ 
     24;  IDL> sst_validation_scatter_2000_2009_v50, date1, date2 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
     36pro sst_validation_scatter_2000_2009_v50, date1, date2 
    337@common 
    4 ;------------------------------------------------------------ 
    538 
    639reinitplt, /z,/invert 
    … …  
    841 
    942openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1143;; Before running this program, you have to compile the following subrutines 
    1244;;   1.  .r read_era_total 
    … …  
    3163 
    3264nsmooth=1.    ;; statistics are with 7 day smoothed 
    33 ;------------------------------------------------------------------------------------------------------------------------ 
    3465;;   This program will create the following text files with statistics of respective variables 
    35 ;------------------------------------------------------------------------------------------------------------------------ 
    3666close,/all 
    3767 
    … …  
    5888printf,6, 'x     y      cor    bias     std     rmsd      mean_tao' 
    5989 
    60 ;------------------------------------------------------------------------------------------------------------------------ 
    6190;; first reading the whole ERAI uncorrected and corrected data 
    6291 
    6392file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_sst_19890101_20091231_oafluxgrid.nc' 
    6493initncdf, file 
    65 unc=read_ncdf('sst',date1,date2,file=file,/nostr) & unc=unc-273.15  
     94unc=read_ncdf('sst',date1,date2,file=file,/nostr) & unc=unc-273.15 
    6695help, unc 
    6796 
    … …  
    78107fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/sst_ncep2_oafluxgrid_19890101_20091231.nc' 
    79108initncdf, fi 
    80 nce=read_ncdf("sst", date1-1, date2, file=fi,/nostr)  
     109nce=read_ncdf("sst", date1-1, date2, file=fi,/nostr) 
    81110nce=nce-273.15 
    82111help, nce 
    … …  
    93122help, nce1 
    94123 
    95 ;------------------------------------------------------------------------------------------------------------------------ 
    96124nn=n_elements(sitelist) 
    97125date1=date1 
    … …  
    122150      uncr=tropflux 
    123151 
    124       extract_flux_tropflux,cor,box, $   
     152      extract_flux_tropflux,cor,box, $ 
    125153        tropflux 
    126154      corr=tropflux 
    127      
     155 
    128156      extract_flux_tropflux,oaf,box, $ 
    129157        tropflux 
    … …  
    143171 
    144172 
    145       ind=where(finite(sst)) & sst=sst(ind) & uncr_sst=uncr(ind)  & corr_sst=corr(ind)  
     173      ind=where(finite(sst)) & sst=sst(ind) & uncr_sst=uncr(ind)  & corr_sst=corr(ind) 
    146174      oafl=oafl(ind) & ncep=ncep(ind) & tmi=tmi(ind) & ncep1=ncep1(ind) 
    147175      mean_tao=total(sst,/nan)/n_elements(ind) 
    148       
     176 
    149177      statistics_3var_v1, sst, uncr_sst, corr_sst, $ 
    150178         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    151   
     179 
    152180      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)' 
    153181      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)' 
    … …  
    166194 
    167195endfor 
    168 ;----------------------------------------------------------------------------------------- 
    169 close,/all   
     196close,/all 
    170197 
    171198fi_sst_erai='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_v50_erai.txt' 
    … …  
    333360oplot, mean_tao, yfit, color=250, thick=2 
    334361 
    335 ;----------------------------------------------------------- 
    336362closeps 
    337363 
    … …  
    341367return 
    342368end 
    343 ;-------------------------------------------------------------------------- 
    344369function x_site_location, site 
    345370    n1=strpos(site, 's') 
    … …  
    356381return, float(x) 
    357382end 
    358 ;-------------------------------------------------------------------------- 
    359383function y_site_location, site 
    360384    n1=strpos(site, 'e') 
    … …  
    382406return,float(y) 
    383407end 
    384  
    385 ;-------------------------------------------------------------------------- 
    386  

  • trunk/src/paper01/fig3/statistics_3var_v1.pro

    r41 r43  
     1;+ 
     2; .. _statistics_3var_v1.pro 
     3; 
     4; ====================== 
     5; statistics_3var_v1.pro 
     6; ====================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> tao=++ 
     23;  IDL> ++ 
     24;  IDL> statistics_3var_v1, tao, var1, var2, cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
     36 
     37 
    138pro statistics_3var_v1, tao,var1,var2, $ 
    239    cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    340 
    4   
     41 
    542@common 
    6 ;-------------------------------------------- 
    7 ;-------------------------------------------- 
    843x=tao & y=var1 & z=var2 
    944 
    … …  
    1247 
    1348if (ind ne -1.) then begin 
    14    x(ind)=!Values.f_nan & y(ind)=!Values.f_nan & z(ind)=!Values.f_nan  
     49   x(ind)=!Values.f_nan & y(ind)=!Values.f_nan & z(ind)=!Values.f_nan 
    1550endif 
    1651 
    … …  
    4277 
    4378end 
    44 ;--------------------------------------------    
    45  

  • trunk/src/paper01/fig3/ws_validation_scatter_2000_2009_v50.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
     1;+ 
     2; .. _ws_validation_scatter_2000_2009_v50.pro 
     3; 
     4; ======================================= 
     5; ws_validation_scatter_2000_2009_v50.pro 
     6; ======================================= 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> date1=++ 
     23;  IDL> date2=++ 
     24;  IDL> ws_validation_scatter_2000_2009_v50, date1, date2 
     25; 
     26; 
     27; EVOLUTIONS 
     28; ========== 
     29; 
     30; $Id$ 
     31; 
     32; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     33; 
     34;   * minimal header 
     35; 
     36;- 
    237pro ws_validation_scatter_2000_2009_v50,date1,date2 
    338@common 
    4 ;------------------------------------------------------------ 
    539 
    640reinitplt, /z,/invert 
    … …  
    842 
    943openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1144;; Before running this program, you have to compile the following subrutines 
    1245;;   1.  .r read_era_total 
    … …  
    3265da1=10000101 & da2=10081231 
    3366nsmooth=1.    ;; statistics are with 7 day smoothed 
    34 ;------------------------------------------------------------------------------------------------------------------------ 
    3567;;   This program will create the following text files with statistics of respective variables 
    36 ;------------------------------------------------------------------------------------------------------------------------ 
    3768close,/all 
    3869 
    … …  
    5889printf,6, 'x     y      cor    bias     std     rmsd      mean_tao' 
    5990 
    60 ;------------------------------------------------------------------------------------------------------------------------ 
    6191;; first reading the whole ERAI uncorrected and corrected data 
    6292 
    6393file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_ws_19890101_20091231_oafluxgrid.nc' 
    6494initncdf, file 
    65 u=read_ncdf('u10',date1,date2,file=file,/nostr)  
    66 v=read_ncdf('v10',date1,date2,file=file,/nostr)      
    67 unc=sqrt(u*u+v*v)  
     95u=read_ncdf('u10',date1,date2,file=file,/nostr) 
     96v=read_ncdf('v10',date1,date2,file=file,/nostr) 
     97unc=sqrt(u*u+v*v) 
    6898help, unc 
    6999 
    70100file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_ws_19890101_20091231_v20.nc' 
    71101initncdf, file 
    72 cor=read_ncdf('ws',date1,date2,file=file,/nostr)  
     102cor=read_ncdf('ws',date1,date2,file=file,/nostr) 
    73103help, cor 
    74104 
    … …  
    80110fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/uwind_ncep2_oafluxgrid_19890101_20091231.nc' 
    81111initncdf, fi 
    82 u=read_ncdf("u", date1-1, date2, file=fi,/nostr)  
     112u=read_ncdf("u", date1-1, date2, file=fi,/nostr) 
    83113fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/vwind_ncep2_oafluxgrid_19890101_20091231.nc' 
    84114initncdf, fi 
    85 v=read_ncdf("v", date1-1, date2, file=fi,/nostr)    
     115v=read_ncdf("v", date1-1, date2, file=fi,/nostr) 
    86116nce=sqrt(u*u+v*v) 
    87117help, nce 
    … …  
    104134help, nce1 
    105135 
    106 ;------------------------------------------------------------------------------------------------------------------------ 
    107136nn=n_elements(sitelist) 
    108137date1=date1 
    … …  
    130159    uncr=tropflux 
    131160 
    132      extract_flux_tropflux,cor,box, $   
     161     extract_flux_tropflux,cor,box, $ 
    133162        tropflux 
    134163     corr=tropflux 
    135      
     164 
    136165     extract_flux_tropflux,oaf,box, $ 
    137166        tropflux 
    … …  
    150179     ncep1=tropflux 
    151180 
    152     ind=where(finite(ws)) & ws=ws(ind) & uncr_ws=uncr(ind)  & corr_ws=corr(ind)  
     181    ind=where(finite(ws)) & ws=ws(ind) & uncr_ws=uncr(ind)  & corr_ws=corr(ind) 
    153182    oafl=oafl(ind) & ncep=ncep(ind) & tmi=tmi(ind) & ncep1=ncep1(ind) 
    154183 
    155184    mean_tao=total(ws)/n_elements(ws) 
    156   
     185 
    157186    statistics_3var_v1, ws, uncr_ws, corr_ws, $ 
    158187         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    159   
     188 
    160189    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,f5.2)' 
    161190    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,f5.2)' 
    … …  
    172201 
    173202endfor 
    174 ;----------------------------------------------------------------------------------------- 
    175 close,/all   
     203close,/all 
    176204fi_ws_erai='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_erai_v50.txt' 
    177205res=read_ascii(fi_ws_erai,data_start=1) 
    … …  
    338366oplot, mean_tao, yfit, color=250, thick=2 
    339367 
    340 ;----------------------------------------------------------- 
    341368closeps 
    342369 
    … …  
    346373return 
    347374end 
    348 ;-------------------------------------------------------------------------- 
    349375function x_site_location, site 
    350376    n1=strpos(site, 's') 
    … …  
    361387return, float(x) 
    362388end 
    363 ;-------------------------------------------------------------------------- 
    364389function y_site_location, site 
    365390    n1=strpos(site, 'e') 
    … …  
    387412return,float(y) 
    388413end 
    389  
    390 ;-------------------------------------------------------------------------- 
    391  
Note: See TracChangeset for help on using the changeset viewer.