Changeset 85 for trunk/src/paper01/fig13


Ignore:
Timestamp:
08/17/11 14:45:28 (13 years ago)
Author:
pinsard
Message:

restart consolidation of paper01 software

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/paper01/fig13/net_flux_validation_scatter_2000_2007.pro

    • Property svn:keywords set to URL
    r50 r85  
    99; =========== 
    1010; 
    11 ; This program will create the following text files with statistics  
    12 ; of respective variables 
     11; This program will create the text files with statistics of respective variables 
     12; 
     13; :file:`${PROJECT_OD}/TropFlux_swr_19890101_20071231_DT.nc` 
     14; containing ++ 
     15; has been produced by 
     16; :ref:`TropFlux_swr_DT_19890101_20071231.pro`. 
     17; 
     18; :file:`${PROJECT_OD}/TropFlux_swr_19890101_20091231_NRT.nc` 
     19; containing ++ 
     20; has been produced by 
     21; :ref:`TropFlux_swr_NRT_19890101_20091231.pro`. 
     22; 
     23; :file:`${PROJECT_OD}/TropFlux_19890101_20091231.nc` 
     24; ; containing ++ 
     25; has been produced by 
     26; :ref:`TropFlux_19890101_20091231.pro`. 
    1327; 
    1428; .. graphviz:: 
    1529; 
    1630;    digraph net_flux_validation_scatter_2000_2007 { 
    17 ;       graph [ 
    18 ;       rankdir="LR", 
    19 ;       ] 
    20 ;       tropflux_general [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/TropFlux_19890101_20091231.nc"]; 
    21 ;       tropflux_dt_swr [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/TropFlux_swr_19890101_20071231_DT_v50.nc"]; 
    22 ;       tropflux_nrt_swr [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/TropFlux_swr_19890101_20091231_NRT_v50.nc"]; 
    23 ;       ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/NCEP2_flux_19890101_20090729.nc"]; 
    24 ;       erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_fluxes_19930101_20090801_TROP_oafluxgrid.nc"]; 
    25 ;       shf_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/shf_oafluxgrid_1985_2008.nc"]; 
    26 ;       lhf_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lhf_oafluxgrid_1985_2008.nc"]; 
    27 ;       lwr_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lwr_oafluxgrid_1985_2007.nc"]; 
    28 ;       ncep1 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/fluxe_ncep1_19890101_20091231.nc"]; 
    29 ;       net_nrt_erai [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/net_2000_2007_erai_v50.txt"]; 
    30 ;       net_nrt_trop [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/net_2000_2007_trop_v50.txt"]; 
    31 ;       net_nrt_ncep [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/net_2000_2007_ncep_v50.txt"]; 
    32 ;       net_nrt_ncep1 [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/net_2000_2007_ncep1_v50.txt"]; 
    33 ;       net_nrt_tropd [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/net_2000_2007_tropd_v50.txt"]; 
     31; 
     32;       tropflux_general [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_19890101_20091231.nc"]; 
     33;       tropflux_dt_swr [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_swr_19890101_20071231_DT.nc"]; 
     34;       tropflux_nrt_swr [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_swr_19890101_20091231_NRT.nc"]; 
     35;       ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/NCEP2_flux_19890101_20090729.nc"]; 
     36;       erai [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_fluxes_19930101_20090801_TROP_oafluxgrid.nc"]; 
     37;       shf_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/shf_oafluxgrid_1985_2008.nc"]; 
     38;       lhf_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/lhf_oafluxgrid_1985_2008.nc"]; 
     39;       lwr_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/lwr_oafluxgrid_1985_2007.nc"]; 
     40;       ncep1 [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/fluxe_ncep1_19890101_20091231.nc"]; 
     41;       net_nrt_erai [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/net_2000_2007_erai.txt"]; 
     42;       net_nrt_trop [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/net_2000_2007_trop.txt"]; 
     43;       net_nrt_ncep [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/net_2000_2007_ncep.txt"]; 
     44;       net_nrt_ncep1 [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/net_2000_2007_ncep1.txt"]; 
     45;       net_nrt_tropd [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/net_2000_2007_tropd.txt"]; 
    3446; 
    3547;       figure [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/net_flux_validation_scatter_2000_2007.ps"]; 
    … …  
    6880; ==== 
    6981; 
    70 ; make it work on cratos : missing data TropFlux_swr_19890101_20071231_DT_v50.nc, etc. 
     82; make it work on cratos : missing data 
    7183; 
    7284; coding rules 
    73 ; 
    74 ; get rid of spwan 
    7585; 
    7686; complete description 
    … …  
    8292; 
    8393; $Id$ 
     94; 
     95; $URL$ 
     96; 
     97; - fplod 20110817T101928Z aedon.locean-ipsl.upmc.fr (Darwin) 
     98; 
     99;   * remove v50 from filenames 
     100;   * replace PROJECT_ID by PROJECT_OD 
     101;   * complete description (not finish) 
    84102; 
    85103; - fplod 20110420T114953Z aedon.locean-ipsl.upmc.fr (Darwin) 
    … …  
    100118key_portrait = 1 
    101119 
    102 openps, FILENAME = 'idl.ps' 
     120openps, FILENAME = project_od_env+'net_flux_validation_scatter_2000_2007.ps' 
    103121 
    104122sitelist=['8s67e','12s55e', '8s55e', '8s80.5e', '1.5s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $ 
    … …  
    120138ocean='net_flux' 
    121139nsmooth=1. 
    122 file_trop=project_id_env+'TropFlux_19890101_20091231.nc' 
     140file_trop=project_od_env+'TropFlux_19890101_20091231.nc' 
    123141initncdf, file_trop 
    124142shf_trop=read_ncdf("shf", date1, date2, file=file_trop,/nostr) 
    … …  
    126144lwr_trop=read_ncdf("lwr", date1, date2, file=file_trop,/nostr) 
    127145 
    128 file_trop=project_id_env+'TropFlux_swr_19890101_20071231_DT_v50.nc' 
     146file_trop=project_od_env+'TropFlux_swr_19890101_20071231_DT.nc' 
    129147initncdf, file_trop 
    130148swr_dt=read_ncdf("swr", date1, date2, file=file_trop,/nostr) 
    … …  
    132150help, trop 
    133151 
    134 file_trop=project_id_env+'TropFlux_swr_19890101_20091231_NRT_v50.nc' 
     152file_trop=project_od_env+'TropFlux_swr_19890101_20091231_NRT.nc' 
    135153initncdf, file_trop 
    136154swr_nrt=read_ncdf("sw", date1-1, date2, file=file_trop,/nostr) 
    … …  
    138156help, tropd 
    139157 
    140 file_ncep=project_id_env+'NCEP2_flux_19890101_20090729.nc' 
     158file_ncep=project_od_env+'NCEP2_flux_19890101_20090729.nc' 
    141159initncdf, file_ncep 
    142160shf=read_ncdf("shf", date1, date2, file=file_ncep,/nostr) 
    … …  
    149167help, ncep 
    150168 
    151 file_erai=project_id_env+'erai_fluxes_19930101_20090801_TROP_oafluxgrid.nc' 
     169file_erai=project_od_env+'erai_fluxes_19930101_20090801_TROP_oafluxgrid.nc' 
    152170initncdf, file_erai 
    153171shf=read_ncdf("shf", date1, date2, file=file_erai,/nostr) 
    … …  
    158176help, erai 
    159177 
    160 file_oaf=project_id_env+'shf_oafluxgrid_1985_2008.nc' 
     178file_oaf=project_od_env+'shf_oafluxgrid_1985_2008.nc' 
    161179initncdf, file_oaf 
    162180shf=read_ncdf("shf", date1, date2, file=file_oaf,/nostr) 
    163 file_oaf=project_id_env+'lhf_oafluxgrid_1985_2008.nc' 
     181file_oaf=project_od_env+'lhf_oafluxgrid_1985_2008.nc' 
    164182initncdf, file_oaf 
    165183lhf=read_ncdf("lhf", date1, date2, file=file_oaf,/nostr) 
    166 file_oaf=project_id_env+'lwr_oafluxgrid_1985_2007.nc' 
     184file_oaf=project_od_env+'lwr_oafluxgrid_1985_2007.nc' 
    167185initncdf, file_oaf 
    168186lwr=read_ncdf("lwr", date1, date2, file=file_oaf,/nostr) 
    169 file_oaf=project_id_env+'swr_oafluxgrid_1985_2007.nc' 
     187file_oaf=project_od_env+'swr_oafluxgrid_1985_2007.nc' 
    170188initncdf, file_oaf 
    171189swr=read_ncdf("swr", date1, date2, file=file_oaf,/nostr) & swr_trop=swr 
    … …  
    173191help, oafl 
    174192 
    175 file=project_id_env+'fluxe_ncep1_19890101_20091231.nc' 
     193file=project_od_env+'fluxe_ncep1_19890101_20091231.nc' 
    176194initncdf, file 
    177195lwr=read_ncdf("lwr", date1, date2, file=file,/nostr) 
    … …  
    185203close,/all 
    186204 
    187 fi_net_NRT_erai=project_id_env+'net_2000_2007_erai_v50.txt' 
     205fi_net_NRT_erai=project_od_env+'net_2000_2007_erai.txt' 
    188206openw,1,fi_net_NRT_erai 
    189 fi_net_NRT_trop=project_id_env+'net_2000_2007_trop_v50.txt' 
     207fi_net_NRT_trop=project_od_env+'net_2000_2007_trop.txt' 
    190208openw,2,fi_net_NRT_trop 
    191 fi_net_NRT_oaflx=project_id_env+'net_2000_2007_oaflx_v50.txt' 
     209fi_net_NRT_oaflx=project_od_env+'net_2000_2007_oaflx.txt' 
    192210openw,3,fi_net_NRT_oaflx 
    193 fi_net_NRT_ncep=project_id_env+'net_2000_2007_ncep_v50.txt' 
     211fi_net_NRT_ncep=project_od_env+'net_2000_2007_ncep.txt' 
    194212openw,4,fi_net_NRT_ncep 
    195 fi_net_NRT_ncep1=project_id_env+'net_2000_2007_ncep1_v50.txt' 
     213fi_net_NRT_ncep1=project_od_env+'net_2000_2007_ncep1.txt' 
    196214openw,5,fi_net_NRT_ncep1 
    197 fi_net_NRT_tropd=project_id_env+'net_2000_2007_tropd_v50.txt' 
     215fi_net_NRT_tropd=project_od_env+'net_2000_2007_tropd.txt' 
    198216openw,6,fi_net_NRT_tropd 
    199217 
    … …  
    276294close,/all 
    277295 
    278 fi_net_NRT_erai=project_id_env+'net_2000_2007_erai_v50.txt' 
     296fi_net_NRT_erai=project_od_env+'net_2000_2007_erai.txt' 
    279297res=read_ascii(fi_net_NRT_erai,data_start=1) 
    280298ff=res.field1 
    … …  
    308326 
    309327 
    310 fi_net_NRT_trop=project_id_env+'net_2000_2007_trop_v50.txt' 
     328fi_net_NRT_trop=project_od_env+'net_2000_2007_trop.txt' 
    311329res=read_ascii(fi_net_NRT_trop,data_start=1) 
    312330ff=res.field1 
    … …  
    340358 
    341359 
    342 fi_net_NRT_oaflx=project_id_env+'net_2000_2007_oaflx_v50.txt' 
     360fi_net_NRT_oaflx=project_od_env+'net_2000_2007_oaflx.txt' 
    343361res=read_ascii(fi_net_NRT_oaflx,data_start=1) 
    344362ff=res.field1 
    … …  
    373391 
    374392 
    375 fi_net_NRT_ncep=project_id_env+'net_2000_2007_ncep_v50.txt' 
     393fi_net_NRT_ncep=project_od_env+'net_2000_2007_ncep.txt' 
    376394res=read_ascii(fi_net_NRT_ncep,data_start=1) 
    377395ff=res.field1 
    … …  
    405423 
    406424 
    407 fi_net_NRT_ncep1=project_id_env+'net_2000_2007_ncep1_v50.txt' 
     425fi_net_NRT_ncep1=project_od_env+'net_2000_2007_ncep1.txt' 
    408426res=read_ascii(fi_net_NRT_ncep1,data_start=1) 
    409427ff=res.field1 
    … …  
    437455 
    438456 
    439 fi_net_NRT_tropd=project_id_env+'net_2000_2007_tropd_v50.txt' 
     457fi_net_NRT_tropd=project_od_env+'net_2000_2007_tropd.txt' 
    440458res=read_ascii(fi_net_NRT_tropd,data_start=1) 
    441459ff=res.field1 
    … …  
    470488closeps 
    471489 
    472 fig=project_od_env+'net_flux_validation_scatter_2000_2007.ps' 
    473 spawn, 'mv '+psdir+'idl.ps '+fig 
    474 spawn, 'gv '+fig 
    475490end 
Note: See TracChangeset for help on using the changeset viewer.