Notes on Mars topography
The topography software has been altered to support the generation of Mars topography out-of-the-box. A namelist variable was added for the user to specify Mars gravity so that PHIS is computed correctly (--gravity 3.72076). Since the Mars dataset is not high resolution a namelist variable was added to compute SGH30 from SGH (--compute_sgh30_from_sgh_fac 0.5) by scaling SGH with 0.5.
The command used to produce Mars topography on Casper is:
qcmd -l walltime=24:00:00 -A P03010039 -l select=1:ncpus=4:mem=200GB -- ./cube_to_target --compute_sgh30_from_sgh_fac 0.50000000000000000 --jmax_segments 5000000 --smoothing_over_ocean --grid_descriptor_file /glade/campaign/cesm/cesmdata/inputdata/share/meshes/ne16np4_ESMFmesh_cdf5_c20211018.nc --source_data_identifier mars --intermediate_cs_name /glade/campaign/cgd/amp/pel/topo/cubedata/mars-ncube3000.nc --output_grid ne16np4 --smoothing_scale 200.0 --name_email_of_creator Lauritzen --gravity 3.72076000 >& output.mars.ne16pg3.txt &
When running on Derecho the model stopped with error message "Killed". It looks like it ran out of memory. Problem fixed by logging into casper and run on a node with:
execcasper -A P03010039 -l select=1:ncpus=4:mem=100GB