Changeset 50 for trunk/src/paper01/fig3


Ignore:
Timestamp:
04/29/11 11:09:12 (13 years ago)
Author:
pinsard
Message:

replace TROPFLUX by PROJECT

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

Legend:

Unmodified
Added
Removed

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

    r46 r50  
    1717;       ] 
    1818; 
    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_OD}/air_2000_2009_trop_v50.txt"]; 
    22 ;       air_oaflux [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/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_OD}/air_2000_2009_ncep1_v50.txt"]; 
    25 ;       air_ncep2 [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/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"]; 
     19;       air_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/air_2000_2009_erai_v50.txt"]; 
     20; 
     21;       air_trop [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/air_2000_2009_trop_v50.txt"]; 
     22;       air_oaflux [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/air_2000_2009_oaflx_v50.txt"]; 
     23;       air_ncep  [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/air_2000_2009_ncep_v50.txt"]; 
     24;       air_ncep1  [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/air_2000_2009_ncep1_v50.txt"]; 
     25;       air_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/air_2000_2009_ncep2_v50.txt"]; 
     26;       air_ncep2_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/air_2m_ncep2_oafluxgrid_19890101_20091231.nc"]; 
     27; 
     28;       t2m_ncep1 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/t2m_ncep1_19890101_20091231.nc"]; 
     29;       figure [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/air_validation_scatter_2000_2009_v50.ps"]; 
     30; 
     31;       erai_t2m [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_t2m_19890101_20091231_oafluxgrid.nc"]; 
     32; 
     33;       tropflux_t2m [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/TropFlux_t2m_19890101_20091231_v50.nc"]; 
     34; 
     35;       oaflux_basic [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/OAFlux_basic_variables_1985_2009.nc"]; 
    3636; 
    3737;       air_validation_scatter_2000_2009_v50 [shape=box, 
    … …  
    4848; ======== 
    4949; 
    50 ; :ref:`tropflux_profile.sh` 
    51 ; :ref:`tropflux_init.pro` 
    52 ; :ref:`cm_tropflux.pro` 
     50; :ref:`project_profile.sh` 
     51; :ref:`project_init.pro` 
     52; :ref:`cm_project.pro` 
    5353; 
    5454; :func:`x_site_location` 
    … …  
    6363; :: 
    6464; 
    65 ;  IDL> @tropflux_init 
    6665;  IDL> date1=20001001L 
    6766;  IDL> date2=20091231L 
    … …  
    102101pro air_validation_scatter_2000_2009_v50, date1, date2 
    103102 
    104 @cm_tropflux 
     103@cm_project 
    105104 
    106105reinitplt, /z,/invert 
    … …  
    129128close,/all 
    130129 
    131 fi_air_erai=tropflux_id_env+'air_2000_2009_erai_v50.txt' 
     130fi_air_erai=project_id_env+'air_2000_2009_erai_v50.txt' 
    132131openw,1,fi_air_erai 
    133 fi_air_trop=tropflux_id_env+'air_2000_2009_trop_v50.txt' 
     132fi_air_trop=project_id_env+'air_2000_2009_trop_v50.txt' 
    134133openw,2,fi_air_trop 
    135 fi_air_oaflx=tropflux_id_env+'air_2000_2009_oaflx_v50.txt' 
     134fi_air_oaflx=project_id_env+'air_2000_2009_oaflx_v50.txt' 
    136135openw,3,fi_air_oaflx 
    137 fi_air_ncep2=tropflux_id_env+'air_2000_2009_ncep2_v50.txt' 
     136fi_air_ncep2=project_id_env+'air_2000_2009_ncep2_v50.txt' 
    138137openw,4,fi_air_ncep2 
    139 fi_air_ncep1=tropflux_id_env+'air_2000_2009_ncep1_v50.txt' 
     138fi_air_ncep1=project_id_env+'air_2000_2009_ncep1_v50.txt' 
    140139openw,5,fi_air_ncep1 
    141140 
    … …  
    149148;; first reading the whole ERAI uncorrected and corrected data 
    150149 
    151 file=tropflux_id_env+'erai_t2m_19890101_20091231_oafluxgrid.nc' 
     150file=project_id_env+'erai_t2m_19890101_20091231_oafluxgrid.nc' 
    152151initncdf, file 
    153152unc=read_ncdf('t2m',date1,date2,file=file,/nostr) & unc=unc-273.15 
    154153help, unc 
    155154 
    156 file=tropflux_id_env+'TropFlux_t2m_19890101_20091231_v50.nc' 
     155file=project_id_env+'TropFlux_t2m_19890101_20091231_v50.nc' 
    157156initncdf, file 
    158157cor=read_ncdf('t2m',date1,date2,file=file,/nostr) & cor=cor-273.15 
    159158help, cor 
    160159 
    161 file=tropflux_id_env+'OAFlux_basic_variables_1985_2009.nc' 
     160file=project_id_env+'OAFlux_basic_variables_1985_2009.nc' 
    162161initncdf, file 
    163162oaf=read_ncdf("air", date1, date2, file=file,/nostr) 
    164163help, oaf 
    165164 
    166 fi=tropflux_id_env+'air_2m_ncep2_oafluxgrid_19890101_20091231.nc' 
     165fi=project_id_env+'air_2m_ncep2_oafluxgrid_19890101_20091231.nc' 
    167166initncdf, fi 
    168167nce=read_ncdf("air", date1-1, date2, file=fi,/nostr) 
    … …  
    170169help, nce 
    171170 
    172 file=tropflux_id_env+'t2m_ncep1_19890101_20091231.nc' 
     171file=project_id_env+'t2m_ncep1_19890101_20091231.nc' 
    173172initncdf, file 
    174173nce1=read_ncdf("t2m", date1, date2, file=file,/nostr) 
    … …  
    245244close,/all 
    246245 
    247 fi_air_erai=tropflux_id_env+'air_2000_2009_erai_v50.txt' 
     246fi_air_erai=project_id_env+'air_2000_2009_erai_v50.txt' 
    248247res=read_ascii(fi_air_erai,data_start=1) 
    249248ff=res.field1 
    … …  
    271270oplot, mean_tao, yfit, color=250, thick=2 
    272271 
    273 fi_air_trop=tropflux_id_env+'air_2000_2009_trop_v50.txt' 
     272fi_air_trop=project_id_env+'air_2000_2009_trop_v50.txt' 
    274273res=read_ascii(fi_air_trop,data_start=1) 
    275274ff=res.field1 
    … …  
    298297oplot, mean_tao, yfit, color=250, thick=2 
    299298 
    300 fi_air_oaflx=tropflux_id_env+'air_2000_2009_oaflx_v50.txt' 
     299fi_air_oaflx=project_id_env+'air_2000_2009_oaflx_v50.txt' 
    301300res=read_ascii(fi_air_oaflx,data_start=1) 
    302301ff=res.field1 
    … …  
    327326 
    328327 
    329 fi_air_ncep=tropflux_id_env+'air_2000_2009_ncep_v50.txt' 
     328fi_air_ncep=project_id_env+'air_2000_2009_ncep_v50.txt' 
    330329res=read_ascii(fi_air_ncep,data_start=1) 
    331330ff=res.field1 
    … …  
    354353oplot, mean_tao, yfit, color=250, thick=2 
    355354 
    356 fi_air_ncep1=tropflux_id_env+'air_2000_2009_ncep1_v50.txt' 
     355fi_air_ncep1=project_id_env+'air_2000_2009_ncep1_v50.txt' 
    357356res=read_ascii(fi_air_ncep1,data_start=1) 
    358357ff=res.field1 
    … …  
    383382closeps 
    384383 
    385 fig=tropflux_od_env+'air_validation_scatter_2000_2009_v50.ps' 
     384fig=project_od_env+'air_validation_scatter_2000_2009_v50.ps' 
    386385spawn, 'mv '+psdir+'idl.ps '+fig 
    387386spawn, 'gv '+fig 

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

    r47 r50  
    1717;       ] 
    1818; 
    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_OD}/q2m_2000_2009_trop_v50.txt"]; 
    22 ;       q2m_oaflux [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/q2m_2000_2009_oaflx_v50.txt"]; 
    23 ;       q2m_ncep  [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/q2m_2000_2009_ncep_v50.txt"]; 
    24 ;       q2m_ncep1  [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/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"]; 
     19;       q2m_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/q2m_2000_2009_erai_v50.txt"]; 
     20; 
     21;       q2m_trop [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/q2m_2000_2009_trop_v50.txt"]; 
     22;       q2m_oaflux [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/q2m_2000_2009_oaflx_v50.txt"]; 
     23;       q2m_ncep  [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/q2m_2000_2009_ncep_v50.txt"]; 
     24;       q2m_ncep1  [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/q2m_2000_2009_ncep1_v50.txt"]; 
     25; 
     26;       erai_q2m_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_t2m_19890101_20091231_oafluxgrid.nc"]; 
     27; 
     28;       tropflux_q2m [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/TropFlux_q2m_19890101_20091231_v50.nc"]; 
     29; 
     30;       oaflux_basic [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/OAFlux_basic_variables_1985_2009.nc"]; 
     31; 
     32;       sphum_ncep2_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/sphum_ncep2_oafluxgrid_19890101_20091231.nc"]; 
    3333; 
    3434;       q2m_validation_scatter_2000_2009_v50 [shape=box, 
    … …  
    3838;       label="${TROPFLUX}/src/paper01/fig3/q2m_validation_scatter_2000_2009_v50.pro"]; 
    3939; 
    40 ;       figure [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/q2m_validation_scatter_2000_2009_v50.ps"]; 
     40;       figure [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/q2m_validation_scatter_2000_2009_v50.ps"]; 
    4141; 
    4242;       {q2m_erai erai_q2m_oafluxgrid tropflux_q2m oaflux_basic sphum_ncep2_oafluxgrid} -> {q2m_validation_scatter_2000_2009_v50} -> {q2m_trop q2m_oaflux q2m_ncep q2m_ncep1 figure} 
    … …  
    4747; ======== 
    4848; 
    49 ; :ref:`tropflux_profile.sh` 
    50 ; :ref:`tropflux_init.pro` 
    51 ; :ref:`cm_tropflux.pro` 
     49; :ref:`project_profile.sh` 
     50; :ref:`project_init.pro` 
     51; :ref:`cm_project.pro` 
    5252; 
    5353; :func:`x_site_location` 
    … …  
    6262; :: 
    6363; 
    64 ;  IDL> @tropflux_init 
    6564;  IDL> date1=20001001L 
    6665;  IDL> date2=20091231L 
    … …  
    104103pro q2m_validation_scatter_2000_2009_v50, date1, date2 
    105104 
    106 @cm_tropflux 
     105@cm_project 
    107106 
    108107reinitplt, /z,/invert 
    … …  
    131130close,/all 
    132131 
    133 fi_q2m_erai=tropflux_id_env+'q2m_2000_2009_erai_v50.txt' 
     132fi_q2m_erai=project_id_env+'q2m_2000_2009_erai_v50.txt' 
    134133openw,1,fi_q2m_erai 
    135 fi_q2m_trop=tropflux_id_env+'q2m_2000_2009_trop_v50.txt' 
     134fi_q2m_trop=project_id_env+'q2m_2000_2009_trop_v50.txt' 
    136135openw,2,fi_q2m_trop 
    137 fi_q2m_oaflx=tropflux_id_env+'q2m_2000_2009_oaflx_v50.txt' 
     136fi_q2m_oaflx=project_id_env+'q2m_2000_2009_oaflx_v50.txt' 
    138137openw,3,fi_q2m_oaflx 
    139 fi_q2m_ncep=tropflux_id_env+'q2m_2000_2009_ncep_v50.txt' 
     138fi_q2m_ncep=project_id_env+'q2m_2000_2009_ncep_v50.txt' 
    140139openw,4,fi_q2m_ncep 
    141 fi_q2m_ncep1=tropflux_id_env+'q2m_2000_2009_ncep1_v50.txt' 
     140fi_q2m_ncep1=project_id_env+'q2m_2000_2009_ncep1_v50.txt' 
    142141openw,5,fi_q2m_ncep1 
    143142 
    … …  
    151150;; first reading the whole ERAI uncorrected and corrected data 
    152151 
    153 file=tropflux_id_env+'erai_q2m_19890101_20091231_oafluxgrid.nc' 
     152file=project_id_env+'erai_q2m_19890101_20091231_oafluxgrid.nc' 
    154153initncdf, file 
    155154unc=read_ncdf('q2m',date1,date2,file=file,/nostr) 
    156155help, unc 
    157156 
    158 file=tropflux_id_env+'TropFlux_q2m_19890101_20091231_v20.nc' 
     157file=project_id_env+'TropFlux_q2m_19890101_20091231_v20.nc' 
    159158initncdf, file 
    160159cor=read_ncdf('q2m',date1,date2,file=file,/nostr) 
    161160help, cor 
    162161 
    163 file=tropflux_id_env+'OAFlux_basic_variables_1985_2009.nc' 
     162file=project_id_env+'OAFlux_basic_variables_1985_2009.nc' 
    164163initncdf, file 
    165164oaf=read_ncdf("sphum", date1, date2, file=file,/nostr) 
    166165help, oaf 
    167166 
    168 fi=tropflux_id_env+'sphum_ncep2_oafluxgrid_19890101_20091231.nc' 
     167fi=project_id_env+'sphum_ncep2_oafluxgrid_19890101_20091231.nc' 
    169168initncdf, fi 
    170169nce=read_ncdf("shum", date1-1, date2, file=fi,/nostr) 
    … …  
    172171help, nce 
    173172 
    174 file=tropflux_id_env+'q2m_ncep1_19890101_20091231.nc' 
     173file=project_id_env+'q2m_ncep1_19890101_20091231.nc' 
    175174initncdf, file 
    176175nce1=read_ncdf("q2m", date1, date2, file=file,/nostr) 
    … …  
    246245close,/all 
    247246 
    248 fi_q2m_erai=tropflux_id_env+'q2m_2000_2009_erai_v50.txt' 
     247fi_q2m_erai=project_id_env+'q2m_2000_2009_erai_v50.txt' 
    249248res=read_ascii(fi_q2m_erai,data_start=1) 
    250249ff=res.field1 
    … …  
    272271oplot, mean_tao, yfit, color=250, thick=2 
    273272 
    274 fi_q2m_trop=tropflux_id_env+'q2m_2000_2009_trop_v50.txt' 
     273fi_q2m_trop=project_id_env+'q2m_2000_2009_trop_v50.txt' 
    275274res=read_ascii(fi_q2m_trop,data_start=1) 
    276275ff=res.field1 
    … …  
    299298oplot, mean_tao, yfit, color=250, thick=2 
    300299 
    301 fi_q2m_oaflx=tropflux_id_env+'q2m_2000_2009_oaflx_v50.txt' 
     300fi_q2m_oaflx=project_id_env+'q2m_2000_2009_oaflx_v50.txt' 
    302301res=read_ascii(fi_q2m_oaflx,data_start=1) 
    303302ff=res.field1 
    … …  
    327326 
    328327 
    329 fi_q2m_ncep=tropflux_id_env+'q2m_2000_2009_ncep_v50.txt' 
     328fi_q2m_ncep=project_id_env+'q2m_2000_2009_ncep_v50.txt' 
    330329res=read_ascii(fi_q2m_ncep,data_start=1) 
    331330ff=res.field1 
    … …  
    354353oplot, mean_tao, yfit, color=250, thick=2 
    355354 
    356 fi_q2m_ncep1=tropflux_id_env+'q2m_2000_2009_ncep1_v50.txt' 
     355fi_q2m_ncep1=project_id_env+'q2m_2000_2009_ncep1_v50.txt' 
    357356res=read_ascii(fi_q2m_ncep1,data_start=1) 
    358357ff=res.field1 
    … …  
    383382closeps 
    384383 
    385 fig=tropflux_od_env+'q2m_validation_scatter_2000_2009_v50.ps' 
     384fig=project_od_env+'q2m_validation_scatter_2000_2009_v50.ps' 
    386385spawn, 'mv '+psdir+'idl.ps '+fig 
    387386spawn, 'gv '+fig 

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

    r46 r50  
    1616;       rankdir="LR", 
    1717;       ] 
    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"]; 
     18;       rad [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/rad*_dy.cdf"]; 
     19;       qlat[shape=ellipse,fontname=Courier,label="${PROJECT_ID}/qlat*_dy.cdf"]; 
     20;       pos[shape=ellipse,fontname=Courier,label="${PROJECT_ID}/pos*_dy.cdf"]; 
     21;       met[shape=ellipse,fontname=Courier,label="${PROJECT_ID}/met*_dy.cdf"]; 
    2222; 
    2323;       at [shape=tab,fontname=Courier,label="at"] 
    … …  
    4242; ======== 
    4343; 
    44 ; :ref:`tropflux_profile.sh` 
     44; :ref:`project_profile.sh` 
    4545; 
    4646; :ref:`air_validation_scatter_2000_2009_v50.pro` 
    … …  
    6060; :: 
    6161; 
    62 ;  IDL> @tropflux_init 
    6362;  IDL> csite=++ 
    6463;  IDL> ++ 
    … …  
    8988pro read_variables_v2, csite, date1, date2, nsmooth, $ 
    9089                     at, sw, rh, sst, wu, wv, ws, lh 
    91 @cm_tropflux 
     90@cm_project 
    9291; 
    9392;; DEFINE THE OUTPUT TIME AXIS 
    … …  
    104103; 
    105104;;SHORTWAVE 
    106 fi=tropflux_id_env+'rad'+csite+'_dy.cdf' 
     105fi=project_id_env+'rad'+csite+'_dy.cdf' 
    107106f=file_test(fi) 
    108107sw=fltarr(jpt)+!values.f_nan 
    … …  
    126125 
    127126;;LHF 
    128 fi=tropflux_id_env+'qlat'+csite+'_dy.cdf' 
     127fi=project_id_env+'qlat'+csite+'_dy.cdf' 
    129128f=file_test(fi) 
    130129lh=fltarr(jpt)+!values.f_nan 
    … …  
    145144 
    146145;;POSITION 
    147 fi=tropflux_id_env+'pos'+csite+'_dy.cdf' 
     146fi=project_id_env+'pos'+csite+'_dy.cdf' 
    148147f=file_test(fi) 
    149148lat=fltarr(jpt)+!values.f_nan 
    … …  
    165164 
    166165;;MET PARAMETERS 
    167 fi=tropflux_id_env+'met'+csite+'_dy.cdf' 
     166fi=project_id_env+'met'+csite+'_dy.cdf' 
    168167f=file_test(fi) 
    169168wu=fltarr(jpt)+!values.f_nan 

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

    r47 r50  
    1616;       ] 
    1717; 
    18 ;       sst_erai [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/sst_2000_2009_erai_v50.txt"]; 
    19 ;       sst_tropflux [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/sst_2000_2009_trop_v50.txt"]; 
    20 ;       sst_oaflux [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/sst_2000_2009_oaflx_v50.txt"]; 
    21 ;       sst_tmi [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/sst_2000_2009_tmi_v50.txt"]; 
    22 ;       sst_ncep [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/sst_2000_2009_ncep_v50.txt"]; 
    23 ;       sst_ncep1 [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/sst_2000_2009_ncep1_v50.txt"]; 
    24 ;       sst_erai_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_sst_19890101_20091231_oafluxgrid.nc"]; 
    25 ;       sst_tropflux_2 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_sst_19890101_20091231_v20.nc"]; 
    26 ;       oaflux_basic [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/OAFlux_basic_variables_1985_2009.nc"]; 
    27 ;       sst_ncep2 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/sst_NCEP2_global_monthly_1989_2009.nc"]; 
    28 ;       sst_tmi_2 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/tmi_OI_oafluxgrid_30N30S.nc"]; 
    29 ;       sst_ncep1_2 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/sst_ncep1_19890101_20091231.nc"]; 
    30 ; 
    31 ;       figure [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/WHOTS_validation_net_2004_2007_all_v10.ps"]; 
     18;       sst_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/sst_2000_2009_erai_v50.txt"]; 
     19;       sst_tropflux [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/sst_2000_2009_trop_v50.txt"]; 
     20;       sst_oaflux [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/sst_2000_2009_oaflx_v50.txt"]; 
     21;       sst_tmi [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/sst_2000_2009_tmi_v50.txt"]; 
     22;       sst_ncep [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/sst_2000_2009_ncep_v50.txt"]; 
     23;       sst_ncep1 [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/sst_2000_2009_ncep1_v50.txt"]; 
     24;       sst_erai_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_sst_19890101_20091231_oafluxgrid.nc"]; 
     25;       sst_tropflux_2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/TropFlux_sst_19890101_20091231_v20.nc"]; 
     26;       oaflux_basic [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/OAFlux_basic_variables_1985_2009.nc"]; 
     27;       sst_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/sst_NCEP2_global_monthly_1989_2009.nc"]; 
     28;       sst_tmi_2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/tmi_OI_oafluxgrid_30N30S.nc"]; 
     29;       sst_ncep1_2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/sst_ncep1_19890101_20091231.nc"]; 
     30; 
     31;       figure [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/WHOTS_validation_net_2004_2007_all_v10.ps"]; 
    3232; 
    3333;       sst_validation_scatter_2000_2009_v50 [shape=box, 
    … …  
    4343; ======== 
    4444; 
    45 ; :ref:`tropflux_profile.sh` 
    46 ; :ref:`tropflux_init.pro` 
    47 ; :ref:`cm_tropflux.pro` 
     45; :ref:`project_profile.sh` 
     46; :ref:`project_init.pro` 
     47; :ref:`cm_project.pro` 
    4848; 
    4949; :func:`x_site_location` 
    … …  
    5858; :: 
    5959; 
    60 ;  IDL> @tropflux_init 
    6160;  IDL> date1=19890101L 
    6261;  IDL> date2=20091231L 
    … …  
    9695pro sst_validation_scatter_2000_2009_v50, date1, date2 
    9796@cm_general 
    98 @cm_tropflux 
     97@cm_project 
    9998reinitplt, /z,/invert 
    10099key_portrait = 1 
    … …  
    122121close,/all 
    123122 
    124 fi_sst_erai=tropflux_id_env+'sst_2000_2009_erai_v50.txt' 
     123fi_sst_erai=project_id_env+'sst_2000_2009_erai_v50.txt' 
    125124openw,1,fi_sst_erai 
    126 fi_sst_trop=tropflux_id_env+'sst_2000_2009_trop_v50.txt' 
     125fi_sst_trop=project_id_env+'sst_2000_2009_trop_v50.txt' 
    127126openw,2,fi_sst_trop 
    128 fi_sst_oaflx=tropflux_id_env+'sst_2000_2009_oaflx_v50.txt' 
     127fi_sst_oaflx=project_id_env+'sst_2000_2009_oaflx_v50.txt' 
    129128openw,3,fi_sst_oaflx 
    130 fi_sst_ncep=tropflux_id_env+'sst_2000_2009_ncep_v50.txt' 
     129fi_sst_ncep=project_id_env+'sst_2000_2009_ncep_v50.txt' 
    131130openw,4,fi_sst_ncep 
    132 fi_sst_tmi=tropflux_id_env+'sst_2000_2009_tmi_v50.txt' 
     131fi_sst_tmi=project_id_env+'sst_2000_2009_tmi_v50.txt' 
    133132openw,5,fi_sst_tmi 
    134 fi_sst_ncep1=tropflux_id_env+'sst_2000_2009_ncep1_v50.txt' 
     133fi_sst_ncep1=project_id_env+'sst_2000_2009_ncep1_v50.txt' 
    135134openw,6,fi_sst_ncep1 
    136135 
    … …  
    146145;; first reading the whole ERAI uncorrected and corrected data 
    147146 
    148 file=tropflux_id_env+'erai_sst_19890101_20091231_oafluxgrid.nc' 
     147file=project_id_env+'erai_sst_19890101_20091231_oafluxgrid.nc' 
    149148initncdf, file 
    150149unc=read_ncdf('sst',date1,date2,file=file,/nostr) & unc=unc-273.15 
    151150help, unc 
    152151 
    153 file=tropflux_id_env+'TropFlux_sst_19890101_20091231_v20.nc' 
     152file=project_id_env+'TropFlux_sst_19890101_20091231_v20.nc' 
    154153initncdf, file 
    155154cor=read_ncdf('sst',date1,date2,file=file,/nostr) & cor=cor-273.15 
    156155help, cor 
    157156 
    158 file=tropflux_id_env+'OAFlux_basic_variables_1985_2009.nc' 
     157file=project_id_env+'OAFlux_basic_variables_1985_2009.nc' 
    159158initncdf, file 
    160159oaf=read_ncdf("sst", date1, date2, file=file,/nostr) 
    161160help, oaf 
    162161 
    163 fi=tropflux_id_env+'sst_ncep2_oafluxgrid_19890101_20091231.nc' 
     162fi=project_id_env+'sst_ncep2_oafluxgrid_19890101_20091231.nc' 
    164163initncdf, fi 
    165164nce=read_ncdf("sst", date1-1, date2, file=fi,/nostr) 
    … …  
    167166help, nce 
    168167 
    169 fi=tropflux_id_env+'tmi_OI_oafluxgrid_30N30S.nc' 
     168fi=project_id_env+'tmi_OI_oafluxgrid_30N30S.nc' 
    170169initncdf, fi 
    171170sst_tmi=read_ncdf("sst", date1, date2, file=fi,/nostr) 
    172171help, sst_tmi 
    173172 
    174 file=tropflux_id_env+'sst_ncep1_19890101_20091231.nc' 
     173file=project_id_env+'sst_ncep1_19890101_20091231.nc' 
    175174initncdf, file 
    176175nce1=read_ncdf("sst", date1, date2, file=file,/nostr) 
    … …  
    250249close,/all 
    251250 
    252 fi_sst_erai=tropflux_id_env+'sst_2000_2009_v50_erai.txt' 
     251fi_sst_erai=project_id_env+'sst_2000_2009_v50_erai.txt' 
    253252res=read_ascii(fi_sst_erai,data_start=1) 
    254253ff=res.field1 
    … …  
    276275oplot, mean_tao, yfit, color=250, thick=2 
    277276 
    278 fi_sst_trop=tropflux_id_env+'sst_2000_2009_v50_trop.txt' 
     277fi_sst_trop=project_id_env+'sst_2000_2009_v50_trop.txt' 
    279278res=read_ascii(fi_sst_trop,data_start=1) 
    280279ff=res.field1 
    … …  
    303302oplot, mean_tao, yfit, color=250, thick=2 
    304303 
    305 fi_sst_oaflx=tropflux_id_env+'sst_2000_2009_v50_oaflx.txt' 
     304fi_sst_oaflx=project_id_env+'sst_2000_2009_v50_oaflx.txt' 
    306305res=read_ascii(fi_sst_oaflx,data_start=1) 
    307306ff=res.field1 
    … …  
    332331 
    333332 
    334 fi_sst_ncep=tropflux_id_env+'sst_2000_2009_v50_ncep.txt' 
     333fi_sst_ncep=project_id_env+'sst_2000_2009_v50_ncep.txt' 
    335334res=read_ascii(fi_sst_ncep,data_start=1) 
    336335ff=res.field1 
    … …  
    360359 
    361360 
    362 fi_sst_tmi=tropflux_id_env+'sst_2000_2009_v50_tmi.txt' 
     361fi_sst_tmi=project_id_env+'sst_2000_2009_v50_tmi.txt' 
    363362res=read_ascii(fi_sst_tmi,data_start=1) 
    364363ff=res.field1 
    … …  
    387386oplot, mean_tao, yfit, color=250, thick=2 
    388387 
    389 fi_sst_ncep1=tropflux_id_env+'sst_2000_2009_v50_ncep1.txt' 
     388fi_sst_ncep1=project_id_env+'sst_2000_2009_v50_ncep1.txt' 
    390389res=read_ascii(fi_sst_ncep1,data_start=1) 
    391390ff=res.field1 
    … …  
    416415closeps 
    417416 
    418 fig=tropflux_od_env+'sst_validation_scatter_2000_2009_v50.ps' 
     417fig=project_od_env+'sst_validation_scatter_2000_2009_v50.ps' 
    419418spawn, 'mv '+psdir+'idl.ps '+fig 
    420419spawn, 'gv '+fig 

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

    r46 r50  
    1212; ======== 
    1313; 
    14 ; :ref:`tropflux_profile.sh` 
     14; :ref:`project_profile.sh` 
    1515; 
    1616; :ref:`net_flux_validation_scatter_2000_2007.pro` 
    … …  
    3333; :: 
    3434; 
    35 ;  IDL> @tropflux_init 
    3635;  IDL> tao=++ 
    3736;  IDL> ++ 
    … …  
    5150    cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    5251 
    53 @cm_tropflux 
     52@cm_project 
    5453x=tao & y=var1 & z=var2 
    5554 

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

    r47 r50  
    1616;       ] 
    1717; 
    18 ;       ws_erai [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/ws_2000_2009_erai_v50.txt"]; 
    19 ;       ws_tropflux [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/ws_2000_2009_trop_v50.txt"]; 
    20 ;       ws_oaflux [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/ws_2000_2009_oaflx_v50.txt"]; 
    21 ;       ws_ncep [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/ws_2000_2009_ncep_v50.txt"]; 
    22 ;       ws_ncep1 [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/ws_2000_2009_ncep1_v50.txt"]; 
    23 ;       ws_erai_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_ws_19890101_20091231_oafluxgrid.nc"]; 
    24 ;       ws_tropflux_2 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_ws_19890101_20091231_v20.nc"]; 
    25 ;       oaflux_basic [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/OAFlux_basic_variables_1985_2009.nc"]; 
    26 ;       uwind_ncep2 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/uwind_ncep2_oafluxgrid_19890101_20091231.nc"]; 
    27 ;       vwind_ncep2 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/bwind_ncep2_oafluxgrid_19890101_20091231.nc"]; 
    28 ;       ws_tmi [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/zonal_wind_speed_oafluxgrid_30N30S.nc"]; 
    29 ; 
    30 ;       figure [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/ws_validation_scatter_2000_2009_v50.ps"]; 
     18;       ws_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/ws_2000_2009_erai_v50.txt"]; 
     19;       ws_tropflux [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/ws_2000_2009_trop_v50.txt"]; 
     20;       ws_oaflux [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/ws_2000_2009_oaflx_v50.txt"]; 
     21;       ws_ncep [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/ws_2000_2009_ncep_v50.txt"]; 
     22;       ws_ncep1 [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/ws_2000_2009_ncep1_v50.txt"]; 
     23;       ws_erai_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_ws_19890101_20091231_oafluxgrid.nc"]; 
     24;       ws_tropflux_2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/TropFlux_ws_19890101_20091231_v20.nc"]; 
     25;       oaflux_basic [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/OAFlux_basic_variables_1985_2009.nc"]; 
     26;       uwind_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/uwind_ncep2_oafluxgrid_19890101_20091231.nc"]; 
     27;       vwind_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/bwind_ncep2_oafluxgrid_19890101_20091231.nc"]; 
     28;       ws_tmi [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/zonal_wind_speed_oafluxgrid_30N30S.nc"]; 
     29; 
     30;       figure [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/ws_validation_scatter_2000_2009_v50.ps"]; 
    3131; 
    3232;       ws_validation_scatter_2000_2009_v50 [shape=box, 
    … …  
    4242; ======== 
    4343; 
    44 ; :ref:`tropflux_profile.sh` 
    45 ; :ref:`tropflux_init.pro` 
    46 ; :ref:`cm_tropflux.pro` 
     44; :ref:`project_profile.sh` 
     45; :ref:`project_init.pro` 
     46; :ref:`cm_project.pro` 
    4747; 
    4848; :func:`x_site_location` 
    … …  
    5757; :: 
    5858; 
    59 ;  IDL> @tropflux_init 
    6059;  IDL> date1=19890101L 
    6160;  IDL> date2=20091231L 
    … …  
    9695pro ws_validation_scatter_2000_2009_v50, date1, date2 
    9796@cm_general 
    98 @cm_tropflux 
     97@cm_project 
    9998reinitplt, /z,/invert 
    10099key_portrait = 1 
    … …  
    123122close,/all 
    124123 
    125 fi_ws_erai=tropflux_id_env+'ws_2000_2009_erai_v50.txt' 
     124fi_ws_erai=project_id_env+'ws_2000_2009_erai_v50.txt' 
    126125openw,1,fi_ws_erai 
    127 fi_ws_trop=tropflux_id_env+'ws_2000_2009_trop_v50.txt' 
     126fi_ws_trop=project_id_env+'ws_2000_2009_trop_v50.txt' 
    128127openw,2,fi_ws_trop 
    129 fi_ws_oaflx=tropflux_id_env+'ws_2000_2009_oaflx_v50.txt' 
     128fi_ws_oaflx=project_id_env+'ws_2000_2009_oaflx_v50.txt' 
    130129openw,3,fi_ws_oaflx 
    131 fi_ws_ncep=tropflux_id_env+'ws_2000_2009_ncep_v50.txt' 
     130fi_ws_ncep=project_id_env+'ws_2000_2009_ncep_v50.txt' 
    132131openw,4,fi_ws_ncep 
    133 fi_ws_tmi=tropflux_id_env+'ws_2000_2009_tmi_v50.txt' 
     132fi_ws_tmi=project_id_env+'ws_2000_2009_tmi_v50.txt' 
    134133openw,5,fi_ws_tmi 
    135 fi_ws_ncep1=tropflux_id_env+'ws_2000_2009_ncep1_v50.txt' 
     134fi_ws_ncep1=project_id_env+'ws_2000_2009_ncep1_v50.txt' 
    136135openw,6,fi_ws_ncep1 
    137136 
    … …  
    146145;; first reading the whole ERAI uncorrected and corrected data 
    147146 
    148 file=tropflux_id_env+'erai_ws_19890101_20091231_oafluxgrid.nc' 
     147file=project_id_env+'erai_ws_19890101_20091231_oafluxgrid.nc' 
    149148initncdf, file 
    150149u=read_ncdf('u10',date1,date2,file=file,/nostr) 
    … …  
    153152help, unc 
    154153 
    155 file=tropflux_id_env+'TropFlux_ws_19890101_20091231_v20.nc' 
     154file=project_id_env+'TropFlux_ws_19890101_20091231_v20.nc' 
    156155initncdf, file 
    157156cor=read_ncdf('ws',date1,date2,file=file,/nostr) 
    158157help, cor 
    159158 
    160 file=tropflux_id_env+'OAFlux_basic_variables_1985_2009.nc' 
     159file=project_id_env+'OAFlux_basic_variables_1985_2009.nc' 
    161160initncdf, file 
    162161oaf=read_ncdf("wind", date1, date2, file=file,/nostr) 
    163162help, oaf 
    164163 
    165 fi=tropflux_id_env+'uwind_ncep2_oafluxgrid_19890101_20091231.nc' 
     164fi=project_id_env+'uwind_ncep2_oafluxgrid_19890101_20091231.nc' 
    166165initncdf, fi 
    167166u=read_ncdf("u", date1-1, date2, file=fi,/nostr) 
    168 fi=tropflux_id_env+'vwind_ncep2_oafluxgrid_19890101_20091231.nc' 
     167fi=project_id_env+'vwind_ncep2_oafluxgrid_19890101_20091231.nc' 
    169168initncdf, fi 
    170169v=read_ncdf("v", date1-1, date2, file=fi,/nostr) 
    … …  
    172171help, nce 
    173172 
    174 fi=tropflux_id_env+'zonal_wind_speed_oafluxgrid_30N30S.nc' 
     173fi=project_id_env+'zonal_wind_speed_oafluxgrid_30N30S.nc' 
    175174initncdf, fi 
    176175u=read_ncdf("u", date1, date2, file=fi,/nostr) 
    177 fi=tropflux_id_env+'meridional_wind_speed_oafluxgrid_30N30S.nc' 
     176fi=project_id_env+'meridional_wind_speed_oafluxgrid_30N30S.nc' 
    178177initncdf, fi 
    179178v=read_ncdf("v", date1, date2, file=fi,/nostr) 
    … …  
    181180help, ws_tmi 
    182181 
    183 file=tropflux_id_env+'wind_ncep1_19890101_20091231.nc' 
     182file=project_id_env+'wind_ncep1_19890101_20091231.nc' 
    184183initncdf, file 
    185184u=read_ncdf("u", date1, date2, file=file,/nostr) 
    … …  
    254253endfor 
    255254close,/all 
    256 fi_ws_erai=tropflux_id_env+'ws_2000_2009_erai_v50.txt' 
     255fi_ws_erai=project_id_env+'ws_2000_2009_erai_v50.txt' 
    257256res=read_ascii(fi_ws_erai,data_start=1) 
    258257ff=res.field1 
    … …  
    280279oplot, mean_tao, yfit, color=250, thick=2 
    281280 
    282 fi_ws_trop=tropflux_id_env+'ws_2000_2009_trop_v50.txt' 
     281fi_ws_trop=project_id_env+'ws_2000_2009_trop_v50.txt' 
    283282res=read_ascii(fi_ws_trop,data_start=1) 
    284283ff=res.field1 
    … …  
    307306oplot, mean_tao, yfit, color=250, thick=2 
    308307 
    309 fi_ws_oaflx=tropflux_id_env+'ws_2000_2009_oaflx_v50.txt' 
     308fi_ws_oaflx=project_id_env+'ws_2000_2009_oaflx_v50.txt' 
    310309res=read_ascii(fi_ws_oaflx,data_start=1) 
    311310ff=res.field1 
    … …  
    336335 
    337336 
    338 fi_ws_ncep=tropflux_id_env+'ws_2000_2009_ncep_v50.txt' 
     337fi_ws_ncep=project_id_env+'ws_2000_2009_ncep_v50.txt' 
    339338res=read_ascii(fi_ws_ncep,data_start=1) 
    340339ff=res.field1 
    … …  
    363362oplot, mean_tao, yfit, color=250, thick=2 
    364363 
    365 fi_ws_tmi=tropflux_id_env+'ws_2000_2009_tmi_v50.txt' 
     364fi_ws_tmi=project_id_env+'ws_2000_2009_tmi_v50.txt' 
    366365res=read_ascii(fi_ws_tmi,data_start=1) 
    367366ff=res.field1 
    … …  
    391390 
    392391 
    393 fi_ws_ncep1=tropflux_id_env+'ws_2000_2009_ncep1_v50.txt' 
     392fi_ws_ncep1=project_id_env+'ws_2000_2009_ncep1_v50.txt' 
    394393res=read_ascii(fi_ws_ncep1,data_start=1) 
    395394ff=res.field1 
    … …  
    420419closeps 
    421420 
    422 fig=tropflux_od_env+'ws_validation_scatter_2000_2009_v50.ps' 
     421fig=project_od_env+'ws_validation_scatter_2000_2009_v50.ps' 
    423422spawn, 'mv '+psdir+'idl.ps '+fig 
    424423spawn, 'gv '+fig 
Note: See TracChangeset for help on using the changeset viewer.