|
--- 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.
|
