Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/callphys.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/callphys.def (revision 223) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/callphys.def (revision 223) @@ -0,0 +1,92 @@ +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" ?) +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) +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 ?) +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +10 +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 nb q(2) mr 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/namelist.input_full =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/namelist.input_full (revision 217) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/namelist.input_full (revision 223) @@ -17,11 +17,8 @@ io_form_boundary = 2 !! (*) Choice of NETCDF for ouputs debug_level = 0 !! (*) Verbose level - !! - !! OPTIONAL - !! + !!!!! OPTIONAL !!!!!!!!!!!!!!! interval_seconds = 3700 !! (p2) Frequency of large-scale fields update (s) input_from_file = T, !! (n)(p2) Initialize a given domain with an input file / - &domains @@ -33,7 +30,5 @@ e_vert = 61, !! (r)(p2) Number of vertical levels p_top_requested = 5 !! (p3) Chosen value of pressure at the top of the model - !! - !! OPTIONAL - !! + !!!!! OPTIONAL !!!!!!!!!!!!!!! time_step_fract_num = 0 !! Additional fraction to time_step: numerator time_step_fract_den = 1 !! Additional fraction to time_step: denominator @@ -53,7 +48,5 @@ &physics - !! - !! OPTIONAL - !! + !!!!! OPTIONAL !!!!!!!!!!!!!!! radt = 1, !! Ratio between physical and dynamical time step mars = 0, !! (r)(p2) Configuration of tracers: @@ -71,11 +64,9 @@ &dynamics - !! - !! OPTIONAL - !! + !!!!! OPTIONAL !!!!!!!!!!!!!!! time_step_sound = 6, !! Ratio of time step dynamic/acoustic integration !! NB: an increase could help solve instabilities non_hydrostatic = T, !! Integrate in non-hydrostatic/hydrostatic mode - pd_scalar = F, !! Positive-definite advection scheme for tracers + pd_scalar = T, !! Positive-definite advection scheme for tracers !! diff_opt = 1 !! (*d) Diffusion option [set to 0 if LES or GCM] @@ -93,7 +84,5 @@ &bdy_control - !! - !! OPTIONAL - !! + !!!!! OPTIONAL !!!!!!!!!!!!!!! specified = T, !! (n)(p3) Boundary conditions specified by GCM nested = F, !! (n)(p3) Boundary conditions from parent domain @@ -109,4 +98,5 @@ / + !!!!! DO NOT MODIFY !!!!!!!!!! &grib2 / Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/namelist.input_minim =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/namelist.input_minim (revision 217) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/namelist.input_minim (revision 223) @@ -1,31 +1,30 @@ &time_control - start_year = 2024, !! (p1) Start Martian Year (20XX for MY XX) - start_month = 07, !! (p1) Start Martian Month - start_day = 01, !! (p1) Start Martian Day - start_hour = 06, !! (p1) Start Martian Hour (at longitude 0) - end_year = 2024, !! (p1) End Martian Year (20XX for MY XX) - end_month = 07, !! (p1) End Martian Month - end_day = 02, !! (p1) End Martian Day - end_hour = 06, !! (p1) End Martian Hour (at longitude 0) - history_interval = 37, !! Frequency of outputs (37 --> 3700s = 1 Martian hour) - frames_per_outfile = 24, !! Size of time dimension in files - restart = .false. !! (*) Output restart files ? - restart_interval = 8880 !! (*) Frequency of output restart files ? - io_form_history = 2 !! (*) Choice of NETCDF for ouputs - io_form_restart = 2 !! (*) Choice of NETCDF for ouputs - io_form_input = 2 !! (*) Choice of NETCDF for ouputs - io_form_boundary = 2 !! (*) Choice of NETCDF for ouputs - debug_level = 0 !! (*) Verbose level + 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, + 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 = 50 !! Dynamical timestep - dx = 20000, !! (p2) Horizontal resolution - dy = 20000, !! (p2) Horizontal resolution (should be equal to dx) - e_we = 61, !! (r)(p2) Number of longitude grid points - e_sn = 61, !! (r)(p2) Number of latitude grid points - e_vert = 61, !! (r)(p2) Number of vertical levels - p_top_requested = 5 !! (p3) Chosen value of pressure at the top of the model + time_step = 50 + dx = 20000, + dy = 20000, + e_we = 61, + e_sn = 61, + e_vert = 61, + p_top_requested = 5 / @@ -45,6 +44,6 @@ / - &namelist_quilt !! (*) - nio_tasks_per_group = 0, !! (*) - nio_groups = 1, !! (*) - / !! (*) + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/WPS/wps_mars/namelist.wps_TEST =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/WPS/wps_mars/namelist.wps_TEST (revision 217) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WPS/wps_mars/namelist.wps_TEST (revision 223) @@ -1,38 +1,35 @@ &share - wrf_core = 'ARW', !! [do not modify: choice of dynamical core] - max_dom = 1, !! number of simulation domains - start_date = '0000-00-00_00:00:00' !! YYYY-MM-DD_HH:mm:ss start date - end_date = '1111-11-11_11:11:11' !! YYYY-MM-DD_HH:mm:ss end date - interval_seconds = 3700 !! frequency of GCM updates [1 Mars hour = 3700 s] - io_form_geogrid = 2, !! [do not modify: choice of NETCDF outputs] - debug_level = 0, !! verbose level of the programs - opt_output_from_geogrid_path='./' !! location of the geogrid outputs + wrf_core = 'ARW', !! [do not modify: choice of dynamical core] + max_dom = 1, !! number of simulation domains + start_date = '0000-00-00_00:00:00' !! YYYY-MM-DD_HH:mm:ss start date + end_date = '1111-11-11_11:11:11' !! YYYY-MM-DD_HH:mm:ss end date + interval_seconds = 3700 !! frequency of GCM updates [1 Mars hour = 3700 s] + io_form_geogrid = 2, !! [do not modify: choice of NETCDF outputs] + debug_level = 0, !! verbose level of the programs + opt_output_from_geogrid_path='./' !! location of the geogrid outputs / - &geogrid - parent_id = 1, !! number identifying the related parent domain - parent_grid_ratio = 1, !! ratio between parent and nested domains - i_parent_start = 1, !! x-position of the southwest corner of nest - j_parent_start = 1, !! y-position of the southwest corner of nest - e_we = 61, !! number of longitude grid points - e_sn = 61, !! number of latitude grid points - geog_data_res = 'gcm' !! choice of static data sources - !! NB: possible: '64ppd', '32ppd', ... - !! NB: please glance at geogrid/GEOGRID.TBL - dx = 20000, !! resolution (meters) in the x-dimension - dy = 20000, !! resolution (meters) in the y-dimension - map_proj = 'mercator', !! map projection: 'mercator', 'lambert' or 'polar' - ref_lat = -12., !! north latitude of the center of the domain - ref_lon = 239., !! east longitude of the center of the domain - truelat1 = 0.0, !! (lambert or polar) lat position of projection cone - truelat2 = 0.0, !! [do not modify] - stand_lon = 0.0, !! (lambert or polar) lon position of projection cone - geog_data_path = './WPS_GEOG', !! [do not modify: symbolic link in the WPS folder] + parent_id = 1, !! number identifying the related parent domain + parent_grid_ratio = 1, !! ratio between parent and nested domains + i_parent_start = 1, !! x-position of the southwest corner of nest + j_parent_start = 1, !! y-position of the southwest corner of nest + e_we = 61, !! number of longitude grid points + e_sn = 61, !! number of latitude grid points + geog_data_res = 'gcm' !! static data sources: '64ppd','32ppd',... cf.GEOGRID.TBL + dx = 20000, !! resolution (meters) in the x-dimension + dy = 20000, !! resolution (meters) in the y-dimension + map_proj = 'mercator', !! map projection: 'mercator', 'lambert' or 'polar' + ref_lat = -12., !! north latitude of the center of the domain + ref_lon = 239., !! east longitude of the center of the domain + truelat1 = 0.0, !! (lambert or polar) lat position of projection cone + truelat2 = 0.0, !! [do not modify] + stand_lon = 0.0, !! (lambert or polar) lon position of projection cone + geog_data_path = './WPS_GEOG', !! [do not modify: symbolic link in the WPS folder] / &metgrid - fg_name = './WPSFEED/LMD' !! [do not modify: symbolic link in the WPS folder] - io_form_metgrid = 2, !! [do not modify: choice of NETCDF outputs] + fg_name = './WPSFEED/LMD' !! [do not modify: symbolic link in the WPS folder] + io_form_metgrid = 2, !! [do not modify: choice of NETCDF outputs] opt_output_from_metgrid_path='./WRFFEED/current' !! [do not modify: symbolic link] /