[MITgcm-support] obcs_apply_ptracer

[Martin Losch]

Taka,

the tracers need to be defined at the tracer points, but whether you  
have to do anything at an open boundary or is determined by S/W  
masks. For example:
> C Northern boundary
> #ifdef ALLOW_OBCS_NORTH
>        J_obc = OB_Jn(I,bi,bj)
>        IF (J_obc.NE.0) THEN
>         obc_mask = _maskS(I,J_obc,K,bi,bj)
>         DO J = J_obc, J_obc+Oly
>          pFld(I,J)=OBNptr(I,K,bi,bj,iTracer)*obc_mask
> CML        pFld(I,J)=pFld(I,J_obc-1)*obc_mask
>         ENDDO
>        ENDIF
> #endif
if (i,j_obs) is an open boundary point, then (i,j_obs-1) is the first  
interior point. Only if both points are "wet" then maskS(i,j_obs) >  
0; you can do this check too with maskC, but you'd have to check if  
maskC(i,j_obs) and maskC(i,j_obs-1) are > 0.

Martin



On 9 Apr 2008, at 05:30, ito at atmos.colostate.edu wrote:
> Hello,
>
> I've been doing offline tracer simulations with OBCS and realized  
> that S/R
> obcs_apply_ptracer (similar to obcs_apply_ts) uses obc_mask based  
> on maskW
> and maskS.  Should the open boundary tracer properties be defined  
> at the
> grid faces?  Or should it be defined at the tracer point in the  
> center of
> grid cell (in which case we would have to use maskC)?
>
> Taka
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support