[MITgcm-support] Trying to relax ptracer using RBCS

[kunal madkaiker]

Hi Martin,
Sorry for the confusion. I was trying with a different file but the error I
get is the same.

 RBCS_FIELDS_LOAD,   1051200 : iP,iLd,i0,i1=  365    0  365    1 ; Wght=
 0.5000000000  0.5000000000
 RBCS_FIELDS_LOAD, it=   1051200 : Reading new data, i0,i1=  365    1
(prev=  365    0 )
(PID.TID 0000.0001)  MDS_READ_FIELD: opening global file:
nitrate_surf_J2D.bin

I get the error after this. I am referring to
https://github.com/MITgcm/MITgcm/blob/9af873c53257354cb8e18bb8591a3ee68c6ae679/verification/exp4/input/data.rbcs
Can you help finding the dimensions size of 'relaxMaskFile' and
'relaxPtracerFile'. It'll be easier to feed the files accordingly.

Thanks,
Kunal



On Mon, 18 Dec 2023, 21:40 Martin Losch, <Martin.Losch at awi.de> wrote:

> It looks like, it doesn’t have anything to do with pkg/rbcs, as
> ‘nitrate_surf_J2D.bin’ does not appear in your data.rbcs.
>
> It’s a different file in a different package (pkg/bling?)
>
> M.
>
> On 18. Dec 2023, at 15:27, kunal madkaiker <kunal.madkaiker02 at gmail.com>
> wrote:
>
> 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
>>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20231219/45067439/attachment.html>