pkg/mnc/mnc_init.F

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