Changeset 44 for trunk/src/paper01/fig3


Ignore:
Timestamp:
04/12/11 18:14:21 (13 years ago)
Author:
pinsard
Message:

start consolidation of paper01 materials

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

Legend:

Unmodified
Added
Removed

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

    r43 r44  
    11;+ 
    2 ; _air_validation_scatter_2000_2009_v50.pro: 
     2; 
     3; .. _air_validation_scatter_2000_2009_v50.pro: 
    34; 
    45; ======================================== 
    … …  
    910; =========== 
    1011; 
     12; .. graphviz:: 
     13; 
     14;    digraph air_validation_scatter_2000_2009_v50 { 
     15;       graph [ 
     16;       rankdir="LR", 
     17;       ] 
     18; 
     19;       air_erai [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/air_2000_2009_erai_v50.txt"]; 
     20; 
     21;       air_trop [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/air_2000_2009_trop_v50.txt"]; 
     22;       air_oaflux [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/air_2000_2009_oaflx_v50.txt"]; 
     23;       air_ncep  [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/air_2000_2009_ncep_v50.txt"]; 
     24;       air_ncep1  [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/air_2000_2009_ncep1_v50.txt"]; 
     25;       air_ncep2 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/air_2000_2009_ncep2_v50.txt"]; 
     26;       air_ncep2_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/air_2m_ncep2_oafluxgrid_19890101_20091231.nc"]; 
     27; 
     28;       t2m_ncep1 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/t2m_ncep1_19890101_20091231.nc"]; 
     29;       figure [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/air_validation_scatter_2000_2009_v50.ps"]; 
     30; 
     31;       erai_t2m [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_t2m_19890101_20091231_oafluxgrid.nc"]; 
     32; 
     33;       tropflux_t2m [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_t2m_19890101_20091231_v50.nc"]; 
     34;      
     35;       oaflux_basic [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/OAFlux_basic_variables_1985_2009.nc"]; 
     36; 
     37;       air_validation_scatter_2000_2009_v50 [shape=box, 
     38;       fontname=Courier, 
     39;       color=blue, 
     40;       URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/paper01/fig3/air_validation_scatter_2000_2009_v50.pro", 
     41;       label="${TROPFLUX}/src/paper01/fig3/air_validation_scatter_2000_2009_v50.pro"]; 
     42; 
     43;       {air_erai air_trop air_oaflux air_ncep air_ncep1 air_ncep2 erai_t2m tropflux_t2m oaflux_basic air_ncep2_oafluxgrid t2m_ncep1} -> {air_validation_scatter_2000_2009_v50} -> {figure} 
     44; 
     45;    } 
     46; 
    1147; SEE ALSO 
    1248; ======== 
    … …  
    1450; :ref:`tropflux_profile.sh` 
    1551; 
     52; :ref:`read_variables_v2.pro` 
     53; 
    1654; EXAMPLES 
    1755; ======== 
    … …  
    2058; 
    2159;  IDL> @tropflux_init 
    22 ;  IDL> date1=++ 
    23 ;  IDL> date2=++ 
     60;  IDL> date1=20001001L 
     61;  IDL> date2=20091231L 
    2462;  IDL> air_validation_scatter_2000_2009_v50, date1, date2 
    2563; 
     64; TODO 
     65; ==== 
     66; 
     67; make it work (missing data) 
     68; 
     69; Praveen confirm date selection 
     70; 
     71; coding rules 
     72; 
     73; get rid of spwan 
     74; 
     75; complete description 
     76; 
     77; externalize function 
     78; 
     79; handle IO error 
    2680; 
    2781; EVOLUTIONS 
    … …  
    3084; $Id$ 
    3185; 
     86; - fplod 20110412T101323Z aedon.locean-ipsl.upmc.fr (Darwin) 
     87; 
     88;   * remove hard coding path 
     89;   * add graphviz 
     90; 
    3291; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
    3392; 
    … …  
    3594; 
    3695;- 
    37  
    3896pro air_validation_scatter_2000_2009_v50, date1, date2 
    39 @common 
     97 
     98@cm_tropflux 
    4099 
    41100reinitplt, /z,/invert 
    … …  
    68127close,/all 
    69128 
    70 fi_air_erai='/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_erai_v50.txt' 
     129fi_air_erai=tropflux_id_env+'air_2000_2009_erai_v50.txt' 
    71130openw,1,fi_air_erai 
    72 fi_air_trop='/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_trop_v50.txt' 
     131fi_air_trop=tropflux_id_env+'air_2000_2009_trop_v50.txt' 
    73132openw,2,fi_air_trop 
    74 fi_air_oaflx='/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_oaflx_v50.txt' 
     133fi_air_oaflx=tropflux_id_env+'air_2000_2009_oaflx_v50.txt' 
    75134openw,3,fi_air_oaflx 
    76 fi_air_ncep2='/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_ncep2_v50.txt' 
     135fi_air_ncep2=tropflux_id_env+'air_2000_2009_ncep2_v50.txt' 
    77136openw,4,fi_air_ncep2 
    78 fi_air_ncep1='/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_ncep1_v50.txt' 
     137fi_air_ncep1=tropflux_id_env+'air_2000_2009_ncep1_v50.txt' 
    79138openw,5,fi_air_ncep1 
    80139 
    … …  
    88147;; first reading the whole ERAI uncorrected and corrected data 
    89148 
    90 file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_t2m_19890101_20091231_oafluxgrid.nc' 
     149file=tropflux_id_env+'erai_t2m_19890101_20091231_oafluxgrid.nc' 
    91150initncdf, file 
    92151unc=read_ncdf('t2m',date1,date2,file=file,/nostr) & unc=unc-273.15 
    … …  
    94153 
    95154;file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_t2m_19890101_20091231_v3.nc' 
    96 file='/Users/pkb/data/TropFlux/TropFlux_t2m_19890101_20091231_v50.nc' 
     155file=tropflux_id_env+'TropFlux_t2m_19890101_20091231_v50.nc' 
    97156initncdf, file 
    98157cor=read_ncdf('t2m',date1,date2,file=file,/nostr) & cor=cor-273.15 
    99158help, cor 
    100159 
    101 file='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/OAFlux_basic_variables_1985_2009.nc' 
     160file=tropflux_id_env+'OAFlux_basic_variables_1985_2009.nc' 
    102161initncdf, file 
    103162oaf=read_ncdf("air", date1, date2, file=file,/nostr) 
    104163help, oaf 
    105164 
    106 fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/air_2m_ncep2_oafluxgrid_19890101_20091231.nc' 
     165fi=tropflux_id_env+'air_2m_ncep2_oafluxgrid_19890101_20091231.nc' 
    107166initncdf, fi 
    108167nce=read_ncdf("air", date1-1, date2, file=fi,/nostr) 
    … …  
    110169help, nce 
    111170 
    112 file='/Volumes/Iomega_HDD/flux_reconstruction/ncep/t2m_ncep1_19890101_20091231.nc' 
     171file=tropflux_id_env+'t2m_ncep1_19890101_20091231.nc' 
    113172initncdf, file 
    114173nce1=read_ncdf("t2m", date1, date2, file=file,/nostr) 
    … …  
    187246close,/all 
    188247 
    189 fi_air_erai='/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_erai_v50.txt' 
     248fi_air_erai=tropflux_id_env+'air_2000_2009_erai_v50.txt' 
    190249res=read_ascii(fi_air_erai,data_start=1) 
    191250ff=res.field1 
    … …  
    213272oplot, mean_tao, yfit, color=250, thick=2 
    214273 
    215 fi_air_trop='/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_trop_v50.txt' 
     274fi_air_trop=tropflux_id_env+'air_2000_2009_trop_v50.txt' 
    216275res=read_ascii(fi_air_trop,data_start=1) 
    217276ff=res.field1 
    … …  
    240299oplot, mean_tao, yfit, color=250, thick=2 
    241300 
    242 fi_air_oaflx='/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_oaflx_v50.txt' 
     301fi_air_oaflx=tropflux_id_env+'air_2000_2009_oaflx_v50.txt' 
    243302res=read_ascii(fi_air_oaflx,data_start=1) 
    244303ff=res.field1 
    … …  
    269328 
    270329 
    271 fi_air_ncep='/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_ncep_v50.txt' 
     330fi_air_ncep=tropflux_id_env+'air_2000_2009_ncep_v50.txt' 
    272331res=read_ascii(fi_air_ncep,data_start=1) 
    273332ff=res.field1 
    … …  
    296355oplot, mean_tao, yfit, color=250, thick=2 
    297356 
    298 fi_air_ncep1='/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_ncep1_v50.txt' 
     357fi_air_ncep1=tropflux_id_env+'air_2000_2009_ncep1_v50.txt' 
    299358res=read_ascii(fi_air_ncep1,data_start=1) 
    300359ff=res.field1 
    … …  
    325384closeps 
    326385 
    327 fig='air_validation_scatter_2000_2009_v50.ps' 
    328 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 
    329 spawn, 'gv '+cpsdir+fig 
     386fig=tropflux_od_env+'air_validation_scatter_2000_2009_v50.ps' 
     387spawn, 'mv '+psdir+'idl.ps '+fig 
     388spawn, 'gv '+fig 
    330389return 
    331390end 

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

    r43 r44  
    11;+ 
    2 ; .. _q2m_validation_scatter_2000_2009_v50.pro 
     2; 
     3; .. _q2m_validation_scatter_2000_2009_v50.pro: 
    34; 
    45; ======================================== 
    … …  
    910; =========== 
    1011; 
     12; .. graphviz:: 
     13; 
     14;    digraph q2m_validation_scatter_2000_2009_v50 { 
     15;       graph [ 
     16;       rankdir="LR", 
     17;       ] 
     18; 
     19;       q2m_erai [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/q2m_2000_2009_erai_v50.txt"]; 
     20; 
     21;       q2m_trop [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/q2m_2000_2009_trop_v50.txt"]; 
     22;       q2m_oaflux [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/q2m_2000_2009_oaflx_v50.txt"]; 
     23;       q2m_ncep  [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/q2m_2000_2009_ncep_v50.txt"]; 
     24;       q2m_ncep1  [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/q2m_2000_2009_ncep1_v50.txt"]; 
     25; 
     26;       erai_q2m_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_t2m_19890101_20091231_oafluxgrid.nc"]; 
     27; 
     28;       tropflux_q2m [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_q2m_19890101_20091231_v50.nc"]; 
     29; 
     30;       oaflux_basic [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/OAFlux_basic_variables_1985_2009.nc"]; 
     31; 
     32;       sphum_ncep2_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/sphum_ncep2_oafluxgrid_19890101_20091231.nc"]; 
     33; 
     34;       q2m_validation_scatter_2000_2009_v50 [shape=box, 
     35;       fontname=Courier, 
     36;       color=blue, 
     37;       URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/paper01/fig3/q2m_validation_scatter_2000_2009_v50.pro", 
     38;       label="${TROPFLUX}/src/paper01/fig3/q2m_validation_scatter_2000_2009_v50.pro"]; 
     39; 
     40;       figure [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/q2m_validation_scatter_2000_2009_v50.ps"]; 
     41; 
     42;       {q2m_erai q2m_trop q2m_oaflux q2m_ncep q2m_ncep1 erai_q2m_oafluxgrid tropflux_q2m oaflux_basic sphum_ncep2_oafluxgrid} -> {q2m_validation_scatter_2000_2009_v50} -> {figure} 
     43; 
     44;     } 
     45; 
    1146; SEE ALSO 
    1247; ======== 
    … …  
    2055; 
    2156;  IDL> @tropflux_init 
    22 ;  IDL> date1=++ 
    23 ;  IDL> date2=++ 
     57;  IDL> date1=20001001L 
     58;  IDL> date2=20091231L 
    2459;  IDL> q2m_validation_scatter_2000_2009_v50, date1, date2 
     60; 
     61; TODO 
     62; ==== 
     63; 
     64; make it work (missing data) 
     65; 
     66; Praveen confirm date selection 
     67; 
     68; coding rules 
     69; 
     70; get rid of spwan 
     71; 
     72; complete description 
     73; 
     74; externalize function 
     75; 
     76; handle IO error 
    2577; 
    2678; EVOLUTIONS 
    … …  
    2981; $Id$ 
    3082; 
     83; - fplod 20110412T104125Z aedon.locean-ipsl.upmc.fr (Darwin) 
     84; 
     85;   * remove hard coding path 
     86;   * add graphviz 
     87; 
    3188; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
    3289; 
    … …  
    3592;- 
    3693pro q2m_validation_scatter_2000_2009_v50, date1, date2 
    37 @common 
     94 
     95@cm_tropflux 
    3896 
    3997reinitplt, /z,/invert 
    … …  
    66124close,/all 
    67125 
    68 fi_q2m_erai='/Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_erai_v50.txt' 
     126fi_q2m_erai='tropflux_id_env+'q2m_2000_2009_erai_v50.txt' 
    69127openw,1,fi_q2m_erai 
    70 fi_q2m_trop='/Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_trop_v50.txt' 
     128fi_q2m_trop='tropflux_id_env+'q2m_2000_2009_trop_v50.txt' 
    71129openw,2,fi_q2m_trop 
    72 fi_q2m_oaflx='/Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_oaflx_v50.txt' 
     130fi_q2m_oaflx='tropflux_id_env+'q2m_2000_2009_oaflx_v50.txt' 
    73131openw,3,fi_q2m_oaflx 
    74 fi_q2m_ncep='/Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_ncep_v50.txt' 
     132fi_q2m_ncep='tropflux_id_env+'q2m_2000_2009_ncep_v50.txt' 
    75133openw,4,fi_q2m_ncep 
    76 fi_q2m_ncep1='/Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_ncep1_v50.txt' 
     134fi_q2m_ncep1='tropflux_id_env+'q2m_2000_2009_ncep1_v50.txt' 
    77135openw,5,fi_q2m_ncep1 
    78136 
    … …  
    86144;; first reading the whole ERAI uncorrected and corrected data 
    87145 
    88 file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_q2m_19890101_20091231_oafluxgrid.nc' 
     146file=tropflux_id_env+'erai_q2m_19890101_20091231_oafluxgrid.nc' 
    89147initncdf, file 
    90148unc=read_ncdf('q2m',date1,date2,file=file,/nostr) 
    91149help, unc 
    92150 
    93 file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_q2m_19890101_20091231_v20.nc' 
     151file=tropflux_id_env+'TropFlux_q2m_19890101_20091231_v20.nc' 
    94152initncdf, file 
    95153cor=read_ncdf('q2m',date1,date2,file=file,/nostr) 
    96154help, cor 
    97155 
    98 file='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/OAFlux_basic_variables_1985_2009.nc' 
     156file=tropflux_id_env+'OAFlux_basic_variables_1985_2009.nc' 
    99157initncdf, file 
    100158oaf=read_ncdf("sphum", date1, date2, file=file,/nostr) 
    101159help, oaf 
    102160 
    103 fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/sphum_ncep2_oafluxgrid_19890101_20091231.nc' 
     161fi=tropflux_id_env+'sphum_ncep2_oafluxgrid_19890101_20091231.nc' 
    104162initncdf, fi 
    105163nce=read_ncdf("shum", date1-1, date2, file=fi,/nostr) 
    … …  
    320378closeps 
    321379 
    322 fig='q2m_validation_scatter_2000_2009_v50.ps' 
    323 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 
    324 spawn, 'gv '+cpsdir+fig 
    325 return 
     380fig=tropflux_od_env+'q2m_validation_scatter_2000_2009_v50.ps' 
     381spawn, 'mv '+psdir+'idl.ps '+fig 
     382spawn, 'gv '+fig 
    326383end 
     384 
    327385function x_site_location, site 
    328386    n1=strpos(site, 's') 

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

    r43 r44  
    11;+ 
    2 ; .. _read_variables_v2.pro 
     2; 
     3; .. _read_variables_v2.pro: 
    34; 
    45; ===================== 
    … …  
    910; =========== 
    1011; 
     12; .. graphviz:: 
     13; 
     14;    digraph read_variables_v2 { 
     15;       graph [ 
     16;       rankdir="LR", 
     17;       ] 
     18;       rad [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/rad*_dy.cdf"]; 
     19;       qlat[shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/qlat*_dy.cdf"]; 
     20;       pos[shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/pos*_dy.cdf"]; 
     21;       met[shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/met*_dy.cdf"]; 
     22; 
     23;       at [shape=tab,fontname=Courier,label="at"] 
     24;       sw [shape=tab,fontname=Courier,label="sw"] 
     25;       rh [shape=tab,fontname=Courier,label="rh"] 
     26;       sst [shape=tab,fontname=Courier,label="sst"] 
     27;       wu [shape=tab,fontname=Courier,label="wu"] 
     28;       wv [shape=tab,fontname=Courier,label="wv"] 
     29;       ws [shape=tab,fontname=Courier,label="ws"] 
     30;       lh [shape=tab,fontname=Courier,label="lh"] 
     31; 
     32;       read_variables_v2 [shape=box, 
     33;       fontname=Courier, 
     34;       color=blue, 
     35;       URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/paper01/fig2/read_variables_v2.pro", 
     36;       label="${TROPFLUX}/src/paper01/fig2/read_variables_v2.pro"]; 
     37; 
     38;       {qlat rad pos met} -> {read_variables_v2} -> {at sw rh sst wu wv ws lh} 
     39;    } 
     40; 
    1141; SEE ALSO 
    1242; ======== 
    1343; 
    1444; :ref:`tropflux_profile.sh` 
     45; 
     46; :ref:`air_validation_scatter_2000_2009_v50.pro` 
     47; :ref:`q2m_validation_scatter_2000_2009_v50.pro` 
     48; :ref:`sst_validation_scatter_2000_2009_v50.pro` 
     49; :ref:`ws_validation_scatter_2000_2009_v50.pro` 
    1550; 
    1651; EXAMPLES 
    … …  
    2257;  IDL> csite=++ 
    2358;  IDL> ++ 
    24 ;  IDL> read_variables_v2, csite,date1,date2,nsmooth, at, sw,rh,sst,wu,wv,ws, lh 
     59;  IDL> read_variables_v2, csite, date1, date2, nsmooth, at, sw, rh, sst, wu, wv, ws, lh 
     60; 
     61; TODO 
     62; ==== 
     63;  
     64; externalize function 
     65; 
     66; make it work (missing data) 
    2567; 
    2668; EVOLUTIONS 
    … …  
    2971; $Id$ 
    3072; 
     73; - fplod 20110412T110534Z aedon.locean-ipsl.upmc.fr (Darwin) 
     74; 
     75;   * remove hard coding path 
     76;   * add graphviz 
     77; 
    3178; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
    3279; 
    … …  
    3481; 
    3582;- 
    36 pro read_variables_v2, csite,date1,date2,nsmooth, $ 
    37                      at, sw,rh,sst,wu,wv,ws, lh 
    38  
    39 ; 
    40 @common 
    41  
    42 dir='/Volumes/Iomega_HDD/work/flux_reconstruction/mooriing_all/' 
     83pro read_variables_v2, csite, date1, date2, nsmooth, $ 
     84                     at, sw, rh, sst, wu, wv, ws, lh 
     85@cm_tropflux 
    4386; 
    4487;; DEFINE THE OUTPUT TIME AXIS 
    … …  
    5497;;   - if not, fill variable with missing values 
    5598; 
    56  
    57  
    5899;;SHORTWAVE 
    59 fi=dir+'rad'+csite+'_dy.cdf' 
     100fi=tropflux_id_env+'rad'+csite+'_dy.cdf' 
    60101f=file_test(fi) 
    61102sw=fltarr(jpt)+!values.f_nan 
    … …  
    69110  ind=where((tt0 ge -1e-5) and (tt0 le jpt-1+1e-5)) 
    70111  if (ind(0) ne -1) then begin 
    71   sw(tt0(ind))=sw0(ind) 
    72   sw_q(tt0(ind))=sw_q0(ind) 
    73   endif 
    74 endif 
     112     sw(tt0(ind))=sw0(ind) 
     113     sw_q(tt0(ind))=sw_q0(ind) 
     114  endif 
     115endif else begin 
     116  print, fi, ' not found' 
     117endelse 
    75118 
    76119ind_sw=where(sw_q ne 1. and sw_q ne 2.) 
    77120 
    78121;;LHF 
    79 fi=dir+'qlat'+csite+'_dy.cdf' 
     122fi=tropflux_id_env+'qlat'+csite+'_dy.cdf' 
    80123f=file_test(fi) 
    81124lh=fltarr(jpt)+!values.f_nan 
    … …  
    88131  ind=where((tt0 ge -1e-5) and (tt0 le jpt-1+1e-5)) 
    89132  if (ind(0) ne -1) then begin 
    90   lh(tt0(ind))=lh0(ind) 
    91   lh_q(tt0(ind))=lh_q0(ind) 
    92   endif 
    93 endif 
     133     lh(tt0(ind))=lh0(ind) 
     134     lh_q(tt0(ind))=lh_q0(ind) 
     135  endif 
     136endif else begin 
     137  print, fi, ' not found' 
     138endelse 
    94139 
    95140;;POSITION 
    96 fi=dir+'pos'+csite+'_dy.cdf' 
     141fi=tropflux_id_env+'pos'+csite+'_dy.cdf' 
    97142f=file_test(fi) 
    98143lat=fltarr(jpt)+!values.f_nan 
    … …  
    108153    lon(tt0(ind))=lon0(ind) 
    109154  endif 
    110 endif 
     155endif else begin 
     156  print, fi, ' not found' 
     157endelse 
    111158 
    112159 
    113160;;MET PARAMETERS 
    114 fi=dir+'met'+csite+'_dy.cdf' 
     161fi=tropflux_id_env+'met'+csite+'_dy.cdf' 
    115162f=file_test(fi) 
    116163wu=fltarr(jpt)+!values.f_nan 
    … …  
    151198    sst_q(tt0(ind))=sst_q0(ind) 
    152199  endif 
    153 endif 
     200endif else begin 
     201  print, fi, ' not found' 
     202endelse 
    154203 
    155204ind_at=where(at_q ne 1 and at_q ne 2)   &  ind_ws=where(ws_q ne 1 and ws_q ne 2) 

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

    r43 r44  
    11;+ 
    2 ; .. _sst_validation_scatter_2000_2009_v50.pro 
     2; .. _sst_validation_scatter_2000_2009_v50.pro: 
    33; 
    44; ======================================== 
    … …  
    362362closeps 
    363363 
    364 fig='sst_validation_scatter_2000_2009_v50.ps' 
    365 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 
    366 spawn, 'gv '+cpsdir+fig 
     364fig=tropflux_od_env+'sst_validation_scatter_2000_2009_v50.ps' 
     365spawn, 'mv '+psdir+'idl.ps '+fig 
     366spawn, 'gv '+fig 
    367367return 
    368368end 

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

    r43 r44  
    11;+ 
    2 ; .. _statistics_3var_v1.pro 
     2; .. _statistics_3var_v1.pro: 
    33; 
    44; ====================== 

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

    r43 r44  
    11;+ 
    2 ; .. _ws_validation_scatter_2000_2009_v50.pro 
     2; .. _ws_validation_scatter_2000_2009_v50.pro: 
    33; 
    44; ======================================= 
    … …  
    117117help, nce 
    118118 
    119 dir='/Volumes/Iomega_HDD/work/flux_reconstruction/gridded_data/' 
    120 fi=dir+'zonal_wind_speed_oafluxgrid_30N30S.nc' 
     119fi=tropflux_id_env+'zonal_wind_speed_oafluxgrid_30N30S.nc' 
    121120initncdf, fi 
    122121u=read_ncdf("u", date1, date2, file=fi,/nostr) 
    123 fi=dir+'meridional_wind_speed_oafluxgrid_30N30S.nc' 
     122fi=tropflux_id_env+'meridional_wind_speed_oafluxgrid_30N30S.nc' 
    124123initncdf, fi 
    125124v=read_ncdf("v", date1, date2, file=fi,/nostr) 
    … …  
    368367closeps 
    369368 
    370 fig='ws_validation_scatter_2000_2009_v50.ps' 
    371 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 
    372 spawn, 'gv '+cpsdir+fig 
     369fig=tropflux_od_env+'ws_validation_scatter_2000_2009_v50.ps' 
     370spawn, 'mv '+psdir+'idl.ps '+fig 
     371spawn, 'gv '+fig 
    373372return 
    374373end 
Note: See TracChangeset for help on using the changeset viewer.