Difference between revisions of "Radiative Generic Condensable Specie"
Noe clement (talk | contribs) (→What to add it in my set up) |
Noe clement (talk | contribs) (→What to add it in my set up) |
||
Line 26: | Line 26: | ||
- in datadir/ (datagcm), add a file « table_tacers_condensable » with the lines : | - in datadir/ (datagcm), add a file « table_tacers_condensable » with the lines : | ||
<syntaxhighlight lang="python"> | <syntaxhighlight lang="python"> | ||
+ | </syntaxhighlight> | ||
+ | |||
+ | - callphys.def | ||
+ | <syntaxhighlight lang="python"> | ||
+ | # Generic Condensable Tracer | ||
+ | # ALREADY AVAILABLE | ||
+ | generic_condensation = .true. # to activate condensation | ||
+ | qvap_deep = 0.20 # kg/kg fixed quantity in the deep atmosphere - works as a sink/source | ||
+ | qvap_top = 1.5e-3 # kg/kg fixed quantity at the top of atmosphere - works as a sink/source | ||
+ | perfect_vap_profile = .true. # to put the cold trap value in the whole stratosphere | ||
+ | generic_rain = .true. # to make condensates precipitating | ||
+ | evap_prec_generic = .true. # to re-evaporate precipitations | ||
+ | rainthreshold_generic = 1e-10 # kg/kg - default value is 0 – how much condensates are kept | ||
+ | evap_coeff_generic = 0.5e5 # evaporation efficiency - default value is 1. | ||
+ | |||
+ | # SOON AVAILABLE | ||
+ | virtual_convadj = .true. # in convadj, virtual potential temperature instead of potential temperature | ||
+ | moistadjustment = .true. # moist adjustment for generic tracer | ||
</syntaxhighlight> | </syntaxhighlight> | ||
Revision as of 10:41, 20 October 2023
What to add it in my set up
- profile_TRACERNAME_vap & profile_TRACERNAME_ice files (similar to « profile » for temperature) to initialize
- in « traceur.def » (flag is « tracer » but name file is « traceur »…), add the lines :
#ModernTrac-v1
2
TRACERNAME_vap is_condensable=1
TRACERNAME_ice is_condensable=1
- in « gases.def » :
# gases
3
H2_ (or other main gaz)
He_ (or other main gaz)
TRACERNAME
0.8
0.15
-1
- in datadir/ (datagcm), add a file « table_tacers_condensable » with the lines :
- callphys.def
# Generic Condensable Tracer
# ALREADY AVAILABLE
generic_condensation = .true. # to activate condensation
qvap_deep = 0.20 # kg/kg fixed quantity in the deep atmosphere - works as a sink/source
qvap_top = 1.5e-3 # kg/kg fixed quantity at the top of atmosphere - works as a sink/source
perfect_vap_profile = .true. # to put the cold trap value in the whole stratosphere
generic_rain = .true. # to make condensates precipitating
evap_prec_generic = .true. # to re-evaporate precipitations
rainthreshold_generic = 1e-10 # kg/kg - default value is 0 – how much condensates are kept
evap_coeff_generic = 0.5e5 # evaporation efficiency - default value is 1.
# SOON AVAILABLE
virtual_convadj = .true. # in convadj, virtual potential temperature instead of potential temperature
moistadjustment = .true. # moist adjustment for generic tracer
Radiative considerations
The model can take into account the radiative effects of Generic condensable specie in the radiative transfer. This is only tested for Hot Jupiter simulations for now. Here is just a copy of the merge request to include the code into the master branch. This will need further editing:
We can activate the scheme by putting aerogeneric = # of aerosols in callphys.def. This is the only needed thing for activating the radiative effects. They must be tracer in modern tracer.def
Commented out the abort if we use more than 4 aerosols
Added reading of optprop files for Radiative Generic Condensable Aerosols
We use the same file for IR and VI channel. For now, only MnS, Cr,Fe and Mg2SiO4 can be read. If you want to add another specie, check the code, it is explained how to quickly do that (right above the Initializations)
Added radii calculation for Radiative Generic Condensable aerosols
Changed the hardcoded size of the totalemit array
The hardcoded size is now 1900 instead of 100 so we don't exceed the array size when working at high spectral resolution (very rare case)
Added opacity computation for Radiative Generic Condensable aerosols
We do this computation in the same fashion as what's performed on water and dust.
switch iniaerosol and initracer order, to prepare for the RGCS scheme
Needed to switch the order of initialization so we can use the RGCS scheme without the assumption that ice and vap tracer of the same specie are following each other in the traceur.def file All my commits regarding the radiative effects of generic tracers. A lot of changes have been made. I've rebased the branch on the latest revision of master + squashed a lot of commits together to avoid the "avalanche de commits". By default, the scheme is NOT activated. Read above on how to activate it
Adding stuff on how to really add a specie (because I'm once again stuck even though I coded that)
go to suaer_corrk. There, need to add the reading of the optical properties file.
Check that in your datadir, the table_tracers_condensable containe the species you wan. You need to have delta_vapH Tref, Pref, mass and metallicity_coeff. In you have all that, you need your traceur.def with the mmol (actually we do't since it's in the table ? check that because weird), the radius and is_condensable=1, is_Rgcs=1