[Mitgcm-support] Re: free slip/no slip

mitgcm-support at dev.mitgcm.org mitgcm-support at dev.mitgcm.org
Wed Jul 9 15:48:14 EDT 2003


Hi Alistair,

thank for the fast answer. We will make another test with free/no slip
on Monday. 

Another question: Is the shaved cell option still in the mitgcmuv code?

Arne

Alistair Adcroft wrote:
> 
> Hi Arne,
> 
> C     Eddy component of vertical flux
>       IF (no_slip_bottom) THEN
>        DO j=jMin,jMax
>         DO i=iMin,iMax
>          vf(i,j) =
>      &   -KappaRU(i,j,k)
>      & *_RA_AT_W
>      &   *( ( 2.-_maskW(i,j,k,bi,bj) )*uVel(i,j,kM1,bi,bj)
>      &     -     _maskW(i,j,k,bi,bj)  *uVel(i,j,k,bi,bj)
>      &    )*rkFac*recip_drC(k)
>         ENDDO
>        ENDDO
>       ELSE
>        DO j=jMin,jMax
>         DO i=iMin,iMax
>          vf(i,j) =
>      &   -KappaRU(i,j,k)
>      & *_RA_AT_W
>      &   *( uVel(i,j,km1,bi,bj)-uVel(i,j,k,bi,bj)
>      &    )*rkFac*recip_drC(k)
>      &    *_maskW(i,j,k,bi,bj)*rVelMaskOverride
>         ENDDO
>        ENDDO
>       ENDIF
> 
> no_slip_bottom=.TRUE., in interior maskW(k)=1
> so  ( 2-maskW(k) ) U(k-1) - maskW(k) U(k)
>     =  1 U(k-1) - 1 U(k)
> 
> no_slip_bottom=.TRUE., at bottom maskW(k)=0
> so  ( 2-maskW(k) ) U(k-1) - maskW(k) U(k)
>     =  2 U(k-1) - 0 U(k)
> 
> no_slip_bottom=.FALSE., in interior maskW(k)=1
> so  maskW(k) [ U(k-1) - maskW(k) U(k) ]
>     =  1 [ U(k-1) - 1 U(k) ]
> 
> no_slip_bottom=.FALSE., at bottom maskW(k)=0
> so  maskW(k) [ U(k-1) - maskW(k) U(k) ]
>     =  0 [ U(k-1) - 1 U(k) ]
> 
> Looks right to me. This is the sort of code
> that doesn't get tested properly (beyond me thinking
> about it) so if you're convinced it's wrong
> show me a test.
> 
> Regards,
> 
> Alistair.
> PS> I would dearly love to tidy up calc_mom_rhs(). It
> takes much more CPU than it should mostly because
> of using these intermediate arrays like vf, af and df.
> One day ...

-- 

Dr. Arne Biastoch
Scripps Institution of Oceanography        phone: +1-858-822-3787
Physical Oceanography Research Division    fax  : +1-858-534-9820
MS: 0230                                   email: abiastoch at ucsd.edu
8605 La Jolla Shores Dr.            
La Jolla, CA 92093-0230, U.S.A.   http://www.ecco.ucsd.edu/~biastoch




More information about the MITgcm-support mailing list