Changes between Version 10 and Version 11 of Doc/Chap2


Ignore:
Timestamp:
01/22/08 13:03:24 (16 years ago)
Author:
kolasinski
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • Doc/Chap2

    v10 v11  
    8383This in line command allows to see 2 plots on the same window ('Paysage' format). These are the seasonnal cycle of precipitation (called precip) and OLR (called topl) for the '2L24' experiment along the equator ('x' -> longitude, 't' -> time). The variable is averaged zonally on the box 'pac_eq' which is defined in 'domain_boxes.def' file. The grid type is 'lmdzl'. The file should contains 12 mean months from January to December and those mean months are carried out on 100 years from 1860 to 1959. The result is the following picture.[[BR]] 
    8484 
    85  
    8685[[Image(2L24_PRECIP_OLR_XT_PACEQ_1mm_100y.gif, 60%)]]  
    8786 
     87At the end of the 'post_it.pro' file, the program is launched through the following command : 
     88{{{ 
     89def_work, data_base_list, out_ps, cmdline, out_all, other_file,  spec_base_list 
     90}}} 
    8891 
     92=== Usage of command line ? === 
    8993 
     94{{{ 
     95cmdline = [ $ 
     96; var        on   exp     grid  plt               timeave   date1                spec                    disp   proj   out 
     97'lastline 0' ] 
     98}}} 
    9099 
     100 * 'var' : The name of the field is usually the long_name stored in the netcdf file. If you add '@@' at the beginning, it means that post_it will read 'fld_macros.def' to get the right macro to launch (for instance, compute the standard deviation of the sosstsst field through @@sosstdev)  
     101{{{ 
     102; var : <name> of field 
     103;         @@<name> for macro defined in Defaults/fld_macros.def 
     104;         <name1>=f(<name2>) for scatter plot y=f(x) or 
     105;         <name1>=f(next) to use next line as 'x' (uses time interval of name1) 
     106;         <name>@s<sigma> to plot <name> on isopycnal <sigma> 
    91107 
     108  '@@sosstdev  1  CDT3   T      xyt   1m@t412  187001  196912  1     1   v', $ 
     109}}} 
    92110 
     111 * 'on' : If you want to see a plot, set this variable to '1'. When you want several plots on the same window, you may hide one of these by setting 'on' to '2'. It will leave a blank space where the plot should have been drawn. 
     112{{{ 
     113; on    : 0/1 (2 = empty window in multi-window plot) 
    93114 
     115  '@@sosstdev    1   CDT3   T      xyt   1m@t412  187001  196912  2x2     1   v', $ 
     116  '@@sotoxdev   2   CDT3   U      xyt   1m@t412  187001  196912  1     1   v', $ 
     117  '@@sotoxdev   0   CDT3   U      xyt   1m@t412  187001  196912  1     1   v', $ 
     118  '@@sotoydev   0   CDT3   V      xyt   1m@t412  187001  196912  1     1   v', $ 
     119}}} 
     120 
     121 * 'exp' : If you want to make a difference between 2 experiments, you 2 have options. This is the first one. Pay attention to the fact that the 2 files should have the same 'var' name, the same 'grid', the same 'timeave', the same 'plt', the same 'date1'...  
     122{{{ 
     123; exp   : name of exp. For difference of 2 exp : <exp1>-<exp2> 
     124;                                 For division of 2 exp : <exp1>/<exp2> 
     125 
     126  'sozotaux     1   CD2-2L24  U     xt_pac_eq   1mm   01_1860-1959  12_1860-1959   2P   1  v', $ 
     127}}} 
     128 
     129 * 'grid' : The grid type defines the grids (regular or irregular like the ORCA grid) to read. For instance, post_it will read grids_ncpt62.nc which is stored in 'IDL/Defaults/Grids'. This grid is considered as regular (because you specified it in plt_def.pro in the variable 'nc_grid_list'). This file contains 3 variables : the land sea mask fraction, latitudes as a 1D variable and longitudes as a 1D variables. And  post_it through saxo routines will build the grid with the help of those variables : gphit (eg latitudes), glamt (longitudes), tmask (masks on T-points), e1t (scale factor along x), e2t (scale factor along y) as 2D variables... Usually the atmospheric grids are regular. You can get easily the 'grids_...' files through the IPCC database1 with the variable 'sftlf' and build a new “grid”. At the moment, only 2D field can be read from these grid. [[BR]] 
     130If the variable is stored in an irregular grid, things are a bit more complicated to build a new grid readable by  post_it and saxo. You will need to build a file where the variables gphit, glamt, glamu, glamv, tmask... are already stored. For the moment, the ORCA grid (NEMO, IPSL) and the MICOM grid (BCM model, NERSC) can be read (even 3D fields). These grids are Arakawa C-type with T, U, V and F points. 
     131{{{ 
     132; grid  : grid name or @<grid> (to read grid from data file) 
     133 
     134  'precip   1   NCEP  ncpt62  xy  58y    1948  -     1    1   v', $ 
     135}}}