/[MITgcm]/MITgcm/pkg/timeave/timeave_statvars.F
ViewVC logotype

Diff of /MITgcm/pkg/timeave/timeave_statvars.F

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

revision 1.13 by dimitri, Thu Nov 13 06:35:15 2003 UTC revision 1.14 by jmc, Fri Dec 5 02:33:56 2003 UTC
# Line 17  C     == Global variables === Line 17  C     == Global variables ===
17  #include "PARAMS.h"  #include "PARAMS.h"
18  #include "DYNVARS.h"  #include "DYNVARS.h"
19  #include "GRID.h"  #include "GRID.h"
 #include "FFIELDS.h"  
20  #include "TIMEAVE_STATV.h"  #include "TIMEAVE_STATV.h"
21    
22        LOGICAL  DIFFERENT_MULTIPLE        LOGICAL  DIFFERENT_MULTIPLE
# Line 157  C-    Time Averages of "intermediate" fi Line 156  C-    Time Averages of "intermediate" fi
156        IF ( myIter .NE. nIter0 ) THEN        IF ( myIter .NE. nIter0 ) THEN
157    
158  C-    Time Averages of surface fluxes  C-    Time Averages of surface fluxes
159         IF ( buoyancyRelation .EQ. 'OCEANICP' ) THEN  C     <- moved to external_forcing_surf
         k=Nr  
        ELSE  
         k=1  
        ENDIF  
   
 C     uFlux  
        DO j=1,sNy  
         DO i=1,sNx  
          tempArray(i,j,bi,bj)=fu(i,j,bi,bj)*foFacMom*_maskW(i,j,k,bi,bj)  
         ENDDO  
        ENDDO  
        CALL TIMEAVE_CUMULATE(uFluxtave,tempArray,1,  
      &                                   deltaTclock, bi, bj, myThid)  
   
 C     vFlux  
        DO j=1,sNy  
         DO i=1,sNx  
          tempArray(i,j,bi,bj)=fv(i,j,bi,bj)*foFacMom*_maskS(i,j,k,bi,bj)  
         ENDDO  
        ENDDO  
        CALL TIMEAVE_CUMULATE(vFluxtave,tempArray,1,  
      &                                   deltaTclock, bi, bj, myThid)  
   
 C     tFlux  
        DO j=1,sNy  
         DO i=1,sNx  
          tempArray(i,j,bi,bj)=maskC(i,j,k,bi,bj)*(  
 #ifdef SHORTWAVE_HEATING  
      &    -Qsw(i,j,bi,bj)+  
 #endif  
      &    (surfaceTendencyT(i,j,bi,bj)+surfaceTendencyTice(I,J,bi,bj))*  
      &    HeatCapacity_Cp*  
      &    recip_horiVertRatio*rhoConst*drF(k)*hFacC(i,j,k,bi,bj))  
         ENDDO  
        ENDDO  
        CALL TIMEAVE_CUMULATE(tFluxtave,tempArray,1,  
      &                                   deltaTclock, bi, bj, myThid)  
   
 C     sFlux  
        DO j=1,sNy  
         DO i=1,sNx  
          tempArray(i,j,bi,bj)=maskC(i,j,k,bi,bj)*  
      &    surfaceTendencyS(i,j,bi,bj)*  
      &    recip_horiVertRatio*rhoConst*drF(k)*hFacC(i,j,k,bi,bj)  
         ENDDO  
        ENDDO  
        CALL TIMEAVE_CUMULATE(sFluxtave,tempArray,1,  
      &                                   deltaTclock, bi, bj, myThid)  
160    
161          CALL TIMEAVE_CUMULATE(phiHydtave, totPhihyd,  Nr,          CALL TIMEAVE_CUMULATE(phiHydtave, totPhihyd,  Nr,
162       &                                   deltaTclock, bi, bj, myThid)       &                                   deltaTclock, bi, bj, myThid)
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22