source: trunk/LMDZ.VENUS/deftank/physiq-96x96x90-chemistry-IONOSPHERE.def @ 3094

Last change on this file since 3094 was 2836, checked in by abierjon, 2 years ago

VENUS GCM:

INCLUDING THE IONOSPHERE CODE IN VENUS GCM


ATTENTION: INCREASED MEMORY DEMAND

NEEDS AT LEAST 135 GB of allocated memory

==============================================================
===> LIST OF APPENDED FILES AND INTERNAL ADDITIONS <==========
==============================================================

NEW MODEL SPECIES

  • deftank/traceur-chemistry-IONOSPHERE.def

Neutrals: 4 Species: N, NO, NO2, N(2D)

Ions: 15 species:

CO2+, CO+, O+, O2+, H+,

N2+, H2O+, OH+, C+, HCO+,

H3O+, HCO2+, N+, NO+, elec

FROM 36 to 55 chemical species

NEW KEYWORD OF PHYSIQ.DEF

  • deftank/physiq-96x96x90-chemistry-IONOSPHERE.def
  • ok_ionchem: keyword supposed to activate ion chemistry.

(be careful that n, no, no2, n2d and the ion species are in the deftank/traceur-chemistry-IONOSPHERE.def)

  • ok_jonline: keyword supposed to activate the online photochemistry

(be careful that n, no, no2 and n2d are in the traceur-chemistry-IONOSPHERE.def)

==============================================================
===> LIST OF MODIFIED PROGRAMS <=========================
==============================================================

nonoro_gwd_ran_mod.F90

  • Change EPFLUXMAX value from 5.E-3 to 1.E-3

photochemistry_venus.F90 (krates; photolysis_online; indices;)

  • Import of keywords: ok_ionchem, tuneupperatm & ok_jonline
  • addition of ion species in order
  • Forcing electroneutrality
  • Update of the reactions a001 and a002 with taking into account the other species
  • Change of formula for a002 with the formulation of Baulch et al., 1976 (confirmed by Smith and Robertson, 2008)

photolysis_mod.F90

  • Modification of the rdsolarflux subroutine to include interpolation with Atlas1 and Atlas3 in connection with E10.7

concentration2.F90

  • Added chemical species n2d, no, no2 and n in the conductivity calculation.

The ions have been excluded because their sum is 105 times less dense than the neutrals and
their thermal conductivity is unknown

iono_h.F90

  • Addition of the phdisrate routine
  • replace the electronic temperature of Mars by that of

origin = 1: Theis et al. 1980 (Venus) with bilinear interpolation altitude/cos(SZA)
origin = 2: Theis et al. 1984 (Venus) with the formula of the electronic temperature at high solar activity

  • addition of an ion temperature model based on VIRA

cleshphys.h

  • added ok_ionchem & ok_jonline in COMMON/clesphys_l/

conf_phys.f90

  • add ok_ionchem & ok_jonline as parameters to read from physiq.def file set to .false. by default

chemparam_mod.F90

  • Add species in M_tr and corresponding i_X. Set all i_X to zero before reading traceur-chemistry-IONOSPHERE.def
  • Added Type_tr table to differentiate species: 1 == neutral, 2 == ION, 3 == liquid, 10 == others


euvheat.F90; hrtherm.F; jthermcalc_e107.F; param_read_e107.F

  • Normalization with Mars

A.M

File size: 3.5 KB
RevLine 
[2836]1## $Header$
2#
3##  Include diurnal cycle or not
4cycle_diurne=y
5##  Include soil model or not
6soil_model=y
7## Use orodr or not for orography
8ok_orodr=n
9## Use orolf or not for orography
10ok_orolf=n
11## Use non-orographic Gravity Waves of not
12ok_gw_nonoro=y
13## Number of calls to the radiative transfer (per day)
14nbapp_rad=24000
15## Number of calls to the chemistry routines (per day)
16nbapp_chim=1
17##  Flag for convection : 1 pour LMD, 2 pour Tiedtke, 3 KE(new version JYG), 30 KE(version IPCC AR4), 4 KE vect
18iflag_con=0
19#
20# orbital parameters
21##~~~~~~~~~~~~~~~~~~
22#
23# VENUS ##
24# R_ecc = 0.006787
25R_ecc = 0.
26R_peri = 0.
27# R_incl = 177.4
28R_incl = 0.
29# solar: effective flux, given at 1 UA
30solaire = 328.
31#
32# parameters for the tracers
33##~~~~~~~~~~~~~~~~~~~~~~~~~
34#
35# 0: nothing (passive tracers)
36# 1: pseudo-chemistry relaxation (phytrac_relax)
37# 2: surface emission (phytrac_emiss)
38# 3: full chemistry (phytrac_chem)
39tr_scheme = 3
40#
41# Reinitialization of tracer abundances
42reinit_trac=n
43#
44# parameters for chemistry and microphysics
45##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
46#
47# use chemistry ?
48ok_chem=y
49# use clouds ? (needed for chemistry)
50ok_cloud=y
51# if yes, which scheme:
52#  1 => simple scheme (Aurelien Stolzenbach)
53#  2 => full microphysical scheme (Sabrina Guilbon)
54cl_scheme=1
55# use sedimentation (goes with clouds)
56ok_sedim=y
57nb_mode=3
58#
59# parameters for the boundary layer
60##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
61#
62iflag_pbl=8
63z0 = 0.01
64lmixmin = 35.
65ksta = 1.e-7
66ok_kzmin=n
67#
68# Surface thermal inertia
69inertie=2000
70#
71# dry convection parameters
72##~~~~~~~~~~~~~~~~~~~~~~~~~
73#
74iflag_ajs = 1
75
76# Solar radiation module
77##~~~~~~~~~~~~~~~~~~~~~~
78## 1 = Rainer Haus Tables
79## 2 = Generic solar module
80#
81solarchoice=1
82
83# Thermosphere and nlte parameters
84##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
85#
86## Compute non-LTE or not
87callnlte=y
88## choice of species profile to use
89nltemodel=2
90## Compute CO2 IR absorption or not
91callnirco2=y
92##CO2 IR Absorption model
93nircorr=0
94##include thermosphere or not
95callthermos=y
96
97## Thermospheric options
98##~~~~~~~~~~~~~~~~~~~~~~
99## Tuning of photochemistry for oxygen production ?
100tuneupperatm = n
101
102## include photolysis j online calculated or not
103ok_jonline = y
104
105## include ion chemistry or not
106ok_ionchem = y
107
108#Method to include solar variability?
109#0-> Old method   1-> Variability with E10.7 as observed
110solvarmod=0
111
112## (Solar min=70 ave=140 max=300)
113fixed_euv_value = 200.
114
115# value for the UV heating efficiency
116##(experimental values between 0.19 and 0.23, lower values may
117## be used to compensate for low 15 um cooling)
118euveff = 0.195       
119#
120#
121#
122#
123# Parameters for IOIPSL output files
124##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
125## OLD. Now we use XIOS => see context_lmdz_physics.xml to taylor the output files
126#
127### OK_journe= y for daily output file histday.nc, =n no histday.nc output
128### Meaningless for Venus
129OK_journe=n
130### OK_mensuel= y for monthly output file histmth.nc, =n no histmth.nc
131### For Venus, only these averaged outputs
132OK_mensuel=n
133## rate (in days) at which the Venus histmth file is to be written               
134# sets the output rate in histmth and/or histins
135ecritphy=0.1
136### OK_instan=y, make some "instantaneous" outputs (same rate as histmth)
137OK_instan=n
138#
139# Output levels for the various output files
140#
141# output level for  "day" lev_histday
142# - lev_hist*=1 => baseline 2D fields
143# - lev_hist*=2 => baseline 3D fields (default)
144# - lev_hist*=3 => radiative transfert
145# - lev_hist*=4 => 3D tendencies
146# - lev_hist*=5 => tracers and others
147lev_histday=2
148#output level for "mth" lev_histmth
149lev_histmth=2
150#output level for "ins" lev_histins
151lev_histins=2
152
Note: See TracBrowser for help on using the repository browser.