Changeset 43 for trunk/src/paper01/fig7


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/fig7
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/paper01/fig7/air_validation_scatter_2000_2009_basin.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
     1;+ 
     2; .. _air_validation_scatter_2000_2009.pro 
     3; 
     4; ==================================== 
     5; air_validation_scatter_2000_2009.pro 
     6; ==================================== 
     7; DESCRIPTION 
     8; =========== 
     9; 
     10; SEE ALSO 
     11; ======== 
     12; 
     13; :ref:`tropflux_profile.sh` 
     14; 
     15; EXAMPLES 
     16; ======== 
     17; 
     18; :: 
     19; 
     20;  IDL> @tropflux_init 
     21;  IDL> .r read_era_total 
     22;  IDL> .r read_variables_v2 
     23;  IDL> .r read_era_box 
     24;  IDL> date1=++ 
     25;  IDL> date2=++ 
     26;  IDL> air_validation_scatter_2000_2009, date1, date2 
     27; 
     28; EVOLUTIONS 
     29; ========== 
     30; 
     31; $Id$ 
     32; 
     33; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     34; 
     35;   * minimal header 
     36; 
     37;- 
    238pro air_validation_scatter_2000_2009,date1,date2 
    339@common 
    4 ;------------------------------------------------------------ 
    540 
    641reinitplt, /z,/invert 
    … …  
    843 
    944openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    11 ;; Before running this program, you have to compile the following subrutines 
    12 ;;   1.  .r read_era_total 
    13 ;;   2.  .r read_variables_v2 
    14 ;;   3.  .r read_era_box 
    1545 
    1646;; Give the location of mooring for validation of basic meteorological variables 
    … …  
    6595 
    6696nsmooth=1.    ;; statistics are with 7 day smoothed 
    67 ;------------------------------------------------------------------------------------------------------------------------ 
    6897;;   This program will create the following text files with statistics of respective variables 
    69 ;------------------------------------------------------------------------------------------------------------------------ 
    7098close,/all 
    7199 
    … …  
    88116printf,5, 'x     y      cor    bias     std     rmsd    mean_tao' 
    89117 
    90 ;------------------------------------------------------------------------------------------------------------------------ 
    91118;; first reading the whole ERAI uncorrected and corrected data 
    92119 
    93120file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_t2m_19890101_20091231_oafluxgrid.nc' 
    94121initncdf, file 
    95 unc=read_ncdf('t2m',date1,date2,file=file,/nostr) & unc=unc-273.15  
     122unc=read_ncdf('t2m',date1,date2,file=file,/nostr) & unc=unc-273.15 
    96123help, unc 
    97124 
    … …  
    109136;fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/air_2m_ncep2_oafluxgrid_19890101_20091231.nc' 
    110137;initncdf, fi 
    111 ;nce=read_ncdf("air", date1-1, date2, file=fi,/nostr)  
     138;nce=read_ncdf("air", date1-1, date2, file=fi,/nostr) 
    112139;nce=nce-273.15 
    113140;help, nce 
    … …  
    115142;file='/Volumes/Iomega_HDD/flux_reconstruction/ncep/t2m_ncep1_19890101_20091231.nc' 
    116143;initncdf, file 
    117 ;nce1=read_ncdf("t2m", date1, date2, file=file,/nostr)  
     144;nce1=read_ncdf("t2m", date1, date2, file=file,/nostr) 
    118145;nce1=nce1-273.15 
    119146;help, nce1 
    … …  
    123150help, unc, cor, oaf, nce, nce1 
    124151 
    125 ;------------------------------------------------------------------------------------------------------------------------ 
    126152nn=n_elements(sitelist) 
    127153date1=date1 
    … …  
    145171                                                ;;        lh  ->  latent heat flux 
    146172                                                ;;        rh  ->  relative humidity 
    147                                                 ;;   wu,wv,ws ->  wind speed     
    148 ;;     
     173                                                ;;   wu,wv,ws ->  wind speed 
     174;; 
    149175 
    150176 
    … …  
    156182    uncr=tropflux 
    157183 
    158      extract_flux_tropflux,cor,box, $   
     184     extract_flux_tropflux,cor,box, $ 
    159185        tropflux 
    160186     corr=tropflux 
    161      
     187 
    162188     extract_flux_tropflux,oaf,box, $ 
    163189        tropflux 
    … …  
    173199 
    174200 
    175     ind=where(finite(air)) & air=air(ind) & uncr_air=uncr(ind)  & corr_air=corr(ind)  
     201    ind=where(finite(air)) & air=air(ind) & uncr_air=uncr(ind)  & corr_air=corr(ind) 
    176202    oafl=oafl(ind) & ncep=ncep(ind) & ncep1=ncep1(ind) 
    177203    mean_tao=total(air,/nan)/n_elements(ind) 
    178       
     204 
    179205    statistics_3var_v1, air, uncr_air, corr_air, $ 
    180206         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    181   
     207 
    182208    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)' 
    183209    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)' 
    … …  
    193219 
    194220endfor 
    195 ;----------------------------------------------------------------------------------------- 
    196 close,/all   
     221close,/all 
    197222 
    198223fi_air_erai='/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_erai_'+ocean+'.txt' 
    … …  
    332357oplot, mean_tao, yfit, color=250, thick=2 
    333358 
    334 ;----------------------------------------------------------- 
    335359closeps 
    336360 
    … …  
    340364return 
    341365end 
    342 ;-------------------------------------------------------------------------- 
    343366function x_site_location, site 
    344367    n1=strpos(site, 's') 
    … …  
    355378return, float(x) 
    356379end 
    357 ;-------------------------------------------------------------------------- 
    358380function y_site_location, site 
    359381    n1=strpos(site, 'e') 
    … …  
    381403return,float(y) 
    382404end 
    383  
    384 ;-------------------------------------------------------------------------- 
    385  

  • trunk/src/paper01/fig7/q2m_validation_scatter_2000_2009_basin.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
    2 pro q2m_validation_scatter_2000_2009_basin,date1,date2 
     1;+ 
     2; .. _q2m_validation_scatter_2000_2009_basin.pro 
     3; 
     4; ========================================== 
     5; q2m_validation_scatter_2000_2009_basin.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_basin, 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;- 
     36 
     37pro q2m_validation_scatter_2000_2009_basin, 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 
    … …  
    6598 
    6699nsmooth=1.    ;; statistics are with 7 day smoothed 
    67 ;------------------------------------------------------------------------------------------------------------------------ 
    68100;;   This program will create the following text files with statistics of respective variables 
    69 ;------------------------------------------------------------------------------------------------------------------------ 
    70101close,/all 
    71102 
    … …  
    88119printf,5, 'x     y      cor    bias     std     rmsd    mean_tao' 
    89120 
    90 ;------------------------------------------------------------------------------------------------------------------------ 
    91121;; first reading the whole ERAI uncorrected and corrected data 
    92122 
    … …  
    98128;file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_q2m_19890101_20091231_v20.nc' 
    99129;initncdf, file 
    100 ;cor=read_ncdf('q2m',date1,date2,file=file,/nostr)  
     130;cor=read_ncdf('q2m',date1,date2,file=file,/nostr) 
    101131;help, cor 
    102132 
    … …  
    108138;fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/sphum_ncep2_oafluxgrid_19890101_20091231.nc' 
    109139;initncdf, fi 
    110 ;nce=read_ncdf("shum", date1-1, date2, file=fi,/nostr)  
     140;nce=read_ncdf("shum", date1-1, date2, file=fi,/nostr) 
    111141;nce=nce*1000 
    112142;help, nce 
    … …  
    122152help, unc, cor, oaf, nce, nce1 
    123153 
    124 ;------------------------------------------------------------------------------------------------------------------------ 
    125154nn=n_elements(sitelist) 
    126155date1=date1 
    … …  
    144173                                                ;;        lh  ->  latent heat flux 
    145174                                                ;;        rh  ->  relative humidity 
    146                                                 ;;   wu,wv,ws ->  wind speed     
    147 ;;     
     175                                                ;;   wu,wv,ws ->  wind speed 
     176;; 
    148177 
    149178 
    … …  
    155184    uncr=tropflux 
    156185 
    157      extract_flux_tropflux,cor,box, $   
     186     extract_flux_tropflux,cor,box, $ 
    158187        tropflux 
    159188     corr=tropflux 
    160      
     189 
    161190     extract_flux_tropflux,oaf,box, $ 
    162191        tropflux 
    … …  
    171200     ncep1=tropflux 
    172201 
    173     ind=where(finite(q2m)) & q2m=q2m(ind) & uncr_q2m=uncr(ind)  & corr_q2m=corr(ind)  
    174     oafl=oafl(ind) & ncep=ncep(ind) & ncep1=ncep1(ind)  
     202    ind=where(finite(q2m)) & q2m=q2m(ind) & uncr_q2m=uncr(ind)  & corr_q2m=corr(ind) 
     203    oafl=oafl(ind) & ncep=ncep(ind) & ncep1=ncep1(ind) 
    175204    mean_tao=total(q2m,/nan)/n_elements(ind) 
    176       
     205 
    177206    statistics_3var_v1, q2m, uncr_q2m, corr_q2m, $ 
    178207         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    179   
     208 
    180209    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)' 
    181210    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)' 
    … …  
    191220 
    192221endfor 
    193 ;----------------------------------------------------------------------------------------- 
    194 close,/all   
     222close,/all 
    195223 
    196224fi_q2m_erai='/Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_erai_'+ocean+'.txt' 
    … …  
    329357oplot, mean_tao, yfit, color=250, thick=2 
    330358 
    331 ;----------------------------------------------------------- 
    332359closeps 
    333360 
    … …  
    337364return 
    338365end 
    339 ;-------------------------------------------------------------------------- 
    340366function x_site_location, site 
    341367    n1=strpos(site, 's') 
    … …  
    352378return, float(x) 
    353379end 
    354 ;-------------------------------------------------------------------------- 
    355380function y_site_location, site 
    356381    n1=strpos(site, 'e') 
    … …  
    378403return,float(y) 
    379404end 
    380  
    381 ;-------------------------------------------------------------------------- 
    382  

  • trunk/src/paper01/fig7/sst_validation_scatter_2000_2009_basin.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
     1;+ 
     2; .. _sst_validation_scatter_2000_2009_basin.pro 
     3; 
     4; ========================================== 
     5; sst_validation_scatter_2000_2009_basin.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_basin, 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 sst_validation_scatter_2000_2009_basin,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 
    … …  
    4376          '21n23w', '4n23w', '4n38w', '6s10w', '8n38w', '8s30w'] 
    4477 
    45 nsmooth=1.     
    46 ;------------------------------------------------------------------------------------------------------------------------ 
     78nsmooth=1. 
    4779;;   This program will create the following text files with statistics of respective variables 
    48 ;------------------------------------------------------------------------------------------------------------------------ 
    4980close,/all 
    5081 
    … …  
    71102printf,6, 'x     y      cor    bias     std     rmsd      mean_tao' 
    72103 
    73 ;------------------------------------------------------------------------------------------------------------------------ 
    74104;; first reading the whole ERAI uncorrected and corrected data 
    75105 
    76106;file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_sst_19890101_20091231_oafluxgrid.nc' 
    77107;initncdf, file 
    78 ;unc=read_ncdf('sst',date1,date2,file=file,/nostr) & unc=unc-273.15  
     108;unc=read_ncdf('sst',date1,date2,file=file,/nostr) & unc=unc-273.15 
    79109;help, unc 
    80110 
    … …  
    91121;fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/sst_ncep2_oafluxgrid_19890101_20091231.nc' 
    92122;initncdf, fi 
    93 ;nce=read_ncdf("sst", date1-1, date2, file=fi,/nostr)  
     123;nce=read_ncdf("sst", date1-1, date2, file=fi,/nostr) 
    94124;nce=nce-273.15 
    95125;help, nce 
    … …  
    112142help, unc, cor, oaf, nce, sst_tmi, nce1 
    113143;; data from 20000101, 20090931 
    114 ;------------------------------------------------------------------------------------------------------------------------ 
    115144nn=n_elements(sitelist) 
    116145date1=date1 
    … …  
    141170      uncr=tropflux 
    142171 
    143       extract_flux_tropflux,cor,box, $   
     172      extract_flux_tropflux,cor,box, $ 
    144173        tropflux 
    145174      corr=tropflux 
    146      
     175 
    147176      extract_flux_tropflux,oaf,box, $ 
    148177        tropflux 
    … …  
    162191 
    163192 
    164       ind=where(finite(sst)) & sst=sst(ind) & uncr_sst=uncr(ind)  & corr_sst=corr(ind)  
     193      ind=where(finite(sst)) & sst=sst(ind) & uncr_sst=uncr(ind)  & corr_sst=corr(ind) 
    165194      oafl=oafl(ind) & ncep=ncep(ind) & tmi=tmi(ind) & ncep1=ncep1(ind) 
    166195      mean_tao=total(sst,/nan)/n_elements(ind) 
    167       
     196 
    168197      statistics_3var_v1, sst, uncr_sst, corr_sst, $ 
    169198         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    170   
     199 
    171200      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)' 
    172201      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)' 
    … …  
    185214 
    186215endfor 
    187 ;----------------------------------------------------------------------------------------- 
    188 close,/all   
     216close,/all 
    189217 
    190218fi_sst_erai='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_erai_'+ocean+'.txt' 
    … …  
    352380oplot, mean_tao, yfit, color=250, thick=2 
    353381 
    354 ;----------------------------------------------------------- 
    355382closeps 
    356383 
    … …  
    360387return 
    361388end 
    362 ;-------------------------------------------------------------------------- 
    363389function x_site_location, site 
    364390    n1=strpos(site, 's') 
    … …  
    375401return, float(x) 
    376402end 
    377 ;-------------------------------------------------------------------------- 
    378403function y_site_location, site 
    379404    n1=strpos(site, 'e') 
    … …  
    401426return,float(y) 
    402427end 
    403  
    404 ;-------------------------------------------------------------------------- 
    405  

  • trunk/src/paper01/fig7/ws_validation_scatter_2000_2009_basin.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
     1;+ 
     2; 
     3; .. _ws_validation_scatter_2000_2009_basin.pro 
     4; 
     5; ========================================= 
     6; ws_validation_scatter_2000_2009_basin.pro 
     7; ========================================= 
     8; 
     9; DESCRIPTION 
     10; =========== 
     11; 
     12; SEE ALSO 
     13; ======== 
     14; 
     15; :ref:`tropflux_profile.sh` 
     16; 
     17; EXAMPLES 
     18; ======== 
     19; 
     20; :: 
     21; 
     22;  IDL> @tropflux_init 
     23;  IDL> date1=++ 
     24;  IDL> date2=++ 
     25;  IDL> ws_validation_scatter_2000_2009_basin, date1, date2 
     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_basin,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 
    … …  
    4982sitelist=['8s125w', '8s110w', '8s95w','5s155w', '5s140w', '5s125w', '5s110w', '5s95w','2s155w', '2s140w', $ 
    5083          '2s125w', '2s110w', '2s95w', '0n155w', '0n140w', '0n125w', '0n110w', '0n95w', '2n155w', '2n140w', $ 
    51           '2n125w', '2n110w', '2n95w',  '5n155w', '5n140w', '5n125w', '5n110w', '5n95w']  
     84          '2n125w', '2n110w', '2n95w',  '5n155w', '5n140w', '5n125w', '5n110w', '5n95w'] 
    5285ocean='Pac_coldtongue' 
    5386 
    … …  
    6598 
    6699nsmooth=1 
    67 ;------------------------------------------------------------------------------------------------------------------------ 
    68100;;   This program will create the following text files with statistics of respective variables 
    69 ;------------------------------------------------------------------------------------------------------------------------ 
    70101close,/all 
    71102 
    … …  
    91122printf,6, 'x     y      cor    bias     std     rmsd      mean_tao' 
    92123 
    93 ;------------------------------------------------------------------------------------------------------------------------ 
    94124;; first reading the whole ERAI uncorrected and corrected data 
    95125;; date1=20000101 & date2=20090931 
    … …  
    97127file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_ws_19890101_20091231_oafluxgrid.nc' 
    98128initncdf, file 
    99 u=read_ncdf('u10',date1,date2,file=file,/nostr)  
    100 ;v=read_ncdf('v10',date1,date2,file=file,/nostr)      
    101 ;unc=sqrt(u*u+v*v)  
     129u=read_ncdf('u10',date1,date2,file=file,/nostr) 
     130;v=read_ncdf('v10',date1,date2,file=file,/nostr) 
     131;unc=sqrt(u*u+v*v) 
    102132;help, unc 
    103133 
    104134;file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_ws_19890101_20091231_v20.nc' 
    105135;initncdf, file 
    106 ;cor=read_ncdf('ws',date1,date2,file=file,/nostr)  
     136;cor=read_ncdf('ws',date1,date2,file=file,/nostr) 
    107137;help, cor 
    108138 
    … …  
    114144;fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/uwind_ncep2_oafluxgrid_19890101_20091231.nc' 
    115145;initncdf, fi 
    116 ;u=read_ncdf("u", date1-1, date2, file=fi,/nostr)  
     146;u=read_ncdf("u", date1-1, date2, file=fi,/nostr) 
    117147;fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/vwind_ncep2_oafluxgrid_19890101_20091231.nc' 
    118148;initncdf, fi 
    119 ;v=read_ncdf("v", date1-1, date2, file=fi,/nostr)    
     149;v=read_ncdf("v", date1-1, date2, file=fi,/nostr) 
    120150;nce=sqrt(u*u+v*v) 
    121151;help, nce 
    … …  
    142172help, unc, cor, oaf, nce, ws_tmi, nce1 
    143173 
    144 ;------------------------------------------------------------------------------------------------------------------------ 
    145174nn=n_elements(sitelist) 
    146175date1=date1 
    … …  
    168197    uncr=tropflux 
    169198 
    170      extract_flux_tropflux,cor,box, $   
     199     extract_flux_tropflux,cor,box, $ 
    171200        tropflux 
    172201     corr=tropflux 
    173      
     202 
    174203     extract_flux_tropflux,oaf,box, $ 
    175204        tropflux 
    … …  
    188217     ncep1=tropflux 
    189218 
    190     ind=where(finite(ws)) & ws=ws(ind) & uncr_ws=uncr(ind)  & corr_ws=corr(ind)  
     219    ind=where(finite(ws)) & ws=ws(ind) & uncr_ws=uncr(ind)  & corr_ws=corr(ind) 
    191220    oafl=oafl(ind) & ncep=ncep(ind) & tmi=tmi(ind) & ncep1=ncep1(ind) 
    192221 
    193222    mean_tao=total(ws)/n_elements(ws) 
    194   
     223 
    195224    statistics_3var_v1, ws, uncr_ws, corr_ws, $ 
    196225         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    197   
     226 
    198227    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)' 
    199228    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)' 
    … …  
    210239 
    211240endfor 
    212 ;----------------------------------------------------------------------------------------- 
    213 close,/all   
     241close,/all 
    214242fi_ws_erai='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_erai_'+ocean+'.txt' 
    215243res=read_ascii(fi_ws_erai,data_start=1) 
    … …  
    376404oplot, mean_tao, yfit, color=250, thick=2 
    377405 
    378 ;----------------------------------------------------------- 
    379406closeps 
    380407 
    … …  
    384411return 
    385412end 
    386 ;-------------------------------------------------------------------------- 
    387413function x_site_location, site 
    388414    n1=strpos(site, 's') 
    … …  
    399425return, float(x) 
    400426end 
    401 ;-------------------------------------------------------------------------- 
    402427function y_site_location, site 
    403428    n1=strpos(site, 'e') 
    … …  
    425450return,float(y) 
    426451end 
    427  
    428 ;-------------------------------------------------------------------------- 
    429  
Note: See TracChangeset for help on using the changeset viewer.