/[MITgcm]/MITgcm/diags/src/timeaver_final_xyz.F
ViewVC logotype

Contents of /MITgcm/diags/src/timeaver_final_xyz.F

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

Revision 1.1 - (show annotations) (download)
Wed Jul 1 19:49:36 1998 UTC (25 years, 10 months ago) by adcroft
Branch: MAIN
CVS Tags: checkpoint11, checkpoint13, checkpoint12, checkpoint10, branch-point-rdot
Branch point for: branch-rdot 
The time-averaging diagnostics has been modified to include
averages of "intermediate" quantities that exist only on a tile.
We currently only time-average K13, K23 and wVel but these
should be usable as a template for other diagnostics.
We have now split off these routines into a separate module
called diags/. Use of these routines is enabled by the
CPP flag ALLOW_DIAGNOSTICS set in model/inc/CPP_OPTIONS.h
1 #include "CPP_EEOPTIONS.h"
2
3 CStartofinterface
4 SUBROUTINE TIMEAVER_FINAL_XYZ( fldtave, timeave, bi, bj, K, myThid )
5 C /==========================================================\
6 C | SUBROUTINE TIMEAVER_FINAL_XY |
7 C | o Get final average of field |
8 C \==========================================================/
9
10 C == Global variables ===
11 #include "SIZE.h"
12 #include "EEPARAMS.h"
13
14 C == Routine arguments ==
15 C myThid - Thread number for this instance of the routine.
16 C timeave - time averaged
17 C fldtave - time averaged Field
18 INTEGER bi, bj, K, myThid
19 _RL fldtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy)
20 REAL timeave(Nz,nSx,nSy)
21
22 CEndofinterface
23
24 C == Local variables ==
25 C i,j,k,bi,bj - Loop counters
26 INTEGER i, j
27
28 C DO bj = myByLo(myThid), myByHi(myThid)
29 C DO bi = myBxLo(myThid), myBxHi(myThid)
30 C DO k=1,Nz
31 DO j=1,sNy
32 DO i=1,sNx
33 fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj) / timeave(k,bi,bj)
34 ENDDO
35 ENDDO
36 C ENDDO
37 C ENDDO
38 C ENDDO
39
40 RETURN
41 END
  ViewVC Help
Powered by ViewVC 1.1.22