[MITgcm-support] RBCS relax u velocities problem

[Martin Losch]

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