/[MITgcm]/MITgcm/pkg/diagnostics/diagnostics_init_varia.F
ViewVC logotype

Contents of /MITgcm/pkg/diagnostics/diagnostics_init_varia.F

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


Revision 1.10 - (show annotations) (download)
Mon Jun 5 18:15:53 2006 UTC (17 years, 11 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint58l_post, checkpoint58n_post, checkpoint58h_post, checkpoint58q_post, checkpoint58j_post, checkpoint58i_post, checkpoint58g_post, checkpoint58o_post, checkpoint58k_post, checkpoint58p_post, checkpoint58m_post
Changes since 1.9: +5 -1 lines
Implement periodic averaging diagnostics (e.g., mean seasonal cycle,
 mean diurnal cycle)

1 C $Header: /u/gcmpack/MITgcm/pkg/diagnostics/diagnostics_init_varia.F,v 1.9 2005/06/26 16:51:49 jmc Exp $
2 C $Name: $
3
4 #include "DIAG_OPTIONS.h"
5
6 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7 CBOP 0
8 C !ROUTINE: DIAGNOSTICS_INIT_VARIA
9
10 C !INTERFACE:
11 SUBROUTINE DIAGNOSTICS_INIT_VARIA(
12 I myThid )
13
14 C !DESCRIPTION:
15 C Initialize the qdiag array which accumulates during integration
16
17 C !USES:
18 IMPLICIT NONE
19 #include "SIZE.h"
20 #include "EEPARAMS.h"
21 #include "PARAMS.h"
22 #include "DIAGNOSTICS_SIZE.h"
23 #include "DIAGNOSTICS.h"
24
25 C !INPUT PARAMETERS:
26 INTEGER myThid
27 CEOP
28
29 C !LOCAL VARIABLES:
30 INTEGER i,j,k,bi,bj
31
32 C Zero out the qdiag array which accumulates during integration
33 DO bj = myByLo(myThid), myByHi(myThid)
34 DO bi = myBxLo(myThid), myBxHi(myThid)
35 DO k = 1,numdiags
36 DO j = 1-Oly,sNy+Oly
37 DO i = 1-Olx,sNx+Olx
38 qdiag(i,j,k,bi,bj) = 0. _d 0
39 ENDDO
40 ENDDO
41 C Zero out the counters for the qdiag array
42 ndiag(k,bi,bj) = 0
43 ENDDO
44 DO k = 1,numlists
45 C Zero out the index array for periodic averaging diagnostic
46 pdiag(k,bi,bj) = 0
47 ENDDO
48 ENDDO
49 ENDDO
50
51 C Zero out the qSdiag array (statistics) which accumulates during integration
52 DO bj = myByLo(myThid), myByHi(myThid)
53 DO bi = myBxLo(myThid), myBxHi(myThid)
54 DO k = 1,diagSt_size
55 DO j = 0,nRegions
56 DO i = 0,nStats
57 qSdiag(i,j,k,bi,bj) = 0. _d 0
58 ENDDO
59 ENDDO
60 ENDDO
61 ENDDO
62 ENDDO
63
64 CALL DIAGNOSTICS_READ_PICKUP( myThid )
65
66 CALL DIAGSTATS_INI_IO( myThid )
67
68 CALL DIAGNOSTICS_SUMMARY( startTime, nIter0, myThid )
69
70 RETURN
71 END

  ViewVC Help
Powered by ViewVC 1.1.22