Use of 'define key_sediment' when using PISCES

ze miscellaneous
1

Hi everyone,

I am trying to vary some of the PISCES settings to see how the results will change for my model setup. In cppdefs.h, I have previously set ‘undef key_sediment’ in the PISCES options and my model is able to run fine. However, when I tried setting ‘define key_sediment’, the model encounters errors at the start and is unable to continue running. Looking at the log file, the errors first appear in the following lines:

   r p4z_rem_init : Initialization of remineralization parameters

ctl_opn bad opening
ctl_opn bad opening
   Namelist parameters for remineralization, nampisrem
      remineralization rate of DOC              xremik    =
 0.300000000000000
      remineralization rate of Si               xsirem    =
 3.000000000000000E-003
      fast remineralization rate of Si          xsiremlab =
 3.000000000000000E-002
      fraction of labile biogenic silica        xsilab    =
 0.500000000000000
      NH4 nitrification rate                    nitrif    =
 5.000000000000000E-002
      Bacterial Fe/C ratio                      feratb    =
 1.000000000000000E-005
      Half-saturation constant for bact. Fe/C   xkferb    =
 3.000000000000000E-010

p4z_poc_init : Initialization of remineralization parameters

ctl_opn bad opening
ctl_opn bad opening
ctl_opn bad opening
ctl_opn bad opening
ctl_opn bad opening
ctl_opn bad opening
Namelist : nampispoc
remineralisation rate of POC xremip =
3.500000000000000E-002
Number of lability classes for POC jcpoc = 15
Shape factor of the gamma distribution rshape =
1.00000000000000
ctl_opn bad opening
ctl_opn bad opening
ate constant for nanoparticle formation kfep =
1.000000000000000E-002

I have not seen the ‘ctl_opn bad opening’ error message before. If anyone is familiar with this issue, I would greatly appreciate any help. The CROCO version I am currently using is 1.3.1.

Thanks,
Lawrence

2

Did you check that you have the namelist_sediment_ref and cfg files in your run directory ? Be aware that you can not run PISCES with its sediment module without XIOS, otherwise you will not get outputs from the sediment module

3

Thank you so much for your reply Renaud. I have never seen the namelist_sediment_ref and cfg files in my run folder before. I will try to look at all the PISCES-related folders for a template. I have never used XIOS before either, so I will see if I am able to run it.