Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/TESTCASE/EXEC_TO_GET_data.sh =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/TESTCASE/EXEC_TO_GET_data.sh (revision 1854) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/TESTCASE/EXEC_TO_GET_data.sh (revision 1854) @@ -0,0 +1,6 @@ +#! /bin/bash + +www=http://www.lmd.jussieu.fr/~aslmd/mesoscale_model/TESTCASE/ + +wget $www/wrfinput_d01 . +wget $www/wrfbdy_d01 . Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/TESTCASE/callphys.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/TESTCASE/callphys.def (revision 1854) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/TESTCASE/callphys.def (revision 1854) @@ -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 ?) +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) +37 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)) +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/TESTCASE/dustopacity.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/TESTCASE/dustopacity.def (revision 1854) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/TESTCASE/dustopacity.def (revision 1854) @@ -0,0 +1,1 @@ +0.3 Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/TESTCASE/namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/TESTCASE/namelist.input (revision 1854) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/TESTCASE/namelist.input (revision 1854) @@ -0,0 +1,70 @@ + &time_control + start_year = 2024, !! (p1) Start Martian Year (20XX for MY XX) + start_month = 01, !! (p1) Start Martian Month + start_day = 17, !! (p1) Start Martian Day + start_hour = 02, !! (p1) Start Martian Hour (at longitude 0) + end_year = 2024, !! (p1) End Martian Year (20XX for MY XX) + end_month = 01, !! (p1) End Martian Month + end_day = 17, !! (p1) End Martian Day + end_hour = 12, !! (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 + / + + + &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 + !! + !! OPTIONAL + !! + num_metgrid_levels = 33 !! (p1) number of vertical levels in GCM inputs (+1) + / + + &physics + !! + !! OPTIONAL + !! + init_TI = 100., !! Define constant thermal inertia value + init_AL = 0.27, !! Define constant albedo value + / + + &dynamics + !! + !! OPTIONAL + !! + pd_scalar = F + non_hydrostatic = F, !! Integrate in non-hydrostatic/hydrostatic mode + / + + &bdy_control + !! + !! OPTIONAL + !! + spec_bdy_width = 3 !! (p3) Width of transition zone with specified=T + !! NB: spec_bdy_width must equal relax_zone+1 + relax_zone = 2 !! (p3) Width of relaxation zone with specified=T + / + + &grib2 + / + + &fdda + / + + &namelist_quilt !! (*) + nio_tasks_per_group = 0, !! (*) + nio_groups = 1, !! (*) + / !! (*) Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/TESTCASE/namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/TESTCASE/namelist.wps (revision 1854) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/TESTCASE/namelist.wps (revision 1854) @@ -0,0 +1,35 @@ +&share + wrf_core = 'ARW', + max_dom = 1, + start_date = '2024-01-17_02:00:00', + end_date = '2024-01-17_12: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 = 61, + e_sn = 61, + geog_data_res = '64ppd', + dx = 20000, + dy = 20000, + map_proj = 'mercator', + ref_lat = -12., + ref_lon = 239., + 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/TESTCASE/wrf.exe =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/TESTCASE/wrf.exe (revision 1854) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/TESTCASE/wrf.exe (revision 1854) @@ -0,0 +1,1 @@ +link ../../../LATEST/wrf_x61_y61_z61_d1_t1_p1.exe Index: /trunk/MESOSCALE/MANUAL/SRC/compile_exec.tex =================================================================== --- /trunk/MESOSCALE/MANUAL/SRC/compile_exec.tex (revision 1853) +++ /trunk/MESOSCALE/MANUAL/SRC/compile_exec.tex (revision 1854) @@ -161,14 +161,13 @@ \sk -In order to test the compiled executables, a ready-to-use test case (with pre-generated initial and boundary conditions) is proposed in the \ttt{LMD\_MM\_MARS\_TESTCASE.tar.gz} archive that you can download in the following FTP site \url{ftp://ftp.lmd.jussieu.fr/pub/aslmd/LMD_MM_MARS_TESTCASE.tar.gz}. This test case simulates the hydrostatic atmospheric flow around Arsia Mons (Figure~\ref{arsia}) during half a sol in springtime with constant thermal inertia, albedo and dust opacity\footnote{Though the simulation reproduces some reasonable features of the mesoscale circulation around Arsia Mons (e.g. slope winds), it should not be used for scientific purpose, for the number of grid points is unsufficient for single-domain simulation and the integration time is below the necessary spin-up time.}. - -\sk -To launch the test simulation, please type the following commands, replacing the \ttt{g95\_32\_single} directory with its corresponding value on your system. In the end, the model should run and output the computed meteorological fields in netCDF files named \ttt{wrfout*}. Feel free to browse those files with \ttt{ncview} or your favorite graphical tool to check if the simulated fields look reasonable. +In order to test the compiled executables, a ready-to-use test case (with pre-generated initial and boundary conditions) is proposed in the \ttt{LMD\_MM\_MARS/SIMU/DEF/TESTCASE} folder (use the \ttt{EXEC\_TO\_GET\_data.sh} command to get the initial and boundary condition files). This test case simulates the hydrostatic atmospheric flow around Arsia Mons (Figure~\ref{arsia}) during half a sol in springtime with constant thermal inertia, albedo and dust opacity\footnote{Though the simulation reproduces some reasonable features of the mesoscale circulation around Arsia Mons (e.g. slope winds), it should not be used for scientific purpose, for the number of grid points is unsufficient for single-domain simulation and the integration time is below the necessary spin-up time.}. + +\sk +To launch the test simulation, please type the following commands, replacing if needed the \ttt{LATEST} directory with its corresponding value on your system. In the end, the model should run and output the computed meteorological fields in netCDF files named \ttt{wrfout*}. Feel free to browse those files with \ttt{ncview} or your favorite graphical tool to check if the simulated fields look reasonable. % \begin{verbatim} -cp LMD_MM_MARS_TESTCASE.tar.gz $MMM -tar xzvf LMD_MM_MARS_TESTCASE.tar.gz -cd TESTCASE -ln -sf ../g95_32_single/wrf_x61_y61_z61_d1_t1_p1.exe wrf.exe +cd $MMM/SIMU/DEF/TESTCASE +./EXEC_TO_GET_data.sh (only once) +ln -sf ../../../LATEST/wrf_x61_y61_z61_d1_t1_p1.exe wrf.exe nohup wrf.exe > log_wrf & \end{verbatim} Index: /trunk/MESOSCALE/MANUAL/SRC/installation.tex =================================================================== --- /trunk/MESOSCALE/MANUAL/SRC/installation.tex (revision 1853) +++ /trunk/MESOSCALE/MANUAL/SRC/installation.tex (revision 1854) @@ -22,4 +22,5 @@ \item G95 free compiler\footnote{Sources and binaries available on \url{http://www.g95.org}} \ttt{g95} \item Intel commercial compiler \ttt{ifort} +\item GNU Fortran compiler \ttt{gfortran} \end{itemize} \item your C compiler is \ttt{gcc} and C development libraries are included; @@ -97,5 +98,6 @@ svn checkout the_link -N the_name_of_your_local_destination_folder cd the_name_of_your_local_destination_folder -svn update LMDZ.MARS MESOSCALE +svn update MESOSCALE +[svn update LMDZ.MARS] ## only if you use the new physics (not recommended for a starter) cd MESOSCALE declare -x MESO=$PWD ## put absolute link in your .bashrc @@ -103,5 +105,5 @@ ## to get latest updates later on cd the_name_of_your_local_destination_folder - svn update LMDZ.MARS MESOSCALE + svn update svn log | more \end{verbatim}