Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/DUST_STORM_bak/callphys.def_dust
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/DUST_STORM_bak/callphys.def_dust	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/DUST_STORM_bak/callphys.def_dust	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+T
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+4        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+T
+callcond  (call CO2 condensation ?)
+T
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+F
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+75 
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+1 0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+T F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+T F
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+F
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+F
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/DUST_STORM_bak/namelist.input_dust
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/DUST_STORM_bak/namelist.input_dust	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/DUST_STORM_bak/namelist.input_dust	(revision 348)
@@ -0,0 +1,66 @@
+ &time_control
+ start_year       = 2024,
+ start_month      = 07, 
+ start_day        = 01,
+ start_hour       = 06,
+ end_year         = 2024,
+ end_month        = 07,
+ end_day          = 02, 
+ end_hour         = 06, 
+ history_interval    = 37,   ! *100s : 74 is each 2 MH, 37 is each MH  
+ frames_per_outfile  = 24, 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ interval_seconds = 7400  !! OPT OPT
+ /
+
+ &domains
+ time_step   = 25,
+ dx = 20000,
+ dy = 20000,
+ e_we   = 51,                ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   = 51,                ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert = 61,                ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 5
+ /
+
+ &physics
+ radt = 1,	! call each half hour with radt = 1850/dt
+ mars = 2,      ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF)
+ init_TI = 100.,	!! OPT OPT
+ init_AL = 0.27, 	!! OPT OPT
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ diff_6th_opt = 2,         !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES
+ diff_6th_factor = 0.2,    !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES
+ diff_opt = 1,
+ km_opt = 4,
+ pd_scalar = F,		!! OPT OPT
+ non_hydrostatic = F	!! OPT OPT
+ !!time_step_sound = 10
+ /
+
+ &grib2
+ /
+
+ &bdy_control
+ spec_bdy_width                      = 3,     !! OPT OPT
+ spec_zone                           = 1,     !! OPT OPT  1+1 est instable
+ relax_zone                          = 2,     !! OPT OPT
+ /
+
+!! DO NOT MODIFY
+ &namelist_quilt
+ nio_tasks_per_group = 0,
+ nio_groups = 1,
+ /
+!! DO NOT MODIFY
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/DUST_STORM_bak/stress.def
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/DUST_STORM_bak/stress.def	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/DUST_STORM_bak/stress.def	(revision 348)
@@ -0,0 +1,2 @@
+0.01
+1.
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUb_namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUb_namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUb_namelist.input	(revision 348)
@@ -0,0 +1,54 @@
+ &time_control
+ start_year       = 2024,
+ start_month      = 09, 
+ start_day        = 23,
+ start_hour       = 00,
+ end_year         = 2024,
+ end_month        = 09,
+ end_day          = 26, 
+ end_hour         = 00, 
+ history_interval    = 74,   ! *100s : 74 is each 2 MH, 37 is each MH  
+ frames_per_outfile  = 12, 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ /
+
+ &domains
+ time_step   = 10,  !5
+ dx = 2500,
+ dy = 2500,
+ e_we   = 150,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   = 150,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert = 31,                ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 50
+ /
+
+ &physics
+ radt = 2, !5	! call each half hour with radt = 1850/dt
+ mars = 0,      ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ pd_scalar = F   !! no water cycle (mars = 0)
+ /
+
+ &grib2
+ /
+
+ &bdy_control
+ /
+
+!! DO NOT MODIFY
+ &namelist_quilt
+ nio_tasks_per_group = 0,
+ nio_groups = 1,
+ /
+!! DO NOT MODIFY
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUb_namelist.wps
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUb_namelist.wps	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUb_namelist.wps	(revision 348)
@@ -0,0 +1,36 @@
+&share
+ wrf_core = 'ARW',
+ max_dom = 1,
+ start_date = '2024-09-23_00:00:00', '2026-01-47_00:00:00'
+ end_date   = '2024-09-26_00:00:00', '2026-01-47_00:00:00'
+ interval_seconds = 3700
+ io_form_geogrid = 2,
+ debug_level = 0,
+ opt_output_from_geogrid_path='./'
+/
+
+
+&geogrid
+ parent_id         =   1, 1
+ parent_grid_ratio =   1, 3
+ i_parent_start    =   1, 12
+ j_parent_start    =   1, 12
+ e_we              =  150  !275, 25
+ e_sn              =  150  !275, 25
+ geog_data_res     =  '64ppd_noHRti'     ! 'custom', 'custom'	!'32ppd','32ppd'
+ dx = 2500, !1500,
+ dy = 2500, !1500,
+ map_proj = 'mercator',
+ ref_lat   =  11 !14.76,
+ ref_lon   =  159.34,
+ truelat1  =  0.0,
+ truelat2  =  0.0,
+ stand_lon =  00.00,
+ geog_data_path = './WPS_GEOG',
+/
+
+&metgrid
+ fg_name = './LMD/LMD'
+ io_form_metgrid = 2,
+ opt_output_from_metgrid_path='/tmp15/aslmd/wrf/MET_EM'
+/
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUs_finer_namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUs_finer_namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUs_finer_namelist.input	(revision 348)
@@ -0,0 +1,57 @@
+ &time_control
+ start_year       = 2024,
+ start_month      = 09, 
+ start_day        = 23,
+ start_hour       = 00,
+ end_year         = 2024,
+ end_month        = 09,
+ end_day          = 26, 
+ end_hour         = 00, 
+ history_interval    = 74,   ! *100s : 74 is each 2 MH, 37 is each MH  
+ frames_per_outfile  = 12, 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ /
+
+ &domains
+!!! time_step   = 1,    ! 2 crashe, meme en baissant le top, 1 aussi au bout d'un moment
+ time_step   = 0
+ time_step_fract_num                 = 1,
+ time_step_fract_den                 = 2,
+ dx = 500, !1000,
+ dy = 500, !1000,
+ e_we   = 181,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   = 181,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert = 31,                ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 50
+ /
+
+ &physics
+ radt = 20, !5	! call each half hour with radt = 1850/dt
+ mars = 0,      ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ pd_scalar = F   !! no water cycle (mars = 0)
+ /
+
+ &grib2
+ /
+
+ &bdy_control
+ /
+
+!! DO NOT MODIFY
+ &namelist_quilt
+ nio_tasks_per_group = 0,
+ nio_groups = 1,
+ /
+!! DO NOT MODIFY
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUs_finer_namelist.wps
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUs_finer_namelist.wps	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUs_finer_namelist.wps	(revision 348)
@@ -0,0 +1,38 @@
+&share
+ wrf_core = 'ARW',
+ max_dom = 1,
+ start_date = '2024-09-23_00:00:00', '2026-01-47_00:00:00'
+ end_date   = '2024-09-26_00:00:00', '2026-01-47_00:00:00'
+ interval_seconds = 3700
+ io_form_geogrid = 2,
+ debug_level = 0,
+ opt_output_from_geogrid_path='./'
+/
+
+
+&geogrid
+ parent_id         =   1, 1
+ parent_grid_ratio =   1, 3
+ i_parent_start    =   1, 12
+ j_parent_start    =   1, 12
+ e_we              =  181 !150 !150  !275, 25
+ e_sn              =  181 !150 !150  !275, 25
+ geog_data_res     =  '64ppd_noHRti'     ! 'custom', 'custom'	!'32ppd','32ppd'
+ dx = 500, !1500,
+ dy = 500, !1500,
+ map_proj = 'mercator',
+! ref_lat   =  13.4  !13.04,    
+! ref_lon   =  158.6 !158.8,
+ ref_lat   =  13.04,
+ ref_lon   =  158.95,
+ truelat1  =  0.0,
+ truelat2  =  0.0,
+ stand_lon =  00.00,
+ geog_data_path = './WPS_GEOG',
+/
+
+&metgrid
+ fg_name = './WPSFEED/LMD'
+ io_form_metgrid = 2,
+ opt_output_from_metgrid_path='./WRFFEED'
+/
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUs_namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUs_namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUs_namelist.input	(revision 348)
@@ -0,0 +1,54 @@
+ &time_control
+ start_year       = 2024,
+ start_month      = 09, 
+ start_day        = 23,
+ start_hour       = 00,
+ end_year         = 2024,
+ end_month        = 09,
+ end_day          = 26, 
+ end_hour         = 00, 
+ history_interval    = 74,   ! *100s : 74 is each 2 MH, 37 is each MH  
+ frames_per_outfile  = 12, 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ /
+
+ &domains
+ time_step   = 5
+ dx = 1000,
+ dy = 1000,
+ e_we   = 150,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   = 150,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert = 31,                ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 50
+ /
+
+ &physics
+ radt = 4, !5	! call each half hour with radt = 1850/dt
+ mars = 0,      ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ pd_scalar = F   !! no water cycle (mars = 0)
+ /
+
+ &grib2
+ /
+
+ &bdy_control
+ /
+
+!! DO NOT MODIFY
+ &namelist_quilt
+ nio_tasks_per_group = 0,
+ nio_groups = 1,
+ /
+!! DO NOT MODIFY
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUs_namelist.wps
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUs_namelist.wps	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/DUs_namelist.wps	(revision 348)
@@ -0,0 +1,38 @@
+&share
+ wrf_core = 'ARW',
+ max_dom = 1,
+ start_date = '2024-09-23_00:00:00', '2026-01-47_00:00:00'
+ end_date   = '2024-09-26_00:00:00', '2026-01-47_00:00:00'
+ interval_seconds = 3700
+ io_form_geogrid = 2,
+ debug_level = 0,
+ opt_output_from_geogrid_path='./'
+/
+
+
+&geogrid
+ parent_id         =   1, 1
+ parent_grid_ratio =   1, 3
+ i_parent_start    =   1, 12
+ j_parent_start    =   1, 12
+ e_we              =  150 !150 !150  !275, 25
+ e_sn              =  150 !150 !150  !275, 25
+ geog_data_res     =  '64ppd_noHRti'     ! 'custom', 'custom'	!'32ppd','32ppd'
+ dx = 1000, !1500,
+ dy = 1000, !1500,
+ map_proj = 'mercator',
+ ref_lat   =  13.4  !13.04,    
+ ref_lon   =  158.6 !158.8,
+! ref_lat   =  -4.87,
+! ref_lon   =  170.14,
+ truelat1  =  0.0,
+ truelat2  =  0.0,
+ stand_lon =  00.00,
+ geog_data_path = './WPS_GEOG',
+/
+
+&metgrid
+ fg_name = './WPSFEED/LMD'
+ io_form_metgrid = 2,
+ opt_output_from_metgrid_path='./WRFFEED'
+/
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_3NEST_callphys.def
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_3NEST_callphys.def	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_3NEST_callphys.def	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+F
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+4        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+F
+callcond  (call CO2 condensation ?)
+T
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+F
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+20 
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+F
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+F
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+F
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_3NEST_callphysd2_.def
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_3NEST_callphysd2_.def	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_3NEST_callphysd2_.def	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+F
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+4        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+F
+callcond  (call CO2 condensation ?)
+T
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+F
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+40 
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+F
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+F
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+F
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_3NEST_callphysd3_.def
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_3NEST_callphysd3_.def	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_3NEST_callphysd3_.def	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+F
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+4        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+F
+callcond  (call CO2 condensation ?)
+T
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+F
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+80 
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+F
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+F
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+F
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_3NEST_launch_simu
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_3NEST_launch_simu	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_3NEST_launch_simu	(revision 348)
@@ -0,0 +1,141 @@
+#! /bin/bash
+
+
+#
+# USE: launch_simu &
+#
+
+#
+# must set env. var. 	WRF_current 
+#			
+
+# PRELIMINARY STEP
+# - model must be compiled and linked (wrf.exe, real.exe)
+# - domain must be correctly set with geogrid.exe
+# - put namelist.input, namelist.wps, callphys*.def
+
+# 'ls' must yield:
+# callphys*.def  launch_simu  namelist.input  namelist.wps  real.exe  wrf.exe
+#
+
+
+###############################################################
+gcmstep=0  #1
+wpsstep=0  #1
+realstep=1
+wrfstep=1
+###############################################################
+# GCM STEP
+   daygcm=371
+# WPS STEP
+   name='em2'     # un nom de dossier non pris dans WRFFEED
+# WRF STEP
+   proc=4
+###############################################################
+# FOLDERS
+   wheregcm=$GCM
+   whereprep=$WRF_current/PREP_MARS
+   wherewps=$WRF_current/WPS
+###############################################################
+
+
+#
+# LOG FILE HEADER
+#
+
+where=$PWD
+echo "check what's happening in " ${where}/log_launch
+
+rm -rf ${where}/log_launch 2> /dev/null
+echo SIMU folder is ${where} > ${where}/log_launch
+echo SIMU WILL BE LAUNCHED on the machine >> ${where}/log_launch
+hostname >> ${where}/log_launch
+
+#
+# 0. GCM step
+#
+if [ ${gcmstep} -eq 1 ]
+then
+
+   echo "0. GCM run and output preparation for WPS" >> ${where}/log_launch
+
+   cd ${wheregcm}
+   ln -sf ../STARTBASE_64_48_25_t2/start/start${daygcm}.nc start.nc
+   ln -sf ../STARTBASE_64_48_25_t2/startfi/startfi${daygcm}.nc startfi.nc
+
+   rm -rf restart.nc
+   rm -rf restartfi.nc
+   rm -rf log_gcm
+
+   echo running the GCM - check in ${wheregcm}/log_gcm >> ${where}/log_launch
+   nohup gcm.e > log_gcm
+
+   echo preparing outputs for WPS >> ${where}/log_launch
+   cd ${whereprep}
+   compile_and_exec >> ${where}/log_launch
+
+fi
+
+#
+# 1. WPS step
+#
+if [ ${wpsstep} -eq 1 ]
+then
+
+   echo "1. WPS: horizontal interpolations of static and meteorological data" >> ${where}/log_launch
+
+#   rm -rf ${where}/met_em 
+   mkdir ${where}/met_em
+   cd ${where}/met_em
+   ln -sf ../namelist.wps .
+
+   cd ${wherewps}
+   rm -rf WHERE_SIMU_IS
+   ln -sf ${where}/met_em WHERE_SIMU_IS
+   mkdir ./WRFFEED/current
+   nohup geogrid.exe >> ${where}/log_launch
+   nohup metgrid.exe >> ${where}/log_launch 
+   mv ./WRFFEED/current ./WRFFEED/${name}
+
+   ln -sf $FEED/WRFFEED/${name}/* ${where}/met_em/
+
+fi
+
+#
+# 2. REAL.EXE
+#
+if [ ${realstep} -eq 1 ]
+then
+
+   echo "2. REAL.EXE: vertical interpolation and preparation of boundary conditions" >> ${where}/log_launch
+
+   cd ${where}/met_em
+   rm -rf rsl.* 2> /dev/null
+   rm -rf wrfbdy* 2> /dev/null
+   rm -rf wrfinput* 2> /dev/null
+
+   ln -sf ../namelist.input .
+   ln -sf ../real.exe .
+   real.exe 
+   cd .. 
+   ln -sf ./met_em/wrf*d0* .
+
+fi
+
+
+#
+# 3. WRF.EXE
+#
+if [ ${wrfstep} -eq 1 ]
+then
+
+   echo "3. WRF.EXE: launch the mesoscale simulation" >> ${where}/log_launch
+
+   rm -rf rsl.*
+   rm -rf wrfout*
+
+   $WHERE_MPI/mpirun -np ${proc} wrf.exe < /dev/null &
+   echo OK, simulation launched !
+
+fi
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_3NEST_namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_3NEST_namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_3NEST_namelist.input	(revision 348)
@@ -0,0 +1,138 @@
+ &time_control
+ start_year       = 2024, 2024, 2024,
+ start_month      =   07,   07,   07,
+ start_day        =   01,   01,   01,
+ start_hour       =   01,   01,   01,
+ end_year         = 2024, 2024, 2024,
+ end_month        =   07,   07,   07,
+ end_day          =   04,   04,   04,
+ end_hour         =   01,   01,   01,
+ history_interval   = 37,   37,   37,    
+ frames_per_outfile = 24,   24,   24, 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ input_from_file  =    T,    T,    T,     ! F en second pour que d1 init. d2 (smoother)
+ /
+
+ &domains
+ time_step   = 30 
+ dx     = 30000, 10000,  3333.333, 
+ dy     = 30000, 10000,  3333.333, 
+ e_we   =   153,   157,   157,            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   =   153,   157,   157,            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert =    61,    61,    61,            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 3
+!!!!! NESTING STUFF
+ max_dom = 3	!!!! PLEASE DO NOT FORGET callphysd*_.def	
+ grid_id                = 1,   2,   3,
+ parent_id              = 0,   1,   2,
+ i_parent_start         = 0,  51,  60,		! attention au raccord avec WPS !
+ j_parent_start         = 0,  55,  45,		! attention au raccord avec WPS !
+ parent_grid_ratio      = 1,   3,   3,
+ parent_time_step_ratio = 1,   3,   3,		
+ feedback               = 0  
+ smooth_option          = 2
+ blend_width            = 5
+!!!!! NESTING STUFF
+! p_top_requested = 0.5
+! eta_levels =       1.00000,
+!     0.999500,
+!     0.998000,
+!     0.995000,
+!     0.985000,
+!     0.970000,
+!     0.940000,
+!     0.900000,
+!     0.800000,
+!     0.700000,
+!     0.592538,
+!     0.501572,
+!     0.424571,
+!     0.359392,
+!     0.304219,
+!     0.257516,
+!     0.217982,
+!     0.184518,
+!     0.156191,
+!     0.132213,
+!     0.111916,
+!    0.0947347,
+!    0.0801912,
+!    0.0678804,    
+!    0.0574595,
+!    0.0486384,
+!    0.0411715,
+!    0.0348509,
+!    0.0295007,
+!    0.0249718,
+!    0.0211382,
+!    0.0178931,
+!    0.0151462,
+!    0.0128209,
+!    0.0108527,
+!   0.00918661,
+!   0.00777630,
+!   0.00658249,
+!   0.00557196,
+!   0.00471656,
+!   0.00399248,
+!   0.00337956,
+!   0.00286074,
+!   0.00242156,
+!   0.00204981,
+!   0.00173513,
+!   0.00146875,
+!   0.00124327,
+!   0.00105241,
+!  0.000890844,
+!  0.000754083,
+!  0.000638317,
+!  0.000540324,
+!  0.000457374,
+!  0.000387159,
+!  0.000327723,
+!  0.000277412,
+!  0.000234824,
+!  0.000198774,
+!  0.000168259,
+!      0.00000
+ /
+
+ &physics
+ radt =  1, 1, 1     ! 1 is the best (change iradia in callphys)
+ mars =  0, 0, 0     ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF)
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ diff_6th_opt =      2,   2,   2,     !! NECESSAIRE 
+ diff_6th_factor =   0.2, 0.2, 0.2    !! NECESSAIRE 
+ diff_opt = 1,
+ km_opt = 4,
+ pd_scalar = F   !! no water cycle (mars = 0)
+ time_step_sound =    10,  10,  10
+ !!damp_opt = 1
+ !!dampcoef = 0.1, 0.1, 0.1
+ /
+
+ &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/REF_ARTICLE/EM_3NEST_namelist.wps
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_3NEST_namelist.wps	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_3NEST_namelist.wps	(revision 348)
@@ -0,0 +1,35 @@
+&share
+ wrf_core = 'ARW',
+ max_dom = 3,
+ start_date = '2024-07-01_01:00:00', '2024-07-01_01:00:00', '2024-07-01_01:00:00'
+ end_date   = '2024-07-04_01:00:00', '2024-07-01_01:00:00', '2024-07-01_01:00:00'
+ interval_seconds = 3700
+ io_form_geogrid = 2,
+ debug_level = 0,
+ opt_output_from_geogrid_path='./'
+/
+
+&geogrid
+ parent_id         =        0,           1,          2  
+ parent_grid_ratio =        1,           3,          3
+ i_parent_start    =        0,          51,         60 
+ j_parent_start    =        0,          55,         45
+ e_we              =      153,         157,        157      !!! ATTENTION
+ e_sn              =      153,         157,        157      !!! ATTENTION
+ geog_data_res     =  '32ppd_HRalb',     '32ppd_HRalb',    '32ppd_HRalb'      
+ dx = 30000
+ dy = 30000
+ map_proj = 'mercator',
+ ref_lat   =  24.66,
+ ref_lon   =  146.76,
+ truelat1  =  0.0,
+ truelat2  =  0.0,
+ stand_lon =  00.00,
+ geog_data_path = './WPS_GEOG',
+/
+
+&metgrid
+ fg_name = './WPSFEED/LMD'
+ io_form_metgrid = 2,
+ opt_output_from_metgrid_path='./WRFFEED/current'
+/
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_namelist.input	(revision 348)
@@ -0,0 +1,59 @@
+ &time_control
+ start_year       = 2024,
+ start_month      = 07, 
+ start_day        = 01,
+ start_hour       = 01,
+ end_year         = 2024,
+ end_month        = 07,
+ end_day          = 05, 
+ end_hour         = 01, 
+ history_interval    = 37    ! *100s : 74 is each 2 MH, 37 is each MH  
+ frames_per_outfile  = 24, 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ /
+
+ &domains
+ time_step   = 10, 
+ dx = 7000,
+ dy = 7000,
+ e_we   = 201,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   = 201,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert = 61,                ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 7
+ /
+
+ &physics
+ radt =  1,      ! 1 is the best (change iradia in callphys)
+ mars =  0,      ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF)
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ non_hydrostatic = F
+ diff_6th_opt =      2,      !! NECESSAIRE
+ diff_6th_factor =   0.2,    !! NECESSAIRE
+ diff_opt = 1,
+ km_opt = 4,
+ pd_scalar = F   !! no water cycle (mars = 0)
+ /
+
+ &grib2
+ /
+
+ &bdy_control
+ /
+
+!! DO NOT MODIFY
+ &namelist_quilt
+ nio_tasks_per_group = 0,
+ nio_groups = 1,
+ /
+!! DO NOT MODIFY
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_namelist.wps
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_namelist.wps	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/EM_namelist.wps	(revision 348)
@@ -0,0 +1,36 @@
+&share
+ wrf_core = 'ARW',
+ max_dom = 1,
+ start_date = '2024-07-01_01:00:00'
+ end_date = '2024-07-05_01:00:00'
+ interval_seconds = 3700
+ io_form_geogrid = 2,
+ debug_level = 0,
+ opt_output_from_geogrid_path='./'
+/
+
+
+&geogrid
+ parent_id         =   1,
+ parent_grid_ratio =   1,
+ i_parent_start    =   1,
+ j_parent_start    =   1,
+ e_we              =  201,
+ e_sn              =  201,
+ geog_data_res     = '64ppd_noHRti'
+ dx = 7000,
+ dy = 7000,
+ map_proj = 'mercator',
+ ref_lat   =  24.8,
+ ref_lon   =  146.9,
+ truelat1  =  0.0,
+ truelat2  =  0.0,
+ stand_lon =  00.00,
+ geog_data_path = './WPS_GEOG',
+/
+
+&metgrid
+ fg_name = './WPSFEED/LMD'
+ io_form_metgrid = 2,
+ opt_output_from_metgrid_path='./WRFFEED/current'
+/
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/GCM_namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/GCM_namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/GCM_namelist.input	(revision 348)
@@ -0,0 +1,78 @@
+ &time_control
+ start_year       = 2024, 2024
+ start_month      =   07,   04 
+ start_day        =   01,   01
+ start_hour       =   01,   06
+ end_year         = 2024, 2024
+ end_month        =   07,   04
+ end_day          =   11,   05
+ end_hour         =   00,   06
+ history_interval   = 37,   74   ! *100s : 74 is each 2 MH, 37 is each MH  !! attention, compatibilite avec time_step 
+ frames_per_outfile = 24,   12 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ input_from_file  = T, F     ! F en second pour que d1 init. d2 (smoother)
+ /
+
+ &domains
+ time_step   = 185 
+ max_dz = 4000
+ dx = 296461.63, 20000
+ dy = 296461.63, 20000 
+ e_we   =  74, 211            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   =  50, 211            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert =  26,  61            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 1.
+!!!!!! NESTING STUFF
+! max_dom = 2	
+! grid_id                = 1,  2
+! parent_id              = 0,  1
+! i_parent_start         = 0, 69      ! attention au raccord avec WPS !
+! j_parent_start         = 0, 69      ! attention au raccord avec WPS !
+! parent_grid_ratio      = 1,  3 
+! parent_time_step_ratio = 1,  5  !3  !4 !5
+! feedback               = 1
+! smooth_option          = 2
+! blend_width            = 5
+!!!!!! NESTING STUFF
+ /
+
+ &physics
+ radt =  1, !10    ! call each half hour with radt = 1850/dt
+ mars =  1, 1   ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF)
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ non_hydrostatic = F, T
+ diff_6th_opt = 2,	   !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES !! OPTION 2 SINON PB QH2O
+ diff_6th_factor = 0.5,    !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES !! FAIBLE SINON PB QH2O_ICE (si pas BDY)
+ !!diff_opt = 1,  !! non sinon crash
+ !!km_opt = 4,
+ smdiv = 0.,
+ emdiv = 0.,
+ epssm = 0.,
+ /
+
+ &grib2
+ /
+
+ &bdy_control
+ specified = T, F
+ periodic_x = T, F
+ nested    = F, 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/REF_ARTICLE/GCM_namelist.wps
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/GCM_namelist.wps	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/GCM_namelist.wps	(revision 348)
@@ -0,0 +1,40 @@
+&share
+ wrf_core = 'ARW',
+ max_dom = 1   !!2,
+ start_date = '2024-10-01_01:00:00', '2024-10-01_01:00:00'
+ end_date   = '2024-10-11_00:00:00', '2024-10-01_01:00:00'
+! start_date = '2024-01-01_01:00:00', '2024-10-01_01:00:00'
+! end_date   = '2024-01-11_00:00:00', '2024-10-01_01:00:00'
+! start_date = '2024-07-01_01:00:00', '2024-10-01_01:00:00'
+! end_date   = '2024-07-11_00:00:00', '2024-10-01_01:00:00'
+ interval_seconds = 3700
+ io_form_geogrid = 2,
+ debug_level = 0,
+ opt_output_from_geogrid_path='./'
+/
+
+
+&geogrid
+ parent_id         =   1,    1
+ parent_grid_ratio =   1,    3     
+ i_parent_start    =   1,   69    
+ j_parent_start    =   1,   69    
+ e_we              =  74,  211  !!! ATTENTION
+ e_sn              =  50,  211  !!! ATTENTION
+ geog_data_res     = '32ppd','64ppd_noHRti'
+ dx = 296461.63   !!1deglat 59292.325
+ dy = 296461.63   !!1deglat 59292.325
+ map_proj = 'mercator',
+ ref_lat   =  0.00,
+ ref_lon   =  175.,
+ truelat1  =  0.0,
+ truelat2  =  0.0,
+ stand_lon =  0.0,
+ geog_data_path = './WPS_GEOG',
+/
+
+&metgrid
+ fg_name = './WPSFEED/LMD'
+ io_form_metgrid = 2,
+ opt_output_from_metgrid_path='./WRFFEED'
+/
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/HE_namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/HE_namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/HE_namelist.input	(revision 348)
@@ -0,0 +1,59 @@
+ &time_control
+ start_year       = 2025,
+ start_month      = 07, 
+ start_day        = 02,
+ start_hour       = 00
+ end_year         = 2025,
+ end_month        = 07,
+ end_day          = 03 !!15, 
+ end_hour         = 00, 
+ history_interval    = 1 !!74,   ! *100s : 74 is each 2 MH, 37 is each MH  
+ frames_per_outfile  = 12, 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ interval_seconds = 7400     ! assim assim
+ /
+
+ &domains
+ time_step   = 18 
+ dx = 50000,
+ dy = 50000,
+ e_we   = 150,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   = 150,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert = 31,                ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 2
+ max_dz = 2000
+ /
+
+ &physics
+ radt = 2,	! call each half hour with radt = 1850/dt
+ mars = 0,      ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF)
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ pd_scalar = F   !! no water cycle (mars = 0)      
+ h_mom_adv_order  = 3      !! for better stability
+ h_sca_adv_order  = 3      !! for better stability 
+ /
+
+ &grib2
+ /
+
+ &bdy_control
+ /
+
+!! DO NOT MODIFY
+ &namelist_quilt
+ nio_tasks_per_group = 0,
+ nio_groups = 1,
+ /
+!! DO NOT MODIFY
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/HE_namelist.wps
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/HE_namelist.wps	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/HE_namelist.wps	(revision 348)
@@ -0,0 +1,37 @@
+&share
+ wrf_core = 'ARW',
+ max_dom = 2,
+ start_date = '2025-07-02_00:00:00'  '2026-01-47_00:00:00'
+! end_date   = '2025-07-15_00:00:00', '2026-01-47_00:00:00'
+ end_date   = '2025-07-03_00:00:00'
+ interval_seconds = 7400
+ io_form_geogrid = 2,
+ debug_level = 0,
+ opt_output_from_geogrid_path='./'
+/
+
+
+&geogrid
+ parent_id         =   1, 1
+ parent_grid_ratio =   1, 3
+ i_parent_start    =   1, 40
+ j_parent_start    =   1, 40
+ e_we              =  150, 181
+ e_sn              =  150, 181
+ geog_data_res     =  '32ppd','32ppd'
+ dx = 50000,
+ dy = 50000,
+ map_proj = 'mercator',
+ ref_lat   =  -10.,
+ ref_lon   =  100.,
+ truelat1  =  0.0,
+ truelat2  =  0.0,
+ stand_lon =  00.00,
+ geog_data_path = './WPS_GEOG',
+/
+
+&metgrid
+ fg_name = './WPSFEED/OXF'
+ io_form_metgrid = 2,
+ opt_output_from_metgrid_path='./WRFFEED'
+/
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/LES_namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/LES_namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/LES_namelist.input	(revision 348)
@@ -0,0 +1,63 @@
+ &time_control
+ start_year       = 2024,
+ start_month      = 01, 
+ start_day        = 05,
+ start_hour       = 21,
+ end_year         = 2024,
+ end_month        = 01,
+ end_day          = 06, 
+ end_hour         = 06, 
+ history_interval    = 1,   ! *100s : 74 is each 2 MH, 37 is each MH  
+ frames_per_outfile  = 37, 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ /
+
+ &domains
+ time_step   = 0
+ time_step_fract_num                 = 1,
+ time_step_fract_den                 = 2,
+ dx = 100,
+ dy = 100,
+ e_we   = 150,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   = 150,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert = 71,                ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 250
+ force_sfc_in_vinterp = 5,
+ /
+
+ &physics
+ radt = 60,	! call each half hour with radt = 1850/dt
+ mars = 0,      ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ init_MU=0.01	! ... quasi no-wind
+ init_MV=0.01	! ... quasi no-wind
+ init_LES=T	! LES MODE
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ pd_scalar = F   !! no water cycle (mars = 0)
+ /
+
+ &grib2
+ /
+
+ &bdy_control
+ specified = F
+ periodic_x = T 
+ periodic_y = 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/REF_ARTICLE/LES_namelist.wps
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/LES_namelist.wps	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/LES_namelist.wps	(revision 348)
@@ -0,0 +1,36 @@
+&share
+ wrf_core = 'ARW',
+ max_dom = 1,
+ start_date = '2024-01-05_18:00:00'
+ end_date   = '2024-01-06_06:00:00'
+ interval_seconds = 3700
+ io_form_geogrid = 2,
+ debug_level = 0,
+ opt_output_from_geogrid_path='./'
+/
+
+
+&geogrid
+ parent_id         =   1,
+ parent_grid_ratio =   1,
+ i_parent_start    =   1,
+ j_parent_start    =   1,
+ e_we              =  150,
+ e_sn              =  150,
+ geog_data_res     = '32ppd'
+ dx = 100,
+ dy = 100,
+ map_proj = 'mercator',
+ ref_lat   =  -14.5  !!-14.57,
+ ref_lon   =  175.4  !!175.47,
+ truelat1  =  0.0,
+ truelat2  =  0.0,
+ stand_lon =  00.00,
+ geog_data_path = './WPS_GEOG',
+/
+
+&metgrid
+ fg_name = './WPSFEED/LMD'
+ io_form_metgrid = 2,
+ opt_output_from_metgrid_path='./WRFFEED/current'
+/
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/OMG_namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/OMG_namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/OMG_namelist.input	(revision 348)
@@ -0,0 +1,75 @@
+ &time_control
+ start_year       = 2024, 2024, 2024,
+ start_month      =   01,   01,   01,
+ start_day        =   49,   49,   49,
+ start_hour       =   06,   06,   06,
+ end_year         = 2024, 2024, 2024,
+ end_month        =   01,   01,   01,
+ end_day          =   54,   54,   54,
+ end_hour         =   06,   06,   06,
+ history_interval   = 74,   74,   74,    
+ frames_per_outfile = 12,   12,   12, 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ input_from_file  =    T,    T,    T,     ! F en second pour que d1 init. d2 (smoother)
+ /
+
+ &domains
+ time_step   = 40 
+ dx     = 36000,  9000,  2250,   !!alt 48000 12000 3000
+ dy     = 36000,  9000,  2250, 
+ e_we   =   153,   157,   157,            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   =   153,   157,   157,            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert =    61,    61,    61,            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 3
+!!!!! NESTING STUFF
+ max_dom = 3	!!!! PLEASE DO NOT FORGET callphysd*_.def	
+ grid_id                = 1,   2,   3,
+ parent_id              = 0,   1,   2,
+ i_parent_start         = 0,  57,  57,      ! attention au raccord avec WPS !
+ j_parent_start         = 0,  57,  57,      ! attention au raccord avec WPS !
+ parent_grid_ratio      = 1,   4,   4,
+ parent_time_step_ratio = 1,   2,   2,
+ feedback               = 0  
+ smooth_option          = 2
+ blend_width            = 5
+!!!!! NESTING STUFF
+ /
+
+ &physics
+ radt =  1, 1, 1     ! 1 is the best (change iradia in callphys)
+ mars =  0, 0, 0     ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF)
+!!! init_MU=3.   ! ... quasi no-wind
+!!! init_MV=3.   ! ... quasi no-wind
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ diff_6th_opt =      2,   2,   2,     !! NECESSAIRE 
+ diff_6th_factor =   0.2, 0.2, 0.2    !! NECESSAIRE 
+ diff_opt = 1,
+ km_opt = 4,
+ pd_scalar = F   !! no water cycle (mars = 0)
+ /
+
+ &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/REF_ARTICLE/OMG_namelist.wps
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/OMG_namelist.wps	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/OMG_namelist.wps	(revision 348)
@@ -0,0 +1,35 @@
+&share
+ wrf_core = 'ARW',
+ max_dom = 3,
+ start_date = '2024-01-49_06:00:00', '2024-01-49_06:00:00', '2024-01-49_06:00:00'
+ end_date   = '2024-01-54_06:00:00', '2024-01-49_06:00:00', '2024-01-49_06:00:00'
+ interval_seconds = 3700
+ io_form_geogrid = 2,
+ debug_level = 0,
+ opt_output_from_geogrid_path='./'
+/
+
+&geogrid
+ parent_id         =        0,           1,          2  
+ parent_grid_ratio =        1,           4,          4
+ i_parent_start    =        0,          57,         57 
+ j_parent_start    =        0,          57,         57
+ e_we              =      153,         157,        157      !!! ATTENTION
+ e_sn              =      153,         157,        157      !!! ATTENTION
+ geog_data_res     =  '64ppd',     '64ppd',    '64ppd'      
+ dx = 36000
+ dy = 36000
+ map_proj = 'mercator',
+ ref_lat   =  11.5,
+ ref_lon   =  15.5,
+ truelat1  =  0.0,
+ truelat2  =  0.0,
+ stand_lon =  00.00,
+ geog_data_path = './WPS_GEOG',
+/
+
+&metgrid
+ fg_name = './WPSFEED/LMD'
+ io_form_metgrid = 2,
+ opt_output_from_metgrid_path='./WRFFEED/current'
+/
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/PATH_namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/PATH_namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/PATH_namelist.input	(revision 348)
@@ -0,0 +1,58 @@
+ &time_control
+ start_year       = 2024,
+ start_month      = 06, 
+ start_day        = 02,
+ start_hour       = 00,
+ end_year         = 2024,
+ end_month        = 06,
+ end_day          = 07, 
+ end_hour         = 00, 
+ history_interval    = 37,   ! *100s : 74 is each 2 MH, 37 is each MH  
+ frames_per_outfile  = 24, 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ /
+
+ &domains
+ time_step   = 24,
+ time_step_fract_num = 2,
+ time_step_fract_den = 3,
+ dx = 20000,
+ dy = 20000,
+ e_we   = 181,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   = 181,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert = 61,                ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 5
+ /
+
+ &physics
+ radt = 1,	! call each half hour with radt = 1850/dt
+ mars = 0,      ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ pd_scalar = F   !! no water cycle (mars = 0)
+ diff_6th_opt = 2,         !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES
+ diff_6th_factor = 0.2,    !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES
+ /
+
+ &grib2
+ /
+
+ &bdy_control
+ /
+
+!! DO NOT MODIFY
+ &namelist_quilt
+ nio_tasks_per_group = 0,
+ nio_groups = 1,
+ /
+!! DO NOT MODIFY
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/PATH_namelist.wps
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/PATH_namelist.wps	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/PATH_namelist.wps	(revision 348)
@@ -0,0 +1,36 @@
+&share
+ wrf_core = 'ARW',
+ max_dom = 1,
+ start_date = '2024-06-02_00:00:00', '2026-01-47_00:00:00'
+ end_date   = '2024-06-07_00:00:00', '2026-01-47_00:00:00'
+ interval_seconds = 3700
+ io_form_geogrid = 2,
+ debug_level = 0,
+ opt_output_from_geogrid_path='./'
+/
+
+
+&geogrid
+ parent_id         =   1, 1
+ parent_grid_ratio =   1, 3
+ i_parent_start    =   1, 12
+ j_parent_start    =   1, 12
+ e_we              =  181, 25
+ e_sn              =  181, 25
+ geog_data_res     =  '64ppd_noHRti'     !'32ppd','32ppd'
+ dx = 20000,
+ dy = 20000, 
+ map_proj = 'lambert' !!'mercator',
+ ref_lat   =  32. !!19.0996,
+ ref_lon   =  -40. !!-33.2526,
+ truelat1  =  30.,
+ truelat2  =  0.0,
+ stand_lon =  -33.,
+ geog_data_path = './WPS_GEOG',
+/
+
+&metgrid
+ fg_name = './WPSFEED/LMD'
+ io_form_metgrid = 2,
+ opt_output_from_metgrid_path='./WRFFEED'
+/
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/PO/callphys.def
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/PO/callphys.def	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/PO/callphys.def	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+F
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+4        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+T
+callcond  (call CO2 condensation ?)
+T
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+F
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+50 
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+F
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+F
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+F
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/PO/namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/PO/namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/PO/namelist.input	(revision 348)
@@ -0,0 +1,60 @@
+ &time_control
+ start_year       = 2024,
+ start_month      = 04, 
+ start_day        = 01,
+ start_hour       = 02,
+ end_year         = 2024,
+ end_month        = 04,
+ end_day          = 12, 
+ end_hour         = 00, 
+ history_interval    = 74,   ! *100s : 74 is each 2 MH, 37 is each MH  
+ frames_per_outfile  = 12, 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ interval_seconds = 7400
+ /
+
+ &domains
+ time_step   = 37, 
+ dx = 40000,
+ dy = 40000,
+ e_we   = 201,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   = 201,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert = 61,                ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 10
+ /
+
+ &physics
+ radt =  1,      ! 1 is the best (change iradia in callphys)
+ mars =  0,      ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF)
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ non_hydrostatic = F
+ diff_6th_opt =      2,      !! NECESSAIRE
+ diff_6th_factor =   0.2,    !! NECESSAIRE
+ diff_opt = 1,
+ km_opt = 4,
+ pd_scalar = F   !! no water cycle (mars = 0)
+ /
+
+ &grib2
+ /
+
+ &bdy_control
+ /
+
+!! DO NOT MODIFY
+ &namelist_quilt
+ nio_tasks_per_group = 0,
+ nio_groups = 1,
+ /
+!! DO NOT MODIFY
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/PO/namelist.wps
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/PO/namelist.wps	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/PO/namelist.wps	(revision 348)
@@ -0,0 +1,36 @@
+&share
+ wrf_core = 'ARW',
+ max_dom = 1,
+ start_date = '2024-04-01_02:00:00'
+ end_date   = '2024-04-12_00:00:00'
+ interval_seconds = 7400
+ io_form_geogrid = 2,
+ debug_level = 0,
+ opt_output_from_geogrid_path='./'
+/
+
+
+&geogrid
+ parent_id         =   1, 
+ parent_grid_ratio =   1, 
+ i_parent_start    =   1, 
+ j_parent_start    =   1, 
+ e_we              =  201,
+ e_sn              =  201,
+ geog_data_res     = '32ppd'
+ dx = 40000
+ dy = 40000
+ map_proj = 'polar',
+ ref_lat   =  -90,
+ ref_lon   =  0,
+ truelat1  =  -90.0,
+ truelat2  =  0.0,
+ stand_lon =  180.00,
+ geog_data_path = './WPS_GEOG',
+/
+
+&metgrid
+ fg_name = './WPSFEED/LMD'
+ io_form_metgrid = 2,
+ opt_output_from_metgrid_path='./WRFFEED/current'
+/
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TESTARSIA_namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TESTARSIA_namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TESTARSIA_namelist.input	(revision 348)
@@ -0,0 +1,66 @@
+ &time_control
+ start_year       = 2024,
+ start_month      = 07, 
+ start_day        = 01,
+ start_hour       = 06,
+ end_year         = 2024,
+ end_month        = 07,
+ end_day          = 02, 
+ end_hour         = 06, 
+ history_interval    = 37,   ! *100s : 74 is each 2 MH, 37 is each MH  
+ frames_per_outfile  = 24, 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ interval_seconds = 7400  !! OPT OPT
+ /
+
+ &domains
+ time_step   = 25,
+ dx = 20000,
+ dy = 20000,
+ e_we   = 51,                ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   = 51,                ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert = 61,                ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 5
+ /
+
+ &physics
+ radt = 1,	! call each half hour with radt = 1850/dt
+ mars = 0,      ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF)
+ init_TI = 100.,	!! OPT OPT
+ init_AL = 0.27, 	!! OPT OPT
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ diff_6th_opt = 2,         !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES
+ diff_6th_factor = 0.2,    !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES
+ diff_opt = 1,
+ km_opt = 4,
+ pd_scalar = F,		!! OPT OPT
+ non_hydrostatic = F	!! OPT OPT
+ !!time_step_sound = 10
+ /
+
+ &grib2
+ /
+
+ &bdy_control
+ spec_bdy_width                      = 3,     !! OPT OPT
+ spec_zone                           = 1,     !! OPT OPT  1+1 est instable
+ relax_zone                          = 2,     !! OPT OPT
+ /
+
+!! DO NOT MODIFY
+ &namelist_quilt
+ nio_tasks_per_group = 0,
+ nio_groups = 1,
+ /
+!! DO NOT MODIFY
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TESTARSIA_namelist.wps
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TESTARSIA_namelist.wps	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TESTARSIA_namelist.wps	(revision 348)
@@ -0,0 +1,36 @@
+&share
+ wrf_core = 'ARW',
+ max_dom = 1,
+ start_date = '2024-07-01_06:00:00'
+ end_date   = '2024-07-02_06:00:00'
+ interval_seconds = 7400
+ io_form_geogrid = 2,
+ debug_level = 0,
+ opt_output_from_geogrid_path='./'
+/
+
+
+&geogrid
+ parent_id         =   1,
+ parent_grid_ratio =   1,
+ i_parent_start    =   1,
+ j_parent_start    =   1,
+ e_we              =  51,
+ e_sn              =  51,
+ geog_data_res     = '64ppd_noHRti'
+ dx = 20000,
+ dy = 20000,
+ map_proj = 'mercator',
+ ref_lat   =  -9.5, 
+ ref_lon   =  239.5,
+ truelat1  =  0.0,
+ truelat2  =  0.0,
+ stand_lon =  0.0,
+ geog_data_path = './WPS_GEOG',
+/
+
+&metgrid
+ fg_name = './WPSFEED/LMD'
+ io_form_metgrid = 2,
+ opt_output_from_metgrid_path='./WRFFEED/current'
+/
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TH_ARSIA_namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TH_ARSIA_namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TH_ARSIA_namelist.input	(revision 348)
@@ -0,0 +1,73 @@
+ &time_control
+ start_year       = 2024, 2024
+ start_month      =   05,   05,  !!06,   06 
+ start_day        =   01,   01,  !!37,   37
+ start_hour       =   06,   06
+ end_year         = 2024, 2024
+ end_month        =   05,   05,  !!06,   06
+ end_day          =   05,   05,  !!43,   43
+ end_hour         =   06,   06
+ history_interval   = 74,   74   ! *100s : 74 is each 2 MH, 37 is each MH  !! attention, compatibilite avec time_step 
+ frames_per_outfile = 12,   12 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ input_from_file  = T, T     ! F en second pour que d1 init. d2 (smoother)
+ /
+ fine_input_stream = 2, 2
+
+ &domains
+ time_step   = 20 
+ dx = 20000, 5000
+ dy = 20000, 5000 
+ e_we   = 205, 209            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   = 205, 209            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert =  61,  61            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 3
+!!!!! NESTING STUFF
+ max_dom = 2	!! PLEASE DO NOT FORGET callphysd2_.def	
+ grid_id                = 1,  2
+ parent_id              = 0,  1
+ i_parent_start         = 0, 75      ! attention au raccord avec WPS !
+ j_parent_start         = 0, 33      ! attention au raccord avec WPS !
+ parent_grid_ratio      = 1,  4 
+ parent_time_step_ratio = 1,  2  !!4
+ feedback               = 0  !!1
+ smooth_option          = 2
+ blend_width            = 5
+!!!!! NESTING STUFF
+ /
+
+ &physics
+ radt =  1, 1     ! call each half hour with radt = 1850/dt
+ mars =  1, 1     ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF)
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ diff_6th_opt =      2,   2,       !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES
+ diff_6th_factor =   0.2, 0.2,     !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES
+ diff_opt = 1,
+ km_opt = 4,
+ /
+
+ &grib2
+ /
+
+ &bdy_control
+ specified = T, F
+ nested    = F, 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/REF_ARTICLE/TH_ARSIA_namelist.wps
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TH_ARSIA_namelist.wps	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TH_ARSIA_namelist.wps	(revision 348)
@@ -0,0 +1,36 @@
+&share
+ wrf_core = 'ARW',
+ max_dom = 2,
+ start_date = '2024-05-01_06:00:00','2024-05-01_06:00:00'
+ end_date   = '2024-05-05_06:00:00','2024-05-01_06:00:00'  
+ interval_seconds = 3700
+ io_form_geogrid = 2,
+ debug_level = 0,
+ opt_output_from_geogrid_path='./'
+/
+
+
+&geogrid
+ parent_id         =   1,    1
+ parent_grid_ratio =   0,    4  
+ i_parent_start    =   0,   75    
+ j_parent_start    =   0,   33    
+ e_we              =  205, 209   !!! ATTENTION
+ e_sn              =  205, 209   !!! ATTENTION
+ geog_data_res     = '32ppd_HRalb','32ppd_HRalb'
+ dx = 20000,
+ dy = 20000,
+ map_proj = 'mercator',
+ ref_lat   =  5.00,
+ ref_lon   =  -120.00,
+ truelat1  =  0.0,
+ truelat2  =  0.0,
+ stand_lon =  0.0,
+ geog_data_path = './WPS_GEOG',
+/
+
+&metgrid
+ fg_name = './WPSFEED/LMD'
+ io_form_metgrid = 2,
+ opt_output_from_metgrid_path='./WRFFEED/current'
+/
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TH_NEST_namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TH_NEST_namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TH_NEST_namelist.input	(revision 348)
@@ -0,0 +1,72 @@
+ &time_control
+ start_year       = 2024, 2024
+ start_month      =   04,   04 
+ start_day        =   01,   01
+ start_hour       =   06,   06
+ end_year         = 2024, 2024
+ end_month        =   04,   04
+ end_day          =   01,   01   !!05,   05
+ end_hour         =   12,   12   !!06,   06
+ history_interval   = 74,   74   ! *100s : 74 is each 2 MH, 37 is each MH  !! attention, compatibilite avec time_step 
+ frames_per_outfile = 12,   12 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ input_from_file  = T, T   !!F     ! F en second pour que d1 init. d2 (smoother)
+ /
+ fine_input_stream = 2, 2
+
+ &domains
+ time_step   = 100  !!60 
+ dx = 60000, 20000
+ dy = 60000, 20000 
+ e_we   = 207, 211            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   = 207, 211            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert =  61,  61            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 3
+!!!!! NESTING STUFF
+ max_dom = 2	!! PLEASE DO NOT FORGET callphysd2_.def	
+ grid_id                = 1,  2
+ parent_id              = 0,  1
+ i_parent_start         = 0, 69      ! attention au raccord avec WPS !
+ j_parent_start         = 0, 69      ! attention au raccord avec WPS !
+ parent_grid_ratio      = 1,  3 
+ parent_time_step_ratio = 1,  5  !3  !4 !5
+ feedback               = 1
+ smooth_option          = 2
+ blend_width            = 5
+!!!!! NESTING STUFF
+ /
+
+ &physics
+ radt =  5,  5 	!10, 10	! call each half hour with radt = 1850/dt
+ mars =  1,  1  ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF)
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ non_hydrostatic = F, T
+ diff_6th_opt = 1,         !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES
+ diff_6th_factor = 0.8,    !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES
+ /
+
+ &grib2
+ /
+
+ &bdy_control
+ specified = T, F
+ nested    = F, 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/REF_ARTICLE/TH_NEST_namelist.wps
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TH_NEST_namelist.wps	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TH_NEST_namelist.wps	(revision 348)
@@ -0,0 +1,36 @@
+&share
+ wrf_core = 'ARW',
+ max_dom = 2,
+ start_date = '2024-04-01_06:00:00','2024-04-01_06:00:00'
+ end_date   = '2024-04-05_06:00:00','2024-04-01_06:00:00'  
+ interval_seconds = 3700
+ io_form_geogrid = 2,
+ debug_level = 0,
+ opt_output_from_geogrid_path='./'
+/
+
+
+&geogrid
+ parent_id         =   1,    1
+ parent_grid_ratio =   1,    3     
+ i_parent_start    =   1,   69    
+ j_parent_start    =   1,   69    
+ e_we              =  207,211  !!! ATTENTION
+ e_sn              =  207,211  !!! ATTENTION
+ geog_data_res     = '32ppd','64ppd_noHRti'
+ dx = 60000, !40000,
+ dy = 60000, !40000,
+ map_proj = 'mercator',
+ ref_lat   =  5.00,
+ ref_lon   =  -120.00,
+ truelat1  =  0.0,
+ truelat2  =  0.0,
+ stand_lon =  0.0,
+ geog_data_path = './WPS_GEOG',
+/
+
+&metgrid
+ fg_name = './WPSFEED/LMD'
+ io_form_metgrid = 2,
+ opt_output_from_metgrid_path='./WRFFEED'
+/
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TH_namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TH_namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TH_namelist.input	(revision 348)
@@ -0,0 +1,57 @@
+ &time_control
+ start_year       = 2024,
+ start_month      = 04, 
+ start_day        = 01,
+ start_hour       = 06,
+ end_year         = 2024,
+ end_month        = 04,
+ end_day          = 05, 
+ end_hour         = 06, 
+ history_interval    = 74,   ! *100s : 74 is each 2 MH, 37 is each MH  
+ frames_per_outfile  = 12, 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ /
+
+ &domains
+ time_step = 20	    
+ dx = 20000,
+ dy = 20000,
+ e_we   = 201,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   = 201,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert = 61,                ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 3
+ /
+
+ &physics
+ radt = 1,   ! call each half hour with radt = 1850/dt
+ mars = 1,   ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF)
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ diff_6th_opt = 2          !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES
+ diff_6th_factor = 0.2,    !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES
+ diff_opt = 1,    
+ km_opt = 4,
+ /
+
+ &grib2
+ /
+
+ &bdy_control
+ /
+
+!! DO NOT MODIFY
+ &namelist_quilt
+ nio_tasks_per_group = 0,
+ nio_groups = 1,
+ /
+!! DO NOT MODIFY
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TH_namelist.wps
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TH_namelist.wps	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/TH_namelist.wps	(revision 348)
@@ -0,0 +1,36 @@
+&share
+ wrf_core = 'ARW',
+ max_dom = 1,
+ start_date = '2024-04-01_06:00:00'
+ end_date   = '2024-04-05_06:00:00'
+ interval_seconds = 3700
+ io_form_geogrid = 2,
+ debug_level = 0,
+ opt_output_from_geogrid_path='./'
+/
+
+
+&geogrid
+ parent_id         =   1,
+ parent_grid_ratio =   1,
+ i_parent_start    =   1,
+ j_parent_start    =   1,
+ e_we              =  201,
+ e_sn              =  201,
+ geog_data_res     = '64ppd_noHRti'
+ dx = 20000,
+ dy = 20000,
+ map_proj = 'mercator',
+ ref_lat   =  5.00,
+ ref_lon   =  -120.00,
+ truelat1  =  0.0,
+ truelat2  =  0.0,
+ stand_lon =  0.0,
+ geog_data_path = './WPS_GEOG',
+/
+
+&metgrid
+ fg_name = './WPSFEED/LMD'
+ io_form_metgrid = 2,
+ opt_output_from_metgrid_path='./WRFFEED/current'
+/
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/VM_callphys.def
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/VM_callphys.def	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/VM_callphys.def	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+F
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+4        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+F
+callcond  (call CO2 condensation ?)
+T
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+F
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+100 
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+F
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+F
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+F
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/VM_launch_simu
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/VM_launch_simu	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/VM_launch_simu	(revision 348)
@@ -0,0 +1,141 @@
+#! /bin/bash
+
+
+#
+# USE: launch_simu &
+#
+
+#
+# must set env. var. 	WRF_current 
+#			
+
+# PRELIMINARY STEP
+# - model must be compiled and linked (wrf.exe, real.exe)
+# - domain must be correctly set with geogrid.exe
+# - put namelist.input, namelist.wps, callphys*.def
+
+# 'ls' must yield:
+# callphys*.def  launch_simu  namelist.input  namelist.wps  real.exe  wrf.exe
+#
+
+
+###############################################################
+gcmstep=0
+wpsstep=1
+realstep=1
+wrfstep=1
+###############################################################
+# GCM STEP
+   daygcm=593
+# WPS STEP
+   name='vm'     # un nom de dossier non pris dans WRFFEED
+# WRF STEP
+   proc=4
+###############################################################
+# FOLDERS
+   wheregcm=$GCM
+   whereprep=$WRF_current/PREP_MARS
+   wherewps=$WRF_current/WPS
+###############################################################
+
+
+#
+# LOG FILE HEADER
+#
+
+where=$PWD
+echo "check what's happening in " ${where}/log_launch
+
+rm -rf ${where}/log_launch 2> /dev/null
+echo SIMU folder is ${where} > ${where}/log_launch
+echo SIMU WILL BE LAUNCHED on the machine >> ${where}/log_launch
+hostname >> ${where}/log_launch
+
+#
+# 0. GCM step
+#
+if [ ${gcmstep} -eq 1 ]
+then
+
+   echo "0. GCM run and output preparation for WPS" >> ${where}/log_launch
+
+   cd ${wheregcm}
+   ln -sf ../STARTBASE_64_48_25_t2/start/start${daygcm}.nc start.nc
+   ln -sf ../STARTBASE_64_48_25_t2/startfi/startfi${daygcm}.nc startfi.nc
+
+   rm -rf restart.nc
+   rm -rf restartfi.nc
+   rm -rf log_gcm
+
+   echo running the GCM - check in ${wheregcm}/log_gcm >> ${where}/log_launch
+   nohup gcm.e > log_gcm
+
+   echo preparing outputs for WPS >> ${where}/log_launch
+   cd ${whereprep}
+   compile_and_exec >> ${where}/log_launch
+
+fi
+
+#
+# 1. WPS step
+#
+if [ ${wpsstep} -eq 1 ]
+then
+
+   echo "1. WPS: horizontal interpolations of static and meteorological data" >> ${where}/log_launch
+
+#   rm -rf ${where}/met_em 
+   mkdir ${where}/met_em
+   cd ${where}/met_em
+   ln -sf ../namelist.wps .
+
+   cd ${wherewps}
+   rm -rf WHERE_SIMU_IS
+   ln -sf ${where}/met_em WHERE_SIMU_IS
+   mkdir ./WRFFEED/current
+   nohup geogrid.exe >> ${where}/log_launch
+   nohup metgrid.exe >> ${where}/log_launch 
+   mv ./WRFFEED/current ./WRFFEED/${name}
+
+   ln -sf $FEED/WRFFEED/${name}/* ${where}/met_em/
+
+fi
+
+#
+# 2. REAL.EXE
+#
+if [ ${realstep} -eq 1 ]
+then
+
+   echo "2. REAL.EXE: vertical interpolation and preparation of boundary conditions" >> ${where}/log_launch
+
+   cd ${where}/met_em
+   rm -rf rsl.* 2> /dev/null
+   rm -rf wrfbdy* 2> /dev/null
+   rm -rf wrfinput* 2> /dev/null
+
+   ln -sf ../namelist.input .
+   ln -sf ../real.exe .
+   real.exe 
+   cd .. 
+   ln -sf ./met_em/wrf*d0* .
+
+fi
+
+
+#
+# 3. WRF.EXE
+#
+if [ ${wrfstep} -eq 1 ]
+then
+
+   echo "3. WRF.EXE: launch the mesoscale simulation" >> ${where}/log_launch
+
+   rm -rf rsl.*
+   rm -rf wrfout*
+
+   $WHERE_MPI/mpirun -np ${proc} wrf.exe < /dev/null &
+   echo OK, simulation launched !
+
+fi
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/VM_namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/VM_namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/VM_namelist.input	(revision 348)
@@ -0,0 +1,58 @@
+ &time_control
+ start_year       = 2024,
+ start_month      = 11, 
+ start_day        = 34,
+ start_hour       = 00,
+ end_year         = 2024,
+ end_month        = 11,
+ end_day          = 38, 
+ end_hour         = 00, 
+ history_interval    = 74,   ! *100s : 74 is each 2 MH, 37 is each MH  
+ frames_per_outfile  = 12, 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ /
+
+ &domains
+ time_step   = 10, 
+ dx = 12000,
+ dy = 12000,
+ e_we   = 401,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   = 121,               ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert = 61,                ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 3	 !! 1 et 2 ne passent pas
+ /
+
+ &physics
+ radt = 1,  	! call each half hour with radt = 1850/dt
+ mars = 0,      ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF)
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ diff_6th_opt = 2,         !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES
+ diff_6th_factor = 0.2,    !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES
+ diff_opt = 1,
+ km_opt = 4,
+ pd_scalar = F,
+ /
+
+ &grib2
+ /
+
+ &bdy_control
+ /
+
+!! DO NOT MODIFY
+ &namelist_quilt
+ nio_tasks_per_group = 0,
+ nio_groups = 1,
+ /
+!! DO NOT MODIFY
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/VM_namelist.wps
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/VM_namelist.wps	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/VM_namelist.wps	(revision 348)
@@ -0,0 +1,36 @@
+&share
+ wrf_core = 'ARW',
+ max_dom = 1,
+ start_date = '2024-11-34_00:00:00'
+ end_date = '2024-11-38_00:00:00'
+ interval_seconds = 3700
+ io_form_geogrid = 2,
+ debug_level = 0,
+ opt_output_from_geogrid_path='./'
+/
+
+
+&geogrid
+ parent_id         =   1,
+ parent_grid_ratio =   1,
+ i_parent_start    =   1,
+ j_parent_start    =   1,
+ e_we              =  401,
+ e_sn              =  121,
+ geog_data_res     = '32ppd_HRalb',
+ dx = 12000,
+ dy = 12000,
+ map_proj = 'mercator',
+ ref_lat   =  -8.00,
+ ref_lon   =  -68.00,
+ truelat1  =  0.0,
+ truelat2  =  0.0,
+ stand_lon =  0.0,
+ geog_data_path = './WPS_GEOG',
+/
+
+&metgrid
+ fg_name = './WPSFEED/LMD'
+ io_form_metgrid = 2,
+ opt_output_from_metgrid_path='./WRFFEED/current'
+/
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.assim
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.assim	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.assim	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+F
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+1352        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+T
+callcond  (call CO2 condensation ?)
+T
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+F
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+1
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+F
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+F
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+F
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.gcm.mars1
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.gcm.mars1	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.gcm.mars1	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+T
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+4        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+T
+callcond  (call CO2 condensation ?)
+T
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+T
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+10 1 
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+T
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+T
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+T
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars0
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars0	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars0	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+F
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+4        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+T
+callcond  (call CO2 condensation ?)
+T
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+F
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+1 
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+F
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+F
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+F
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars0_PATH
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars0_PATH	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars0_PATH	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+F
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+4        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+T F
+callcond  (call CO2 condensation ?)
+T
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+F
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+15 
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+F
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+F
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+F
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars0_noconvadj
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars0_noconvadj	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars0_noconvadj	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+F
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+4        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+F
+callcond  (call CO2 condensation ?)
+T
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+F
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+1 
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+F
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+F
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+F
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars1
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars1	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars1	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+T
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+4        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+T
+callcond  (call CO2 condensation ?)
+T
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+F
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+1 
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+T
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+T
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+T
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars1.MGSdust
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars1.MGSdust	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars1.MGSdust	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+T
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+3        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+T
+callcond  (call CO2 condensation ?)
+T
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+F
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+1 
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+T
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+T
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+T
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars1_noconvadj
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars1_noconvadj	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.mars1_noconvadj	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+T
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+4        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+F
+callcond  (call CO2 condensation ?)
+T
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+F
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+1 
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+T
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+T
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+T
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.microscale
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.microscale	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/callphys.def.microscale	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+F
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+1        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+F
+callcond  (call CO2 condensation ?)
+F
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+F
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+1 
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+F
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+F
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+F
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/dustopacity.def
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/dustopacity.def	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/dustopacity.def	(revision 348)
@@ -0,0 +1,1 @@
+0.3
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/old_namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/old_namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/old_namelist.input	(revision 348)
@@ -0,0 +1,133 @@
+ &time_control
+ run_days                            = 0,
+ run_hours                           = 0,
+ run_minutes                         = 0,
+ run_seconds                         = 0,
+ start_year                          = 2025,	2026,
+ start_month                         = 07 !06,	01, 
+ start_day                           = 02 !49,	47, 
+ start_hour                          = 00,	00, 
+ start_minute                        = 00,	00, 
+ start_second                        = 00,	00, 
+ end_year                            = 2025,	2026,
+ end_month                           = 07, 	07,
+ end_day                             = 15, 	09,
+ end_hour                            = 00, 	22,
+ end_minute                          = 00, 	00,
+ end_second                          = 00, 	00,
+ interval_seconds		     = 7400
+ input_from_file                     = .true., .false.		! .false. en second pour que d1 init. d2
+ fine_input_stream		     = 0
+ history_interval                    = 74,	74,		! *100s : 74 is each 2 MH, 37 is each MH  
+ frames_per_outfile                  = 12,	12, 
+ restart                             = .false.
+ restart_interval                    = 8880,
+ io_form_history                     = 2
+ io_form_restart                     = 2
+ io_form_input                       = 2
+ io_form_boundary                    = 2
+ debug_level                         = 0	!200
+ /
+
+
+ &domains
+ time_step                           = 18 !37, 
+ time_step_fract_num                 = 0,
+ time_step_fract_den                 = 1,
+max_dom = 1
+s_we    = 1,1,
+e_we    = 150,49,
+s_sn    = 1,1,
+e_sn    = 150,49,
+s_vert  = 1,1,
+e_vert  = 31,50,
+ num_metgrid_levels                  = 26, 
+ max_dz                              = 2000.
+ sfcp_to_sfcp                        = .true.
+ lowest_lev_from_sfc                 = .false.
+ p_top_requested                     = 2,
+ dx                                  = 50000,	5000,
+ dy                                  = 50000,	5000,
+ grid_id                             = 1,	2,
+ parent_id                           = 0,	1,
+ i_parent_start                      = 0,	12,	! attention au raccord avec WPS !	
+ j_parent_start                      = 0,	12,	! attention au raccord avec WPS !
+ parent_grid_ratio                   = 1,	3,
+ parent_time_step_ratio              = 1,	2,
+ feedback                            = 1,
+ smooth_option                       = 2
+ blend_width			     = 5
+ interp_type                         = 2,       ! pas d'influence sur la stabilite
+ zap_close_levels                    = 0,       ! set it to >0 when mountaneous domain
+ force_sfc_in_vinterp                = 9        ! default: 8 (number of first layers prescribed for PBL)
+ /
+ eta_levels		= 1.0,0.95,0.9,0.85,0.8,0.75,0.7,0.65,0.6,0.55,0.5,0.45,0.4,0.35,0.3,0.25,0.2,0.15,0.1,0.05,0
+ lagrange_order                      = 1
+ interp_type			     = 2
+ adjust_heights                      = .false.
+ ztop                                = 60000.
+
+
+ 
+
+ &physics
+ mp_physics                          = 0,	0,
+ ra_lw_physics                       = 0,	0,
+ ra_sw_physics                       = 0,	0,
+ radt                                = 2,	10,	
+ sf_sfclay_physics                   = 0,	0,
+ sf_surface_physics                  = 0,	0,
+ bl_pbl_physics                      = 0,	0,
+ bldt                                = 0,	0,
+ cu_physics                          = 0,	0,
+ cudt                                = 400 !200,       1,		
+ isfflx                              = 1,	
+ ifsnow                              = 0,	
+ icloud                              = 1,	
+ num_soil_layers                     = 5,	
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ dyn_opt                             = 2,
+ rk_ord                              = 3,
+ diff_opt                            = 0,
+ km_opt                              = 4,
+ diff_6th_opt                        = 0,
+ diff_6th_factor                     = 0.12,
+ base_temp                           = 220.
+ base_pres			     = 610.
+ base_lapse			     = 30.
+ damp_opt                            = 0,
+ zdamp                               = 3000,
+ dampcoef                            = 0.05,
+ w_damping			     = 0 	
+ smdiv                               = 0.1,	0.1,
+ emdiv                               = 0.01,	0.01,
+ epssm                               = 0.1,	0.1,
+ non_hydrostatic                     = .true.,	.true.,
+ time_step_sound                     = 4,	4,
+ h_mom_adv_order                     = 3, !5,	5,
+ v_mom_adv_order                     = 3,	3,
+ h_sca_adv_order                     = 3, !5,	5,
+ v_sca_adv_order                     = 3,	3,
+
+
+ &grib2
+ /
+
+ &bdy_control
+ spec_bdy_width                      = 5,
+ spec_zone                           = 1,
+ relax_zone                          = 4,
+ specified                           = .true.,	.false.,
+ nested                              = .false., .true.,
+ /
+
+ &namelist_quilt
+ nio_tasks_per_group = 0,
+ nio_groups = 1,
+ /
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/yeyeye/callphys.def
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/yeyeye/callphys.def	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/yeyeye/callphys.def	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+F
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+4        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+T F
+callcond  (call CO2 condensation ?)
+T
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+F
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+20 
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+F
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+F
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+F
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/yeyeye/callphysd2_.def
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/yeyeye/callphysd2_.def	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/yeyeye/callphysd2_.def	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+F
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+4        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+T F
+callcond  (call CO2 condensation ?)
+T
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+F
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+40 
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+F
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+F
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+F
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/yeyeye/callphysd3_.def
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/yeyeye/callphysd3_.def	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/yeyeye/callphysd3_.def	(revision 348)
@@ -0,0 +1,94 @@
+General options
+~~~~~~~~~~~~~~~
+tracer    (Run with or without tracer transport ?)
+F
+diurnal   (Diurnal cycle ?  if diurnal=F, diurnal averaged solar heating)
+T
+season    (Seasonal cycle ? if season=F, Ls stays constant like in "start")
+T
+lwrite    (want some more output on the screen ?) 
+F
+stats     (Saving statistics in file "cumul" ?) -- noWRF
+F
+calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF
+F
+Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
+~~~~~~~~~~~~~
+iaervar  (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario 
+4        (=4 Mars Year 24 from TES assimilation)
+iddist  (Dust vertical distribution: =0: old distrib. (Pollack90) 
+3       (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario )
+topdustref (Dust top altitude (km). Matter only if iddist=1)
+55.
+Physical Parameterizations :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+callrad   (call radiative transfer ?)
+T
+callnlte (call NLTE radiative schemes ?   matter only if callrad=T)
+F
+callnirco2 (call CO2 NIR absorption ?   matter only if callrad=T)
+T
+calldifv  (call turbulent vertical diffusion ?)
+T
+calladj   (call convective adjustment ?)
+T F
+callcond  (call CO2 condensation ?)
+T
+callsoil  (call thermal conduction in the soil ?)
+T
+calllott  (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF
+F
+Radiative transfer options :
+~~~~~~~~~~~~~~~~~~~~~~~~~~
+iradia    (the rad.transfer is computed every "iradia" physical timestep)
+80 
+callg2d   (Output of the exchange coefficient mattrix ? for diagnostic only)
+F
+rayleigh  (Rayleigh scattering : should be =F for now)
+F
+Tracer (dust water, ice and/or chemical species) options (use if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+dustbin    (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin))
+0
+active     (DUST: Radiatively active dust ? (uses q(1) to q(dustbin))
+F
+doubleq    (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?)
+F
+lifting    (DUST: lifted by GCM surface winds ?)
+F
+dustdevil  (DUST: lifted by dust devils ?)
+F
+scavenging (DUST: Scavenging by CO2 snowfall ?)
+F
+sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?)
+F
+iceparty   (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1))
+F
+activice   (WATERICE: Radiatively active transported atmospheric water ice ?)
+F
+water      (WATER: Compute water cycle using q(nqmx) )
+F
+caps       (WATER: put the current permanent caps at both poles)
+F
+photochem  (PHOTOCHEMISTRY: chemical species included)
+F
+Thermospheric options (relevant if tracer=T) :
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+callthermos  (call thermosphere ?)
+F
+thermoswater  (WATER: included without cycle	only if water=F)
+F
+callconduct  (call thermal conduction ?     matter only if callthermos=T)
+F
+calleuv  (call EUV heating ?                matter only if callthermos=T)
+F
+callmolvis  (call molecular viscosity ?     matter only if callthermos=T)
+F
+callmoldiff  (call molecular diffusion ?    matter only if callthermos=T)
+F
+thermochem  (call thermospheric photochemistry ?  matter only if callthermos=T)
+F
+solarcondate (date for solar flux calculation: 1985 < date < 2002))
+1993.4       (Solar min=1996.4 ave=1993.4 max=1990.6)
+
+
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/yeyeye/makemeso_input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/yeyeye/makemeso_input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/yeyeye/makemeso_input	(revision 348)
@@ -0,0 +1,6 @@
+1
+153
+153
+61
+3
+1
Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/yeyeye/namelist.input
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/yeyeye/namelist.input	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/yeyeye/namelist.input	(revision 348)
@@ -0,0 +1,76 @@
+ &time_control
+ start_year       = 2024, 2024, 2024,
+ start_month      =   01,   01,   01,
+ start_day        =   49,   49,   49,
+ start_hour       =   06,   06,   06,
+ end_year         = 2024, 2024, 2024,
+ end_month        =   01,   01,   01,
+ end_day          =   54,   54,   54,
+ end_hour         =   06,   06,   06,
+ history_interval   = 37,   37,   37,    
+! history_interval_s   = 1,1,1
+ frames_per_outfile = 24,   24,   24, 
+ restart          = .false.
+ restart_interval = 8880,
+ io_form_history  = 2
+ io_form_restart  = 2
+ io_form_input    = 2
+ io_form_boundary = 2
+ debug_level      = 0
+ input_from_file  =    T,    T,    T,     ! F en second pour que d1 init. d2 (smoother)
+ /
+
+ &domains
+ time_step   = 40 
+ dx     = 36000,  9000,  2250,   !!alt 48000 12000 3000
+ dy     = 36000,  9000,  2250, 
+ e_we   =   153,   157,   157,            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_sn   =   153,   157,   157,            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ e_vert =    61,    61,    61,            ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS
+ p_top_requested = 3
+!!!!! NESTING STUFF
+ max_dom = 3	!!!! PLEASE DO NOT FORGET callphysd*_.def	
+ grid_id                = 1,   2,   3,
+ parent_id              = 0,   1,   2,
+ i_parent_start         = 0,  57,  57,      ! attention au raccord avec WPS !
+ j_parent_start         = 0,  57,  57,      ! attention au raccord avec WPS !
+ parent_grid_ratio      = 1,   4,   4,
+ parent_time_step_ratio = 1,   2,   2,
+ feedback               = 0  
+ smooth_option          = 2
+ blend_width            = 5
+!!!!! NESTING STUFF
+ /
+
+ &physics
+ radt =  1, 1, 1     ! 1 is the best (change iradia in callphys)
+ mars =  0, 0, 0     ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF)
+!!! init_MU=3.   ! ... quasi no-wind
+!!! init_MV=3.   ! ... quasi no-wind
+ /
+
+ &fdda
+ /
+
+ &dynamics
+ diff_6th_opt =      2,   2,   2,     !! NECESSAIRE 
+ diff_6th_factor =   0.2, 0.2, 0.2    !! NECESSAIRE 
+ diff_opt = 1,
+ km_opt = 4,
+ pd_scalar = F   !! no water cycle (mars = 0)
+ /
+
+ &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/REF_ARTICLE/yeyeye/namelist.wps
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/yeyeye/namelist.wps	(revision 348)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/REF_ARTICLE/yeyeye/namelist.wps	(revision 348)
@@ -0,0 +1,35 @@
+&share
+ wrf_core = 'ARW',
+ max_dom = 3,
+ start_date = '2024-01-49_06:00:00', '2024-01-49_06:00:00', '2024-01-49_06:00:00'
+ end_date   = '2024-01-54_06:00:00', '2024-01-49_06:00:00', '2024-01-49_06:00:00'
+ interval_seconds = 3700
+ io_form_geogrid = 2,
+ debug_level = 0,
+ opt_output_from_geogrid_path='./'
+/
+
+&geogrid
+ parent_id         =        0,           1,          2  
+ parent_grid_ratio =        1,           4,          4
+ i_parent_start    =        0,          57,         57 
+ j_parent_start    =        0,          57,         57
+ e_we              =      153,         157,        157      !!! ATTENTION
+ e_sn              =      153,         157,        157      !!! ATTENTION
+ geog_data_res     =  '64ppd',     '64ppd',    '64ppd'      
+ dx = 36000
+ dy = 36000
+ map_proj = 'mercator',
+ ref_lat   =  11.5,
+ ref_lon   =  15.5,
+ truelat1  =  0.0,
+ truelat2  =  0.0,
+ stand_lon =  00.00,
+ geog_data_path = './WPS_GEOG',
+/
+
+&metgrid
+ fg_name = './WPSFEED/LMD'
+ io_form_metgrid = 2,
+ opt_output_from_metgrid_path='./WRFFEED/current'
+/
