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

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

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

revision 1.18 by cnh, Sun Sep 6 14:45:11 1998 UTC revision 1.19 by cnh, Sun Sep 6 17:35:20 1998 UTC
# Line 82  C     aS2d: integral in Z Ay/dY Line 82  C     aS2d: integral in Z Ay/dY
82         cg2dNorm = myNorm         cg2dNorm = myNorm
83        _BEGIN_MASTER( myThid )        _BEGIN_MASTER( myThid )
84  CcnhDebugStarts  CcnhDebugStarts
85         WRITE(msgBuf,'(A,F)') '// CG2D normalisation factor = ', cg2dNorm         WRITE(msgBuf,'(A,E40.25)') '// CG2D normalisation factor = ', cg2dNorm
86         CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)         CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
87         WRITE(msgBuf,*) '                               '         WRITE(msgBuf,*) '                               '
88         CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)         CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
# Line 133  C           defaults to 0.51 but can be Line 133  C           defaults to 0.51 but can be
133            aC = -(            aC = -(
134       &     aW2d(I,J,bi,bj) + aW2d(I+1,J  ,bi,bj)       &     aW2d(I,J,bi,bj) + aW2d(I+1,J  ,bi,bj)
135       &    +aS2d(I,J,bi,bj) + aS2D(I  ,J+1,bi,bj)       &    +aS2d(I,J,bi,bj) + aS2D(I  ,J+1,bi,bj)
136       &    +freeSurfFac*myNorm*  Gravity*rhoConst*       &    +freeSurfFac*myNorm* horiVertRatio*
137       &     rA(I,J,bi,bj)/deltaTMom/deltaTMom       &     rA(I,J,bi,bj)/deltaTMom/deltaTMom
138       &    )       &    )
139            aCs = -(            aCs = -(
140       &     aW2d(I,J-1,bi,bj) + aW2d(I+1,J-1,bi,bj)       &     aW2d(I,J-1,bi,bj) + aW2d(I+1,J-1,bi,bj)
141       &    +aS2d(I,J-1,bi,bj) + aS2d(I  ,J  ,bi,bj)       &    +aS2d(I,J-1,bi,bj) + aS2d(I  ,J  ,bi,bj)
142       &    +freeSurfFac*myNorm*  Gravity*rhoConst*       &    +freeSurfFac*myNorm* horiVertRatio*
143       &     rA(I,J-1,bi,bj)/deltaTMom/deltaTMom       &     rA(I,J-1,bi,bj)/deltaTMom/deltaTMom
144       &    )       &    )
145            aCw = -(            aCw = -(
146       &     aW2d(I-1,J,bi,bj) + aW2d(I  ,J  ,bi,bj)       &     aW2d(I-1,J,bi,bj) + aW2d(I  ,J  ,bi,bj)
147       &    +aS2d(I-1,J,bi,bj) + aS2d(I-1,J+1,bi,bj)       &    +aS2d(I-1,J,bi,bj) + aS2d(I-1,J+1,bi,bj)
148       &    +freeSurfFac*myNorm*  Gravity*rhoConst*       &    +freeSurfFac*myNorm* horiVertRatio*
149       &     rA(I-1,J,bi,bj)/deltaTMom/deltaTMom       &     rA(I-1,J,bi,bj)/deltaTMom/deltaTMom
150       &    )       &    )
151            IF ( aC .EQ. 0. ) THEN            IF ( aC .EQ. 0. ) THEN
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22