pkg/my82/my82_readparms.F

C $Header: /u/gcmpack/MITgcm/pkg/my82/my82_readparms.F,v 1.1 2004/09/02 09:11:54 mlosch Exp $
C $Name:  $

#include "MY82_OPTIONS.h"

CBOP
C !ROUTINE: MY82_READPARMS

C !INTERFACE: ==========================================================
      SUBROUTINE MY82_READPARMS( myThid )

C !DESCRIPTION:
C     Initialize MY82 parameters, read in data.pp

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

C !INPUT PARAMETERS: ===================================================
C  myThid       :: thread number
      INTEGER myThid

C !OUTPUT PARAMETERS: ==================================================
C  none

#ifdef ALLOW_MY82

C !LOCAL VARIABLES: ====================================================
C  iUnit        :: unit number for I/O
C  msgBuf       :: message buffer
      INTEGER iUnit
      CHARACTER*(MAX_LEN_MBUF) msgBuf
CEOP

      NAMELIST /MY_PARM01/
     &     MYviscMax,
     &     MYdiffMax,
     &     MYhblScale,
     &     MYdumpFreq,
     &     MYMixingMaps,
     &     MYwriteState,
     &     MYtaveFreq

C This routine has been called by the main model so we set our
C internal flag to indicate we are in business
      MYisON=.TRUE.

C Set defaults values for parameters in MY82.h
      MYviscMax    = 1. _d 0
      MYdiffMax    = 1. _d 0
      MYhblScale   = 0.1 _d 0
      RiMax        = 0.1950 _d 0
      MYdumpFreq   = dumpFreq
      MYtaveFreq   = taveFreq
      MYMixingMaps = .FALSE.
      MYwriteState = .FALSE.

C Open and read the data.pp file
      _BEGIN_MASTER(myThid)
      WRITE(msgBuf,'(A)') ' MY82_READPARMS: opening data.my82'
      CALL PRINT_MESSAGE(msgBuf, standardMessageUnit,
     &                   SQUEEZE_RIGHT , 1)
      CALL OPEN_COPY_DATA_FILE(
     I                   'data.my82', 'MY82_READPARMS',
     O                   iUnit,
     I                   myThid )
      READ(UNIT=iUnit,NML=MY_PARM01)
      WRITE(msgBuf,'(A)')
     &  ' MY82_READPARMS: finished reading data.my82'
      CALL PRINT_MESSAGE(msgBuf, standardMessageUnit,
     &                   SQUEEZE_RIGHT , 1)

C Close the open data file
      CLOSE(iUnit)
      _END_MASTER(myThid)

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

C Now set-up any remaining parameters that result from the input parameters
      IF ( MYviscMax .LE. 0. ) THEN
       WRITE(msgBuf,'(A)') 'MYviscMax must be greater than zero'
       CALL PRINT_ERROR( msgBuf , 1)
       STOP 'ABNORMAL END: S/R MY82_READPARMS'
      ENDIF
      IF ( MYdiffMax .LE. 0. ) THEN
       WRITE(msgBuf,'(A)') 'MYdiffMax must be greater than zero'
       CALL PRINT_ERROR( msgBuf , 1)
       STOP 'ABNORMAL END: S/R MY82_READPARMS'
      ENDIF
#endif /* ALLOW_MY82 */

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/my82/my82_readparms.F,v 1.1 2004/09/02 09:11:54 mlosch Exp $
2	C $Name: $
3
4	#include "MY82_OPTIONS.h"
5
6	CBOP
7	C !ROUTINE: MY82_READPARMS
8
9	C !INTERFACE: ==========================================================
10	SUBROUTINE MY82_READPARMS( myThid )
11
12	C !DESCRIPTION:
13	C Initialize MY82 parameters, read in data.pp
14
15	C !USES: ===============================================================
16	IMPLICIT NONE
17	#include "SIZE.h"
18	#include "EEPARAMS.h"
19	#include "MY82.h"
20	#include "PARAMS.h"
21
22	C !INPUT PARAMETERS: ===================================================
23	C myThid :: thread number
24	INTEGER myThid
25
26	C !OUTPUT PARAMETERS: ==================================================
27	C none
28
29	#ifdef ALLOW_MY82
30
31	C !LOCAL VARIABLES: ====================================================
32	C iUnit :: unit number for I/O
33	C msgBuf :: message buffer
34	INTEGER iUnit
35	CHARACTER*(MAX_LEN_MBUF) msgBuf
36	CEOP
37
38	NAMELIST /MY_PARM01/
39	& MYviscMax,
40	& MYdiffMax,
41	& MYhblScale,
42	& MYdumpFreq,
43	& MYMixingMaps,
44	& MYwriteState,
45	& MYtaveFreq
46
47	C This routine has been called by the main model so we set our
48	C internal flag to indicate we are in business
49	MYisON=.TRUE.
50
51	C Set defaults values for parameters in MY82.h
52	MYviscMax = 1. _d 0
53	MYdiffMax = 1. _d 0
54	MYhblScale = 0.1 _d 0
55	RiMax = 0.1950 _d 0
56	MYdumpFreq = dumpFreq
57	MYtaveFreq = taveFreq
58	MYMixingMaps = .FALSE.
59	MYwriteState = .FALSE.
60
61	C Open and read the data.pp file
62	_BEGIN_MASTER(myThid)
63	WRITE(msgBuf,'(A)') ' MY82_READPARMS: opening data.my82'
64	CALL PRINT_MESSAGE(msgBuf, standardMessageUnit,
65	& SQUEEZE_RIGHT , 1)
66	CALL OPEN_COPY_DATA_FILE(
67	I 'data.my82', 'MY82_READPARMS',
68	O iUnit,
69	I myThid )
70	READ(UNIT=iUnit,NML=MY_PARM01)
71	WRITE(msgBuf,'(A)')
72	& ' MY82_READPARMS: finished reading data.my82'
73	CALL PRINT_MESSAGE(msgBuf, standardMessageUnit,
74	& SQUEEZE_RIGHT , 1)
75
76	C Close the open data file
77	CLOSE(iUnit)
78	_END_MASTER(myThid)
79
80	C Everyone else must wait for the parameters to be loaded
81	_BARRIER
82
83	C Now set-up any remaining parameters that result from the input parameters
84	IF ( MYviscMax .LE. 0. ) THEN
85	WRITE(msgBuf,'(A)') 'MYviscMax must be greater than zero'
86	CALL PRINT_ERROR( msgBuf , 1)
87	STOP 'ABNORMAL END: S/R MY82_READPARMS'
88	ENDIF
89	IF ( MYdiffMax .LE. 0. ) THEN
90	WRITE(msgBuf,'(A)') 'MYdiffMax must be greater than zero'
91	CALL PRINT_ERROR( msgBuf , 1)
92	STOP 'ABNORMAL END: S/R MY82_READPARMS'
93	ENDIF
94	#endif /* ALLOW_MY82 */
95
96	RETURN
97	END