Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/callphys.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/callphys.def (revision 752) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/callphys.def (revision 752) @@ -0,0 +1,122 @@ +## Thermals Options +##~~~~~~~~~~~~~~~~~ +nircorr=0 +# Activate thermals : +calltherm = .false. +callrichsl = .true. +yamada4_flag=.true. + +##General options +##~~~~~~~~~~~~~~~ +#Run with or without tracer transport ? +tracer=.true. + +#Diurnal cycle ? if diurnal=False, diurnal averaged solar heating +diurnal=.true. + +#Seasonal cycle ? if season=False, Ls stays constant, to value set in "start" +season = .true. + +#write some more output on the screen ? +lwrite = .false. + +#Save statistics in file "stats.nc" ? +stats =.true. + +#Save EOF profiles in file "profiles" for Climate Database? +calleofdump = .false. + +## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +## ~~~~~~~~~~~~~ +# =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario, +# =4 Mars Year 24 from TES assimilation (same as =24 for now) +# =24 Mars Year 24 from TES assimilation (ie: MCD reference case) +# =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm) +# =26 Mars Year 26 from TES assimilation +iaervar = 1 +# Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case) +tauvis=0.5 +# Dust vertical distribution: +# (=0: old distrib. (Pollack90), =1: top set by "topdustref", +# =2: Viking scenario; =3 MGS scenario) +iddist = 1 +# Dust top altitude (km). (Matters only if iddist=1) +topdustref = 55. + +## Physical Parameterizations : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call radiative transfer ? +callrad = .true. +# call NLTE radiative schemes ? matters only if callrad=T +callnlte = .false. +# call CO2 NIR absorption ? matters only if callrad=T +callnirco2 = .false. +# call turbulent vertical diffusion ? +calldifv = .true. +# call convective adjustment ? +calladj = .false. +# call CO2 condensation ? +callcond =.true. +# call thermal conduction in the soil ? +callsoil = .true. +# call Lott's gravity wave/subgrid topography scheme ?' +calllott = .true. + +## Radiative transfer options : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# the rad.transfer is computed every "iradia" physical timestep +iradia = 1 +# Output of the exchange coefficient mattrix ? for diagnostic only +callg2d = .false. +# Rayleigh scattering : (should be .false. for now) +rayleigh = .false. + +## Tracer (dust water, ice and/or chemical species) options (used if tracer=T): +## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# DUST: Transported dust ? (if >0, use 'dustbin' dust bins) +dustbin = 0 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .false. +# DUST: use mass and number mixing ratios to predict dust size ? +doubleq = .false. +# DUST: lifted by GCM surface winds ? +lifting = .false. +# DUST: lifted by dust devils ? +callddevil = .false. +# DUST: Scavenging by CO2 snowfall ? +scavenging = .false. +# DUST/WATERICE: Gravitationnal sedimentation ? +sedimentation = .false. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .false. +# WATER: Compute water cycle +water = .false. +# WATER: current permanent caps at both poles. True IS RECOMMENDED +# (with .true., North cap is a source of water and South pole +# is a cold trap) +caps = .true. +# PHOTOCHEMISTRY: include chemical species +photochem = .false. + +## Thermospheric options (relevant if tracer=T) : +##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call thermosphere ? +callthermos = .false. +# WATER: included without cycle (only if water=.false.) +thermoswater = .false. +# call thermal conduction ? (only if callthermos=.true.) +callconduct = .false. +# call EUV heating ? (only if callthermos=.true.) +calleuv=.false. +# call molecular viscosity ? (only if callthermos=.true.) +callmolvis = .false. +# call molecular diffusion ? (only if callthermos=.true.) +callmoldiff = .false. +# call thermospheric photochemistry ? (only if callthermos=.true.) +thermochem = .false. +# date for solar flux calculation: (1985 < date < 2002) +## (Solar min=1996.4 ave=1993.4 max=1990.6) +solarcondate = 1993.4 + + + Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/callphysd2_.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/callphysd2_.def (revision 752) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/callphysd2_.def (revision 752) @@ -0,0 +1,122 @@ +## Thermals Options +##~~~~~~~~~~~~~~~~~ +nircorr=0 +# Activate thermals : +calltherm = .false. +callrichsl = .true. +yamada4_flag=.true. + +##General options +##~~~~~~~~~~~~~~~ +#Run with or without tracer transport ? +tracer=.true. + +#Diurnal cycle ? if diurnal=False, diurnal averaged solar heating +diurnal=.true. + +#Seasonal cycle ? if season=False, Ls stays constant, to value set in "start" +season = .true. + +#write some more output on the screen ? +lwrite = .false. + +#Save statistics in file "stats.nc" ? +stats =.true. + +#Save EOF profiles in file "profiles" for Climate Database? +calleofdump = .false. + +## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +## ~~~~~~~~~~~~~ +# =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario, +# =4 Mars Year 24 from TES assimilation (same as =24 for now) +# =24 Mars Year 24 from TES assimilation (ie: MCD reference case) +# =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm) +# =26 Mars Year 26 from TES assimilation +iaervar = 1 +# Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case) +tauvis=0.5 +# Dust vertical distribution: +# (=0: old distrib. (Pollack90), =1: top set by "topdustref", +# =2: Viking scenario; =3 MGS scenario) +iddist = 1 +# Dust top altitude (km). (Matters only if iddist=1) +topdustref = 55. + +## Physical Parameterizations : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call radiative transfer ? +callrad = .true. +# call NLTE radiative schemes ? matters only if callrad=T +callnlte = .false. +# call CO2 NIR absorption ? matters only if callrad=T +callnirco2 = .false. +# call turbulent vertical diffusion ? +calldifv = .true. +# call convective adjustment ? +calladj = .false. +# call CO2 condensation ? +callcond =.true. +# call thermal conduction in the soil ? +callsoil = .true. +# call Lott's gravity wave/subgrid topography scheme ?' +calllott = .true. + +## Radiative transfer options : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# the rad.transfer is computed every "iradia" physical timestep +iradia = 1 +# Output of the exchange coefficient mattrix ? for diagnostic only +callg2d = .false. +# Rayleigh scattering : (should be .false. for now) +rayleigh = .false. + +## Tracer (dust water, ice and/or chemical species) options (used if tracer=T): +## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# DUST: Transported dust ? (if >0, use 'dustbin' dust bins) +dustbin = 0 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .false. +# DUST: use mass and number mixing ratios to predict dust size ? +doubleq = .false. +# DUST: lifted by GCM surface winds ? +lifting = .false. +# DUST: lifted by dust devils ? +callddevil = .false. +# DUST: Scavenging by CO2 snowfall ? +scavenging = .false. +# DUST/WATERICE: Gravitationnal sedimentation ? +sedimentation = .false. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .false. +# WATER: Compute water cycle +water = .false. +# WATER: current permanent caps at both poles. True IS RECOMMENDED +# (with .true., North cap is a source of water and South pole +# is a cold trap) +caps = .true. +# PHOTOCHEMISTRY: include chemical species +photochem = .false. + +## Thermospheric options (relevant if tracer=T) : +##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call thermosphere ? +callthermos = .false. +# WATER: included without cycle (only if water=.false.) +thermoswater = .false. +# call thermal conduction ? (only if callthermos=.true.) +callconduct = .false. +# call EUV heating ? (only if callthermos=.true.) +calleuv=.false. +# call molecular viscosity ? (only if callthermos=.true.) +callmolvis = .false. +# call molecular diffusion ? (only if callthermos=.true.) +callmoldiff = .false. +# call thermospheric photochemistry ? (only if callthermos=.true.) +thermochem = .false. +# date for solar flux calculation: (1985 < date < 2002) +## (Solar min=1996.4 ave=1993.4 max=1990.6) +solarcondate = 1993.4 + + + Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/callphysd3_.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/callphysd3_.def (revision 752) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/callphysd3_.def (revision 752) @@ -0,0 +1,122 @@ +## Thermals Options +##~~~~~~~~~~~~~~~~~ +nircorr=0 +# Activate thermals : +calltherm = .false. +callrichsl = .true. +yamada4_flag=.true. + +##General options +##~~~~~~~~~~~~~~~ +#Run with or without tracer transport ? +tracer=.true. + +#Diurnal cycle ? if diurnal=False, diurnal averaged solar heating +diurnal=.true. + +#Seasonal cycle ? if season=False, Ls stays constant, to value set in "start" +season = .true. + +#write some more output on the screen ? +lwrite = .false. + +#Save statistics in file "stats.nc" ? +stats =.true. + +#Save EOF profiles in file "profiles" for Climate Database? +calleofdump = .false. + +## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +## ~~~~~~~~~~~~~ +# =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario, +# =4 Mars Year 24 from TES assimilation (same as =24 for now) +# =24 Mars Year 24 from TES assimilation (ie: MCD reference case) +# =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm) +# =26 Mars Year 26 from TES assimilation +iaervar = 1 +# Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case) +tauvis=0.5 +# Dust vertical distribution: +# (=0: old distrib. (Pollack90), =1: top set by "topdustref", +# =2: Viking scenario; =3 MGS scenario) +iddist = 1 +# Dust top altitude (km). (Matters only if iddist=1) +topdustref = 55. + +## Physical Parameterizations : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call radiative transfer ? +callrad = .true. +# call NLTE radiative schemes ? matters only if callrad=T +callnlte = .false. +# call CO2 NIR absorption ? matters only if callrad=T +callnirco2 = .false. +# call turbulent vertical diffusion ? +calldifv = .true. +# call convective adjustment ? +calladj = .false. +# call CO2 condensation ? +callcond =.true. +# call thermal conduction in the soil ? +callsoil = .true. +# call Lott's gravity wave/subgrid topography scheme ?' +calllott = .true. + +## Radiative transfer options : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# the rad.transfer is computed every "iradia" physical timestep +iradia = 1 +# Output of the exchange coefficient mattrix ? for diagnostic only +callg2d = .false. +# Rayleigh scattering : (should be .false. for now) +rayleigh = .false. + +## Tracer (dust water, ice and/or chemical species) options (used if tracer=T): +## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# DUST: Transported dust ? (if >0, use 'dustbin' dust bins) +dustbin = 0 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .false. +# DUST: use mass and number mixing ratios to predict dust size ? +doubleq = .false. +# DUST: lifted by GCM surface winds ? +lifting = .false. +# DUST: lifted by dust devils ? +callddevil = .false. +# DUST: Scavenging by CO2 snowfall ? +scavenging = .false. +# DUST/WATERICE: Gravitationnal sedimentation ? +sedimentation = .false. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .false. +# WATER: Compute water cycle +water = .false. +# WATER: current permanent caps at both poles. True IS RECOMMENDED +# (with .true., North cap is a source of water and South pole +# is a cold trap) +caps = .true. +# PHOTOCHEMISTRY: include chemical species +photochem = .false. + +## Thermospheric options (relevant if tracer=T) : +##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call thermosphere ? +callthermos = .false. +# WATER: included without cycle (only if water=.false.) +thermoswater = .false. +# call thermal conduction ? (only if callthermos=.true.) +callconduct = .false. +# call EUV heating ? (only if callthermos=.true.) +calleuv=.false. +# call molecular viscosity ? (only if callthermos=.true.) +callmolvis = .false. +# call molecular diffusion ? (only if callthermos=.true.) +callmoldiff = .false. +# call thermospheric photochemistry ? (only if callthermos=.true.) +thermochem = .false. +# date for solar flux calculation: (1985 < date < 2002) +## (Solar min=1996.4 ave=1993.4 max=1990.6) +solarcondate = 1993.4 + + + Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/namelist.input =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/namelist.input (revision 752) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/namelist.input (revision 752) @@ -0,0 +1,72 @@ + &time_control + start_year = 2024, 2024, 2024 + start_month = 09, 09, 09 + start_day = 05, 05, 05 + start_hour = 00, 00, 00 + start_minute = 00, 00, 00 + start_second = 00, 00, 00 + end_year = 2024, 2024, 2024 + end_month = 09, 09, 09 + end_day = 10, 10, 10 + end_hour = 00, 00, 00 + history_interval_s = 1850, 1850, 1850 + frames_per_outfile = 48, 48, 48 + restart = .false. + restart_interval = 360 + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + !!!! OPT + interval_seconds = 3700 !Freq of GCM update + + / + + &domains + time_step = 40 + dx = 135000, 45000, 15000, + dy = 135000, 45000, 15000, + e_we = 57, 61, 61 + e_sn = 57, 61, 61 + e_vert = 61, 61, 61 + p_top_requested = 5 + max_dom = 3 + grid_id = 1, 2, 3, + parent_id = 0, 1, 2, + i_parent_start = 0, 19, 20, + j_parent_start = 0, 19, 25, + parent_grid_ratio = 1, 3, 3, + parent_time_step_ratio = 1, 3, 3, ! ok in non_hydrostatic/hydrostatic with time_step = 20, and feedback=0 ! essayer 1,2,2 sans feedback pour voir ? + feedback = 1 !voir aussi smooth_option et blend_width + smooth_option = 2 ! =0 nan au bout de quelques heures +/ + + &physics + radt = 1, ! call each half hour with radt = 1850/dt + mars = 0, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + num_soil_layers = 18 !!! NOUVELLE PHYSIQUE + / + + &fdda + / + + &dynamics + diff_6th_factor = 0.2, 0.2, 0.2 + non_hydrostatic = T, T, T + / + + &grib2 + / + + &bdy_control + specified = T, F, F + nested = F, T, T + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/namelist.wps =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/namelist.wps (revision 752) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/namelist.wps (revision 752) @@ -0,0 +1,38 @@ +&share + wrf_core = 'ARW', !! [do not modify: choice of dynamical core] + max_dom = 3 + start_date = '2024-09-05_00:00:00','2024-09-05_00:00:00','2024-09-05_00:00:00' !! YYYY-MM-DD_HH:mm:ss start date + end_date = '2024-09-10_00:00:00','2024-09-05_00:00:00','2024-09-05_00:00:00' !! YYYY-MM-DD_HH:mm:ss end date + interval_seconds = 3700 !Freq of GCM update + io_form_geogrid = 2, !! [do not modify: choice of NETCDF outputs] + debug_level = 0, !! verbose level of the programs + opt_output_from_geogrid_path='./' !! location of the geogrid outputs +/ + + +&geogrid + parent_id = 0, 1, 2, + parent_grid_ratio = 1, 3, 3, + i_parent_start = 0, 19, 20, + j_parent_start = 0, 19, 25, + e_we = 57, 61, 61 + e_sn = 57, 61, 61 + dx = 135000, + dy = 135000, + geog_data_res = '64ppd', '64ppd', '64ppd' !! choice of static data sources + !! NB: possible: '64ppd', '32ppd', ... + !! NB: please glance at geogrid/GEOGRID.TBL + map_proj = 'mercator', !! map projection: 'mercator', 'lambert' or 'polar' + ref_lat = -2, !! north latitude of the center of the domain + ref_lon = -6, !! east longitude of the center of the domain + truelat1 = -2, !! (lambert or polar) lat position of projection cone + truelat2 = 0.0, !! [do not modify] + stand_lon = -6, !! (lambert or polar) lon position of projection cone + geog_data_path = './WPS_GEOG', !! [do not modify: symbolic link in the WPS folder] +/ + +&metgrid + fg_name = './WPSFEED/LMD' !! [do not modify: symbolic link in the WPS folder] + io_form_metgrid = 2, !! [do not modify: choice of NETCDF outputs] + opt_output_from_metgrid_path='./WRFFEED/current' !! [do not modify: symbolic link] +/ Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/run.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/run.def (revision 752) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_EXOMARS_twowaynests/run.def (revision 752) @@ -0,0 +1,110 @@ +# +#----------------------------------------------------------------------- +#Parametres de controle du run: +#------------------------------ +# Nombre de jours d'integration + nday=1 + +# nombre de pas par jour (multiple de iperiod) ( ici pour dt = 1 min ) + day_step = 480 + +# periode pour le pas Matsuno (en pas) + iperiod=5 + +# periode de sortie des variables de controle (en pas) + iconser=120 + +# periode d'ecriture du fichier histoire (en jour) + iecri=100 + +# periode de stockage fichier histmoy (en jour) + periodav=60. + +# periode de la dissipation (en pas) + idissip=1 + +# choix de l'operateur de dissipation (star ou non star ) + lstardis=.true. + +# avec ou sans coordonnee hybrides + hybrid=.true. + +# nombre d'iterations de l'operateur de dissipation gradiv +nitergdiv=1 + +# nombre d'iterations de l'operateur de dissipation nxgradrot +nitergrot=2 + +# nombre d'iterations de l'operateur de dissipation divgrad + niterh=2 + +# temps de dissipation des plus petites long.d ondes pour u,v (gradiv) + tetagdiv= 4000. + +# temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot) + tetagrot= 5000. + +# temps de dissipation des plus petites long.d ondes pour h ( divgrad) + tetatemp= 5000. + +# coefficient pour gamdissip + coefdis=0. + +# choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog) + purmats=.false. + +# avec ou sans physique + physic=.true. + +# periode de la physique (en pas) + iphysiq= 10 +# choix d'une grille reguliere + grireg=.true. + +# frequence (en pas) de l'ecriture du fichier diagfi + ecritphy= 60 +# ecritphy= 60 +# longitude en degres du centre du zoom + clon=63. + +# latitude en degres du centre du zoom + clat=0. + +# facteur de grossissement du zoom,selon longitude + grossismx=1. + +# facteur de grossissement du zoom ,selon latitude + grossismy=1. + +# Fonction f(y) hyperbolique si = .true. , sinon sinusoidale + fxyhypb=.false. + +# extension en longitude de la zone du zoom ( fraction de la zone totale) + dzoomx= 0. + +# extension en latitude de la zone du zoom ( fraction de la zone totale) + dzoomy=0. + +# raideur du zoom en X + taux=2. + +# raideur du zoom en Y + tauy=2. + +# Fonction f(y) avec y = Sin(latit.) si = .TRUE. , Sinon y = latit. + ysinus= .false. + +# Avec sponge layer + callsponge = .true. + +# Sponge: mode0(u=v=0), mode1(u=umoy,v=0), mode2(u=umoy,v=vmoy) + mode_sponge= 2 + +# Sponge: hauteur de sponge (km) + hsponge= 130 + +# Sponge: tetasponge (secondes) + tetasponge = 50000 + +# some definitions for the physics, in file 'callphys.def' +INCLUDEDEF=callphys.def