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

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

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

revision 1.50 by jmc, Tue Dec 8 21:49:23 2009 UTC revision 1.51 by jmc, Sat Jan 23 00:04:03 2010 UTC
# Line 334  C-- explicit part of the NH pressure gra Line 334  C-- explicit part of the NH pressure gra
334        IF ( use3Dsolver .AND. implicitNHPress.NE.1. _d 0 ) THEN        IF ( use3Dsolver .AND. implicitNHPress.NE.1. _d 0 ) THEN
335          nhFac = pfFacMom*(1. _d 0 - implicitNHPress)          nhFac = pfFacMom*(1. _d 0 - implicitNHPress)
336       &         *recip_deepFacC(k)*recip_rhoFacC(k)       &         *recip_deepFacC(k)*recip_rhoFacC(k)
337           IF ( exactConserv ) THEN
338          DO j=jMin,jMax          DO j=jMin,jMax
339           DO i=iMin,iMax           DO i=iMin,iMax
340            gUtmp(i,j) = gUtmp(i,j)            gUtmp(i,j) = gUtmp(i,j)
341       &               - nhFac*_recip_dxC(i,j,bi,bj)*       &               - nhFac*_recip_dxC(i,j,bi,bj)*
342       &                (phi_nh(i,j,k,bi,bj)-phi_nh(i-1,j,k,bi,bj))       &                ( (phi_nh(i,j,k,bi,bj)-phi_nh(i-1,j,k,bi,bj))
343         &                 -( dPhiNH(i,j,bi,bj) - dPhiNH(i-1,j,bi,bj) )
344         &                )
345            gVtmp(i,j) = gVtmp(i,j)            gVtmp(i,j) = gVtmp(i,j)
346       &               - nhFac*_recip_dyC(i,j,bi,bj)*       &               - nhFac*_recip_dyC(i,j,bi,bj)*
347       &                (phi_nh(i,j,k,bi,bj)-phi_nh(i,j-1,k,bi,bj))       &                ( (phi_nh(i,j,k,bi,bj)-phi_nh(i,j-1,k,bi,bj))
348         &                 -( dPhiNH(i,j,bi,bj) - dPhiNH(i,j-1,bi,bj) )
349         &                )
350           ENDDO           ENDDO
351          ENDDO          ENDDO
352           ELSE
353            DO j=jMin,jMax
354             DO i=iMin,iMax
355              gUtmp(i,j) = gUtmp(i,j)
356         &               - nhFac*_recip_dxC(i,j,bi,bj)*
357         &                  (phi_nh(i,j,k,bi,bj)-phi_nh(i-1,j,k,bi,bj))
358              gVtmp(i,j) = gVtmp(i,j)
359         &               - nhFac*_recip_dyC(i,j,bi,bj)*
360         &                  (phi_nh(i,j,k,bi,bj)-phi_nh(i,j-1,k,bi,bj))
361             ENDDO
362            ENDDO
363           ENDIF
364        ENDIF        ENDIF
365  #endif /* ALLOW_NONHYDROSTATIC */  #endif /* ALLOW_NONHYDROSTATIC */
366    
Legend:
Removed from v.1.50  
changed lines
  Added in v.1.51
  ViewVC Help
Powered by ViewVC 1.1.22