/[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.117 by molod, Mon May 23 20:49:37 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 161  C                      index into fVerTe Line 160  C                      index into fVerTe
160        INTEGER i, j        INTEGER i, j
161        INTEGER k, km1, kp1, kup, kDown        INTEGER k, km1, kp1, kup, kDown
162    
163        LOGICAL  DIFFERENT_MULTIPLE  #ifdef ALLOW_DIAGNOSTICS
164        EXTERNAL DIFFERENT_MULTIPLE        _RL tmpFld  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
165          LOGICAL  DIAGNOSTICS_IS_ON
166          EXTERNAL DIAGNOSTICS_IS_ON
167    #endif /* ALLOW_DIAGNOSTICS */
168    
169    
170  C---    The algorithm...  C---    The algorithm...
171  C  C
# Line 471  C--     End implicit Vertical advection Line 474  C--     End implicit Vertical advection
474        ENDIF        ENDIF
475  #endif  #endif
476    
477    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
478    
479  Cml(  Cml(
480  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
481  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 484  C      CALL REMOVE_MEAN_RL( 1, phiHydLow
484  C     &                'phiHydLow', myThid )  C     &                'phiHydLow', myThid )
485  Cml)  Cml)
486    
487    #ifdef ALLOW_DIAGNOSTICS
488          IF ( usediagnostics ) THEN
489    
490           CALL DIAGNOSTICS_FILL(totPhihyd,'PHIHYD  ',0,Nr,0,1,1,myThid)
491    
492           IF ( DIAGNOSTICS_IS_ON('PHIHYDSQ',myThid) ) THEN
493            DO bj = myByLo(myThid), myByHi(myThid)
494            DO bi = myBxLo(myThid), myBxHi(myThid)
495             DO k = 1,Nr
496              DO j = 1,sNy
497               DO i = 1,sNx
498                 tmpFld(i,j) = totPhihyd(i,j,k,bi,bj)*totPhihyd(i,j,k,bi,bj)
499               ENDDO
500              ENDDO
501              CALL DIAGNOSTICS_FILL(tmpFld,'PHIHYDSQ',k,1,2,bi,bj,myThid)
502             ENDDO
503            ENDDO
504            ENDDO
505           ENDIF
506    
507           CALL DIAGNOSTICS_FILL(phiHydLow,'PHIBOT  ',0,1,0,1,1,myThid)
508    
509           IF ( DIAGNOSTICS_IS_ON('PHIBOTSQ',myThid) ) THEN
510            DO bj = myByLo(myThid), myByHi(myThid)
511             DO bi = myBxLo(myThid), myBxHi(myThid)
512              DO j = 1,sNy
513               DO i = 1,sNx
514                 tmpFld(i,j) = phiHydLow(i,j,bi,bj)*phiHydLow(i,j,bi,bj)
515               ENDDO
516              ENDDO
517              CALL DIAGNOSTICS_FILL(tmpFld,'PHIBOTSQ',0,1,2,bi,bj,myThid)
518             ENDDO
519            ENDDO
520           ENDIF
521    
522          ENDIF
523    #endif /* ALLOW_DIAGNOSTICS */
524          
525  #ifdef ALLOW_DEBUG  #ifdef ALLOW_DEBUG
526        If ( debugLevel .GE. debLevB ) THEN        If ( debugLevel .GE. debLevB ) THEN
527         CALL DEBUG_STATS_RL(1,EtaN,'EtaN (DYNAMICS)',myThid)         CALL DEBUG_STATS_RL(1,EtaN,'EtaN (DYNAMICS)',myThid)
# Line 487  Cml) Line 530  Cml)
530         CALL DEBUG_STATS_RL(Nr,wVel,'Wvel (DYNAMICS)',myThid)         CALL DEBUG_STATS_RL(Nr,wVel,'Wvel (DYNAMICS)',myThid)
531         CALL DEBUG_STATS_RL(Nr,theta,'Theta (DYNAMICS)',myThid)         CALL DEBUG_STATS_RL(Nr,theta,'Theta (DYNAMICS)',myThid)
532         CALL DEBUG_STATS_RL(Nr,salt,'Salt (DYNAMICS)',myThid)         CALL DEBUG_STATS_RL(Nr,salt,'Salt (DYNAMICS)',myThid)
533         CALL DEBUG_STATS_RL(Nr,Gu,'Gu (DYNAMICS)',myThid)         CALL DEBUG_STATS_RL(Nr,gU,'Gu (DYNAMICS)',myThid)
534         CALL DEBUG_STATS_RL(Nr,Gv,'Gv (DYNAMICS)',myThid)         CALL DEBUG_STATS_RL(Nr,gV,'Gv (DYNAMICS)',myThid)
535         CALL DEBUG_STATS_RL(Nr,Gt,'Gt (DYNAMICS)',myThid)         CALL DEBUG_STATS_RL(Nr,gT,'Gt (DYNAMICS)',myThid)
536         CALL DEBUG_STATS_RL(Nr,Gs,'Gs (DYNAMICS)',myThid)         CALL DEBUG_STATS_RL(Nr,gS,'Gs (DYNAMICS)',myThid)
537         CALL DEBUG_STATS_RL(Nr,GuNm1,'GuNm1 (DYNAMICS)',myThid)  #ifndef ALLOW_ADAMSBASHFORTH_3
538         CALL DEBUG_STATS_RL(Nr,GvNm1,'GvNm1 (DYNAMICS)',myThid)         CALL DEBUG_STATS_RL(Nr,guNm1,'GuNm1 (DYNAMICS)',myThid)
539         CALL DEBUG_STATS_RL(Nr,GtNm1,'GtNm1 (DYNAMICS)',myThid)         CALL DEBUG_STATS_RL(Nr,gvNm1,'GvNm1 (DYNAMICS)',myThid)
540         CALL DEBUG_STATS_RL(Nr,GsNm1,'GsNm1 (DYNAMICS)',myThid)         CALL DEBUG_STATS_RL(Nr,gtNm1,'GtNm1 (DYNAMICS)',myThid)
541           CALL DEBUG_STATS_RL(Nr,gsNm1,'GsNm1 (DYNAMICS)',myThid)
542    #endif
543        ENDIF        ENDIF
544  #endif  #endif
545    
Legend:
Removed from v.1.112  
changed lines
  Added in v.1.117
  ViewVC Help
Powered by ViewVC 1.1.22