pkg/diagnostics/diagnostics_init_varia.F

C $Header: /u/gcmpack/MITgcm/pkg/diagnostics/diagnostics_init_varia.F,v 1.1 2004/12/13 21:43:54 jmc Exp $
C $Name:  $

#include "DIAG_OPTIONS.h"
      
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
CBOP 0
C     !ROUTINE: DIAGNOSTICS_INIT_VARIA

C     !INTERFACE:
      SUBROUTINE DIAGNOSTICS_INIT_VARIA(
     I     myThid )

C     !DESCRIPTION:
C     Initialize the qdiag array which accumulates during integration
      
C     !USES:
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "DIAGNOSTICS_SIZE.h"
#include "DIAGNOSTICS.h"

C     !INPUT PARAMETERS:
      INTEGER myThid
CEOP

C     !LOCAL VARIABLES:
      INTEGER i,j,n,bi,bj
      CHARACTER*(80) fn
      LOGICAL glf

C     Zero out the qdiag array which accumulates during integration
      DO bj = myByLo(myThid), myByHi(myThid)
        DO bi = myBxLo(myThid), myBxHi(myThid)
          DO n = 1,numdiags
            DO j = 1-Oly,sNy+Oly
              DO i = 1-Olx,sNx+Olx
                qdiag(i,j,n,bi,bj) = 0. _d 0
              ENDDO
            ENDDO
          ENDDO
        ENDDO
      ENDDO


C     Add pickup capability
      IF (diag_pickup_read) THEN

#ifdef ALLOW_MNC
C       IF (useMNC .AND. diag_pickup_read_mnc) THEN
C       ENDIF
#endif
        
        IF (diag_pickup_read_mdsio) THEN
          DO i = 1,80
            fn(i:i) = ' '
          ENDDO
          write(fn,'(a)') 'pickup_diagnostics'
          glf = globalFiles
          CALL MDSREADFIELD(fn,readBinaryPrec,glf,'RL',
     &         numdiags,qdiag,1,myThid)
        ENDIF

      ENDIF

      RETURN 
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/diagnostics/diagnostics_init_varia.F,v 1.1 2004/12/13 21:43:54 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,n,bi,bj
31	CHARACTER*(80) fn
32	LOGICAL glf
33
34	C Zero out the qdiag array which accumulates during integration
35	DO bj = myByLo(myThid), myByHi(myThid)
36	DO bi = myBxLo(myThid), myBxHi(myThid)
37	DO n = 1,numdiags
38	DO j = 1-Oly,sNy+Oly
39	DO i = 1-Olx,sNx+Olx
40	qdiag(i,j,n,bi,bj) = 0. _d 0
41	ENDDO
42	ENDDO
43	ENDDO
44	ENDDO
45	ENDDO
46
47
48	C Add pickup capability
49	IF (diag_pickup_read) THEN
50
51	#ifdef ALLOW_MNC
52	C IF (useMNC .AND. diag_pickup_read_mnc) THEN
53	C ENDIF
54	#endif
55
56	IF (diag_pickup_read_mdsio) THEN
57	DO i = 1,80
58	fn(i:i) = ' '
59	ENDDO
60	write(fn,'(a)') 'pickup_diagnostics'
61	glf = globalFiles
62	CALL MDSREADFIELD(fn,readBinaryPrec,glf,'RL',
63	& numdiags,qdiag,1,myThid)
64	ENDIF
65
66	ENDIF
67
68	RETURN
69	END