[MITgcm-support] Trying to relax ptracer using RBCS

[kunal madkaiker]

Dear All,

I need some help with the RBCS package. I am trying to restore the
'nitrate' ptracer at the surface (I am using BLING). Since its just one
mask to be read, I have set maskLEN=1 in RBCS_SIZE.h. So, my iTr=4 and irbc
should be 6 (2+iTr) if I am not mistaken.

However, when I run the setup, I get an error: Fortran runtime error:
Non-existing record number.
After keeping maskLEN> irbc, I still get the same error for
'coastal_mask.bin' file.

My data.rbcs looks like this,
# RBCS package parameters:
 &RBCS_PARM01
   relaxMaskFile(6)='coastal_mask.bin',
   rbcsForcingPeriod=432000.,
   rbcsForcingCycle=31536000.,
 &

# RBCS for pTracers (read this namelist only when ptracers pkg is compiled)
 &RBCS_PARM02
  useRBCptrnum(4)    = .TRUE.,
  tauRelaxPTR(4)     = 432000.,
  relaxPtracerFile(4)='nitrate_surf_pentadclim.bin',
 &

Kindly let me know your thoughts and any modifications that I need to make.

Thanks,
Kunal
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20231215/16793d2f/attachment.html>