[MITgcm-support] Trying to relax ptracer using RBCS

[kunal madkaiker]

Dear All,

Just to add more information to my previous query, both the relaxMaskFile
and relaxPtracerFile have dimensions X*Y*73 i.e., 73 timesteps as both are
being repeated every fifth day.
I also tried changing it to daily (X*Y*365) but still get the same error:
At line 2203 of file mdsio_read_field.f (unit = 9, file =
'nitrate_surf_J2D.bin')
Fortran runtime error: Non-existing record number

Do I need to add an extra timestep (73+1)? Can anyone help me find the
actual dimensions required for 'relaxMaskFile' and 'relaxPtracerFile'.

Thanks,
Kunal

On Fri, Dec 15, 2023 at 3:11 PM kunal madkaiker <kunal.madkaiker02 at gmail.com>
wrote:

> 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/20231218/b90b7d8e/attachment.html>