Mars 1D testphys1d program

From Planets
Revision as of 16:41, 22 August 2023 by Emillour (talk | contribs) (Created page with "== Overview == It is possible to run the Mars PCM in a "single-column" configuration (quite useful for some studies and or when developing and testing parametrizations). == C...")

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

Overview

It is possible to run the Mars PCM in a "single-column" configuration (quite useful for some studies and or when developing and testing parametrizations).

Compiling testphys1d

The main program testphys1d is compiled using the same compilation script makelmdz_fcm than the 3D model except that:

  • the -d option requires only one argument, the number of vertical levels
  • the main program to compile is testphys1d rather than gcm

So for instance to compile a case for 54 vertical levels one would run something like

makelmdz_fcm -arch somearch -d 54 -p mars testphys1d

testphys1d inputs

Just like the gcm the 1D model needs some inputs to run. The main difference is that it does not require a start (initial conditions) file. Some def files are however needed:

  • z2sig.def : definition of vertical levels
  • traceur.def : definition of tracers the model will run with. Note that at run time testphys1d looks for optional profile_sometracername files which contain initial values (surface value on the first line and then one value for each layer to describe the profile) of the tracer quantities.
  • callphys.def : definition of parametrizations to switch on (just like in the 3D PCM)
  • run.def : the run configuration file; similar to the one for the 3D PCM described here The run.def Input File, but adapted to the 1D case. Note that there is an example of a run.def file for testphys1d in LMDZ.MARS/deftank (file run.def.1d)

testphys1d outputs

Testphys1d can output diagfi.nc and stats.nc files, and just like for the GCM optional diagfi.def and 'stats.def ASCII files can be used to specify which variable needs be outputted.