[MITgcm-support] How does MITgcm handle the 2D problem

q li qliuri at yahoo.com
Fri Dec 7 20:46:34 EST 2012


I got it. Thank you very much, Jean-Michel.

Li


________________________________
 From: Jean-Michel Campin <jmc at ocean.mit.edu>
To: mitgcm-support at mitgcm.org 
Sent: Friday, December 7, 2012 9:10 PM
Subject: Re: [MITgcm-support] How does MITgcm handle the 2D problem
 
Hi Li,

Yes, no "boundary condition" in the missing direction.
the 1 grid point value is copied over the full over regions (1-Olx:0
and sNx+1:sNx+Olx) so that all gradient along x will vanish.

Cheers,
Jean-Michel

On Fri, Dec 07, 2012 at 01:10:55AM -0800, q li wrote:
> Hi Jean-Michel,
> 
> Thanks for the quick reply.
> 
> Does this imply that no boundary condition  is needed in the missing direction? Or a pseudo-periodic condition?
> 
> Best regards,
> 
> Li
> 
> 
> ________________________________
>  From: Jean-Michel Campin <jmc at ocean.mit.edu>
> To: mitgcm-support at mitgcm.org 
> Sent: Friday, December 7, 2012 1:11 PM
> Subject: Re: [MITgcm-support] How does MITgcm handle the 2D problem
>  
> Hi Li,
> 
> The way a 2-D (x-z or y-z) set-up works is quiet simple, 
> there is just 1 grid point in the missing direction (respectively, y or x),
> still periodic along this direction, therefore all gradients
> are zero; there is no boundary in the missing direction, 
> and the 2 horizontal components of the velocity are stepped forward
> (so that coriolis effects are accounted for).
> 
> Does this answers your question ?
> 
> Cheers,
> Jean-Michel
> 
> On Thu, Dec 06, 2012 at 05:09:34PM -0800, q li wrote:
> > Hi,
> > 
> > I'm confusing about how the MITgcm handles the 2D runs, after receiving a comment from a paper reviewer.
> > 
> > For example, in the "verification/internal_wave" case, tidal forcing is in east-west direction.  delY=5.E3 in the "data" file means only 1 grid in y-direction?
> > On the north/south boundary, how come u and v are set to be zero? Does it mean u and v have to disappear in the whole domain?
> > 
> > Could anyone give me some hints?
> > 
> > Thanks,
> > 
> > Li
> > 
> > "verification/internal_wave/code/obcs_calc.F"
> > 
> > 
> > C         Northern OB, template for forcing
> >       IF (useOrlanskiNorth) THEN
> >         CALL ORLANSKI_NORTH(
> >      &          bi, bj, futureTime,
> >      &          uVel, vVel, wVel, theta, salt,
> >      &          myThid )
> >       ELSE
> >         DO K=1,Nr
> >           DO I=1-Olx,sNx+Olx
> >             OBNv(I,K,bi,bj)=0.
> >             OBNu(I,K,bi,bj)=0.
> >             OBNt(I,K,bi,bj)=tRef(K)
> >             OBNs(I,K,bi,bj)=sRef(K)
> > #ifdef ALLOW_NONHYDROSTATIC
> >             OBNw(I,K,bi,bj)=0.
> > #endif
> >           ENDDO
> >         ENDDO
> >       ENDIF
> 
> > _______________________________________________
> > 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

> _______________________________________________
> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20121207/4f7ff42c/attachment.htm>


More information about the MITgcm-support mailing list