source: trunk/LMDZ.TITAN/deftank/callphys.def.allmuparams

### Orbit / General options
### ~~~~~~~~~~~~~~~~~~~~~~~
# Run with or without tracer transport ?
tracer    = .true.
# Diurnal cycle ?  If diurnal = .false., diurnally averaged solar heating
diurnal   = .true.
# Seasonal cycle ? If season = .false., Ls stays constant, to value set in "start"
season    = .true. 
# Tidally resonant orbit ? Must have diurnal = .false., correct rotation rate in newstart
tlocked   = .false.
# Tidal resonance ratio ? Ratio T_orbit to T_rotation
nres      = 1
# Write some more output on the screen ?
lwrite    = .false.
# Save statistics in file "stats.nc" ?
callstats = .false.
# Test energy conservation of model physics ?
enertest  = .true.
# Check if cpp matches gases.def (should be false in 1d)
check_cpp_match = .false.
# Force cpp and mugaz values from gases.def (should be false in 1d)
force_cpp = .false.
# Compute effective variations of gravity field with altitude in physics ?
# Better be true if you run with chemistry
eff_gz     = .false.

### Directory where external input files are
### ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
datadir   = /data/emoisan/datagcm


### Radiative transfer options
### ~~~~~~~~~~~~~~~~~~~~~~~~~~
# Call radiative transfer?
callrad    = .true.
# The rad. transfer is computed every "iradia" physical timestep
iradia     = 4
# Use blackbody for stellar spectrum ?
stelbbody  = .false.
# Stellar blackbody temperature ?
stelTbb    = 2000.000
# Call multilayer correlated-k radiative transfer ?
corrk      = .true.
# Folder in which correlated-k data is stored ?
corrkdir   = recombin
# Use corr-k recombination instead of pre-mixed ?
corrk_recombin = .true.
# Call visible gaseous absorption in radiative transfer ?
callgasvis = .true.
# Call continuum in radiative transfer ?
Continuum  = .true.
# Allow temperatures above kmatrix tmax
strictboundcorrk = .false.
# Include Rayleigh scattering in the visible ?
rayleigh   = .true.
# Uniform absorption coefficient in radiative transfer?
graybody   = .false.
# Constant absorption coefficient in visible 
#      (in m^2/kg; only if graybody=true):
#      tau_surf= kappa*P/g
kappa_VI   = 0.
# Constant absorption coefficient in IR
#      (in m^2/kg; only if graybody=true):
kappa_IR   = 0.
# Use Newtonian cooling in place of radiative transfer ?
newtonian  = .false.
# Radiative timescale for Newtonian cooling ? [only if newtonian = T]
tau_relax  = 30.00000
# Test physics timescale in 1D ?
testradtimes = .false.
# Characteristic planetary equilibrium (black body) temperature
# This is used only in the aerosol radiative transfer setup. (see aerave.F)
tplanet    = 94.
# Output spectral OLR in 1D/3D?
specOLR    = .false.
# Output global radiative balance in file 'rad_bal.out' - slow for 1D!!
meanOLR    = .false.
# Remove lower boundary (e.g. for gas giant sims)
noradsurf  = .false. 
# Output optical diagnostics (in otpcv & optci) ?
diagdtau   = .false.


### Chemistry
### ~~~~~~~~~
# Call chemistry ? (must have tracer = .true.)
callchim   = .false.
# Activate zonal mean ? 
moyzon_ch  = .true.
# Chemistry is computed every "ichim" phys. timestep
ichim      = 20
# Infinite tank of CH4 ? (.false. == lakes at the poles)
resCH4_inf = .false.


### Microphysics
### ~~~~~~~~~~~~
# Call microphysics ? (must have tracer = .true.)
callmufi   = .true.
# Activate zonal mean ? (must be .true. if callmufi by now ?)
moyzon_mu  = .false.
# If yes, compute clouds ?
callclouds = .true.
# If yes, number of ices ? (must be compatible with traceur.def AND microphysical model)
nices      = 1
# Use new optics for clouds ?
opt4clouds = .true.
# Cloudy fraction in the Clear / Dark column method ?
Fcloudy    = 0.2
# Adjustment of ssca albedo for cloudy dropplet ?
Fssa       = 0.90
# Adjustment of aerosol properties in the VIS ?
FHVIS      = 1.0
# Adjustment of aerosol properties in the IR ?
FHIR       = 1.0
# Desactivate radiative coupling of microphysics (i.e. use mean opacity profile)
uncoupl_optic_haze = .false.
# Use haze seasonal model ? (if you run without coupled microphysics)
seashaze   = .false.

# ~~~~~ Parameters for microphysics ~~~~~
# Microphys. config file ?
config_mufi = config.ne15.cfg
# Pressure level of aer. production (Pa) ?
p_prod     = 1.0
# Aer. production rate (kg.m-2.s-1) ?
tx_prod    = 3.5e-13
# Equivalent radius production (m) ?
rc_prod    = 1.0e-8
# Radius of air (nitrogen) molecule (m) ?
air_rad    = 1.75e-10
# Monomer radius (m)
rm                     = 50e-9
# Fractal dimension
df                     = 2.0
# Aerosol density (kg.m-3)
rho_aer                = 600.
# Enable/disable Haze production process 
haze_production        = T
# Enable/disable Haze coagulation process 
haze_coagulation       = T
# Coagulation interactions, a combination of:
#    0 - no interactions (same as haze_coagulation == F)
#    1 - SS interactions
#    2 - SF interactions
#    4 - FF interactions.
# (for example : 5 = 4+1 --> SS and FF coagulation only)
haze_coag_interactions = 7
# Enable/disable Haze sedimentation process 
haze_sedimentation     = T
# Disable Fiadero correction for sedimentation process
no_fiadero             = T
# Fiadero correction minimum ratio threshold
fiadero_min_ratio      = 0.1
# Fiadero correction maximum ratio threshold
fiadero_max_ratio      = 10.
# Force settling velocity to M0
wsed_m0                = T
# Force settling velocity to M3
wsed_m3                = F
# Enable/disable clouds sedimentation process
# (automatically set to F if clouds microphysics is not enabled)
clouds_sedimentation   = T
# Enable/disable clouds nucleation and condensation processes
# (automatically set to F if clouds microphysics is not enabled)
clouds_nuc_cond        = T
# Condensible species configuration file
# (not needed if clouds microphysics is not enabled)
specie_cfg             = mp2m_species.cfg
# Enable/disable spherical mode transfert probability
transfert_probability = T
# Spherical mode transfert probability look-up tables file
# (optional if 'transfert_probability' is False) 
ps2s_file             = mmp_ps2s_rm50_ne15.nc
# Electric charging coagulation correction
# If set to .false. then no correction is assumed
electric_charging     = T 
# File for the electric charging correction factor.
# (optional if 'electric_charging' is False) 
mq_file               = mmp_qmean_rm50_ne15.nc
### Thresholds related parameters
# Aerosol spherical mode total number min. threshold
m0as_min = 1e-10
# Aerosol spherical mode min. caracteristic radius threshold
rcs_min  = 1e-9
# Aerosol fractal mode total number min. threshold
m0af_min = 1e-10
# Aerosol fractal mode min. caracteristic radius threshold (updated to monomer radius at runtime if needed)
rcf_min  = 1e-9
# cloud drop total number min. threshold
m0n_min = 1e-10


### Star type
### ~~~~~~~~~
startype = 1
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# The choices are:
#
#   startype = 1        Sol        (G2V-class main sequence)
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Stellar flux at 1 AU. Examples:
# 1366.0 W m-2          Sol today
# 1024.5 W m-2          Sol today x 0.75 = weak early Sun
Fat1AU = 1366.0


### Other physics options
### ~~~~~~~~~~~~~~~~~~~~~
# Call turbulent vertical diffusion ?
calldifv    = .true.
# Use turbdiff instead of vdifc ?
UseTurbDiff = .true.
# Call convective adjustment ?
calladj     = .true.
# Call thermal conduction in the soil ?
callsoil    = .true.
# Call nudging of zonal wind ?
nudging_u   = .false.