pkg/mnc/mnc_init.F

C $Header: /u/gcmpack/MITgcm/pkg/mnc/mnc_init.F,v 1.23 2006/03/10 22:01:54 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 "SIZE.h"
#include "EEPARAMS.h"
#include "EESUPPORT.h"
#include "PARAMS.h"
#include "mnc_common.h"
#include "MNC_PARAMS.h"

C     !INPUT PARAMETERS:
      integer myThid
CEOP

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

C     Write blanks or zeros to all the internal names and ID tables
      DO i = 1,MNC_MAX_CHAR
        blank(i:i) = ' '
      ENDDO
      DO i = 1,MNC_MAX_PATH
        bpath(i:i) = ' '
      ENDDO

      mnc_blank_name(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
      DO i = 1,MNC_MAX_FID
        mnc_f_names(i)(1:MNC_MAX_PATH) = bpath(1:MNC_MAX_PATH)
        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
      DO i = 1,MNC_MAX_ID
        mnc_d_size(i) = 0
        mnc_d_ids(i)  = 0
        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)
      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

      DO i = 1,2
       mnc_def_imv(i) = UNSET_I
       mnc_def_rmv(i) = UNSET_FLOAT4
       mnc_def_dmv(i) = UNSET_FLOAT8
      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	mlosch	1.24	C $Header: /u/gcmpack/MITgcm/pkg/mnc/mnc_init.F,v 1.23 2006/03/10 22:01:54 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.20	#include "SIZE.h"
20			#include "EEPARAMS.h"
21			#include "EESUPPORT.h"
22			#include "PARAMS.h"
23	mlosch	1.24	#include "mnc_common.h"
24			#include "MNC_PARAMS.h"
25	edhill	1.1
26	edhill	1.13	C !INPUT PARAMETERS:
27	edhill	1.11	integer myThid
28	edhill	1.14	CEOP
29	edhill	1.1
30	edhill	1.13	C !LOCAL VARIABLES:
31	edhill	1.22	integer i,j, g
32	edhill	1.1	character blank*(MNC_MAX_CHAR)
33	edhill	1.22	character bpath*(MNC_MAX_PATH)
34	edhill	1.1
35	edhill	1.8	C Write blanks or zeros to all the internal names and ID tables
36	edhill	1.10	DO i = 1,MNC_MAX_CHAR
37	edhill	1.17	blank(i:i) = ' '
38	edhill	1.1	ENDDO
39	edhill	1.22	DO i = 1,MNC_MAX_PATH
40			bpath(i:i) = ' '
41			ENDDO
42
43	edhill	1.2	mnc_blank_name(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
44	edhill	1.23	DO i = 1,MNC_MAX_FID
45	edhill	1.22	mnc_f_names(i)(1:MNC_MAX_PATH) = bpath(1:MNC_MAX_PATH)
46	edhill	1.3	DO j = 1,MNC_MAX_INFO
47			mnc_f_info(i,j) = 0
48			mnc_fv_ids(i,j) = 0
49			mnc_fd_ind(i,j) = 0
50	edhill	1.4	mnc_f_alld(i,j) = 0
51	edhill	1.3	ENDDO
52			ENDDO
53	edhill	1.23	DO i = 1,MNC_MAX_ID
54			mnc_d_size(i) = 0
55			mnc_d_ids(i) = 0
56			mnc_g_names(i)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
57			mnc_v_names(i)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
58			mnc_d_names(i)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
59			ENDDO
60	edhill	1.3
61	edhill	1.15	C Blank the CW tables
62	edhill	1.6	DO g = 1,MNC_MAX_ID
63	edhill	1.7	mnc_cw_gname(g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
64	edhill	1.8	mnc_cw_vname(g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
65			mnc_cw_vgind(g) = 0
66	edhill	1.21	mnc_cw_vfmv(g) = 0
67	edhill	1.6	mnc_cw_ndim(g) = 0
68	edhill	1.12	mnc_cw_fgnm(g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
69			mnc_cw_fgud(g) = 0
70	edhill	1.16	mnc_cw_fgis(g) = 0
71	edhill	1.20	mnc_cw_fgci(g) = 1
72	edhill	1.19	C mnc_cw_cvnm(g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
73	edhill	1.6	DO i = 1,MNC_CW_MAX_I
74			mnc_cw_dn(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
75	edhill	1.8	mnc_cw_vtnm(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
76			mnc_cw_vinm(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
77			mnc_cw_vdnm(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
78			mnc_cw_vtat(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)
79	edhill	1.6	mnc_cw_dims(i,g) = 0
80			mnc_cw_is(i,g) = 0
81			mnc_cw_ie(i,g) = 0
82	edhill	1.8	mnc_cw_viat(i,g) = 0
83			mnc_cw_vdat(i,g) = 0.0D0
84			ENDDO
85			DO i = 1,3
86			mnc_cw_vnat(i,g) = 0
87	edhill	1.10	ENDDO
88			DO i = 1,2
89			mnc_cw_vbij(i,g) = 0
90	edhill	1.21	mnc_cw_vmvi(i,g) = 0
91			mnc_cw_vmvr(i,g) = 0.0
92			mnc_cw_vmvd(i,g) = 0.0D0
93	edhill	1.19	C mnc_cw_cvse(i,g) = 0
94	edhill	1.6	ENDDO
95			ENDDO
96	edhill	1.20
97			DO i = 1,MNC_MAX_INFO
98			mnc_cw_cit(1,i) = 0
99			mnc_cw_cit(2,i) = 0
100			mnc_cw_cit(3,i) = -1
101			ENDDO
102
103	mlosch	1.24	DO i = 1,2
104			mnc_def_imv(i) = UNSET_I
105			mnc_def_rmv(i) = UNSET_FLOAT4
106			mnc_def_dmv(i) = UNSET_FLOAT8
107			ENDDO
108
109	edhill	1.20	C The default for all file types is to add the iter to the name and
110			C grow them in "lock step" together
111			mnc_cw_cit(1,1) = 1
112			mnc_cw_cit(2,1) = nIter0
113
114			C Here, we do not add the iter to the file name
115			mnc_cw_cit(1,2) = -1
116			mnc_cw_cit(2,2) = -1
117
118			C For checkpoint files, we want to use the current iter but we do
119			C not (by default, anyway) want to update the current iter for
120			C everything else
121			mnc_cw_cit(1,3) = 3
122			mnc_cw_cit(2,3) = nIter0
123
124	edhill	1.19	C DO i = 1,MNC_CW_CVDAT
125			C mnc_cw_cvdt(i) = 0.0D0
126			C ENDDO
127	edhill	1.6
128	edhill	1.3	RETURN
129			END
130
131	edhill	1.6	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|