/[MITgcm]/MITgcm/model/src/timestep_wvel.F
ViewVC logotype

Diff of /MITgcm/model/src/timestep_wvel.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.4 by jmc, Mon Mar 12 23:43:54 2007 UTC revision 1.5 by jmc, Sun Nov 29 03:12:32 2009 UTC
# Line 65  C---+----1----+----2----+----3----+----4 Line 65  C---+----1----+----2----+----3----+----4
65    
66          DO k=2,Nr          DO k=2,Nr
67    
68             IF ( implicitNHPress.NE.1. _d 0 ) THEN
69    C--   add explicit part of NH pressure gradient:
70              tmpFac = pfFacMom*(1. _d 0 - implicitNHPress)
71         &                     *wUnit2rVel(k)*wUnit2rVel(k)
72         &                     *recip_rhoFacF(k)
73              DO j=jMin,jMax
74               DO i=iMin,iMax
75                gW(i,j,k,bi,bj) = gW(i,j,k,bi,bj)
76         &             - tmpFac*rkSign*recip_drC(k)
77         &             *( phi_nh(i,j,k,bi,bj) - phi_nh(i,j,k-1,bi,bj) )
78               ENDDO
79              ENDDO
80             ENDIF
81  C     apply mask to gW and keep a copy of wVel in gW:  C     apply mask to gW and keep a copy of wVel in gW:
82           DO j=1-Oly,sNy+Oly           DO j=1-Oly,sNy+Oly
83            DO i=1-Olx,sNx+Olx            DO i=1-Olx,sNx+Olx
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22