Missing block in PISCESv2 namelist_sediment_cfg and _ref in CROCO v1.3
@x-JMcGov encountered an error at runtime when switching to CROCO v1.3 where the namelist_sediment default file (as included in /PISCES/SED) was missing the following block at the very end (i.e. the last block that was previously available to read from the namelist file was nam_rst)
!----------------------------------------------------------------------- &nam_output ! parameters for outputing the sediment module !----------------------------------------------------------------------- ldefsedpis_avg = .true. ! write averaged output variables cn_sedwri_out = "output_sed.nc" ! name of the input restart file name of the sediment module nrpfsedpis_avg = 0 ! Frequency of the averaged outputs nwrtsedpis_avg = 24 ! Frequency of the averaged outputs ntssedpis_avg = 1 ! ??? /
@x-JMcGov found the missing block correctly formatted here (for NEMO-PISCES): http://forge.ipsl.jussieu.fr/nemo/changeset/14306#file2
The corresponding block of text is read by PISCES/SED/sedini.F90 in line 425
@x-JMcGov became aware of this fault:
- by switching on key_sediment for PISCES, and
- encountering a runtime error as follows: croco1 0000000000C8B0BE setavg_sed_mp_set 2149 setavg_sed_.f90
where the expression (mod(ilc-1,nwrtsedpis_avg).eq.1) was returning a "division by zero" error as nwrtsedpis_avg had not been set and thus had been defaulted to zero.