/[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.64 by jmc, Tue Mar 6 16:59:44 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 539  C       (note: this loop will be replace Line 538  C       (note: this loop will be replace
538            DO j=jMin,jMax            DO j=jMin,jMax
539              DO i=iMin,iMax              DO i=iMin,iMax
540                phiSurfX(i,j) = _recip_dxC(i,j,bi,bj)*gBaro                phiSurfX(i,j) = _recip_dxC(i,j,bi,bj)*gBaro
541       &           *(cg2d_x(i,j,bi,bj)-cg2d_x(i-1,j,bi,bj))       &           *(etaN(i,j,bi,bj)-etaN(i-1,j,bi,bj))
542                phiSurfY(i,j) = _recip_dyC(i,j,bi,bj)*gBaro                phiSurfY(i,j) = _recip_dyC(i,j,bi,bj)*gBaro
543       &           *(cg2d_x(i,j,bi,bj)-cg2d_x(i,j-1,bi,bj))       &           *(etaN(i,j,bi,bj)-etaN(i,j-1,bi,bj))
544              ENDDO              ENDDO
545            ENDDO            ENDDO
546          ENDIF          ENDIF
# Line 670  c         CALL WRITE_FLD_XYZ_RL('PH.',su Line 669  c         CALL WRITE_FLD_XYZ_RL('PH.',su
669  c       ENDIF  c       ENDIF
670  Cjmc(end)  Cjmc(end)
671    
672  #ifdef INCLUDE_DIAGNOSTICS_INTERFACE_CODE  #ifdef ALLOW_TIMEAVE
673          IF (taveFreq.GT.0.) THEN          IF (taveFreq.GT.0.) THEN
674           DO K=1,Nr            CALL TIMEAVE_CUMULATE(phiHydtave, phiHyd, Nr,
675            CALL TIMEAVER_1FLD_XYZ(phiHyd, phiHydtave,       I                              deltaTclock, bi, bj, myThid)
      I                              deltaTclock, bi, bj, K, myThid)  
676            IF (ivdc_kappa.NE.0.) THEN            IF (ivdc_kappa.NE.0.) THEN
677              CALL TIMEAVER_1FLD_XYZ(ConvectCount, ConvectCountTave,              CALL TIMEAVE_CUMULATE(ConvectCountTave, ConvectCount, Nr,
678       I                              deltaTclock, bi, bj, K, myThid)       I                              deltaTclock, bi, bj, myThid)
679            ENDIF            ENDIF
          ENDDO  
680          ENDIF          ENDIF
681  #endif /* INCLUDE_DIAGNOSTICS_INTERFACE_CODE */  #endif /* ALLOW_TIMEAVE */
682    
683         ENDDO         ENDDO
684        ENDDO        ENDDO
Legend:
Removed from v.1.63  
changed lines
  Added in v.1.64
  ViewVC Help
Powered by ViewVC 1.1.22