Ignore:
Timestamp:
Feb 27, 2012, 10:44:32 AM (13 years ago)
Author:
acolaitis
Message:

27/02/12 == AC

Continuation of thermals setting, comparisons with mesoscale results on Case C
Added possibility to call gcm (or 1d) with constant prescribed sensible heat flux, in the spirit of direct comparisons with LES

... This is directly comparable to the variable tke_heat_flux in namelist.input
... Requires the use of yamada4.F from terrestrial GCM (mixes more, seem more numerically stable)
... Usually requires high timesteps (>1000) to avoids crashes. Best approach is to compare
height of first model level z1 and teta_1 between LES and 1D, and increase the timesteps until results
between the two models are comparable (might require a slitghly different tke_heat_flux between the two models
due to difference in vertical diffusion schemes and subgrid effects)
... Syntax for use is to add "tke_heat_flux = ..." in callphys.def

Corrected some stuff with tke transport in thermals (which is anyway deactivated but one day maybe...)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/libf/phymars/callkeys.h

    r485 r544  
    1818     &   ,dustbin,nqchem_min,nltemodel,nircorr
    1919     
    20       COMMON/callkeys_r/topdustref,solarcondate,semi,alphan
     20      COMMON/callkeys_r/topdustref,solarcondate,semi,alphan,            &
     21     &   tke_heat_flux
    2122     
    2223      LOGICAL callrad,calldifv,calladj,callcond,callsoil,               &
     
    3637      real alphan
    3738      real solarcondate
     39      real tke_heat_flux
    3840
    3941      integer iddist
Note: See TracChangeset for help on using the changeset viewer.