Changeset 1042 for trunk/MESOSCALE


Ignore:
Timestamp:
Sep 19, 2013, 2:15:35 AM (11 years ago)
Author:
aslmd
Message:

MESOSCALE. def files closer to callphys.def.MCD5

Location:
trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER
Files:
2 added
4 edited
2 copied

Legend:

Unmodified
Added
Removed

  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars0

    r1041 r1042  
    11##General options
    22##~~~~~~~~~~~~~~~
    3 #Directory where external input files are
    4 datadir = /san/home/spiga/DATAFILE/
    5 datadir = /d5/aslmd/datafile/
    6 datadir=/users/forget/WWW/datagcm/datafile/
     3#Run with or without tracer transport ?
     4tracer=.false.
    75
    8 #Run with or without tracer transport ?
    9 tracer=.true.
     6#Directory where external input files are:
     7datadir=/planeto/emlmd/run_mcd5.v0_new/datadir
     8#datadir = /san/home/spiga/DATAFILE
     9#datadir = /d5/aslmd/datafile
     10#datadir=/users/forget/WWW/datagcm/datafile
    1011
    1112#Diurnal cycle ?  if diurnal=False, diurnal averaged solar heating
    … …  
    2728## ~~~~~~~~~~~~~
    2829#  =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario,
    29 #  =4 Mars Year 24 from TES assimilation (same as =24 for now)
    30 #  =24 Mars Year 24 from TES assimilation (ie: MCD reference case)
    31 #  =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm)
    32 #  =26 Mars Year 26 from TES assimilation
    33 iaervar = 26
    34 # Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case)
    35 tauvis=0.2
     30#  =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file)
     31#  =6 "cold" (low dust) scenario ; =7 "warm" (high dust) scenario
     32#  =24 Mars Year 24  ; =25 Mars Year 25 (year with a global dust storm) ; ...
     33#  =30 Mars Year 30
     34iaervar = 24
     35iaervar = 26
     36# Dust opacity at 610 Pa (when constant, i.e. for the iaervar=1 case)
     37tauvis=0.1
    3638# Dust vertical distribution:
    3739# (=0: old distrib. (Pollack90), =1: top set by "topdustref",
    … …  
    4648callrad   = .true.
    4749# call NLTE radiative schemes ?   matters only if callrad=T
    48 callnlte = .false. #.true.
     50callnlte = .true.
     51callnlte = .false.
     52# NLTE 15um scheme to use.
     53# 0-> Old scheme, static oxygen
     54# 1-> Old scheme, dynamic oxygen
     55# 2-> New scheme
     56nltemodel = 2
     57nltemodel = 0
    4958# call CO2 NIR absorption ?   matters only if callrad=T
    50 callnirco2 = .false. #.true.
     59callnirco2 = .true.
     60# NIR NLTE correction ? matters only if callnirco2=T
     61nircorr=1
     62nircorr=0
    5163# call turbulent vertical diffusion ?
    5264calldifv  = .true.
    5365# call convective adjustment ?
    5466calladj  = .true.
     67# Thermals
     68calltherm = .true.
     69callrichsl = .true.
    5570# call CO2 condensation ?
    5671callcond =.true.
    … …  
    6176# Impose polar cap surface albedos as observed by TES?
    6277TESicealbedo = .true.
     78TESicealbedo = .false.
     79## Coefficient for Northern cap albedoes
     80TESice_Ncoef=1.6
     81## Coefficient for Southern cap albedoes
     82TESice_Scoef=1.6
     83
    6384
    6485## Radiative transfer options :
    … …  
    7495## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    7596# DUST: Transported dust ? (if >0, use 'dustbin' dust bins)
    76 dustbin    = 0
     97dustbin    = 2
     98dustbin    = 0
    7799# DUST: Radiatively active dust ? (matters if dustbin>0)
     100active  = .true.
    78101active  = .false.
    79102# DUST: use mass and number mixing ratios to predict dust size ?
    80 doubleq   = .false.
     103doubleq   = .true.
     104doubleq   = .false.
    81105# DUST: lifted by GCM surface winds ?
     106lifting = .true.
    82107lifting = .false.
    83108# DUST: lifted by dust devils ?
    84109callddevil = .false.
    85 # DUST: Scavenging by CO2 snowfall ?
     110# DUST: Scavenging by H2O/CO2 snowfall ?
     111scavenging = .true.
    86112scavenging = .false.
    87113# DUST/WATERICE: Gravitationnal sedimentation ?
     114sedimentation = .true.
    88115sedimentation = .false.
    89116# WATERICE: Radiatively active transported atmospheric water ice ?
     117activice   = .true.
    90118activice   = .false.
    91119# WATER: Compute water cycle
     120water = .true.
    92121water = .false.
    93122# WATER: Microphysical scheme for water-ice clouds?
     123microphys = .true.
    94124microphys = .false.
    95125# WATER: parameter contact
    96126mteta = 0.95
     127# WATER: Effective variance for sedimentation for the log-normal
     128#        distribution of ice particles ?
     129nuice_sed=0.1
    97130# WATER: current permanent caps at both poles. True IS RECOMMENDED
    98131#        (with .true., North cap is a source of water and South pole
    99132#         is a cold trap)
    100133caps  = .true.
     134caps  = .false.
    101135# WATER: Water ice albedo ?
    102136albedo_h2o_ice = 0.35
    … …  
    105139# WATER: Frost thickness threshold for albedo
    106140frost_albedo_threshold = 0.005
    107 nuice_sed=0.1
    108 icedryness=1
    109 nircorr=0
    110 imicro = 2
    111 imicro = 1
    112141# PHOTOCHEMISTRY: include chemical species
     142photochem  = .true.
    113143photochem  = .false.
    114144
    … …  
    116146##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    117147# call thermosphere ?
     148callthermos = .true.
    118149callthermos = .false.
    119150# WATER: included without cycle (only if water=.false.)
    120151thermoswater = .false.
    121152# call thermal conduction ? (only if callthermos=.true.)
     153callconduct = .true.
    122154callconduct = .false.
    123 # call EUV heating ? (only if callthermos=.true.)
    124 calleuv=.false.
    125155# call molecular viscosity ? (only if callthermos=.true.)
     156callmolvis = .true.
    126157callmolvis = .false.
    127158# call molecular diffusion ? (only if callthermos=.true.)
     159callmoldiff = .true.
    128160callmoldiff = .false.
    129161# call thermospheric photochemistry ? (only if callthermos=.true.)
     162thermochem = .true.
    130163thermochem = .false.
     164# call EUV heating ? (only if callthermos=.true.)
     165calleuv=.true.
     166calleuv=.false.
     167#Method to include solar variability?
     168#0-> Old method (fixed EUV input)  1-> Variability with E10.7 as observed
     169solvarmod=0
    131170# date for solar flux calculation: (1985 < date < 2002)
    132 ## (Solar min=1996.4 ave=1993.4 max=1990.6)
     171## (Solar min=1996.4 ave=1993.4 max=1990.6) ; Only used if solvarmod=0
    133172solarcondate = 1993.4
    134        
     173#Solar variability as observed for MY? (must bebetween MY23 and MY30)
     174# (only matters if solvarmod=1)
     175solvaryear=24
     176# value for the UV heating efficiency
     177##(experimental values between 0.19 and 0.23, lower values may
     178## be used to compensate for low 15 um cooling)
     179euveff = 0.21       
    135180
    136 ## Thermals
    137 ## ~~~~~~~~
    138 calltherm = .true.
    139 callrichsl = .true.
     181
     182

  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars0.nothermals

    r1041 r1042  
    33#Run with or without tracer transport ?
    44tracer=.false.
     5
     6#Directory where external input files are:
     7datadir=/planeto/emlmd/run_mcd5.v0_new/datadir
     8#datadir = /san/home/spiga/DATAFILE
     9#datadir = /d5/aslmd/datafile
     10#datadir=/users/forget/WWW/datagcm/datafile
    511
    612#Diurnal cycle ?  if diurnal=False, diurnal averaged solar heating
    … …  
    2329#  =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario,
    2430#  =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file)
    25 #  =24 Mars Year 24 from TES assimilation (ie: MCD reference case)
    26 #  =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm)
    27 #  =26 Mars Year 26 from TES assimilation
    28 iaervar = 26
    29 # Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case)
    30 tauvis=0.2
     31#  =6 "cold" (low dust) scenario ; =7 "warm" (high dust) scenario
     32#  =24 Mars Year 24  ; =25 Mars Year 25 (year with a global dust storm) ; ...
     33#  =30 Mars Year 30
     34iaervar = 24
     35iaervar = 26
     36# Dust opacity at 610 Pa (when constant, i.e. for the iaervar=1 case)
     37tauvis=0.1
    3138# Dust vertical distribution:
    3239# (=0: old distrib. (Pollack90), =1: top set by "topdustref",
    … …  
    4148callrad   = .true.
    4249# call NLTE radiative schemes ?   matters only if callrad=T
     50callnlte = .true.
    4351callnlte = .false.
     52# NLTE 15um scheme to use.
     53# 0-> Old scheme, static oxygen
     54# 1-> Old scheme, dynamic oxygen
     55# 2-> New scheme
     56nltemodel = 2
     57nltemodel = 0
    4458# call CO2 NIR absorption ?   matters only if callrad=T
    4559callnirco2 = .true.
     60# NIR NLTE correction ? matters only if callnirco2=T
     61nircorr=1
     62nircorr=0
    4663# call turbulent vertical diffusion ?
    4764calldifv  = .true.
    4865# call convective adjustment ?
    4966calladj  = .true.
     67# Thermals
     68calltherm = .true.
     69callrichsl = .true.
     70calltherm = .false.
     71callrichsl = .false.
    5072# call CO2 condensation ?
    5173callcond =.true.
    … …  
    5678# Impose polar cap surface albedos as observed by TES?
    5779TESicealbedo = .true.
     80TESicealbedo = .false.
     81## Coefficient for Northern cap albedoes
     82TESice_Ncoef=1.6
     83## Coefficient for Southern cap albedoes
     84TESice_Scoef=1.6
     85
    5886
    5987## Radiative transfer options :
    … …  
    6997## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    7098# DUST: Transported dust ? (if >0, use 'dustbin' dust bins)
     99dustbin    = 2
    71100dustbin    = 0
    72101# DUST: Radiatively active dust ? (matters if dustbin>0)
     102active  = .true.
    73103active  = .false.
    74104# DUST: use mass and number mixing ratios to predict dust size ?
     105doubleq   = .true.
    75106doubleq   = .false.
    76107# DUST: lifted by GCM surface winds ?
    77 lifting = .false.
     108lifting = .true.
     109lifting = .false.
    78110# DUST: lifted by dust devils ?
    79111callddevil = .false.
    80 # DUST: Scavenging by CO2 snowfall ?
     112# DUST: Scavenging by H2O/CO2 snowfall ?
     113scavenging = .true.
    81114scavenging = .false.
    82115# DUST/WATERICE: Gravitationnal sedimentation ?
     116sedimentation = .true.
    83117sedimentation = .false.
    84118# WATERICE: Radiatively active transported atmospheric water ice ?
     119activice   = .true.
    85120activice   = .false.
    86121# WATER: Compute water cycle
     122water = .true.
    87123water = .false.
     124# WATER: Microphysical scheme for water-ice clouds?
     125microphys = .true.
     126microphys = .false.
     127# WATER: parameter contact
     128mteta = 0.95
     129# WATER: Effective variance for sedimentation for the log-normal
     130#        distribution of ice particles ?
     131nuice_sed=0.1
    88132# WATER: current permanent caps at both poles. True IS RECOMMENDED
    89133#        (with .true., North cap is a source of water and South pole
    90134#         is a cold trap)
     135caps  = .true.
    91136caps  = .false.
     137# WATER: Water ice albedo ?
     138albedo_h2o_ice = 0.35
     139# WATER: Water ice thermal inertia
     140inert_h2o_ice = 800
     141# WATER: Frost thickness threshold for albedo
     142frost_albedo_threshold = 0.005
    92143# PHOTOCHEMISTRY: include chemical species
     144photochem  = .true.
    93145photochem  = .false.
    94146
    … …  
    96148##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    97149# call thermosphere ?
     150callthermos = .true.
    98151callthermos = .false.
    99152# WATER: included without cycle (only if water=.false.)
    100153thermoswater = .false.
    101154# call thermal conduction ? (only if callthermos=.true.)
     155callconduct = .true.
    102156callconduct = .false.
    103 # call EUV heating ? (only if callthermos=.true.)
    104 calleuv=.false.
    105157# call molecular viscosity ? (only if callthermos=.true.)
     158callmolvis = .true.
    106159callmolvis = .false.
    107160# call molecular diffusion ? (only if callthermos=.true.)
     161callmoldiff = .true.
    108162callmoldiff = .false.
    109163# call thermospheric photochemistry ? (only if callthermos=.true.)
     164thermochem = .true.
    110165thermochem = .false.
     166# call EUV heating ? (only if callthermos=.true.)
     167calleuv=.true.
     168calleuv=.false.
     169#Method to include solar variability?
     170#0-> Old method (fixed EUV input)  1-> Variability with E10.7 as observed
     171solvarmod=0
    111172# date for solar flux calculation: (1985 < date < 2002)
    112 ## (Solar min=1996.4 ave=1993.4 max=1990.6)
     173## (Solar min=1996.4 ave=1993.4 max=1990.6) ; Only used if solvarmod=0
    113174solarcondate = 1993.4
    114        
     175#Solar variability as observed for MY? (must bebetween MY23 and MY30)
     176# (only matters if solvarmod=1)
     177solvaryear=24
     178# value for the UV heating efficiency
     179##(experimental values between 0.19 and 0.23, lower values may
     180## be used to compensate for low 15 um cooling)
     181euveff = 0.21       
    115182
    116183
     184

  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars11

    r76 r1042  
    33#Run with or without tracer transport ?
    44tracer=.true.
     5
     6#Directory where external input files are:
     7datadir=/planeto/emlmd/run_mcd5.v0_new/datadir
     8#datadir = /san/home/spiga/DATAFILE
     9#datadir = /d5/aslmd/datafile
     10#datadir=/users/forget/WWW/datagcm/datafile
    511
    612#Diurnal cycle ?  if diurnal=False, diurnal averaged solar heating
    … …  
    2329#  =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario,
    2430#  =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file)
    25 #  =24 Mars Year 24 from TES assimilation (ie: MCD reference case)
    26 #  =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm)
    27 #  =26 Mars Year 26 from TES assimilation
    28 iaervar = 26
    29 # Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case)
    30 tauvis=0.2
     31#  =6 "cold" (low dust) scenario ; =7 "warm" (high dust) scenario
     32#  =24 Mars Year 24  ; =25 Mars Year 25 (year with a global dust storm) ; ...
     33#  =30 Mars Year 30
     34iaervar = 24
     35iaervar = 26
     36# Dust opacity at 610 Pa (when constant, i.e. for the iaervar=1 case)
     37tauvis=0.1
    3138# Dust vertical distribution:
    3239# (=0: old distrib. (Pollack90), =1: top set by "topdustref",
    … …  
    4148callrad   = .true.
    4249# call NLTE radiative schemes ?   matters only if callrad=T
     50callnlte = .true.
    4351callnlte = .false.
     52# NLTE 15um scheme to use.
     53# 0-> Old scheme, static oxygen
     54# 1-> Old scheme, dynamic oxygen
     55# 2-> New scheme
     56nltemodel = 2
     57nltemodel = 0
    4458# call CO2 NIR absorption ?   matters only if callrad=T
    4559callnirco2 = .true.
     60# NIR NLTE correction ? matters only if callnirco2=T
     61nircorr=1
     62nircorr=0
    4663# call turbulent vertical diffusion ?
    4764calldifv  = .true.
    4865# call convective adjustment ?
    4966calladj  = .true.
     67# Thermals
     68calltherm = .true.
     69callrichsl = .true.
    5070# call CO2 condensation ?
    5171callcond =.true.
    … …  
    5676# Impose polar cap surface albedos as observed by TES?
    5777TESicealbedo = .true.
     78TESicealbedo = .false.
     79## Coefficient for Northern cap albedoes
     80TESice_Ncoef=1.6
     81## Coefficient for Southern cap albedoes
     82TESice_Scoef=1.6
     83
    5884
    5985## Radiative transfer options :
    … …  
    7197dustbin    = 2
    7298# DUST: Radiatively active dust ? (matters if dustbin>0)
    73 active  = .false.
     99active  = .true.
    74100# DUST: use mass and number mixing ratios to predict dust size ?
    75 doubleq   = .true. 
     101doubleq   = .true.
    76102# DUST: lifted by GCM surface winds ?
    77 lifting = .true.
     103lifting = .true.
     104lifting = .false.
    78105# DUST: lifted by dust devils ?
    79106callddevil = .false.
    80 # DUST: Scavenging by CO2 snowfall ?
     107# DUST: Scavenging by H2O/CO2 snowfall ?
     108scavenging = .true.
    81109scavenging = .false.
    82110# DUST/WATERICE: Gravitationnal sedimentation ?
    83111sedimentation = .true.
    84112# WATERICE: Radiatively active transported atmospheric water ice ?
     113activice   = .true.
    85114activice   = .false.
    86115# WATER: Compute water cycle
    87116water = .true.
     117# WATER: Microphysical scheme for water-ice clouds?
     118microphys = .true.
     119microphys = .false.
     120# WATER: parameter contact
     121mteta = 0.95
     122# WATER: Effective variance for sedimentation for the log-normal
     123#        distribution of ice particles ?
     124nuice_sed=0.1
    88125# WATER: current permanent caps at both poles. True IS RECOMMENDED
    89126#        (with .true., North cap is a source of water and South pole
    90127#         is a cold trap)
    91 caps  = .false. #.true.
     128caps  = .true.
     129# WATER: Water ice albedo ?
     130albedo_h2o_ice = 0.35
     131# WATER: Water ice thermal inertia
     132inert_h2o_ice = 800
     133# WATER: Frost thickness threshold for albedo
     134frost_albedo_threshold = 0.005
    92135# PHOTOCHEMISTRY: include chemical species
     136photochem  = .true.
    93137photochem  = .false.
    94138
    … …  
    96140##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    97141# call thermosphere ?
     142callthermos = .true.
    98143callthermos = .false.
    99144# WATER: included without cycle (only if water=.false.)
    100145thermoswater = .false.
    101146# call thermal conduction ? (only if callthermos=.true.)
     147callconduct = .true.
    102148callconduct = .false.
    103 # call EUV heating ? (only if callthermos=.true.)
    104 calleuv=.false.
    105149# call molecular viscosity ? (only if callthermos=.true.)
     150callmolvis = .true.
    106151callmolvis = .false.
    107152# call molecular diffusion ? (only if callthermos=.true.)
     153callmoldiff = .true.
    108154callmoldiff = .false.
    109155# call thermospheric photochemistry ? (only if callthermos=.true.)
     156thermochem = .true.
    110157thermochem = .false.
     158# call EUV heating ? (only if callthermos=.true.)
     159calleuv=.true.
     160calleuv=.false.
     161#Method to include solar variability?
     162#0-> Old method (fixed EUV input)  1-> Variability with E10.7 as observed
     163solvarmod=0
    111164# date for solar flux calculation: (1985 < date < 2002)
    112 ## (Solar min=1996.4 ave=1993.4 max=1990.6)
     165## (Solar min=1996.4 ave=1993.4 max=1990.6) ; Only used if solvarmod=0
    113166solarcondate = 1993.4
    114        
     167#Solar variability as observed for MY? (must bebetween MY23 and MY30)
     168# (only matters if solvarmod=1)
     169solvaryear=24
     170# value for the UV heating efficiency
     171##(experimental values between 0.19 and 0.23, lower values may
     172## be used to compensate for low 15 um cooling)
     173euveff = 0.21       
    115174
    116175
     176

  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars11_radact

    r81 r1042  
    11##General options
    22##~~~~~~~~~~~~~~~
    3 #Directory where external input files are
    4 datadir = /san/home/spiga/DATAFILE/
    5 
    63#Run with or without tracer transport ?
    74tracer=.true.
     5
     6#Directory where external input files are:
     7datadir=/planeto/emlmd/run_mcd5.v0_new/datadir
     8#datadir = /san/home/spiga/DATAFILE
     9#datadir = /d5/aslmd/datafile
     10#datadir=/users/forget/WWW/datagcm/datafile
    811
    912#Diurnal cycle ?  if diurnal=False, diurnal averaged solar heating
    … …  
    2629#  =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario,
    2730#  =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file)
    28 #  =24 Mars Year 24 from TES assimilation (ie: MCD reference case)
    29 #  =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm)
    30 #  =26 Mars Year 26 from TES assimilation
    31 iaervar = 26
    32 # Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case)
    33 tauvis=0.2
     31#  =6 "cold" (low dust) scenario ; =7 "warm" (high dust) scenario
     32#  =24 Mars Year 24  ; =25 Mars Year 25 (year with a global dust storm) ; ...
     33#  =30 Mars Year 30
     34iaervar = 24
     35iaervar = 26
     36# Dust opacity at 610 Pa (when constant, i.e. for the iaervar=1 case)
     37tauvis=0.1
    3438# Dust vertical distribution:
    3539# (=0: old distrib. (Pollack90), =1: top set by "topdustref",
    … …  
    4448callrad   = .true.
    4549# call NLTE radiative schemes ?   matters only if callrad=T
     50callnlte = .true.
    4651callnlte = .false.
     52# NLTE 15um scheme to use.
     53# 0-> Old scheme, static oxygen
     54# 1-> Old scheme, dynamic oxygen
     55# 2-> New scheme
     56nltemodel = 2
     57nltemodel = 0
    4758# call CO2 NIR absorption ?   matters only if callrad=T
    4859callnirco2 = .true.
     60# NIR NLTE correction ? matters only if callnirco2=T
     61nircorr=1
     62nircorr=0
    4963# call turbulent vertical diffusion ?
    5064calldifv  = .true.
    5165# call convective adjustment ?
    5266calladj  = .true.
     67# Thermals
     68calltherm = .true.
     69callrichsl = .true.
    5370# call CO2 condensation ?
    5471callcond =.true.
    … …  
    5976# Impose polar cap surface albedos as observed by TES?
    6077TESicealbedo = .true.
     78TESicealbedo = .false.
     79## Coefficient for Northern cap albedoes
     80TESice_Ncoef=1.6
     81## Coefficient for Southern cap albedoes
     82TESice_Scoef=1.6
     83
    6184
    6285## Radiative transfer options :
    … …  
    7699active  = .true.
    77100# DUST: use mass and number mixing ratios to predict dust size ?
    78 doubleq   = .true. 
     101doubleq   = .true.
    79102# DUST: lifted by GCM surface winds ?
    80 lifting = .true.
     103lifting = .true.
     104lifting = .false.
    81105# DUST: lifted by dust devils ?
    82106callddevil = .false.
    83 # DUST: Scavenging by CO2 snowfall ?
     107# DUST: Scavenging by H2O/CO2 snowfall ?
     108scavenging = .true.
    84109scavenging = .false.
    85110# DUST/WATERICE: Gravitationnal sedimentation ?
    … …  
    89114# WATER: Compute water cycle
    90115water = .true.
     116# WATER: Microphysical scheme for water-ice clouds?
     117microphys = .true.
     118microphys = .false.
     119# WATER: parameter contact
     120mteta = 0.95
     121# WATER: Effective variance for sedimentation for the log-normal
     122#        distribution of ice particles ?
     123nuice_sed=0.1
    91124# WATER: current permanent caps at both poles. True IS RECOMMENDED
    92125#        (with .true., North cap is a source of water and South pole
    93126#         is a cold trap)
    94127caps  = .true.
     128# WATER: Water ice albedo ?
     129albedo_h2o_ice = 0.35
     130# WATER: Water ice thermal inertia
     131inert_h2o_ice = 800
     132# WATER: Frost thickness threshold for albedo
     133frost_albedo_threshold = 0.005
    95134# PHOTOCHEMISTRY: include chemical species
     135photochem  = .true.
    96136photochem  = .false.
    97137
    … …  
    99139##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    100140# call thermosphere ?
     141callthermos = .true.
    101142callthermos = .false.
    102143# WATER: included without cycle (only if water=.false.)
    103144thermoswater = .false.
    104145# call thermal conduction ? (only if callthermos=.true.)
     146callconduct = .true.
    105147callconduct = .false.
    106 # call EUV heating ? (only if callthermos=.true.)
    107 calleuv=.false.
    108148# call molecular viscosity ? (only if callthermos=.true.)
     149callmolvis = .true.
    109150callmolvis = .false.
    110151# call molecular diffusion ? (only if callthermos=.true.)
     152callmoldiff = .true.
    111153callmoldiff = .false.
    112154# call thermospheric photochemistry ? (only if callthermos=.true.)
     155thermochem = .true.
    113156thermochem = .false.
     157# call EUV heating ? (only if callthermos=.true.)
     158calleuv=.true.
     159calleuv=.false.
     160#Method to include solar variability?
     161#0-> Old method (fixed EUV input)  1-> Variability with E10.7 as observed
     162solvarmod=0
    114163# date for solar flux calculation: (1985 < date < 2002)
    115 ## (Solar min=1996.4 ave=1993.4 max=1990.6)
     164## (Solar min=1996.4 ave=1993.4 max=1990.6) ; Only used if solvarmod=0
    116165solarcondate = 1993.4
    117        
     166#Solar variability as observed for MY? (must bebetween MY23 and MY30)
     167# (only matters if solvarmod=1)
     168solvaryear=24
     169# value for the UV heating efficiency
     170##(experimental values between 0.19 and 0.23, lower values may
     171## be used to compensate for low 15 um cooling)
     172euveff = 0.21       
    118173
    119174
     175

  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars2

    r76 r1042  
    33#Run with or without tracer transport ?
    44tracer=.true.
     5
     6#Directory where external input files are:
     7datadir=/planeto/emlmd/run_mcd5.v0_new/datadir
     8#datadir = /san/home/spiga/DATAFILE
     9#datadir = /d5/aslmd/datafile
     10#datadir=/users/forget/WWW/datagcm/datafile
    511
    612#Diurnal cycle ?  if diurnal=False, diurnal averaged solar heating
    … …  
    2329#  =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario,
    2430#  =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file)
    25 #  =24 Mars Year 24 from TES assimilation (ie: MCD reference case)
    26 #  =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm)
    27 #  =26 Mars Year 26 from TES assimilation
    28 iaervar = 26
    29 # Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case)
    30 tauvis=0.2
     31#  =6 "cold" (low dust) scenario ; =7 "warm" (high dust) scenario
     32#  =24 Mars Year 24  ; =25 Mars Year 25 (year with a global dust storm) ; ...
     33#  =30 Mars Year 30
     34iaervar = 24
     35iaervar = 26
     36# Dust opacity at 610 Pa (when constant, i.e. for the iaervar=1 case)
     37tauvis=0.1
    3138# Dust vertical distribution:
    3239# (=0: old distrib. (Pollack90), =1: top set by "topdustref",
    … …  
    4148callrad   = .true.
    4249# call NLTE radiative schemes ?   matters only if callrad=T
     50callnlte = .true.
    4351callnlte = .false.
     52# NLTE 15um scheme to use.
     53# 0-> Old scheme, static oxygen
     54# 1-> Old scheme, dynamic oxygen
     55# 2-> New scheme
     56nltemodel = 2
     57nltemodel = 0
    4458# call CO2 NIR absorption ?   matters only if callrad=T
    4559callnirco2 = .true.
     60# NIR NLTE correction ? matters only if callnirco2=T
     61nircorr=1
     62nircorr=0
    4663# call turbulent vertical diffusion ?
    4764calldifv  = .true.
    4865# call convective adjustment ?
    4966calladj  = .true.
     67# Thermals
     68calltherm = .true.
     69callrichsl = .true.
    5070# call CO2 condensation ?
    5171callcond =.true.
    … …  
    5676# Impose polar cap surface albedos as observed by TES?
    5777TESicealbedo = .true.
     78TESicealbedo = .false.
     79## Coefficient for Northern cap albedoes
     80TESice_Ncoef=1.6
     81## Coefficient for Southern cap albedoes
     82TESice_Scoef=1.6
     83
    5884
    5985## Radiative transfer options :
    … …  
    6995## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    7096# DUST: Transported dust ? (if >0, use 'dustbin' dust bins)
     97dustbin    = 2
    7198dustbin    = 1
    7299# DUST: Radiatively active dust ? (matters if dustbin>0)
     100active  = .true.
    73101active  = .false.
    74102# DUST: use mass and number mixing ratios to predict dust size ?
    75 doubleq   = .false.
     103doubleq   = .true.
     104doubleq   = .false.
    76105# DUST: lifted by GCM surface winds ?
    77 lifting = .true. #.false.
     106lifting = .true.
     107lifting = .false.
    78108# DUST: lifted by dust devils ?
    79109callddevil = .false.
    80 # DUST: Scavenging by CO2 snowfall ?
     110# DUST: Scavenging by H2O/CO2 snowfall ?
     111scavenging = .true.
    81112scavenging = .false.
    82113# DUST/WATERICE: Gravitationnal sedimentation ?
    83114sedimentation = .true.
    84115# WATERICE: Radiatively active transported atmospheric water ice ?
     116activice   = .true.
    85117activice   = .false.
    86118# WATER: Compute water cycle
    87 water = .false. #.true.
     119water = .true.
     120water = .false.
     121# WATER: Microphysical scheme for water-ice clouds?
     122microphys = .true.
     123microphys = .false.
     124# WATER: parameter contact
     125mteta = 0.95
     126# WATER: Effective variance for sedimentation for the log-normal
     127#        distribution of ice particles ?
     128nuice_sed=0.1
    88129# WATER: current permanent caps at both poles. True IS RECOMMENDED
    89130#        (with .true., North cap is a source of water and South pole
    90131#         is a cold trap)
    91 caps  = .false. #.true.
     132caps  = .true.
     133caps  = .false.
     134# WATER: Water ice albedo ?
     135albedo_h2o_ice = 0.35
     136# WATER: Water ice thermal inertia
     137inert_h2o_ice = 800
     138# WATER: Frost thickness threshold for albedo
     139frost_albedo_threshold = 0.005
    92140# PHOTOCHEMISTRY: include chemical species
     141photochem  = .true.
    93142photochem  = .false.
    94143
    … …  
    96145##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    97146# call thermosphere ?
     147callthermos = .true.
    98148callthermos = .false.
    99149# WATER: included without cycle (only if water=.false.)
    100150thermoswater = .false.
    101151# call thermal conduction ? (only if callthermos=.true.)
     152callconduct = .true.
    102153callconduct = .false.
    103 # call EUV heating ? (only if callthermos=.true.)
    104 calleuv=.false.
    105154# call molecular viscosity ? (only if callthermos=.true.)
     155callmolvis = .true.
    106156callmolvis = .false.
    107157# call molecular diffusion ? (only if callthermos=.true.)
     158callmoldiff = .true.
    108159callmoldiff = .false.
    109160# call thermospheric photochemistry ? (only if callthermos=.true.)
     161thermochem = .true.
    110162thermochem = .false.
     163# call EUV heating ? (only if callthermos=.true.)
     164calleuv=.true.
     165calleuv=.false.
     166#Method to include solar variability?
     167#0-> Old method (fixed EUV input)  1-> Variability with E10.7 as observed
     168solvarmod=0
    111169# date for solar flux calculation: (1985 < date < 2002)
    112 ## (Solar min=1996.4 ave=1993.4 max=1990.6)
     170## (Solar min=1996.4 ave=1993.4 max=1990.6) ; Only used if solvarmod=0
    113171solarcondate = 1993.4
    114        
     172#Solar variability as observed for MY? (must bebetween MY23 and MY30)
     173# (only matters if solvarmod=1)
     174solvaryear=24
     175# value for the UV heating efficiency
     176##(experimental values between 0.19 and 0.23, lower values may
     177## be used to compensate for low 15 um cooling)
     178euveff = 0.21       
    115179
    116180
     181

  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars3

    r88 r1042  
    11##General options
    22##~~~~~~~~~~~~~~~
    3 #Directory where external input files are
    4 datadir = /san/home/spiga/DATAFILE/
     3#Run with or without tracer transport ?
     4tracer=.false.
    55
    6 #Run with or without tracer transport ?
    7 tracer=.true.
     6#Directory where external input files are:
     7datadir=/planeto/emlmd/run_mcd5.v0_new/datadir
     8#datadir = /san/home/spiga/DATAFILE
     9#datadir = /d5/aslmd/datafile
     10#datadir=/users/forget/WWW/datagcm/datafile
    811
    912#Diurnal cycle ?  if diurnal=False, diurnal averaged solar heating
    … …  
    2629#  =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario,
    2730#  =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file)
    28 #  =24 Mars Year 24 from TES assimilation (ie: MCD reference case)
    29 #  =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm)
    30 #  =26 Mars Year 26 from TES assimilation
    31 iaervar = 26
    32 # Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case)
    33 tauvis=0.2
     31#  =6 "cold" (low dust) scenario ; =7 "warm" (high dust) scenario
     32#  =24 Mars Year 24  ; =25 Mars Year 25 (year with a global dust storm) ; ...
     33#  =30 Mars Year 30
     34iaervar = 24
     35iaervar = 26
     36# Dust opacity at 610 Pa (when constant, i.e. for the iaervar=1 case)
     37tauvis=0.1
    3438# Dust vertical distribution:
    3539# (=0: old distrib. (Pollack90), =1: top set by "topdustref",
    … …  
    4448callrad   = .true.
    4549# call NLTE radiative schemes ?   matters only if callrad=T
     50callnlte = .true.
    4651callnlte = .false.
     52# NLTE 15um scheme to use.
     53# 0-> Old scheme, static oxygen
     54# 1-> Old scheme, dynamic oxygen
     55# 2-> New scheme
     56nltemodel = 2
     57nltemodel = 0
    4758# call CO2 NIR absorption ?   matters only if callrad=T
    4859callnirco2 = .true.
     60# NIR NLTE correction ? matters only if callnirco2=T
     61nircorr=1
     62nircorr=0
    4963# call turbulent vertical diffusion ?
    5064calldifv  = .true.
    5165# call convective adjustment ?
    5266calladj  = .true.
     67# Thermals
     68calltherm = .true.
     69callrichsl = .true.
    5370# call CO2 condensation ?
    5471callcond =.true.
    … …  
    5976# Impose polar cap surface albedos as observed by TES?
    6077TESicealbedo = .true.
     78TESicealbedo = .false.
     79## Coefficient for Northern cap albedoes
     80TESice_Ncoef=1.6
     81## Coefficient for Southern cap albedoes
     82TESice_Scoef=1.6
     83
    6184
    6285## Radiative transfer options :
    … …  
    7497dustbin    = 2
    7598# DUST: Radiatively active dust ? (matters if dustbin>0)
    76 active  = .true. #.false.
     99active  = .true.
    77100# DUST: use mass and number mixing ratios to predict dust size ?
    78 doubleq   = .true. #.false.
     101doubleq   = .true.
    79102# DUST: lifted by GCM surface winds ?
    80 lifting = .true. #.false.
     103lifting = .true.
     104lifting = .false.
    81105# DUST: lifted by dust devils ?
    82106callddevil = .false.
    83 # DUST: Scavenging by CO2 snowfall ?
     107# DUST: Scavenging by H2O/CO2 snowfall ?
     108scavenging = .true.
    84109scavenging = .false.
    85110# DUST/WATERICE: Gravitationnal sedimentation ?
    86111sedimentation = .true.
    87112# WATERICE: Radiatively active transported atmospheric water ice ?
     113activice   = .true.
    88114activice   = .false.
    89115# WATER: Compute water cycle
    90 water = .false. #.true.
     116water = .true.
     117water = .false.
     118# WATER: Microphysical scheme for water-ice clouds?
     119microphys = .true.
     120microphys = .false.
     121# WATER: parameter contact
     122mteta = 0.95
     123# WATER: Effective variance for sedimentation for the log-normal
     124#        distribution of ice particles ?
     125nuice_sed=0.1
    91126# WATER: current permanent caps at both poles. True IS RECOMMENDED
    92127#        (with .true., North cap is a source of water and South pole
    93128#         is a cold trap)
    94 caps  = .false. #.true.
     129caps  = .true.
     130caps  = .false.
     131# WATER: Water ice albedo ?
     132albedo_h2o_ice = 0.35
     133# WATER: Water ice thermal inertia
     134inert_h2o_ice = 800
     135# WATER: Frost thickness threshold for albedo
     136frost_albedo_threshold = 0.005
    95137# PHOTOCHEMISTRY: include chemical species
     138photochem  = .true.
    96139photochem  = .false.
    97140
    … …  
    99142##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    100143# call thermosphere ?
     144callthermos = .true.
    101145callthermos = .false.
    102146# WATER: included without cycle (only if water=.false.)
    103147thermoswater = .false.
    104148# call thermal conduction ? (only if callthermos=.true.)
     149callconduct = .true.
    105150callconduct = .false.
    106 # call EUV heating ? (only if callthermos=.true.)
    107 calleuv=.false.
    108151# call molecular viscosity ? (only if callthermos=.true.)
     152callmolvis = .true.
    109153callmolvis = .false.
    110154# call molecular diffusion ? (only if callthermos=.true.)
     155callmoldiff = .true.
    111156callmoldiff = .false.
    112157# call thermospheric photochemistry ? (only if callthermos=.true.)
     158thermochem = .true.
    113159thermochem = .false.
     160# call EUV heating ? (only if callthermos=.true.)
     161calleuv=.true.
     162calleuv=.false.
     163#Method to include solar variability?
     164#0-> Old method (fixed EUV input)  1-> Variability with E10.7 as observed
     165solvarmod=0
    114166# date for solar flux calculation: (1985 < date < 2002)
    115 ## (Solar min=1996.4 ave=1993.4 max=1990.6)
     167## (Solar min=1996.4 ave=1993.4 max=1990.6) ; Only used if solvarmod=0
    116168solarcondate = 1993.4
    117        
     169#Solar variability as observed for MY? (must bebetween MY23 and MY30)
     170# (only matters if solvarmod=1)
     171solvaryear=24
     172# value for the UV heating efficiency
     173##(experimental values between 0.19 and 0.23, lower values may
     174## be used to compensate for low 15 um cooling)
     175euveff = 0.21       
    118176
    119177
     178
Note: See TracChangeset for help on using the changeset viewer.