/[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.63 by jmc, Tue Feb 20 15:06:21 2001 UTC revision 1.65 by jmc, Thu Mar 8 20:25:01 2001 UTC
# Line 26  C     \================================= Line 26  C     \=================================
26  C     == Global variables ===  C     == Global variables ===
27  #include "SIZE.h"  #include "SIZE.h"
28  #include "EEPARAMS.h"  #include "EEPARAMS.h"
 #include "CG2D.h"  
29  #include "PARAMS.h"  #include "PARAMS.h"
30  #include "DYNVARS.h"  #include "DYNVARS.h"
31  #include "GRID.h"  #include "GRID.h"
# Line 40  C     == Global variables === Line 39  C     == Global variables ===
39  # include "KPP.h"  # include "KPP.h"
40  #endif  #endif
41    
42  #ifdef INCLUDE_DIAGNOSTICS_INTERFACE_CODE  #ifdef ALLOW_TIMEAVE
43  #include "AVER.h"  #include "TIMEAVE_STATV.h"
44  #endif  #endif
45    
46  C     == Routine arguments ==  C     == Routine arguments ==
# Line 68  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 533  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       &           *(cg2d_x(i,j,bi,bj)-cg2d_x(i-1,j,bi,bj))       O         phiSurfX,phiSurfY,
541                phiSurfY(i,j) = _recip_dyC(i,j,bi,bj)*gBaro       I         myThid )                        
      &           *(cg2d_x(i,j,bi,bj)-cg2d_x(i,j-1,bi,bj))  
             ENDDO  
           ENDDO  
542          ENDIF          ENDIF
543    
544  C--     Start of dynamics loop  C--     Start of dynamics loop
# Line 670  c         CALL WRITE_FLD_XYZ_RL('PH.',su Line 665  c         CALL WRITE_FLD_XYZ_RL('PH.',su
665  c       ENDIF  c       ENDIF
666  Cjmc(end)  Cjmc(end)
667    
668  #ifdef INCLUDE_DIAGNOSTICS_INTERFACE_CODE  #ifdef ALLOW_TIMEAVE
669          IF (taveFreq.GT.0.) THEN          IF (taveFreq.GT.0.) THEN
670           DO K=1,Nr            CALL TIMEAVE_CUMULATE(phiHydtave, phiHyd, Nr,
671            CALL TIMEAVER_1FLD_XYZ(phiHyd, phiHydtave,       I                              deltaTclock, bi, bj, myThid)
      I                              deltaTclock, bi, bj, K, myThid)  
672            IF (ivdc_kappa.NE.0.) THEN            IF (ivdc_kappa.NE.0.) THEN
673              CALL TIMEAVER_1FLD_XYZ(ConvectCount, ConvectCountTave,              CALL TIMEAVE_CUMULATE(ConvectCountTave, ConvectCount, Nr,
674       I                              deltaTclock, bi, bj, K, myThid)       I                              deltaTclock, bi, bj, myThid)
675            ENDIF            ENDIF
          ENDDO  
676          ENDIF          ENDIF
677  #endif /* INCLUDE_DIAGNOSTICS_INTERFACE_CODE */  #endif /* ALLOW_TIMEAVE */
678    
679         ENDDO         ENDDO
680        ENDDO        ENDDO
Legend:
Removed from v.1.63  
changed lines
  Added in v.1.65
  ViewVC Help
Powered by ViewVC 1.1.22