source: trunk/LMDZ.MARS/util/solzenangle.def @ 3493

Last change on this file since 3493 was 2434, checked in by abierjon, 4 years ago

Mars GCM:
Update utilities :
concat.F90

  • rewrite and simplify the handling of time and offset so that any file can be concatenated, including files from different years or stats file
  • use modulo to shift starting sols in concat.nc to values between 0 and 669
  • add a "adls" option in addition to "sol" or "ls" to add a Ls 1D variable

while keeping "sol" as the Time axis

  • add conservation of altitude attributes (long_name,units,positive)
  • enable absence of both hybrid (aps,bps) and sigma levels

solzenangle.F90

  • improve calculation for 1 sol file (stats) to use all local time data
  • read the starting sol in the file instead of asking it to the user (except for stats file) ; keep the possibility for user to change it
  • ask mean Ls value for stats file instead of sol number
  • fix crash when using "all variable" option (ticket #66)
  • fix bug on aps,bps by using GCM_layers dimension instead of altitude

localtime.F90

  • fix crash when using "all variable" option (ticket #66)
  • fix bug on aps,bps by using GCM_layers dimension instead of altitude

For all 3 utilities :

  • handle all ndim cases for the variables (ticket #52)
  • change "title" attribute into "long_name" by default (ticket #48)
  • extend size of long_name character string (ticket #57)
  • remove the use of #ifdef NC_DOUBLE (ticket #67)

AB

File size: 807 bytes
RevLine 
[2280]1diagfi.nc
2tsurf
3ps
4temp
5
[2289]6morning
790
[2280]80
9
10-----------------------------------------------------------------------
[2289]11ABOVE is the list of inputs to be fed to "solzenangle.e" if you don t want
[2280]12to reply directly to the program:
13
141) complete name of file to be read
152) List of X variables to be kept (X lines) or 'all'
163) blank line at the end
[2289]174) planet side wrt to noon (morning or evening)
185) solar zenith angle (within [0;180[ deg)
190deg = zenith ; 90deg = terminator at the surface ; >90deg = terminator in altitude
[2434]206) reference time
21    - for a diagfi/concat :
22        - write "n" to keep the date stored in the file (recommended)     
23        - if needed:  sol value at the beginning of the run, wrt Ls=0
24    - for a stats : Ls value at the middle of the run (degree)
[2280]25USE :
[2289]26solzenangle.e < solzenangle.def
[2280]27
28
Note: See TracBrowser for help on using the repository browser.