[MITgcm-support] RBCS relax u velocities problem

[Dustin Carroll]

Hi All,

I am trying to use the RBCS package to create a plume in a 3-D model domain. The T,S, and tracer fields are relaxing properly, but the velocities are not relaxing. 

I looked through the RBCS pkg code and it looks this functionality is supported, are there any special flags / settings I need to get this work?

I am setting:

#undef DISABLE_RBCS_MOM
#define ALLOW_RBCS

in my CPP_OPTIONS.h file.


In my STDOUT.0000 file I see that the rbcsU.bin files are loading. It only shows the loading of two masks though.

(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // End MONITOR EXF statistics
(PID.TID 0000.0001) // =======================================================
RBCS_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=    1    1 (prev=    1    0 )
(PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: rbcsT.bin
(PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: rbcsT.bin
(PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: rbcsS.bin
(PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: rbcsS.bin
(PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: rbcsMask.bin
(PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: rbcsMask.bin
RBCS_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=    1    1 (prev=    1    0 )
(PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: rbcsU.bin
(PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: rbcsU.bin
RBCS_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=    1    1 (prev=    1    0 )
(PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: rbcsU.bin
(PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: rbcsU.bin
RBCS_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=    1    1 (prev=    1    0 )
RBCS_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=    1    1 (prev=    1    0 )

Here is my data.rbcs file, I am using the same mask for T,S, ptracer and u.

&RBCS_PARM01
   tauRelaxT=86400.,
   tauRelaxS=86400.,
   tauRelaxU=86400.,

   rbcsForcingUPeriod=0.,
   rbcsForcingUCycle=0.,
   rbcsForcingUOffset=0.,

   relaxMaskFile(1)='rbcsMask.bin',
   relaxMaskFile(2)='rbcsMask.bin',
   relaxMaskFile(3)='rbcsMask.bin',    

   relaxTFile='rbcsT.bin',
   relaxSFile='rbcsS.bin',
   relaxUFile='rbcsU.bin',

   relaxMaskUFile='rbcsMask.bin',

   useRBCtemp=.TRUE.,
   useRBCsalt=.TRUE.,
   useRBCuVel=.TRUE.,

   useRBCptracers=.TRUE.,
   rbcsIniter=0.,
   rbcsForcingPeriod=0.,
   rbcsForcingCycle=0.,
&

 &RBCS_PARM02

 useRBCptrnum(1)    = .TRUE.,
 tauRelaxPTR(1)     = 1200.,
 relaxPtracerFile(1)='rbcsMask.bin',

 &

Thanks for the help,
Dustin Carroll