/[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.112 by heimbach, Mon Jan 24 17:00:17 2005 UTC revision 1.113 by jmc, Fri Jan 28 01:00:13 2005 UTC
# Line 108  C      |-- IMPLDIFF Line 108  C      |-- IMPLDIFF
108  C      |  C      |
109  C      |-- OBCS_APPLY_UV  C      |-- OBCS_APPLY_UV
110  C      |  C      |
 C      |-- CALL TIMEAVE_CUMUL_1T  
111  C      |-- CALL DEBUG_STATS_RL  C      |-- CALL DEBUG_STATS_RL
112    
113  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
# Line 122  C     myThid - Thread number for this in Line 121  C     myThid - Thread number for this in
121    
122  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
123  C     == Local variables  C     == Local variables
124  C     fVer[STUV]               o fVer: Vertical flux term - note fVer  C     fVer[UV]               o fVer: Vertical flux term - note fVer
125  C                                      is "pipelined" in the vertical  C                                    is "pipelined" in the vertical
126  C                                      so we need an fVer for each  C                                    so we need an fVer for each
127  C                                      variable.  C                                    variable.
128  C     phiHydC    :: hydrostatic potential anomaly at cell center  C     phiHydC    :: hydrostatic potential anomaly at cell center
129  C                   In z coords phiHyd is the hydrostatic potential  C                   In z coords phiHyd is the hydrostatic potential
130  C                      (=pressure/rho0) anomaly  C                      (=pressure/rho0) anomaly
# Line 163  C                      index into fVerTe Line 162  C                      index into fVerTe
162    
163        LOGICAL  DIFFERENT_MULTIPLE        LOGICAL  DIFFERENT_MULTIPLE
164        EXTERNAL DIFFERENT_MULTIPLE        EXTERNAL DIFFERENT_MULTIPLE
165    
166    #ifdef ALLOW_DIAGNOSTICS
167          _RL tmpFld  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
168          LOGICAL  DIAGNOSTICS_IS_ON
169          EXTERNAL DIAGNOSTICS_IS_ON
170    #endif /* ALLOW_DIAGNOSTICS */
171    
172    
173  C---    The algorithm...  C---    The algorithm...
174  C  C
# Line 471  C--     End implicit Vertical advection Line 477  C--     End implicit Vertical advection
477        ENDIF        ENDIF
478  #endif  #endif
479    
480    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
481    
482  Cml(  Cml(
483  C     In order to compare the variance of phiHydLow of a p/z-coordinate  C     In order to compare the variance of phiHydLow of a p/z-coordinate
484  C     run with etaH of a z/p-coordinate run the drift of phiHydLow  C     run with etaH of a z/p-coordinate run the drift of phiHydLow
# Line 479  C      CALL REMOVE_MEAN_RL( 1, phiHydLow Line 487  C      CALL REMOVE_MEAN_RL( 1, phiHydLow
487  C     &                'phiHydLow', myThid )  C     &                'phiHydLow', myThid )
488  Cml)  Cml)
489    
490    #ifdef ALLOW_DIAGNOSTICS
491          IF ( usediagnostics ) THEN
492    
493           CALL DIAGNOSTICS_FILL(totPhihyd,'PHIHYD  ',0,Nr,0,1,1,myThid)
494           CALL DIAGNOSTICS_FILL(phiHydLow,'PHIBOT  ',0,1,0,1,1,myThid)
495    
496           IF ( DIAGNOSTICS_IS_ON('PHIBOTSQ',myThid) ) THEN
497            DO bj = myByLo(myThid), myByHi(myThid)
498             DO bi = myBxLo(myThid), myBxHi(myThid)
499              DO j = 1,sNy
500               DO i = 1,sNx
501                 tmpFld(i,j) = phiHydLow(i,j,bi,bj)*phiHydLow(i,j,bi,bj)
502               ENDDO
503              ENDDO
504              CALL DIAGNOSTICS_FILL(tmpFld,'PHIBOTSQ',0,1,2,bi,bj,myThid)
505             ENDDO
506            ENDDO
507           ENDIF
508    
509          ENDIF
510    #endif /* ALLOW_DIAGNOSTICS */
511          
512  #ifdef ALLOW_DEBUG  #ifdef ALLOW_DEBUG
513        If ( debugLevel .GE. debLevB ) THEN        If ( debugLevel .GE. debLevB ) THEN
514         CALL DEBUG_STATS_RL(1,EtaN,'EtaN (DYNAMICS)',myThid)         CALL DEBUG_STATS_RL(1,EtaN,'EtaN (DYNAMICS)',myThid)
Legend:
Removed from v.1.112  
changed lines
  Added in v.1.113
  ViewVC Help
Powered by ViewVC 1.1.22