/[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.19 by cnh, Sun Sep 6 17:35:20 1998 UTC revision 1.20 by cnh, Wed Oct 28 03:11:37 1998 UTC
# Line 55  C     aS2d: integral in Z Ay/dY Line 55  C     aS2d: integral in Z Ay/dY
55          DO K=1,Nr          DO K=1,Nr
56           DO J=1,sNy           DO J=1,sNy
57            DO I=1,sNx            DO I=1,sNx
58             faceArea = _dyG(I,J,bi,bj)*drF(K)*_hFacW(I,J,K,bi,bj)             faceArea = _dyG(I,J,bi,bj)*drF(K)
59         &               *_hFacW(I,J,K,bi,bj)
60             aW2d(I,J,bi,bj) = aW2d(I,J,bi,bj) +             aW2d(I,J,bi,bj) = aW2d(I,J,bi,bj) +
61       &      gBaro*faceArea*recip_dxC(I,J,bi,bj)       &      gBaro*faceArea*recip_dxC(I,J,bi,bj)
62             faceArea = _dxG(I,J,bi,bj)*drF(K)*_hFacS(I,J,K,bi,bj)             faceArea = _dxG(I,J,bi,bj)*drF(K)
63         &               *_hFacS(I,J,K,bi,bj)
64             aS2d(I,J,bi,bj) = aS2d(I,J,bi,bj) +             aS2d(I,J,bi,bj) = aS2d(I,J,bi,bj) +
65       &      gBaro*faceArea*recip_dyC(I,J,bi,bj)       &      gBaro*faceArea*recip_dyC(I,J,bi,bj)
66            ENDDO            ENDDO
# Line 82  C     aS2d: integral in Z Ay/dY Line 84  C     aS2d: integral in Z Ay/dY
84         cg2dNorm = myNorm         cg2dNorm = myNorm
85        _BEGIN_MASTER( myThid )        _BEGIN_MASTER( myThid )
86  CcnhDebugStarts  CcnhDebugStarts
87         WRITE(msgBuf,'(A,E40.25)') '// CG2D normalisation factor = ', cg2dNorm         WRITE(msgBuf,'(A,E40.25)') '// CG2D normalisation factor = ',
88         &               cg2dNorm
89         CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)         CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
90         WRITE(msgBuf,*) '                               '         WRITE(msgBuf,*) '                               '
91         CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)         CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
Legend:
Removed from v.1.19  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22