Changes between Version 10 and Version 11 of Doc/Chap2
- Timestamp:
- 01/22/08 13:03:24 (16 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
Doc/Chap2
v10 v11 83 83 This 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]] 84 84 85 86 85 [[Image(2L24_PRECIP_OLR_XT_PACEQ_1mm_100y.gif, 60%)]] 87 86 87 At the end of the 'post_it.pro' file, the program is launched through the following command : 88 {{{ 89 def_work, data_base_list, out_ps, cmdline, out_all, other_file, spec_base_list 90 }}} 88 91 92 === Usage of command line ? === 89 93 94 {{{ 95 cmdline = [ $ 96 ; var on exp grid plt timeave date1 spec disp proj out 97 'lastline 0' ] 98 }}} 90 99 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> 91 107 108 '@@sosstdev 1 CDT3 T xyt 1m@t412 187001 196912 1 1 v', $ 109 }}} 92 110 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) 93 114 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]] 130 If 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 }}}