Changeset 49 for trunk


Ignore:
Timestamp:
04/22/11 12:22:17 (13 years ago)
Author:
pinsard
Message:

improve plot_tropflux

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/plot_tropflux.pro

    r39 r49  
    88; 
    99; Produce +todo+ figures with TROPFLUX data 
     10; 
     11; **tropflux_varname* must be one one these strings : 
     12;    q2m t2m sst ws swr lwr lhf shf  
    1013; 
    1114; .. graphviz:: 
     
    1518;       rankdir="LR", 
    1619;       ] 
    17 ;       file_tropflux_q2m [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_q2m_19890101-20091231.nc"]; 
    18 ;       file_tropflux_t2m [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/tropFlux_t2m_19890101-20091231.nc"]; 
    19 ;       file_tropflux_ws [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_ws_19890101-20091231.nc"]; 
    20 ;       file_tropflux_sst [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_sst_19890101-20091231.nc"]; 
    21 ;       file_tropflux [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_19890101-20091231.nc"]; 
     20;       file_tropflux_q2m [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_q2m_19890101_20091231.nc"]; 
     21;       file_tropflux_t2m [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/tropFlux_t2m_19890101_20091231.nc"]; 
     22;       file_tropflux_ws [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_ws_19890101_20091231.nc"]; 
     23;       file_tropflux_sst [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_sst_19890101_20091231.nc"]; 
     24;       file_tropflux [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_19890101_20091231.nc"]; 
    2225;       figure_++_tropflux_++ [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/plot_tropflux_++_year.png"]; 
    2326; 
     
    3538; ======== 
    3639; 
    37 ; To plot ++ of q2m :: 
    38 ; 
    39 ;    IDL> @tropflux_init 
    40 ;    IDL> status = plot_tropflux('q2m') 
     40; To plot ++ of q2m from :file:`${TROPFLUX}/TropFlux_q2m_19890101_20091231.nc`:: 
     41; 
     42;    IDL> @tropflux_init 
     43;    IDL> tropflux_varname='q2m' 
     44;    IDL> status = plot_tropflux(tropflux_varname) 
    4145; 
    4246; **status** is 0 if everything is ok, -1 otherwise 
     
    4448; :file:`${TROPFLUX_OD}/plot_tropflux_q2m_year.png` and 
    4549; :file:`${TROPFLUX_OD}/plot_tropflux_q2m_month.png` have been created. 
     50; 
     51; Same idea for t2m from :file:`${TROPFLUX}/TropFlux_t2m_19890101_20091231.nc`:: 
     52; 
     53;    IDL> @tropflux_init 
     54;    IDL> tropflux_varname='t2m' 
     55;    IDL> status = plot_tropflux(tropflux_varname) 
     56; 
     57; Same idea for sst from :file:`${TROPFLUX}/TropFlux_sst_19890101_20091231.nc`:: 
     58; 
     59;    IDL> @tropflux_init 
     60;    IDL> tropflux_varname='sst' 
     61;    IDL> status = plot_tropflux(tropflux_varname) 
     62; 
     63; Same idea for ws from :file:`${TROPFLUX}/TropFlux_ws_19890101_20091231.nc`:: 
     64; 
     65;    IDL> @tropflux_init 
     66;    IDL> tropflux_varname='ws' 
     67;    IDL> status = plot_tropflux(tropflux_varname) 
     68;  
     69; Same idea for swr from :file:`${TROPFLUX}/TropFlux_19890101_20091231.nc`:: 
     70; 
     71;    IDL> @tropflux_init 
     72;    IDL> tropflux_varname='swr' 
     73;    IDL> status = plot_tropflux(tropflux_varname) 
     74; 
     75; Same idea for lwr from :file:`${TROPFLUX}/TropFlux_19890101_20091231.nc`:: 
     76; 
     77;    IDL> @tropflux_init 
     78;    IDL> tropflux_varname='lwr' 
     79;    IDL> status = plot_tropflux(tropflux_varname) 
     80; 
     81; Same idea for lhf from :file:`${TROPFLUX}/TropFlux_19890101_20091231.nc`:: 
     82; 
     83;    IDL> @tropflux_init 
     84;    IDL> tropflux_varname='lhf' 
     85;    IDL> status = plot_tropflux(tropflux_varname) 
     86; 
     87; Same idea for shf from :file:`${TROPFLUX}/TropFlux_19890101_20091231.nc`:: 
     88; 
     89;    IDL> @tropflux_init 
     90;    IDL> tropflux_varname='shf' 
     91;    IDL> status = plot_tropflux(tropflux_varname) 
    4692; 
    4793; SEE ALSO 
     
    63109; ==== 
    64110; 
    65 ; make it work. when ok set debug to 0 or add a keyword 
    66 ; same for overwrite 
     111; define plots to be put on the website 
     112; 
     113; add a keyword debug 
     114; add a keyworf overwrite 
    67115; 
    68116; lecture du mask oaflux pour initcdf 
     
    70118; coding rules 
    71119; 
    72 ; parametre date debu date fin ? 
     120; parametre date debut date fin ? 
    73121; 
    74122; parametre version ? 
     
    81129; add parameter saveimage 
    82130; 
    83 ; add parameter overwrite 
    84 ; 
    85131; ajuster le format à la dynamique de la variable 
     132; 
    86133; check arg 
    87 ; 
    88 ; pourquoi read_ncdf n'ouvre pas tout de suite le ficher sans passer par la case widget !! 
    89134; 
    90135; loop on month : min and max values for same color scale 
     
    100145; EVOLUTIONS 
    101146; ========== 
     147; 
     148; $Id$ 
     149; 
     150; - fplod 20110422T095159Z aedon.locean-ipsl.upmc.fr (Darwin) 
     151; 
     152;   * fix read_ncdf call 
     153;   * improve example 
     154;   * fix typo 
    102155; 
    103156; - fplod 20110408T101102Z cratos.locean-ipsl.upmc.fr (Linux) 
     
    119172status=-1 
    120173; 
    121 debug=1 
     174debug=0 
    122175overwrite=1 
    123176; 
    124  
    125177; check for tropflux_varname 
    126178; and build input filename 
     
    215267ENDIF 
    216268; 
    217 fullfilename_out_year=iodirout+'plot_tropflux' + tropflux_varname + '_year.png' 
     269fullfilename_out_year=iodirout+'plot_tropflux_' + tropflux_varname + '_year.png' 
    218270print, fullfilename_out_year 
    219271; in order to avoid unexpected overwritten 
     
    226278ENDIF 
    227279; 
    228 fullfilename_out_month=iodirout+'plot_tropflux' + tropflux_varname + '_month.png' 
     280fullfilename_out_month=iodirout+'plot_tropflux_' + tropflux_varname + '_month.png' 
    229281print, fullfilename_out_month 
    230282; in order to avoid unexpected overwritten 
     
    249301; 
    250302; Read data 
    251 tropflux_var=read_ncdf(tropflux_varname,19890101L,19891231L,fullfilename_in[0],/nostruct) 
     303tropflux_var=read_ncdf(tropflux_varname,19890101L,19891231L,filename=fullfilename_in[0],/nostruct) 
    252304; 
    253305IF debug EQ 1 THEN BEGIN 
     
    266318   nl=4 
    267319   nc=2 
    268    ; ++ watinfg for final plot do be done 
     320   ; ++ wating for final plot do be done 
    269321   jpt=nc*nl 
    270322   FOR itime = 0L, jpt - 1 DO BEGIN 
Note: See TracChangeset for help on using the changeset viewer.