pkg/mypackage/mypackage_readparms.F

C $Header: /u/gcmpack/MITgcm/pkg/mypackage/mypackage_readparms.F,v 1.3 2008/01/27 20:52:11 jmc Exp $
C $Name:  $

#include "MYPACKAGE_OPTIONS.h"

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

C !INTERFACE:
      SUBROUTINE MYPACKAGE_READPARMS( myThid )

C     !DESCRIPTION:
C     Initialize MYPACKAGE variables and constants.

C     !USES:
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "MYPACKAGE.h"

C     !INPUT PARAMETERS:
      INTEGER myThid
CEOP

#ifdef ALLOW_MYPACKAGE

C     !LOCAL VARIABLES:
C     msgBuf     :: Informational/error message buffer
C     iUnit      :: Work variable for IO unit number
      CHARACTER*(MAX_LEN_MBUF) msgBuf
      INTEGER iUnit

      NAMELIST /MYPACKAGE_PARM01/
     &       myPa_MNC,
     &       myPa_StaV_Cgrid, myPa_Tend_Cgrid,
     &       myPa_applyTendT, myPa_applyTendS,
     &       myPa_applyTendU, myPa_applyTendV,
     &       myPa_doSwitch1,  myPa_doSwitch2,
     &       myPa_index1,     myPa_index2,
     &       myPa_param1,     myPa_param2,
     &       myPa_string1,    myPa_string2,
     &       myPa_Scal1File,  myPa_Scal2File,
     &       myPa_VelUFile,   myPa_VelVFile,
     &       myPa_Surf1File,  myPa_Surf2File

      _BEGIN_MASTER(myThid)

C--   Default values for MYPACKAGE
      myPa_MNC         = useMNC
      myPa_StaV_Cgrid  = .TRUE.
      myPa_Tend_Cgrid  = .TRUE.
      myPa_applyTendT  = .FALSE.
      myPa_applyTendS  = .FALSE.
      myPa_applyTendU  = .FALSE.
      myPa_applyTendV  = .FALSE.
C-    additional parameters:
      myPa_doSwitch1   = .FALSE.
      myPa_doSwitch2   = .FALSE.
      myPa_index1      = 0
      myPa_index2      = 0
      myPa_param1      = 0. _d 0
      myPa_param2      = 0. _d 0
      myPa_string1     = ' '
      myPa_string2     = ' '
C-    file names for initial conditions:
      myPa_Scal1File   = ' '
      myPa_Scal2File   = ' '
      myPa_VelUFile    = ' '
      myPa_VelVFile    = ' '
      myPa_Surf1File   = ' '
      myPa_Surf2File   = ' '

      WRITE(msgBuf,'(A)') 'MYPACKAGE_READPARMS: opening data.mypackage'
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &                    SQUEEZE_RIGHT , 1)
      CALL OPEN_COPY_DATA_FILE(
     I                     'data.mypackage', 'MYPACKAGE_READPARMS',
     O                     iUnit,
     I                     myThid )

C     Read parameters from open data file
      READ(UNIT=iUnit,NML=MYPACKAGE_PARM01)
      WRITE(msgBuf,'(A)')
     &    'MYPACKAGE_READPARMS: finished reading data.mypackage'
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &    SQUEEZE_RIGHT , 1)
C     Close the open data file
      CLOSE(iUnit)

C     Make sure that we locally honor the global MNC on/off flag
      myPa_MNC = myPa_MNC .AND. useMNC
#ifndef ALLOW_MNC
C     Fix to avoid running without getting any output:
      myPa_MNC = .FALSE.
#endif
      myPa_MDSIO = (.NOT. myPa_MNC) .OR. outputTypesInclusive

      _END_MASTER(myThid)

C--   Everyone else must wait for the parameters to be loaded
      _BARRIER

#endif /* ALLOW_MYPACKAGE */

      RETURN
      END
1	jmc	1.4	C $Header: /u/gcmpack/MITgcm/pkg/mypackage/mypackage_readparms.F,v 1.3 2008/01/27 20:52:11 jmc Exp $
2	heimbach	1.1	C $Name: $
3
4			#include "MYPACKAGE_OPTIONS.h"
5
6	edhill	1.2	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
7			CBOP 0
8			C !ROUTINE: MYPACKAGE_READPARMS
9
10			C !INTERFACE:
11	heimbach	1.1	SUBROUTINE MYPACKAGE_READPARMS( myThid )
12	edhill	1.2
13			C !DESCRIPTION:
14			C Initialize MYPACKAGE variables and constants.
15
16			C !USES:
17	heimbach	1.1	IMPLICIT NONE
18			#include "SIZE.h"
19			#include "EEPARAMS.h"
20			#include "PARAMS.h"
21			#include "MYPACKAGE.h"
22
23	edhill	1.2	C !INPUT PARAMETERS:
24	heimbach	1.1	INTEGER myThid
25	edhill	1.2	CEOP
26	heimbach	1.1
27			#ifdef ALLOW_MYPACKAGE
28
29	jmc	1.4	C !LOCAL VARIABLES:
30			C msgBuf :: Informational/error message buffer
31			C iUnit :: Work variable for IO unit number
32			CHARACTER*(MAX_LEN_MBUF) msgBuf
33			INTEGER iUnit
34
35	heimbach	1.1	NAMELIST /MYPACKAGE_PARM01/
36	jmc	1.3	& myPa_MNC,
37			& myPa_StaV_Cgrid, myPa_Tend_Cgrid,
38			& myPa_applyTendT, myPa_applyTendS,
39			& myPa_applyTendU, myPa_applyTendV,
40			& myPa_doSwitch1, myPa_doSwitch2,
41			& myPa_index1, myPa_index2,
42			& myPa_param1, myPa_param2,
43			& myPa_string1, myPa_string2,
44			& myPa_Scal1File, myPa_Scal2File,
45			& myPa_VelUFile, myPa_VelVFile,
46			& myPa_Surf1File, myPa_Surf2File
47	heimbach	1.1
48			_BEGIN_MASTER(myThid)
49
50			C-- Default values for MYPACKAGE
51	jmc	1.3	myPa_MNC = useMNC
52			myPa_StaV_Cgrid = .TRUE.
53			myPa_Tend_Cgrid = .TRUE.
54			myPa_applyTendT = .FALSE.
55			myPa_applyTendS = .FALSE.
56			myPa_applyTendU = .FALSE.
57			myPa_applyTendV = .FALSE.
58			C- additional parameters:
59			myPa_doSwitch1 = .FALSE.
60			myPa_doSwitch2 = .FALSE.
61			myPa_index1 = 0
62			myPa_index2 = 0
63			myPa_param1 = 0. _d 0
64			myPa_param2 = 0. _d 0
65			myPa_string1 = ' '
66			myPa_string2 = ' '
67			C- file names for initial conditions:
68			myPa_Scal1File = ' '
69			myPa_Scal2File = ' '
70			myPa_VelUFile = ' '
71			myPa_VelVFile = ' '
72			myPa_Surf1File = ' '
73			myPa_Surf2File = ' '
74
75	heimbach	1.1	WRITE(msgBuf,'(A)') 'MYPACKAGE_READPARMS: opening data.mypackage'
76			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
77			& SQUEEZE_RIGHT , 1)
78			CALL OPEN_COPY_DATA_FILE(
79			I 'data.mypackage', 'MYPACKAGE_READPARMS',
80			O iUnit,
81			I myThid )
82
83			C Read parameters from open data file
84			READ(UNIT=iUnit,NML=MYPACKAGE_PARM01)
85	jmc	1.3	WRITE(msgBuf,'(A)')
86	heimbach	1.1	& 'MYPACKAGE_READPARMS: finished reading data.mypackage'
87			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
88			& SQUEEZE_RIGHT , 1)
89			C Close the open data file
90			CLOSE(iUnit)
91
92			C Make sure that we locally honor the global MNC on/off flag
93	jmc	1.3	myPa_MNC = myPa_MNC .AND. useMNC
94	heimbach	1.1	#ifndef ALLOW_MNC
95			C Fix to avoid running without getting any output:
96	jmc	1.3	myPa_MNC = .FALSE.
97	heimbach	1.1	#endif
98	jmc	1.3	myPa_MDSIO = (.NOT. myPa_MNC) .OR. outputTypesInclusive
99	heimbach	1.1
100			_END_MASTER(myThid)
101
102			C-- Everyone else must wait for the parameters to be loaded
103			_BARRIER
104
105			#endif /* ALLOW_MYPACKAGE */
106
107			RETURN
108			END