[MITgcm-support] RBCS forcing period spec

[Jeffery R Scott]

Thanks Martin.

I forgot that I need to specify a 3D file for restoring using RBCS.

Jeff


________________________________________
From: MITgcm-support <mitgcm-support-bounces at mitgcm.org> on behalf of Martin Losch <Martin.Losch at awi.de>
Sent: Tuesday, June 22, 2021 12:31 PM
To: MITgcm Support
Subject: Re: [MITgcm-support] RBCS forcing period spec

Hi Jeff,

from reading the code, this should work.

Did you try to generate debug-messages with debugLevel .GE. debugLevB (2?)? Lines starting with “RBCS_FIELDS_LOAD”? That should give you the record numbers that the model wants to read. Maybe that can give you a clue.

What about the “readBinaryPrec”? Is that being used in the same way as for the files in “data”?

Martin

> On 22. Jun 2021, at 18:20, Jeffery R Scott <jscott at mit.edu> wrote:
>
> Hi,
>
> I'm trying to shift a repeating annual cycle of (monthly mean data) sss restoring as specified in file data (where it works) to data.rbcs.
>
> In data:
>  periodicExternalForcing=.TRUE.,
>  externForcingPeriod=2592000.,
>  externForcingCycle=31104000.,
>  saltClimFile=   'sss25_128x60x12.bin',
>
>
> In rbcs, following the manual example  8.3.2.4.1. Cyclic data in a single file
>
>  rbcsForcingPeriod =2592000.,
>  rbcsForcingCycle=31104000.,
>  relaxSFile='sss25_128x60x12.bin',
>
> (with default settings rbcsSingleTimeFiles FALSE and ForcingOffset=0)
>
> but I get this in STDERR
> (PID.TID 0000.0001) *** ERROR ***  RBCS_READPARAMS: startTime before rbcsForcingOffset+0.5*rbcsForcingPeriod
> (PID.TID 0000.0001) *** ERROR ***  RBCS_READPARAMS: will use last record
>
> and then a crash attempting to access non-existing record. I would think it should average record 12 and record 1 at the start of the simulation
> (nIter0 = 0) which is what I believe is done when restoring through data (?)
>
> Either I'm misunderstanding something, or the manual is not correct, or I have something wrong somewhere (?)
>
> Thanks,
> Jeff
>
>
>
>
>
>
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support