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

[q li]

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