pkg/regrid/regrid_mnc_init.F

C $Header: /u/gcmpack/MITgcm_contrib/eh3/regrid/regrid/regrid_mnc_init.F,v 1.1 2006/08/10 05:00:15 edhill Exp $
C $Name:  $

#include "REGRID_OPTIONS.h"

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
CBOP 0
C     !ROUTINE: REGRID_MNC_INIT

C     !INTERFACE:
      SUBROUTINE REGRID_MNC_INIT( myThid )

C     !DESCRIPTION:
C     Initialize MNC lookup tables for REGRID.

C     !USES:
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "REGRID_SIZE.h"
#include "REGRID.h"
#ifdef ALLOW_MNC
#include "mnc_common.h"
#endif /* ALLOW_MNC */

C     !INPUT PARAMETERS:
      INTEGER myThid
CEOP

C     !LOCAL VARIABLES:
      INTEGER i, k, ndim
#ifdef ALLOW_MNC
      INTEGER CW_MAX_LOC
      PARAMETER ( CW_MAX_LOC = 6 )
      INTEGER dim(CW_MAX_LOC), ib(CW_MAX_LOC), ie(CW_MAX_LOC)
      CHARACTER*(MNC_MAX_CHAR) name
      character*(MNC_MAX_CHAR) dn(CW_MAX_LOC)
      
      IF (useMNC) THEN

C         DO i = 1,CW_MAX_LOC
C           dn(i)(1:MNC_MAX_CHAR)=mnc_blank_name(1:MNC_MAX_CHAR)
C           dim(i) = 0
C           ib(i) = 0
C           ie(i) = 0
C         ENDDO
C         dn(2)(1:1) = 'T'
C         dim(2)     = -1
C         ib(2)      =  1
C         ie(2)      =  1

C         DO k = 1,regrid_ngrids

C           name(1:MNC_MAX_CHAR) = mnc_blank_name(1:MNC_MAX_CHAR)
C           dn(1)(1:MNC_MAX_CHAR) = mnc_blank_name(1:MNC_MAX_CHAR)
C           write(name,'(a,i3)') 'regrid_', k
C           write(dn(1),'(a,i3)') 'regrid_', k
C           dim(1) = REGRID_nout(k)
C           ib(1)  = 1
C           ie(1)  = REGRID_nout(k)
          
C           CALL MNC_CW_ADD_GNAME(name, ndim, dim, dn, ib, ie, myThid)

C         ENDDO
        
      ENDIF

#endif /* ALLOW_MNC */

      RETURN
      END

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
1	C $Header: /u/gcmpack/MITgcm_contrib/eh3/regrid/regrid/regrid_mnc_init.F,v 1.1 2006/08/10 05:00:15 edhill Exp $
2	C $Name: $
3
4	#include "REGRID_OPTIONS.h"
5
6	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
7	CBOP 0
8	C !ROUTINE: REGRID_MNC_INIT
9
10	C !INTERFACE:
11	SUBROUTINE REGRID_MNC_INIT( myThid )
12
13	C !DESCRIPTION:
14	C Initialize MNC lookup tables for REGRID.
15
16	C !USES:
17	IMPLICIT NONE
18	#include "SIZE.h"
19	#include "EEPARAMS.h"
20	#include "PARAMS.h"
21	#include "REGRID_SIZE.h"
22	#include "REGRID.h"
23	#ifdef ALLOW_MNC
24	#include "mnc_common.h"
25	#endif /* ALLOW_MNC */
26
27	C !INPUT PARAMETERS:
28	INTEGER myThid
29	CEOP
30
31	C !LOCAL VARIABLES:
32	INTEGER i, k, ndim
33	#ifdef ALLOW_MNC
34	INTEGER CW_MAX_LOC
35	PARAMETER ( CW_MAX_LOC = 6 )
36	INTEGER dim(CW_MAX_LOC), ib(CW_MAX_LOC), ie(CW_MAX_LOC)
37	CHARACTER*(MNC_MAX_CHAR) name
38	character*(MNC_MAX_CHAR) dn(CW_MAX_LOC)
39
40	IF (useMNC) THEN
41
42	C DO i = 1,CW_MAX_LOC
43	C dn(i)(1:MNC_MAX_CHAR)=mnc_blank_name(1:MNC_MAX_CHAR)
44	C dim(i) = 0
45	C ib(i) = 0
46	C ie(i) = 0
47	C ENDDO
48	C dn(2)(1:1) = 'T'
49	C dim(2) = -1
50	C ib(2) = 1
51	C ie(2) = 1
52
53	C DO k = 1,regrid_ngrids
54
55	C name(1:MNC_MAX_CHAR) = mnc_blank_name(1:MNC_MAX_CHAR)
56	C dn(1)(1:MNC_MAX_CHAR) = mnc_blank_name(1:MNC_MAX_CHAR)
57	C write(name,'(a,i3)') 'regrid_', k
58	C write(dn(1),'(a,i3)') 'regrid_', k
59	C dim(1) = REGRID_nout(k)
60	C ib(1) = 1
61	C ie(1) = REGRID_nout(k)
62
63	C CALL MNC_CW_ADD_GNAME(name, ndim, dim, dn, ib, ie, myThid)
64
65	C ENDDO
66
67	ENDIF
68
69	#endif /* ALLOW_MNC */
70
71	RETURN
72	END
73
74	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|