[MITgcm-support] RBCS relax u velocities problem

[Dustin Carroll]

Hi Martin,

Thanks for the reply. Good question, I am relaxing the plume sufficiently far away from the boundaries or a land mask. 
 
Are there any diagnostics for the rbcs pkg so I can see if the mask are being loaded properly, etc?

Cheers,
Dustin

On Jan 15, 2014, at 9:40 AM, Martin Losch <Martin.Losch at awi.de> wrote:

> Hi Dustin,
> 
> does it make sense to use the same mask for U as for T/S ? The U points are staggered by half a grid cell on a C-grid and if you are only restoring near the boundaries this might lead to unexpected results (i.e. relaxing the velocites on the boundaries, that are always masked to zero)?
> 
> Martin
> 
> On Jan 15, 2014, at 11:50 AM, Dustin Carroll <dcarroll at uoregon.edu> wrote:
> 
>> 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
>> 
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support