Changeset 2432 for trunk/LMDZ.MARS/util/README
- Timestamp:
- Nov 14, 2020, 3:27:44 PM (4 years ago)
- File:
-
- 1 moved
-
trunk/LMDZ.MARS/util/README (moved) (moved from trunk/LMDZ.MARS/util/README.exec) (8 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/util/README
r2431 r2432 1 This directory contains executable files (and their source) that can be1 This directory contains source codes of utilities (programs) that can be 2 2 used to process the LMD Mars GCM output files (like "diagfi" or 3 "stats") , whatever the dimension.4 In addition most output file from one of these program can be5 processe sby another.3 "stats"). 4 In addition most output file from one of these program can then be 5 processed by another. 6 6 7 The executable files should run on any Linux platform. 7 See script "compile" for instructions and recommendations to compile 8 the utilities. Once adequately modified, that script can be used to 9 compile each utility; e.g. to build the zrecast.e utility program: 10 > compile zrecast 8 11 9 Inputs can be provided by (1) replying to questions on screen or10 (2) filling the corresponding *.def and directthe input of these *.def12 Note that inputs can be provided by (1) replying to questions on screen or 13 (2) filling the corresponding *.def and redirecting the input of these *.def 11 14 files instead. For example : 12 15 > concatnc.e < concatnc.def … … 14 17 15 18 -------------------------------------------------------------------- 16 1) concatnc.e 19 1) concatnc.e : concatenate successive GCM output files 17 20 -------------------------------------------------------------------- 18 21 … … 23 26 software like Grads may not be able to use it. To have a linear "Ls" 24 27 timescale, you can use "Ls_Linear.e" (see below). 28 You can also add a "Ls" variable corresponding solar longitude for each 29 timestep by using the "adls" option. 25 30 26 31 Output file is : concat.nc 27 32 28 MODIFICATION:29 07/2008 Utility concatnc.F90 (not used by the gcm): improvement in order30 to 1) concatenate 1D variable and 2) increase the number of input files31 up to 100032 33 34 33 -------------------------------------------------------------------- 35 2) localtime.e34 2) zrecast.e : put GCM data in physical vertical coordinate 36 35 -------------------------------------------------------------------- 37 36 38 Program to redistribute and interpolate the variable a the same 39 local times everywhere (useful to mimic satellite observations, or 40 analyse day to day variations at a given local time). 41 input : diagfi.nc / concat.nc / stats.nc kind of files 42 43 output file is : 44 name_of_input_file_LT.nc with pressure coordinate 45 46 -------------------------------------------------------------------- 47 3) zrecast.e 48 -------------------------------------------------------------------- 37 GCM outputs are in GCM hybrid coordinate which do not correspond to any 38 physical vertical coordinate ! zrecast is NECESSARY to make any publishable 39 scientific figure. 49 40 50 41 This program reads 4D (lon-lat-alt-time) fields from GCM output files … … 65 56 - surface pressure 66 57 - atmospheric temperature 67 - hybrid coordinates aps() and bps(), or sigma levels() (see section 68 1.3.2) 58 - hybrid coordinates aps() and bps(), or sigma levels() 69 59 - ground geopotential (in input file; if not found, it is sought 70 60 in a 'diagfi.nc' file. If not found there, it is then sought in … … 77 67 2.2 in source file) 78 68 69 - Vertical interpolation : note that density (kg.m-3) and species density 70 (e.g. in molecules.cm-3) must be vertically interpolated in log-space. For 71 that purpose, these variables names must start with "rho" (exemple : rho for 72 density, rho_co2, rho_o2, etc....) 73 79 74 output file is : 80 75 name_of_input_file_P.nc with pressure coordinate … … 83 78 name_of_input_file_R.nc with altitude as distance to center of planet 84 79 85 MODIFICATION : 86 01/2010 : correction to interpolate above surface if density is not available. 87 03/2011 : added possibility to have output as distance to center of planet 80 -------------------------------------------------------------------- 81 3) localtime.e : put GCM data in local time coordinate 82 -------------------------------------------------------------------- 83 84 Program to redistribute and interpolate the variable a the same 85 local times everywhere (useful to mimic satellite observations, or 86 analyse day to day variations at a given local time). 87 input : diagfi.nc / concat.nc / stats.nc kind of files 88 89 output file is : name_of_input_file_LT.nc 88 90 89 91 -------------------------------------------------------------------- 90 4) lslin.e 92 4) lslin.e : redistribute and AVERAGE gcm output in Ls coordinate. 91 93 -------------------------------------------------------------------- 92 94 93 95 This program has been designed to interpol data in Solar Longitude (Ls) 94 96 linear time coordinate (usable with grads) from Netcdf diagfi or concatnc 95 files. 96 output file is : lslin.nc 97 lslin also create a lslin.ctl file that can be read 97 files and to AVERAGE the instantaneous data in Ls period (for comparison with 98 binned dataset, for instance). 99 100 lslin also creates a lslin.ctl file that can be read 98 101 directly by grads (>xdfopen lsllin.ctl) to plot in Ls coordinate to 99 102 avoid some problem with grads when grads think that "the time interval 100 103 is too small"... 101 104 102 MODIFICATION 103 10/2007 Utility lslin.F90 (not used by the gcm) 104 changed evaluation of 'start_var' from hard-coded values to a computed value 105 04/2015 Added possibility to bin data (instead of interpolating) over 106 the time intervals 105 Output file is : name_of_input_file_LS.nc 107 106 108 107 -------------------------------------------------------------------- 109 5) hrecast.e108 5) solzenangle.e : select GCM data at a given solar zenith angle 110 109 -------------------------------------------------------------------- 111 110 112 This program can interpolate GCM output on any horizontal grid (regular lat - lon) as long as it cover all the 113 planet. The grid can be given points by points. The best way is to use the redirected input hrecast.def 111 Program to redistribute and interpolate the variable a the same 112 solar zenith angle (notably useful to mimic satellite observations and in particular 113 SOLAR OCCULTATIONS by choosing solar zenith angle = 90°) 114 The user choose between Morning and Evening side. 115 116 input : diagfi.nc / concat.nc / stats.nc kind of files 117 118 On the morning side output file (1 output per sol per grid point) is : 119 name_of_input_file_MO.nc 120 121 On the evening side output file is (1 output per sol per grid point): 122 name_of_input_file_EV.nc 123 124 -------------------------------------------------------------------- 125 6) hrecast.e : interpolate data at another horizontal grid resolution. 126 -------------------------------------------------------------------- 127 128 This program can interpolate GCM output on any horizontal grid (regular 129 lat - lon) as long as it cover all the planet. 130 The grid can be given points by points. The best way is to use the 131 redirected input hrecast.def 114 132 115 133 hrecast.e < hrecast.def … … 117 135 118 136 -------------------------------------------------------------------- 119 6) expandstartfi.e137 7) expandstartfi.e : to plot data in a startfi.nc file 120 138 -------------------------------------------------------------------- 121 139 140 In startfi.nc file, data are not plotable because the horizontal coordinate 141 in only a 1D list of all atmospheric columns on the planets. 122 142 This program takes a physics start file ("startfi.nc") and recasts it 123 143 on the corresponding lonxlat grid (so it contents may easily be displayed 124 using Grads, Ferret, etc.)144 using Python, Grads, Ferret, etc.) 125 145 126 146 Simply run expandstartfi.e as a command line with arguments: … … 132 152 133 153 -------------------------------------------------------------------- 134 6) extract.e 154 8) extract.e : get data at specific coordinates for comparison with observations 135 155 -------------------------------------------------------------------- 136 156
Note: See TracChangeset
for help on using the changeset viewer.
