Changeset 46 for trunk/src/paper01/fig3


Ignore:
Timestamp:
04/20/11 16:45:53 (13 years ago)
Author:
pinsard
Message:

going on 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

    r44 r46  
    1919;       air_erai [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/air_2000_2009_erai_v50.txt"]; 
    2020; 
    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"]; 
     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"]; 
    2323;       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"]; 
     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"]; 
    2626;       air_ncep2_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/air_2m_ncep2_oafluxgrid_19890101_20091231.nc"]; 
    2727; 
    … …  
    3232; 
    3333;       tropflux_t2m [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_t2m_19890101_20091231_v50.nc"]; 
    34 ;      
     34; 
    3535;       oaflux_basic [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/OAFlux_basic_variables_1985_2009.nc"]; 
    3636; 
    … …  
    4141;       label="${TROPFLUX}/src/paper01/fig3/air_validation_scatter_2000_2009_v50.pro"]; 
    4242; 
    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} 
     43;       {air_erai air_ncep erai_t2m tropflux_t2m oaflux_basic air_ncep2_oafluxgrid t2m_ncep1} -> {air_validation_scatter_2000_2009_v50} -> { air_trop air_oaflux air_ncep1 air_ncep2 figure} 
    4444; 
    4545;    } 
    … …  
    4949; 
    5050; :ref:`tropflux_profile.sh` 
     51; :ref:`tropflux_init.pro` 
     52; :ref:`cm_tropflux.pro` 
     53; 
     54; :func:`x_site_location` 
     55; :func:`y_site_location` 
    5156; 
    5257; :ref:`read_variables_v2.pro` 
     58; :ref:`statistics_3var_v1.pro` 
    5359; 
    5460; EXAMPLES 
    … …  
    102108 
    103109openps, FILENAME = 'idl.ps' 
    104 ;; Before running this program, you have to compile the following subrutines 
    105 ;;   1.  .r read_era_total 
    106 ;;   2.  .r read_variables_v2 
    107 ;;   3.  .r read_era_box 
    108110 
    109111;; Give the location of mooring for validation of basic meteorological variables 
    … …  
    152154help, unc 
    153155 
    154 ;file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_t2m_19890101_20091231_v3.nc' 
    155156file=tropflux_id_env+'TropFlux_t2m_19890101_20091231_v50.nc' 
    156157initncdf, file 
    … …  
    176177 
    177178nn=n_elements(sitelist) 
    178 date1=date1 
    179 date2=date2 
    180179for n=0, nn-1 do begin 
    181180 

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

    r44 r46  
    1919;       q2m_erai [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/q2m_2000_2009_erai_v50.txt"]; 
    2020; 
    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"]; 
     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"]; 
    2525; 
    2626;       erai_q2m_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_t2m_19890101_20091231_oafluxgrid.nc"]; 
    … …  
    4040;       figure [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/q2m_validation_scatter_2000_2009_v50.ps"]; 
    4141; 
    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} 
     42;       {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} 
    4343; 
    4444;     } 
    … …  
    4848; 
    4949; :ref:`tropflux_profile.sh` 
     50; :ref:`tropflux_init.pro` 
     51; :ref:`cm_tropflux.pro` 
     52; 
     53; :func:`x_site_location` 
     54; :func:`y_site_location` 
     55; 
     56; :ref:`read_variables_v2.pro` 
     57; :ref:`statistics_3var_v1.pro` 
    5058; 
    5159; EXAMPLES 
    … …  
    99107 
    100108openps, FILENAME = 'idl.ps' 
    101 ;; Before running this program, you have to compile the following subrutines 
    102 ;;   1.  .r read_era_total 
    103 ;;   2.  .r read_variables_v2 
    104 ;;   3.  .r read_era_box 
    105109 
    106110;; Give the location of mooring for validation of basic meteorological variables 
    … …  
    165169help, nce 
    166170 
    167 file='/Volumes/Iomega_HDD/flux_reconstruction/ncep/q2m_ncep1_19890101_20091231.nc' 
     171file=tropflux_id_env+''q2m_ncep1_19890101_20091231.nc' 
    168172initncdf, file 
    169173nce1=read_ncdf("q2m", date1, date2, file=file,/nostr) 
    … …  
    172176 
    173177nn=n_elements(sitelist) 
    174 date1=date1 
    175 date2=date2 
    176178for n=0, nn-1 do begin 
    177179 
    … …  
    241243close,/all 
    242244 
    243 fi_q2m_erai='/Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_erai_v50.txt' 
     245fi_q2m_erai=tropflux_id_env+'q2m_2000_2009_erai_v50.txt' 
    244246res=read_ascii(fi_q2m_erai,data_start=1) 
    245247ff=res.field1 
    … …  
    267269oplot, mean_tao, yfit, color=250, thick=2 
    268270 
    269 fi_q2m_trop='/Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_trop_v50.txt' 
     271fi_q2m_trop=tropflux_id_env+'q2m_2000_2009_trop_v50.txt' 
    270272res=read_ascii(fi_q2m_trop,data_start=1) 
    271273ff=res.field1 
    … …  
    294296oplot, mean_tao, yfit, color=250, thick=2 
    295297 
    296 fi_q2m_oaflx='/Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_oaflx_v50.txt' 
     298fi_q2m_oaflx=tropflux_id_env+'q2m_2000_2009_oaflx_v50.txt' 
    297299res=read_ascii(fi_q2m_oaflx,data_start=1) 
    298300ff=res.field1 
    … …  
    322324 
    323325 
    324 fi_q2m_ncep='/Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_ncep_v50.txt' 
     326fi_q2m_ncep=tropflux_id_env+'q2m_2000_2009_ncep_v50.txt' 
    325327res=read_ascii(fi_q2m_ncep,data_start=1) 
    326328ff=res.field1 
    … …  
    349351oplot, mean_tao, yfit, color=250, thick=2 
    350352 
    351 fi_q2m_ncep1='/Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_ncep1_v50.txt' 
     353fi_q2m_ncep1=tropflux_id_env+'q2m_2000_2009_ncep1_v50.txt' 
    352354res=read_ascii(fi_q2m_ncep1,data_start=1) 
    353355ff=res.field1 

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

    r44 r46  
    4848; :ref:`sst_validation_scatter_2000_2009_v50.pro` 
    4949; :ref:`ws_validation_scatter_2000_2009_v50.pro` 
     50; :ref:`air_validation_scatter_2000_2009_basin.pro` 
     51; :ref:`q2m_validation_scatter_2000_2009_basin.pro` 
     52; :ref:`sst_validation_scatter_2000_2009_basin.pro` 
     53; :ref:`ws_validation_scatter_2000_2009_basin.pro` 
     54; :ref:`lhf_validation_scatter_2000_2009.pro` 
     55; :ref:`shf_validation_scatter_2000_2009.pro` 
    5056; 
    5157; EXAMPLES 
    … …  
    6167; TODO 
    6268; ==== 
    63 ;  
     69; 
    6470; externalize function 
    6571; 

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

    r44 r46  
    99; =========== 
    1010; 
     11; .. graphviz:: 
     12; 
     13;    digraph sst_validation_scatter_2000_2009_v50 { 
     14;       graph [ 
     15;       rankdir="LR", 
     16;       ] 
     17; 
     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"]; 
     32; 
     33;       sst_validation_scatter_2000_2009_v50 [shape=box, 
     34;       fontname=Courier, 
     35;       color=blue, 
     36;       URL="http://forge.ipsl.jussieu.fr/tropflux/broswrer/trunk/src/paper01/fig3/sst_validation_scatter_2000_2009_v50.pro", 
     37;       label="${TROPFLUX}/src/paper01/fig3/sst_validation_scatter_2000_2009_v50.pro"]; 
     38; 
     39;       {sst_erai sst_erai_oafluxgrid sst_tropflux_2 oaflux_basic sst_ncep2 sst_tmi_2 sst_ncep1_2} -> {sst_validation_scatter_2000_2009_v50} -> {sst_tropflux sst_oaflux sst_tmi sst_ncep sst_ncep figure} 
     40;    } 
     41; 
    1142; SEE ALSO 
    1243; ======== 
    1344; 
    1445; :ref:`tropflux_profile.sh` 
     46; :ref:`tropflux_init.pro` 
     47; :ref:`cm_tropflux.pro` 
     48; 
     49; :func:`x_site_location` 
     50; :func:`y_site_location` 
     51; 
     52; :ref:`read_variables_v2.pro` 
     53; :ref:`statistics_3var_v1.pro` 
    1554; 
    1655; EXAMPLES 
    … …  
    2059; 
    2160;  IDL> @tropflux_init 
    22 ;  IDL> date1=++ 
    23 ;  IDL> date2=++ 
     61;  IDL> date1=19890101L 
     62;  IDL> date2=20091231L 
    2463;  IDL> sst_validation_scatter_2000_2009_v50, date1, date2 
     64; 
     65; TODO 
     66; ==== 
     67; 
     68; ++ plus de licences sur cratos 
     69; 
     70; ++ mooring data in graphviz 
     71; 
     72; coding rules 
     73; 
     74; get rid of spwan 
     75; 
     76; complete description 
     77; 
     78; handle IO error 
    2579; 
    2680; EVOLUTIONS 
    … …  
    2983; $Id$ 
    3084; 
     85; - fplod 20110420T131533Z aedon.locean-ipsl.upmc.fr (Darwin) 
     86; 
     87;   * remove hard coding path 
     88;   * add graphviz 
     89;   * externalize functions 
     90; 
    3191; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
    3292; 
    … …  
    3595;- 
    3696pro sst_validation_scatter_2000_2009_v50, date1, date2 
    37 @common 
    38  
     97@cm_general 
     98@cm_tropflux 
    3999reinitplt, /z,/invert 
    40100key_portrait = 1 
    41101 
    42102openps, FILENAME = 'idl.ps' 
    43 ;; Before running this program, you have to compile the following subrutines 
    44 ;;   1.  .r read_era_total 
    45 ;;   2.  .r read_variables_v2 
    46 ;;   3.  .r read_era_box 
    47103 
    48104;; Give the location of mooring for validation of basic meteorological variables 
    … …  
    66122close,/all 
    67123 
    68 fi_sst_erai='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_erai_v50.txt' 
     124fi_sst_erai=tropflux_id_env+'sst_2000_2009_erai_v50.txt' 
    69125openw,1,fi_sst_erai 
    70 fi_sst_trop='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_trop_v50.txt' 
     126fi_sst_trop=tropflux_id_env+'sst_2000_2009_trop_v50.txt' 
    71127openw,2,fi_sst_trop 
    72 fi_sst_oaflx='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_oaflx_v50.txt' 
     128fi_sst_oaflx=tropflux_id_env+'sst_2000_2009_oaflx_v50.txt' 
    73129openw,3,fi_sst_oaflx 
    74 fi_sst_ncep='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_ncep_v50.txt' 
     130fi_sst_ncep=tropflux_id_env+'sst_2000_2009_ncep_v50.txt' 
    75131openw,4,fi_sst_ncep 
    76 fi_sst_tmi='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_tmi_v50.txt' 
     132fi_sst_tmi=tropflux_id_env+'sst_2000_2009_tmi_v50.txt' 
    77133openw,5,fi_sst_tmi 
    78 fi_sst_ncep1='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_ncep1_v50.txt' 
     134fi_sst_ncep1=tropflux_id_env+'sst_2000_2009_ncep1_v50.txt' 
    79135openw,6,fi_sst_ncep1 
    80136 
    … …  
    90146;; first reading the whole ERAI uncorrected and corrected data 
    91147 
    92 file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_sst_19890101_20091231_oafluxgrid.nc' 
     148file=tropflux_id_env+'erai_sst_19890101_20091231_oafluxgrid.nc' 
    93149initncdf, file 
    94150unc=read_ncdf('sst',date1,date2,file=file,/nostr) & unc=unc-273.15 
    95151help, unc 
    96152 
    97 file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_sst_19890101_20091231_v20.nc' 
     153file=tropflux_id_env+'TropFlux_sst_19890101_20091231_v20.nc' 
    98154initncdf, file 
    99155cor=read_ncdf('sst',date1,date2,file=file,/nostr) & cor=cor-273.15 
    100156help, cor 
    101157 
    102 file='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/OAFlux_basic_variables_1985_2009.nc' 
     158file=tropflux_id_env+'OAFlux_basic_variables_1985_2009.nc' 
    103159initncdf, file 
    104160oaf=read_ncdf("sst", date1, date2, file=file,/nostr) 
    105161help, oaf 
    106162 
    107 fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/sst_ncep2_oafluxgrid_19890101_20091231.nc' 
     163fi=tropflux_id_env+'sst_ncep2_oafluxgrid_19890101_20091231.nc' 
    108164initncdf, fi 
    109165nce=read_ncdf("sst", date1-1, date2, file=fi,/nostr) 
    … …  
    111167help, nce 
    112168 
    113 fi='/Volumes/Iomega_HDD/work/flux_reconstruction/gridded_data/tmi_OI_oafluxgrid_30N30S.nc' 
     169fi=tropflux_id_env+'tmi_OI_oafluxgrid_30N30S.nc' 
    114170initncdf, fi 
    115171sst_tmi=read_ncdf("sst", date1, date2, file=fi,/nostr) 
    116172help, sst_tmi 
    117173 
    118 file='/Volumes/Iomega_HDD/flux_reconstruction/ncep/sst_ncep1_19890101_20091231.nc' 
     174file=tropflux_id_env+'sst_ncep1_19890101_20091231.nc' 
    119175initncdf, file 
    120176nce1=read_ncdf("sst", date1, date2, file=file,/nostr) 
    … …  
    123179 
    124180nn=n_elements(sitelist) 
    125 date1=date1 
    126 date2=date2 
    127181for n=0, nn-1 do begin 
    128182 
    … …  
    196250close,/all 
    197251 
    198 fi_sst_erai='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_v50_erai.txt' 
     252fi_sst_erai=tropflux_id_env+'sst_2000_2009_v50_erai.txt' 
    199253res=read_ascii(fi_sst_erai,data_start=1) 
    200254ff=res.field1 
    … …  
    222276oplot, mean_tao, yfit, color=250, thick=2 
    223277 
    224 fi_sst_trop='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_v50_trop.txt' 
     278fi_sst_trop=tropflux_id_env+'sst_2000_2009_v50_trop.txt' 
    225279res=read_ascii(fi_sst_trop,data_start=1) 
    226280ff=res.field1 
    … …  
    249303oplot, mean_tao, yfit, color=250, thick=2 
    250304 
    251 fi_sst_oaflx='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_v50_oaflx.txt' 
     305fi_sst_oaflx=tropflux_id_env+'sst_2000_2009_v50_oaflx.txt' 
    252306res=read_ascii(fi_sst_oaflx,data_start=1) 
    253307ff=res.field1 
    … …  
    278332 
    279333 
    280 fi_sst_ncep='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_v50_ncep.txt' 
     334fi_sst_ncep=tropflux_id_env+'sst_2000_2009_v50_ncep.txt' 
    281335res=read_ascii(fi_sst_ncep,data_start=1) 
    282336ff=res.field1 
    … …  
    306360 
    307361 
    308 fi_sst_tmi='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_v50_tmi.txt' 
     362fi_sst_tmi=tropflux_id_env+'sst_2000_2009_v50_tmi.txt' 
    309363res=read_ascii(fi_sst_tmi,data_start=1) 
    310364ff=res.field1 
    … …  
    333387oplot, mean_tao, yfit, color=250, thick=2 
    334388 
    335 fi_sst_ncep1='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_v50_ncep1.txt' 
     389fi_sst_ncep1=tropflux_id_env+'sst_2000_2009_v50_ncep1.txt' 
    336390res=read_ascii(fi_sst_ncep1,data_start=1) 
    337391ff=res.field1 
    … …  
    365419spawn, 'mv '+psdir+'idl.ps '+fig 
    366420spawn, 'gv '+fig 
    367 return 
    368421end 
    369 function x_site_location, site 
    370     n1=strpos(site, 's') 
    371 if (n1 gt -1) then begin 
    372     ns=-1. 
    373     x=strmid(site, 0, n1) 
    374     x=float(x)*ns 
    375 endif else begin 
    376     n1=strpos(site, 'n') 
    377     x=strmid(site, 0, n1) 
    378     ny=1. 
    379     x=float(x)*ny 
    380 endelse 
    381 return, float(x) 
    382 end 
    383 function y_site_location, site 
    384     n1=strpos(site, 'e') 
    385 if (n1 gt -1) then begin 
    386     n=strpos(site, 's') 
    387     if (n gt -1) then begin 
    388         y=strmid(site, n+1, n1-n-1) 
    389     endif else begin 
    390         n=strpos(site, 'n') 
    391         y=strmid(site, n+1, n1-n-1) 
    392     endelse 
    393  
    394 endif else begin 
    395     n1=strpos(site, 'w') 
    396     n=strpos(site, 's') 
    397     if (n gt -1) then begin 
    398         y=strmid(site, n+1, n1-n-1) 
    399         y=180+(180-float(y)) 
    400     endif else begin 
    401         n=strpos(site, 'n') 
    402         y=strmid(site, n+1, n1-n-1) 
    403         y=180+(180-float(y)) 
    404    endelse 
    405 endelse 
    406 return,float(y) 
    407 end 

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

    r44 r46  
    1313; 
    1414; :ref:`tropflux_profile.sh` 
     15; 
     16; :ref:`net_flux_validation_scatter_2000_2007.pro` 
     17; :ref:`air_validation_scatter_2000_2009_v50.pro` 
     18; :ref:`q2m_validation_scatter_2000_2009_v50.pro` 
     19; :ref:`sst_validation_scatter_2000_2009_v50.pro` 
     20; :ref:`ws_validation_scatter_2000_2009_v50.pro` 
     21; :ref:`air_validation_scatter_2000_2009_basin.pro` 
     22; :ref:`q2m_validation_scatter_2000_2009_basin.pro` 
     23; :ref:`sst_validation_scatter_2000_2009_basin.pro` 
     24; :ref:`ws_validation_scatter_2000_2009_basin.pro` 
     25; :ref:`lhf_validation_scatter_2000_2009.pro` 
     26; :ref:`shf_validation_scatter_2000_2009.pro` 
     27; :ref:`lwr_validation_scatter_2000_2007.pro` 
     28; :ref:`swr_validation_scatter_2000_2007.pro` 
    1529; 
    1630; EXAMPLES 
    … …  
    3448; 
    3549;- 
    36  
    37  
    3850pro statistics_3var_v1, tao,var1,var2, $ 
    3951    cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    4052 
    41  
    42 @common 
     53@cm_tropflux 
    4354x=tao & y=var1 & z=var2 
    4455 

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

    r44 r46  
    99; =========== 
    1010; 
     11; .. graphviz:: 
     12; 
     13;    digraph ws_validation_scatter_2000_2009_v50 { 
     14;       graph [ 
     15;       rankdir="LR", 
     16;       ] 
     17; 
     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"]; 
     31; 
     32;       ws_validation_scatter_2000_2009_v50 [shape=box, 
     33;       fontname=Courier, 
     34;       color=blue, 
     35;       URL="http://forge.ipsl.jussieu.fr/tropflux/broswrer/trunk/src/paper01/fig3/ws_validation_scatter_2000_2009_v50.pro", 
     36;       label="${TROPFLUX}/src/paper01/fig3/ws_validation_scatter_2000_2009_v50.pro"]; 
     37; 
     38;       {ws_erai ws_erai_oafluxgrid ws_tropflux_2 oaflux_basic ws_ncep2 ws_tmi_2 ws_ncep1_2} -> {ws_validation_scatter_2000_2009_v50} -> {ws_tropflux ws_oaflux ws_tmi ws_ncep ws_ncep1 figure} 
     39;    } 
     40; 
    1141; SEE ALSO 
    1242; ======== 
    1343; 
    1444; :ref:`tropflux_profile.sh` 
     45; :ref:`tropflux_init.pro` 
     46; :ref:`cm_tropflux.pro` 
     47; 
     48; :func:`x_site_location` 
     49; :func:`y_site_location` 
     50; 
     51; :ref:`read_variables_v2.pro` 
     52; :ref:`statistics_3var_v1.pro` 
    1553; 
    1654; EXAMPLES 
    … …  
    2058; 
    2159;  IDL> @tropflux_init 
    22 ;  IDL> date1=++ 
    23 ;  IDL> date2=++ 
     60;  IDL> date1=19890101L 
     61;  IDL> date2=20091231L 
    2462;  IDL> ws_validation_scatter_2000_2009_v50, date1, date2 
     63; 
     64; TODO 
     65; ==== 
     66; 
     67; ++ plus de licences sur cratos 
     68; 
     69; ++ mooring data in graphviz 
     70; 
     71; coding rules 
     72; 
     73; get rid of spwan 
     74; 
     75; complete description 
     76; 
     77; handle IO error 
    2578; 
    2679; 
    … …  
    3083; $Id$ 
    3184; 
     85; - fplod 20110420T132401Z aedon.locean-ipsl.upmc.fr (Darwin) 
     86; 
     87;   * remove hard coding path 
     88;   * add graphviz 
     89;   * externalize functions 
     90; 
    3291; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
    3392; 
    … …  
    3594; 
    3695;- 
    37 pro ws_validation_scatter_2000_2009_v50,date1,date2 
    38 @common 
    39  
     96pro ws_validation_scatter_2000_2009_v50, date1, date2 
     97@cm_general 
     98@cm_tropflux 
    4099reinitplt, /z,/invert 
    41100key_portrait = 1 
    42101 
    43102openps, FILENAME = 'idl.ps' 
    44 ;; Before running this program, you have to compile the following subrutines 
    45 ;;   1.  .r read_era_total 
    46 ;;   2.  .r read_variables_v2 
    47 ;;   3.  .r read_era_box 
    48103 
    49104;; Give the location of mooring for validation of basic meteorological variables 
    … …  
    68123close,/all 
    69124 
    70 fi_ws_erai='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_erai_v50.txt' 
     125fi_ws_erai=tropflux_id_env+'ws_2000_2009_erai_v50.txt' 
    71126openw,1,fi_ws_erai 
    72 fi_ws_trop='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_trop_v50.txt' 
     127fi_ws_trop=tropflux_id_env+'ws_2000_2009_trop_v50.txt' 
    73128openw,2,fi_ws_trop 
    74 fi_ws_oaflx='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_oaflx_v50.txt' 
     129fi_ws_oaflx=tropflux_id_env+'ws_2000_2009_oaflx_v50.txt' 
    75130openw,3,fi_ws_oaflx 
    76 fi_ws_ncep='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_ncep_v50.txt' 
     131fi_ws_ncep=tropflux_id_env+'ws_2000_2009_ncep_v50.txt' 
    77132openw,4,fi_ws_ncep 
    78 fi_ws_tmi='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_tmi_v50.txt' 
     133fi_ws_tmi=tropflux_id_env+'ws_2000_2009_tmi_v50.txt' 
    79134openw,5,fi_ws_tmi 
    80 fi_ws_ncep1='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_ncep1_v50.txt' 
     135fi_ws_ncep1=tropflux_id_env+'ws_2000_2009_ncep1_v50.txt' 
    81136openw,6,fi_ws_ncep1 
    82137 
    … …  
    91146;; first reading the whole ERAI uncorrected and corrected data 
    92147 
    93 file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_ws_19890101_20091231_oafluxgrid.nc' 
     148file=tropflux_id_env+'erai_ws_19890101_20091231_oafluxgrid.nc' 
    94149initncdf, file 
    95150u=read_ncdf('u10',date1,date2,file=file,/nostr) 
    … …  
    98153help, unc 
    99154 
    100 file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_ws_19890101_20091231_v20.nc' 
     155file=tropflux_id_env+'TropFlux_ws_19890101_20091231_v20.nc' 
    101156initncdf, file 
    102157cor=read_ncdf('ws',date1,date2,file=file,/nostr) 
    103158help, cor 
    104159 
    105 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' 
    106161initncdf, file 
    107162oaf=read_ncdf("wind", date1, date2, file=file,/nostr) 
    108163help, oaf 
    109164 
    110 fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/uwind_ncep2_oafluxgrid_19890101_20091231.nc' 
     165fi=tropflux_id_env+'uwind_ncep2_oafluxgrid_19890101_20091231.nc' 
    111166initncdf, fi 
    112167u=read_ncdf("u", date1-1, date2, file=fi,/nostr) 
    113 fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/vwind_ncep2_oafluxgrid_19890101_20091231.nc' 
     168fi=tropflux_id_env+'vwind_ncep2_oafluxgrid_19890101_20091231.nc' 
    114169initncdf, fi 
    115170v=read_ncdf("v", date1-1, date2, file=fi,/nostr) 
    … …  
    126181help, ws_tmi 
    127182 
    128 file='/Volumes/Iomega_HDD/flux_reconstruction/ncep/wind_ncep1_19890101_20091231.nc' 
     183file=tropflux_id_env+'wind_ncep1_19890101_20091231.nc' 
    129184initncdf, file 
    130185u=read_ncdf("u", date1, date2, file=file,/nostr) 
    … …  
    134189 
    135190nn=n_elements(sitelist) 
    136 date1=date1 
    137 date2=date2 
    138191for n=0, nn-1 do begin 
    139192 
    … …  
    201254endfor 
    202255close,/all 
    203 fi_ws_erai='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_erai_v50.txt' 
     256fi_ws_erai=tropflux_id_env+'ws_2000_2009_erai_v50.txt' 
    204257res=read_ascii(fi_ws_erai,data_start=1) 
    205258ff=res.field1 
    … …  
    227280oplot, mean_tao, yfit, color=250, thick=2 
    228281 
    229 fi_ws_trop='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_trop_v50.txt' 
     282fi_ws_trop=tropflux_id_env+'ws_2000_2009_trop_v50.txt' 
    230283res=read_ascii(fi_ws_trop,data_start=1) 
    231284ff=res.field1 
    … …  
    254307oplot, mean_tao, yfit, color=250, thick=2 
    255308 
    256 fi_ws_oaflx='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_oaflx_v50.txt' 
     309fi_ws_oaflx=tropflux_id_env+'ws_2000_2009_oaflx_v50.txt' 
    257310res=read_ascii(fi_ws_oaflx,data_start=1) 
    258311ff=res.field1 
    … …  
    283336 
    284337 
    285 fi_ws_ncep='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_ncep_v50.txt' 
     338fi_ws_ncep=tropflux_id_env+'ws_2000_2009_ncep_v50.txt' 
    286339res=read_ascii(fi_ws_ncep,data_start=1) 
    287340ff=res.field1 
    … …  
    310363oplot, mean_tao, yfit, color=250, thick=2 
    311364 
    312 fi_ws_tmi='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_tmi_v50.txt' 
     365fi_ws_tmi=tropflux_id_env+'ws_2000_2009_tmi_v50.txt' 
    313366res=read_ascii(fi_ws_tmi,data_start=1) 
    314367ff=res.field1 
    … …  
    338391 
    339392 
    340 fi_ws_ncep1='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_ncep1_v50.txt' 
     393fi_ws_ncep1=tropflux_id_env+'ws_2000_2009_ncep1_v50.txt' 
    341394res=read_ascii(fi_ws_ncep1,data_start=1) 
    342395ff=res.field1 
    … …  
    370423spawn, 'mv '+psdir+'idl.ps '+fig 
    371424spawn, 'gv '+fig 
    372 return 
    373425end 
    374 function x_site_location, site 
    375     n1=strpos(site, 's') 
    376 if (n1 gt -1) then begin 
    377     ns=-1. 
    378     x=strmid(site, 0, n1) 
    379     x=float(x)*ns 
    380 endif else begin 
    381     n1=strpos(site, 'n') 
    382     x=strmid(site, 0, n1) 
    383     ny=1. 
    384     x=float(x)*ny 
    385 endelse 
    386 return, float(x) 
    387 end 
    388 function y_site_location, site 
    389     n1=strpos(site, 'e') 
    390 if (n1 gt -1) then begin 
    391     n=strpos(site, 's') 
    392     if (n gt -1) then begin 
    393         y=strmid(site, n+1, n1-n-1) 
    394     endif else begin 
    395         n=strpos(site, 'n') 
    396         y=strmid(site, n+1, n1-n-1) 
    397     endelse 
    398  
    399 endif else begin 
    400     n1=strpos(site, 'w') 
    401     n=strpos(site, 's') 
    402     if (n gt -1) then begin 
    403         y=strmid(site, n+1, n1-n-1) 
    404         y=180+(180-float(y)) 
    405     endif else begin 
    406         n=strpos(site, 'n') 
    407         y=strmid(site, n+1, n1-n-1) 
    408         y=180+(180-float(y)) 
    409    endelse 
    410 endelse 
    411 return,float(y) 
    412 end 
Note: See TracChangeset for help on using the changeset viewer.