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

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

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

revision 1.64 by jmc, Tue Mar 6 16:59:44 2001 UTC revision 1.65 by jmc, Thu Mar 8 20:25:01 2001 UTC
# Line 67  C                                      v Line 67  C                                      v
67  C     rhoK, rhoKM1   - Density at current level, and level above  C     rhoK, rhoKM1   - Density at current level, and level above
68  C     phiHyd         - Hydrostatic part of the potential phiHydi.  C     phiHyd         - Hydrostatic part of the potential phiHydi.
69  C                      In z coords phiHydiHyd is the hydrostatic  C                      In z coords phiHydiHyd is the hydrostatic
70  C                      pressure anomaly  C                      Potential (=pressure/rho0) anomaly
71  C                      In p coords phiHydiHyd is the geopotential  C                      In p coords phiHydiHyd is the geopotential
72  C                      surface height  C                      surface height anomaly.
 C                      anomaly.  
73  C     phiSurfX, - gradient of Surface potentiel (Pressure/rho, ocean)  C     phiSurfX, - gradient of Surface potentiel (Pressure/rho, ocean)
74  C     phiSurfY             or geopotentiel (atmos) in X and Y direction  C     phiSurfY             or geopotentiel (atmos) in X and Y direction
75  C     KappaRT,       - Total diffusion in vertical for T and S.  C     KappaRT,       - Total diffusion in vertical for T and S.
# Line 532  C--     Start computation of dynamics Line 531  C--     Start computation of dynamics
531          jMin = 1-OLy+2          jMin = 1-OLy+2
532          jMax = sNy+OLy-1          jMax = sNy+OLy-1
533    
534  C--     Explicit part of the Surface Pressure Gradient (add in TIMESTEP)  C--     Explicit part of the Surface Potentiel Gradient (add in TIMESTEP)
535  C       (note: this loop will be replaced by CALL CALC_GRAD_ETA)  C       (note: this loop will be replaced by CALL CALC_GRAD_ETA)
536          IF (implicSurfPress.NE.1.) THEN          IF (implicSurfPress.NE.1.) THEN
537            DO j=jMin,jMax            CALL CALC_GRAD_PHI_SURF(
538              DO i=iMin,iMax       I         bi,bj,iMin,iMax,jMin,jMax,
539                phiSurfX(i,j) = _recip_dxC(i,j,bi,bj)*gBaro       I         etaN,
540       &           *(etaN(i,j,bi,bj)-etaN(i-1,j,bi,bj))       O         phiSurfX,phiSurfY,
541                phiSurfY(i,j) = _recip_dyC(i,j,bi,bj)*gBaro       I         myThid )                        
      &           *(etaN(i,j,bi,bj)-etaN(i,j-1,bi,bj))  
             ENDDO  
           ENDDO  
542          ENDIF          ENDIF
543    
544  C--     Start of dynamics loop  C--     Start of dynamics loop
Legend:
Removed from v.1.64  
changed lines
  Added in v.1.65
  ViewVC Help
Powered by ViewVC 1.1.22