/[MITgcm]/MITgcm/pkg/mom_common/mom_quasihydrostatic.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_common/mom_quasihydrostatic.F

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

revision 1.3 by jmc, Tue Dec 5 05:27:38 2006 UTC revision 1.4 by jmc, Mon Feb 5 03:16:23 2007 UTC
# Line 47  C  effectiveBuoy :: Density anomaly Line 47  C  effectiveBuoy :: Density anomaly
47  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
48  C  i,j           :: loop indices  C  i,j           :: loop indices
49  C  scalingFactor :: scaling factor (from acceleration to density)  C  scalingFactor :: scaling factor (from acceleration to density)
50          INTEGER iMin,iMax,jMin,jMax
51          PARAMETER( iMin = 0 , iMax = sNx+1 )
52          PARAMETER( jMin = 0 , jMax = sNy+1 )
53        INTEGER i,j        INTEGER i,j
54        _RL scalingFactor        _RL scalingFactor
55  CEOP  CEOP
# Line 56  C-    to reproduce old (wrong) results: Line 59  C-    to reproduce old (wrong) results:
59  c     scalingFactor=gravitySign*recip_gravity  c     scalingFactor=gravitySign*recip_gravity
60    
61        IF ( use3dCoriolis ) THEN        IF ( use3dCoriolis ) THEN
62         DO j=1-Oly,sNy+Oly-1         DO j=jMin,jMax
63          DO i=1-Olx,sNx+Olx-1          DO i=iMin,iMax
64           effectiveBuoy(i,j)=effectiveBuoy(i,j)           effectiveBuoy(i,j)=effectiveBuoy(i,j)
65       &    +scalingFactor*       &    +scalingFactor*
66       &     fCoriCos(i,j,bi,bj)*       &     fCoriCos(i,j,bi,bj)*
# Line 71  c     scalingFactor=gravitySign*recip_gr Line 74  c     scalingFactor=gravitySign*recip_gr
74        ENDIF        ENDIF
75    
76        IF ( useNHMTerms ) THEN        IF ( useNHMTerms ) THEN
77         DO j=1-Oly,sNy+Oly-1         DO j=jMin,jMax
78          DO i=1-Olx,sNx+Olx-1          DO i=iMin,iMax
79           effectiveBuoy(i,j)=effectiveBuoy(i,j)           effectiveBuoy(i,j)=effectiveBuoy(i,j)
80       &    +scalingFactor*       &    +scalingFactor*
81       &     (   (uFld( i ,j,k,bi,bj)*uFld( i ,j,k,bi,bj)       &     (   (uFld( i ,j,k,bi,bj)*uFld( i ,j,k,bi,bj)
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22