source: trunk/LMDZ.MARS/libf/phymars/callkeys.h @ 308

Last change on this file since 308 was 284, checked in by acolaitis, 13 years ago

--- AC 07/09/2011 ---

  • Added new flag for the Richardson-based surface layer :

callrichsl, can be changed in callphys.def

One should always use the thermals model when using this surface layer model.
Somes cases (weakly unstable with low winds), when not using thermals, won't be well represented by the
Richardson surface layer. This stands for Mesoscale and Gcm but not for LES model.

Correct configs :

callrichsl = .true.
calltherm = .true.

callrichsl = .false.
calltherm = .false.

callrichsl = .false.
calltherm = .true.

Previously unstable config :

callrichsl = .true.
calltherm = .false.

  • To be able to run without thermals and with the new surface layer, a modification has been made to

physiq.F to account for gustiness in GCM and MESOSCALE for negative Richardson, so that :

callrichsl = .true.
calltherm = .false.

can now be used without problems, but is not recommended.

  • Consequently, callrichsl = .false. is now the default configuration for thermals.

We recall the available options in callphys.def for thermals :

outptherm = BOOLEAN (.false. by default) : outputs thermals related quantities (lots of diagfi)
nsplit_thermals = INTEGER (50 by default in gcm, 2 in mesoscale) : subtimestep for thermals model.

It is recommended to use at least 40 in the gcm, and at least 2 in the mesoscale.
The user can lower these values but should check it's log for anomalies or errors regarding
tracer transport in the thermals, or "granulosity" in the outputs for wmax, lmax and hfmax.
File size: 2.0 KB
Line 
1!
2! For Fortran 77/Fortran 90 compliance always use line continuation
3! symbols '&' in columns 73 and 6
4!
5! NB: to keep commons aligned, it is better to split them in groups
6!     of given types (logical, integer, real, ...)
7
8      COMMON/callkeys_l/callrad,calldifv,calladj,callcond,callsoil      &
9     &   ,season,diurnal,lwrite,calllott,callstats,calleofdump          &
10     &   ,callnirco2,callnlte,callthermos,callconduct,calleuv           &
11     &   ,callmolvis,callmoldiff,thermochem,thermoswater,callemis       &
12     &   ,callg2d,linear,rayleigh,tracer,active,doubleq,submicron       &
13     &   ,lifting,callddevil,scavenging,sedimentation,activice,water    &
14     &   ,caps,photochem,calltherm,outptherm,callrichsl,callslope
15     
16      COMMON/callkeys_i/iradia,iaervar,iddist,ilwd,ilwb,ilwn,ncouche    &
17     &   ,dustbin,nqchem_min
18     
19      COMMON/callkeys_r/topdustref,solarcondate,semi,alphan
20     
21      LOGICAL callrad,calldifv,calladj,callcond,callsoil,               &
22     &   season,diurnal,lwrite,calllott                                 &
23     &   ,callstats,calleofdump                                         &
24     &   ,callnirco2,callnlte,callthermos,callconduct,                  &
25     &    calleuv,callmolvis,callmoldiff,thermochem,thermoswater        &
26     &   ,calltherm,outptherm,callrichsl,callslope
27
28
29      logical callemis
30      logical callg2d
31      logical linear
32
33      real topdustref
34      real semi
35      real alphan
36      real solarcondate
37
38      integer iddist
39      integer iaervar
40      integer iradia
41      integer ilwd
42      integer ilwb
43      integer ilwn
44      integer ncouche
45
46      logical rayleigh
47      logical tracer
48      integer dustbin
49      logical active,doubleq,submicron,lifting,callddevil,scavenging
50      logical sedimentation,activice,water,caps
51      logical photochem
52      integer nqchem_min
53
54      integer swrtype ! type of short wave (solar wavelength) radiative
55      ! transfer to use 1: Fouquart 2: Toon.
56      parameter (swrtype=2)
57!      parameter (swrtype=2)
Note: See TracBrowser for help on using the repository browser.