Changeset 43 for trunk/src/paper01/fig8


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/fig8
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/paper01/fig8/lhf_validation_scatter_2000_2009.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
     1;+ 
     2; .. _lhf_validation_scatter_2000_2009.pro 
     3; 
     4; ==================================== 
     5; lhf_validation_scatter_2000_2009.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> lhf_validation_scatter_2000_2009, 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;- 
    236pro lhf_validation_scatter_2000_2009,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 
    … …  
    3163ocean='global' 
    3264 
    33 ;------------------------------------------------------------------------------------------------------------------------ 
    3465;;   This program will create the following text files with statistics of respective variables 
    35 ;------------------------------------------------------------------------------------------------------------------------ 
    3666close,/all 
    3767 
    … …  
    5383printf,4, 'x     y      cor    bias     std     rmsd      mean_tao' 
    5484 
    55 ;------------------------------------------------------------------------------------------------------------------------ 
    5685;; first reading the whole ERAI uncorrected and corrected data 
    5786 
    5887file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_lhf_19890101_20091231_oafluxgrid.nc' 
    5988initncdf, file 
    60 unc=read_ncdf('lhf',date1,date2,file=file,/nostr)  
     89unc=read_ncdf('lhf',date1,date2,file=file,/nostr) 
    6190help, unc 
    6291 
    … …  
    74103fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/lhf_ncep2_oafluxgrid_19890101_20091231.nc' 
    75104initncdf, fi 
    76 nce=read_ncdf("lhf", date1, date2, file=fi,/nostr)  
     105nce=read_ncdf("lhf", date1, date2, file=fi,/nostr) 
    77106help, nce 
    78107 
    79108file='/Volumes/Iomega_HDD/flux_reconstruction/ncep/fluxe_ncep1_19890101_20091231.nc' 
    80109initncdf, file 
    81 nce1=-1*read_ncdf("lhf", date1, date2, file=file,/nostr)  
     110nce1=-1*read_ncdf("lhf", date1, date2, file=file,/nostr) 
    82111help, nce1 
    83112 
    84113nsmooth=1. 
    85 ;------------------------------------------------------------------------------------------------------------------------ 
    86114nn=n_elements(sitelist) 
    87115date1=date1 
    … …  
    101129 
    102130    lhf=lh & ind=where(finite(lhf)) & valid=n_elements(ind) 
    103   
     131 
    104132    if (valid ge 180. ) then begin 
    105133 
    … …  
    108136        uncr=tropflux 
    109137 
    110         extract_flux_tropflux,cor,box, $   
     138        extract_flux_tropflux,cor,box, $ 
    111139                 tropflux 
    112140        corr=tropflux 
    113      
     141 
    114142        extract_flux_tropflux,oaf,box, $ 
    115143                tropflux 
    … …  
    130158        statistics_3var_v1, lhf, uncr_lhf, corr_lhf, $ 
    131159                 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    132   
     160 
    133161        printf, 1, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f4.2,3x,f7.2,3x,f4.2,3x,f5.2,3x,f6.2)' 
    134162        printf, 2, x, y, cor2, bias2, std2, rmsd2, mean_tao, format='(f6.2, 3x, f6.2, 3x, f4.2,3x,f7.2,3x,f4.2,3x,f5.2,3x,f6.2)' 
    … …  
    170198 
    171199splot, mean_tao, mean_erai, charsize=1.1, title='LHF - TAO Vs ERAI', $ 
    172      xrange=[20,200], yrange=[20,200], small=[2,3,1], psym=2  
     200     xrange=[20,200], yrange=[20,200], small=[2,3,1], psym=2 
    173201xyouts, 50,25, cstat, charsize=1. 
    174202xyouts, 50,10, 'cor    bias    std    rmsd', charsize=1. 
    … …  
    305333oplot, [20,200], [20,200] 
    306334 
    307 ;----------------------------------------------------------- 
    308335closeps 
    309336 
    … …  
    313340return 
    314341end 
    315 ;-------------------------------------------------------------------------- 
    316342function x_site_location, site 
    317343    n1=strpos(site, 's') 
    … …  
    328354return, float(x) 
    329355end 
    330 ;-------------------------------------------------------------------------- 
    331356function y_site_location, site 
    332357    n1=strpos(site, 'e') 
    … …  
    354379return,float(y) 
    355380end 
    356  
    357 ;-------------------------------------------------------------------------- 
    358  

  • trunk/src/paper01/fig8/shf_validation_scatter_2000_2009.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
     1;+ 
     2; .. _shf_validation_scatter_2000_2009.pro 
     3; 
     4; ==================================== 
     5; shf_validation_scatter_2000_2009.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> shf_validation_scatter_2000_2009, 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;- 
    236pro shf_validation_scatter_2000_2009,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_variables_sh 
    … …  
    3062ocean='global' 
    3163 
    32 nsmooth=1.  
    33 ;------------------------------------------------------------------------------------------------------------------------ 
     64nsmooth=1. 
    3465;;   This program will create the following text files with statistics of respective variables 
    35 ;------------------------------------------------------------------------------------------------------------------------ 
    3666close,/all 
    3767 
    … …  
    5585printf,5, 'x     y      cor    bias     std     rmsd      mean_tao' 
    5686 
    57 ;------------------------------------------------------------------------------------------------------------------------ 
    5887;; first reading the whole ERAI uncorrected and corrected data 
    5988 
    6089file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_shf_19890101_20091231_oafluxgrid.nc' 
    6190initncdf, file 
    62 unc=read_ncdf('shf',date1,date2,file=file,/nostr)  
     91unc=read_ncdf('shf',date1,date2,file=file,/nostr) 
    6392help, unc 
    6493 
    … …  
    76105fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/shf_ncep2_oafluxgrid_19890101_20091231.nc' 
    77106initncdf, fi 
    78 nce=read_ncdf("shf", date1, date2, file=fi,/nostr)  
     107nce=read_ncdf("shf", date1, date2, file=fi,/nostr) 
    79108help, nce 
    80109 
    … …  
    84113help, nce1 
    85114 
    86 ;------------------------------------------------------------------------------------------------------------------------ 
    87115nn=n_elements(sitelist) 
    88116date1=date1 
    … …  
    102130         at, sw,rh,sst,wu,wv,ws,sh 
    103131    shf=sh & ind=where(finite(shf)) & valid=n_elements(ind) 
    104   
     132 
    105133    if (valid ge 180. ) then begin 
    106134 
    … …  
    109137        uncr=tropflux 
    110138 
    111         extract_flux_tropflux,cor,box, $   
     139        extract_flux_tropflux,cor,box, $ 
    112140                 tropflux 
    113141        corr=tropflux 
    114      
     142 
    115143        extract_flux_tropflux,oaf,box, $ 
    116144                tropflux 
    … …  
    130158        statistics_3var_v1, shf, uncr_shf, corr_shf, $ 
    131159                 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    132   
     160 
    133161        printf, 1, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f4.2,3x,f7.2,3x,f4.2,3x,f5.2,3x,f6.2)' 
    134162        printf, 2, x, y, cor2, bias2, std2, rmsd2, mean_tao, format='(f6.2, 3x, f6.2, 3x, f4.2,3x,f7.2,3x,f4.2,3x,f5.2,3x,f6.2)' 
    … …  
    305333oplot, [-5,25], [-5,25] 
    306334 
    307 ;----------------------------------------------------------- 
    308335closeps 
    309336 
    … …  
    313340return 
    314341end 
    315 ;-------------------------------------------------------------------------- 
    316342function x_site_location, site 
    317343    n1=strpos(site, 's') 
    … …  
    328354return, float(x) 
    329355end 
    330 ;-------------------------------------------------------------------------- 
    331356function y_site_location, site 
    332357    n1=strpos(site, 'e') 
    … …  
    354379return,float(y) 
    355380end 
    356  
    357 ;-------------------------------------------------------------------------- 
    358  
Note: See TracChangeset for help on using the changeset viewer.