Index: /trunk/LMDZ.GENERIC/deftank/jupiter1d/callphys.def =================================================================== --- /trunk/LMDZ.GENERIC/deftank/jupiter1d/callphys.def (revision 1678) +++ /trunk/LMDZ.GENERIC/deftank/jupiter1d/callphys.def (revision 1678) @@ -0,0 +1,202 @@ +## Directory where external input files are +levels## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +#datadir = /d5/aslmd/DATAGENERIC +datadir = /planeto/sglmd/datagcm + +## Orbit / general options +## ~~~~~~~~~~~~~~~~~~~~~~~ +# Run with or without tracer transport ? +tracer = .false. +# Diurnal cycle ? if diurnal=false, diurnally averaged solar heating +diurnal = .false. +# 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 +# Planet with rings? +rings_shadow = .false. +# Compute latitude-dependent gravity field?? +oblate = .true. +# Flattening (a-b)/a +flatten = 0.06487 +# Needed if oblate=.true.: J2 +J2 = 0.01470 +# Needed if oblate=.true.: Planet mean radius (m) +Rmean = 69911000. +# Needed if oblate=.true.: Mass of the planet (*1e24 kg) +MassPlanet = 1898.3 +# 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 to see if cpp, mugaz values used match gas mixture defined in gases.def (recommended) ? +check_cpp_match=.false. +# +## Radiative transfer options +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call radiative transfer? +callrad = .true. +# the rad. transfer is computed every "iradia" physical timestep +iradia = 8 +# call multilayer correlated-k radiative transfer ? +corrk = .true. +# Include continuum absorption in radiative transfer (note CO2 is treated separately) ? +continuum = .true. +# folder in which correlated-k data is stored ? +corrkdir = Jupiter_NewMOSES_NH3_HITRAN2012_REY_ISO_ge-12_NoKarko +# call visible gaseous absorption in radiative transfer ? +callgasvis = .true. +# Include Rayleigh scattering in the visible ? +rayleigh = .true. +# Uniform absorption coefficient in IR ? +graybody = .false. +# Characteristic planetary equilibrium (black body) temperature +# This is used only in the aerosol radiative transfer setup. (see aerave.F) +tplanet = 100. +# 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. +# Output spectral OLR in 1D/3D? +specOLR = .true. +# Output global radiative balance in file 'rad_bal.out' - slow for 1D!! +meanOLR = .false. +# Variable gas species: Radiatively active ? +varactive = .false. +# Variable gas species: Fixed vertical distribution ? +varfixed = .false. +# Variable gas species: Saturation percentage value at ground ? +satval = 0.0 +# Use fixed vertical profile, 1 step, no iteration ? +kastprof = .false. +# Use blackbody for stellar spectrum ? +stelbbody = .false. +# Stellar blackbody temperature ? +stelTbb = 5800.000 +# Force specific heat capacity and molecular mass values +force_cpp = .true. +# Specific heat capacity in J K-1 kg-1 [only used if force_cpp = T] +cpp = 11500 +# Molecular mass in g mol-1 [only used if force_cpp = T] +mugaz = 2.343 +#Call conduction in upper atmosphere? +callconduc = .false. +#Flux from thermosphere (if callconduc)? +phitop = 1.e-3 + +## Star type +## ~~~~~~~~~ +startype = 1 +# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# The choices are: +# +# startype = 1 Sol (G2V-class main sequence) +# startype = 2 Ad Leo (M-class, synthetic) +# startype = 3 GJ644 +# startype = 4 HD128167 +# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# Stellar flux at 1 AU. Examples: +# 1366.0 W m-2 Sol today +# 1024.5 W m-2 Sol today x 0.75 = weak Sun 3.8 GYa +# 18.462 W m-2 The feeble GJ581 +# 19.960 W m-2 GJ581 with e=0.38 orbital average +Fat1AU = 1366.0 + +## Tracer and aerosol options +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# Gravitational sedimentation of tracers (KEEP FALSE FOR NOW) ? +sedimentation = .false. + +## 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 = .false. +# Internal heat flux (matters only if callsoil=F) +intheat = 5.44 +# Remove lower boundary (e.g. for gas giant sims) +nosurf = .true. + +######################################################################### +## extra non-standard definitions for Gliese 581d / Early Mars +######################################################################### + +## Tracer and aerosol options +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# Radiatively active CO2 aerosol? +aeroco2 = .false. + +# Ammonia cloud (Saturn/Jupiter)? +aeronh3 = .true. +size_nh3_cloud = 15.e-6 +pres_nh3_cloud = 80000 +tau_nh3_cloud = 7.0 +# Radiatively active aerosol (Saturn/Jupiter)? +aeroback2lay = .true. +obs_tau_col_tropo= 2.0 +size_tropo= 1.e-6 +pres_bottom_tropo=80000. +pres_top_tropo=18000. +obs_tau_col_strato= 0.1D0 + +# Auroral aerosols (Saturn/Jupiter)? +aeroaurora = .true. + +# Fixed CO2 aerosol distribution? +aerofixco2 = .false. +# Radiatively active water aerosol? +aeroh2o = .false. +# Fixed water aerosol distribution? +aerofixh2o = .false. +# Varying H2O cloud fraction? +CLFvarying = .false. +# H2O cloud fraction? +CLFfixval = 0.0 +# number mixing ratio of CO2 ice particles +Nmix_co2 = 100000. +# number mixing ratio of water ice particles +Nmix_h2o = 100000. +# basic dust opacity +dusttau = 0.0 +# water cloud pressure level (norm. by psurf) +cloudlvl = 0.0 + + +## Water options +## ~~~~~~~~~~~~~ +# Model water cycle +water = .false. +# Model water cloud formation +watercond = .false. +# Model water precipitation (including coagulation etc.) +waterrain = .false. +# Precipitation threshold (simple scheme only) ? +rainthreshold = 0.011 +# Surface hydrology? +hydrology = .false. +# H2O snow (and ice) albedo ? +albedosnow = 0.6 +# Maximum sea ice thickness (m) ? +maxicethick = 2.0 +# Freezing point of seawater (degrees C) ? +Tsaldiff = 0.0 +# Evolve surface water sources ? +sourceevol = .false. + +## CO2 options +## ~~~~~~~~~~~ +# call CO2 condensation ? +co2cond = .false. +# Set initial temperature profile to 1 K above CO2 condensation everywhere? +nearco2cond = .false. Index: /trunk/LMDZ.GENERIC/deftank/jupiter1d/gases.def =================================================================== --- /trunk/LMDZ.GENERIC/deftank/jupiter1d/gases.def (revision 1678) +++ /trunk/LMDZ.GENERIC/deftank/jupiter1d/gases.def (revision 1678) @@ -0,0 +1,17 @@ +# gases +5 +H2_ +He_ +CH4 +C2H2 +C2H6 +0.863 +0.134 +0.0018 +1.e-7 +1.e-5 +# First line is number of gases +# Followed by gas names (always 3 characters) +# and then molar mixing ratios. +# mixing ratio -1 means the gas is variable. + Index: /trunk/LMDZ.GENERIC/deftank/jupiter1d/profile =================================================================== --- /trunk/LMDZ.GENERIC/deftank/jupiter1d/profile (revision 1678) +++ /trunk/LMDZ.GENERIC/deftank/jupiter1d/profile (revision 1678) @@ -0,0 +1,65 @@ + 200.0000 + 242.0764 + 238.3593 + 234.7491 + 227.3349 + 214.2193 + 198.9277 + 184.7276 + 171.5412 + 159.2960 + 147.9250 + 137.3656 + 127.9115 + 122.1330 + 118.1073 + 115.3581 + 113.7415 + 113.0411 + 113.1916 + 114.1231 + 115.9067 + 118.6585 + 122.4839 + 126.9294 + 132.2991 + 138.3907 + 142.6255 + 146.7856 + 150.1914 + 152.9115 + 155.0030 + 156.5529 + 157.7320 + 158.5769 + 159.1656 + 159.5710 + 159.7470 + 160.1693 + 160.1162 + 160.1524 + 159.9328 + 159.6921 + 159.4074 + 159.0866 + 158.7370 + 158.3921 + 158.0377 + 157.6633 + 157.2459 + 156.7444 + 156.1365 + 155.4394 + 154.6760 + 153.8250 + 152.8934 + 151.8827 + 150.8160 + 149.7348 + 148.6630 + 147.6131 + 146.6401 + 145.7417 + 144.9731 + 144.2265 + 143.0092 Index: /trunk/LMDZ.GENERIC/deftank/jupiter1d/rcm1d.def =================================================================== --- /trunk/LMDZ.GENERIC/deftank/jupiter1d/rcm1d.def (revision 1678) +++ /trunk/LMDZ.GENERIC/deftank/jupiter1d/rcm1d.def (revision 1678) @@ -0,0 +1,98 @@ +#---------------------------------# +# Run parameters for the 1D model # +#---------------------------------# + +## Planetary constants +## ~~~~~~~~~~~~~~~~~~~~~~~ +## NB: those are mandatory +## ~~~~~~~~~~~~~~~~~~~~~~~ +## LENGTH OF A DAY in s +daysec = 35733. +## GRAVITY in m s-2 +g = 24.79 +##PLANETARY RADIUS in m +rad = 7.1492E+07 +## LENGTH OF A YEAR in days +year_day = 10475.5 +## MIN DIST STAR-PLANET in AU [periastron] +periastr = 4.95 +## MAX DIST STAR-PLANET in AU [apoastron] +apoastr = 5.46 +## DATE OF PERIASTRON in days +peri_day = 1520 +## OBLIQUITY in deg +obliquit = 3.13 +## SURFACE PRESSURE in Pa +psurf = 3.e5 +#psurf = 66000 + +## Time integration parameters +## ~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# Initial date (in solar days,=0 at Ls=0) +day0 = 0 +# Initial local time (in hours, between 0 and 24) +time = 12 +# Number of time steps per sol +day_step = 1 +# Number of sols to run +ndt = 63000 +# ecriture physiq +ecritphy = 200 + +## Vertical levels +## ~~~~~~~~~~~~~~~ +# hybrid vertical coordinate ? (.true. for hybrid and .false. for sigma levels) +hybrid = .false. +# autocompute vertical discretisation? (useful for exoplanet runs) +autozlevs = .false. +# Ceiling pressure (Pa) ? (only used if autozlevs = .true.) +pceil = 0.2 + +## Thermal properties +## ~~~~~~~~~~~~~~~~~~ +# Simulate global averaged conditions ? +global1d = .false. +# Latitude (deg) [only used if global1d = F] +latitude = 13. +# Solar Zenith angle (deg) [only used if global1d = T] +szangle = 60 +# Albedo of bare ground +albedo = 0. +# Emissivity of bare ground +emis = 1.0 +# Soil thermal inertia (SI) +inertia = 0.001 +# Initial CO2 ice on the surface (kg.m-2) +co2ice = 0 + +## Wind profile +## ~~~~~~~~~~~~ +## zonal eastward component of the geostrophic wind (m/s) +u = 0. +# meridional northward component of the geostrophic wind (m/s) +v = 0. + +## Initial atmospheric temperature profile +## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# +# Type of initial temperature profile +# ichoice=1 Constant Temperature: T=tref +# [Mars] ichoice=2 Savidjari profile (as Seiff but with dT/dz=cte) +# [Mars] ichoice=3 Lindner (polar profile) +# [Mars] ichoice=4 inversion +# [Mars] ichoice=5 Seiff (standard profile, based on Viking entry) +# ichoice=6 constant T + gaussian perturbation (levels) +# ichoice=7 constant T + gaussian perturbation (km) +# ichoice=8 Read in an ascii file "profile" +ichoice = 8 +# Reference temperature tref (K) +tref = 170 +# Add a perturbation to profile if isin=1 +isin = 0 +# peak of gaussian perturbation (for ichoice=6 or 7) +pic = 26.522 +# width of the gaussian perturbation (for ichoice=6 or 7) +largeur = 10 +# height of the gaussian perturbation (for ichoice=6 or 7) +hauteur = 30. + Index: /trunk/LMDZ.GENERIC/deftank/jupiter1d/traceur.def =================================================================== --- /trunk/LMDZ.GENERIC/deftank/jupiter1d/traceur.def (revision 1678) +++ /trunk/LMDZ.GENERIC/deftank/jupiter1d/traceur.def (revision 1678) @@ -0,0 +1,2 @@ +0 +dummy Index: /trunk/LMDZ.GENERIC/deftank/jupiter1d/z2sig.def =================================================================== --- /trunk/LMDZ.GENERIC/deftank/jupiter1d/z2sig.def (revision 1678) +++ /trunk/LMDZ.GENERIC/deftank/jupiter1d/z2sig.def (revision 1678) @@ -0,0 +1,65 @@ +25.0 + 1.00000 + 2.00000 + 3.00000 + 5.00000 + 11.0000 + 17.0000 + 23.0000 + 29.0000 + 35.0000 + 41.0000 + 47.0000 + 53.0000 + 59.0000 + 65.0000 + 71.0000 + 77.0000 + 83.0000 + 89.0000 + 95.0000 + 101.000 + 107.000 + 113.000 + 119.000 + 125.000 + 131.000 + 137.000 + 143.000 + 149.000 + 155.000 + 161.000 + 167.000 + 173.000 + 179.000 + 185.000 + 191.000 + 197.000 + 203.000 + 209.000 + 215.000 + 221.000 + 227.000 + 233.000 + 239.000 + 245.000 + 251.000 + 257.000 + 263.000 + 269.000 + 275.000 + 281.000 + 287.000 + 293.000 + 299.000 + 305.000 + 311.000 + 317.000 + 323.000 + 329.000 + 335.000 + 341.000 + 347.000 + 353.000 + 359.000 + 365.000