wiki:MedusaCoupling

Version 3 (modified by smoreira, 7 years ago) (diff)

--

Coupling MEDUSA Sediment Model to iLOVECLIM

Running a simulation with MEDUSA

Running a simulation of the full iLOVECLIM carbon cycle (OCYCC) with MEDUSA requires to configure the original iLOVECLIM carbon cycle model as well as several input files (namelists) for MEDUSA model. The iLOVECLIM version coupled with MEDUSA works in the same way that the original iLOVECLIM for Carbon Cycle except for the new input and configuration files for running MEDUSA.

  • Modify the file: lbm/sources/choixcomposantes.h to have the same options than for the Carbon Cycle procedures and also MEDUSA 1, please check < https://forge.ipsl.jussieu.fr/ludus/wiki/CarbonCyclePage > for more details. Do not forget to establish your initial conditions for restarting using, for example, -s 50800 and -S data/default/ic50800 (to start from the initial conditions)
  • Prepare your restarting files for MEDUSA using reaclay, flx and sedcore files from a previous simulation or in case you don't want to use any initial condition for MEDUSA set the corresponding options of medusa_filelist.nml file to "/dev/null"
  • Configure your MEDUSA input files and namelists in medusa-svn/inputdata-cpl: medusa.rrp, medusa.tsi, medusa_filelist.nml, medusa_seafloor.init.
    • It is specially important to pay attention to medusa_filelist.nml, cfn_ncin_sedcore file is the main file for doing a restart in the MEDUSA simulation. In theory, it is possible to read this input file and create a different output file cfn_ncout_sedcore file. However, there is a small bug and this option does not work in our current version. Therefore, the procedure suggested is to make a copy of the original sedcore input nectdf file, re-name it and use this file as the same file for input and output (that works perfectly!)

Here you can find an example of medusa_filelist.nml

&nml_cfg
cfn_nmlin_init     = "medusa_seafloor_init.nml"
cfn_ncin_init      = "/MEDUSA_RESTART_DIR/mediloveclim_reaclay_file.nc"
cfn_ncin_flx       = "/MEDUSA_RESTART_DIR/mediloveclim_flx_file.nc"
cfn_ncin_sedcore   = "/MEDUSA_RESTART_DIR/mediloveclim_sedcore_file.nc"

cfn_ncout_reaclay  = "outputdata/sediments/mediloveclim_reaclay_newfile.nc"
cfn_ncout_reaction = "outputdata/sediments/mediloveclim_reaction_newfile.nc"
cfn_ncout_bc       = "outputdata/sediments/mediloveclim_bc_newfile.nc"
cfn_ncout_flx      = "outputdata/sediments/mediloveclim_flx_newfile.nc"
cfn_ncout_sedcore  = "/MEDUSA_RESTART_DIR/mediloveclim_sedcore_file.nc"
ctitle_ncfiles     = "Short description"
/
&nml_extra
cfn_ncout_aux      = "outputdata/sediments/mediloveclim_aux.nc"
/