pkg/mnc/mnc_init.F

C $Header: /u/gcmpack/MITgcm/pkg/mnc/mnc_init.F,v 1.20 2005/09/10 18:30:07 edhill Exp $
C $Name:  $
      
#include "MNC_OPTIONS.h"
      
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
CBOP 0
C     !ROUTINE: MNC_INIT

C     !INTERFACE:
      SUBROUTINE MNC_INIT( myThid )

C     !DESCRIPTION:
C     Initialize (zero) the look-up tables.  This routine should
C     \textbf{always} be run before any other MNC subroutines to ensure
C     that the lookup tables start in a well-defined state.
      
C     !USES:
#include "mnc_common.h"
#include "SIZE.h"
#include "EEPARAMS.h"
#include "EESUPPORT.h"
#include "PARAMS.h"

C     !INPUT PARAMETERS:
      integer myThid
CEOP

C     !LOCAL VARIABLES:
      integer i,j, g, v
      character blank*(MNC_MAX_CHAR)

C     Write blanks or zeros to all the internal names and ID tables
      DO i = 1,MNC_MAX_CHAR
        blank(i:i) = ' '
      ENDDO
      mnc_blank_name(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
      DO i = 1,MNC_MAX_ID
        mnc_d_size(i) = 0
        mnc_d_ids(i)  = 0
        mnc_f_names(i)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
        mnc_g_names(i)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
        mnc_v_names(i)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
        mnc_d_names(i)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
        DO j = 1,MNC_MAX_INFO
          mnc_f_info(i,j) = 0
          mnc_fv_ids(i,j) = 0
          mnc_fd_ind(i,j) = 0
          mnc_f_alld(i,j) = 0
        ENDDO
      ENDDO

C     Blank the CW tables
      DO g = 1,MNC_MAX_ID
        mnc_cw_gname(g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
        mnc_cw_vname(g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
        mnc_cw_vgind(g) = 0
        mnc_cw_vfmv(g) = 0
        mnc_cw_ndim(g) = 0
        mnc_cw_fgnm(g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
        mnc_cw_fgud(g) = 0
        mnc_cw_fgis(g) = 0
        mnc_cw_fgci(g) = 1
C       mnc_cw_cvnm(g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
        DO i = 1,MNC_CW_MAX_I
          mnc_cw_dn(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
          mnc_cw_vtnm(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
          mnc_cw_vinm(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
          mnc_cw_vdnm(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
          mnc_cw_vtat(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
          mnc_cw_dims(i,g) = 0
          mnc_cw_is(i,g) = 0
          mnc_cw_ie(i,g) = 0
          mnc_cw_viat(i,g) = 0
          mnc_cw_vdat(i,g) = 0.0D0
        ENDDO
        DO i = 1,3
          mnc_cw_vnat(i,g) = 0
        ENDDO
        DO i = 1,2
           mnc_cw_vbij(i,g) = 0
           mnc_cw_vmvi(i,g) = 0
           mnc_cw_vmvr(i,g) = 0.0
           mnc_cw_vmvd(i,g) = 0.0D0
C          mnc_cw_cvse(i,g) = 0
        ENDDO
      ENDDO

      DO i = 1,MNC_MAX_INFO
        mnc_cw_cit(1,i) = 0
        mnc_cw_cit(2,i) = 0
        mnc_cw_cit(3,i) = -1
      ENDDO

C     The default for all file types is to add the iter to the name and
C     grow them in "lock step" together
      mnc_cw_cit(1,1) = 1
      mnc_cw_cit(2,1) = nIter0

C     Here, we do not add the iter to the file name
      mnc_cw_cit(1,2) = -1
      mnc_cw_cit(2,2) = -1

C     For checkpoint files, we want to use the current iter but we do
C     not (by default, anyway) want to update the current iter for
C     everything else
      mnc_cw_cit(1,3) = 3
      mnc_cw_cit(2,3) = nIter0

C     DO i = 1,MNC_CW_CVDAT
C     mnc_cw_cvdt(i) = 0.0D0
C     ENDDO

      RETURN
      END

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
1	edhill	1.21	C $Header: /u/gcmpack/MITgcm/pkg/mnc/mnc_init.F,v 1.20 2005/09/10 18:30:07 edhill Exp $
2	edhill	1.1	C $Name: $
3
4			#include "MNC_OPTIONS.h"
5
6	edhill	1.6	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
7	edhill	1.14	CBOP 0
8	edhill	1.13	C !ROUTINE: MNC_INIT
9	edhill	1.1
10	edhill	1.13	C !INTERFACE:
11	edhill	1.11	SUBROUTINE MNC_INIT( myThid )
12	edhill	1.1
13	edhill	1.13	C !DESCRIPTION:
14	edhill	1.14	C Initialize (zero) the look-up tables. This routine should
15			C \textbf{always} be run before any other MNC subroutines to ensure
16			C that the lookup tables start in a well-defined state.
17	edhill	1.13
18			C !USES:
19	edhill	1.1	#include "mnc_common.h"
20	edhill	1.20	#include "SIZE.h"
21			#include "EEPARAMS.h"
22			#include "EESUPPORT.h"
23			#include "PARAMS.h"
24	edhill	1.1
25	edhill	1.13	C !INPUT PARAMETERS:
26	edhill	1.11	integer myThid
27	edhill	1.14	CEOP
28	edhill	1.1
29	edhill	1.13	C !LOCAL VARIABLES:
30	edhill	1.8	integer i,j, g, v
31	edhill	1.1	character blank*(MNC_MAX_CHAR)
32
33	edhill	1.8	C Write blanks or zeros to all the internal names and ID tables
34	edhill	1.10	DO i = 1,MNC_MAX_CHAR
35	edhill	1.17	blank(i:i) = ' '
36	edhill	1.1	ENDDO
37	edhill	1.2	mnc_blank_name(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
38	edhill	1.5	DO i = 1,MNC_MAX_ID
39			mnc_d_size(i) = 0
40			mnc_d_ids(i) = 0
41	edhill	1.1	mnc_f_names(i)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
42	edhill	1.2	mnc_g_names(i)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
43			mnc_v_names(i)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
44			mnc_d_names(i)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
45	edhill	1.3	DO j = 1,MNC_MAX_INFO
46			mnc_f_info(i,j) = 0
47			mnc_fv_ids(i,j) = 0
48			mnc_fd_ind(i,j) = 0
49	edhill	1.4	mnc_f_alld(i,j) = 0
50	edhill	1.3	ENDDO
51			ENDDO
52
53	edhill	1.15	C Blank the CW tables
54	edhill	1.6	DO g = 1,MNC_MAX_ID
55	edhill	1.7	mnc_cw_gname(g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
56	edhill	1.8	mnc_cw_vname(g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
57			mnc_cw_vgind(g) = 0
58	edhill	1.21	mnc_cw_vfmv(g) = 0
59	edhill	1.6	mnc_cw_ndim(g) = 0
60	edhill	1.12	mnc_cw_fgnm(g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
61			mnc_cw_fgud(g) = 0
62	edhill	1.16	mnc_cw_fgis(g) = 0
63	edhill	1.20	mnc_cw_fgci(g) = 1
64	edhill	1.19	C mnc_cw_cvnm(g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
65	edhill	1.6	DO i = 1,MNC_CW_MAX_I
66			mnc_cw_dn(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
67	edhill	1.8	mnc_cw_vtnm(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
68			mnc_cw_vinm(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
69			mnc_cw_vdnm(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
70			mnc_cw_vtat(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
71	edhill	1.6	mnc_cw_dims(i,g) = 0
72			mnc_cw_is(i,g) = 0
73			mnc_cw_ie(i,g) = 0
74	edhill	1.8	mnc_cw_viat(i,g) = 0
75			mnc_cw_vdat(i,g) = 0.0D0
76			ENDDO
77			DO i = 1,3
78			mnc_cw_vnat(i,g) = 0
79	edhill	1.10	ENDDO
80			DO i = 1,2
81			mnc_cw_vbij(i,g) = 0
82	edhill	1.21	mnc_cw_vmvi(i,g) = 0
83			mnc_cw_vmvr(i,g) = 0.0
84			mnc_cw_vmvd(i,g) = 0.0D0
85	edhill	1.19	C mnc_cw_cvse(i,g) = 0
86	edhill	1.6	ENDDO
87			ENDDO
88	edhill	1.20
89			DO i = 1,MNC_MAX_INFO
90			mnc_cw_cit(1,i) = 0
91			mnc_cw_cit(2,i) = 0
92			mnc_cw_cit(3,i) = -1
93			ENDDO
94
95			C The default for all file types is to add the iter to the name and
96			C grow them in "lock step" together
97			mnc_cw_cit(1,1) = 1
98			mnc_cw_cit(2,1) = nIter0
99
100			C Here, we do not add the iter to the file name
101			mnc_cw_cit(1,2) = -1
102			mnc_cw_cit(2,2) = -1
103
104			C For checkpoint files, we want to use the current iter but we do
105			C not (by default, anyway) want to update the current iter for
106			C everything else
107			mnc_cw_cit(1,3) = 3
108			mnc_cw_cit(2,3) = nIter0
109
110	edhill	1.19	C DO i = 1,MNC_CW_CVDAT
111			C mnc_cw_cvdt(i) = 0.0D0
112			C ENDDO
113	edhill	1.6
114	edhill	1.3	RETURN
115			END
116
117	edhill	1.6	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|