pkg/shap_filt/shap_filt_readparms.F

C $Header: /u/gcmpack/models/MITgcmUV/pkg/shap_filt/shap_filt_readparms.F,v 1.2 2001/05/29 14:01:40 adcroft Exp $
C $Name:  $

#include "SHAP_FILT_OPTIONS.h"

      SUBROUTINE SHAP_FILT_READPARMS( myThid )
C     /==========================================================\
C     | SUBROUTINE SHAP_FILT_READPARMS                           |
C     | o Routine to initialize Shapiro Filter parameters        |
C     |==========================================================|
C     \==========================================================/
      IMPLICIT NONE

C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "SHAP_FILT.h"

C     === Routine arguments ===
      INTEGER myThid

#ifdef ALLOW_SHAP_FILT

      NAMELIST /SHAP_PARM01/ Shap_funct, 
     &          nShapT,  nShapTrPhys,
     &          nShapUV, nShapUVPhys,
     &          Shap_Trtau, Shap_TrLength,
     &          Shap_uvtau, Shap_uvLength

C     === Local variables ===
C     msgBuf      - Informational/error meesage buffer
C     iUnit       - Work variable for IO unit number
      CHARACTER*(MAX_LEN_MBUF) msgBuf
      INTEGER iUnit

C--   SHAP_FILT_READPARMS has been called so we know that
C     the package is active.
c     SHAPIsOn=.TRUE.

      _BEGIN_MASTER(myThid)

      WRITE(msgBuf,'(A)') ' SHAP_FILT_READPARMS: opening data.shap'
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &                    SQUEEZE_RIGHT , 1)

      CALL OPEN_COPY_DATA_FILE(
     I                          'data.shap', 'SHAP_FILT_READPARMS',
     O                          iUnit,
     I                          myThid )

C--   Default flags and values for Shapiro Filter
      Shap_funct = 2
      nShapT = 0
      nShapUV = 0
      nShapTrPhys = 0
      nShapUVPhys = 0
      Shap_Trtau = deltaTtracer
      Shap_TrLength = 0.
      Shap_uvtau = deltaTMom
      Shap_TrLength = 0.

C--   Read parameters from open data file
      READ(UNIT=iUnit,NML=SHAP_PARM01)

      WRITE(msgBuf,'(A)') 
     &   ' SHAP_FILT_READPARMS: finished reading data.shap'
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &                    SQUEEZE_RIGHT , 1)

      write(0,*) 'Shap_funct, nShap_Tr,UV _Phys=',
     &    Shap_funct, nShapT, nShapUV, nShapTrPhys, nShapUVPhys
      write(0,*) 'Shap_Trtau,Shap_uvtau=',Shap_Trtau,Shap_uvtau
C--   Close the open data file
      CLOSE(iUnit)
      _END_MASTER(myThid)

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

#endif /* ALLOW_SHAP_FILT */
      RETURN
      END
1	C $Header: /u/gcmpack/models/MITgcmUV/pkg/shap_filt/shap_filt_readparms.F,v 1.2 2001/05/29 14:01:40 adcroft Exp $
2	C $Name: $
3
4	#include "SHAP_FILT_OPTIONS.h"
5
6	SUBROUTINE SHAP_FILT_READPARMS( myThid )
7	C /==========================================================\
8	C \| SUBROUTINE SHAP_FILT_READPARMS \|
9	C \| o Routine to initialize Shapiro Filter parameters \|
10	C \|==========================================================\|
11	C \==========================================================/
12	IMPLICIT NONE
13
14	C === Global variables ===
15	#include "SIZE.h"
16	#include "EEPARAMS.h"
17	#include "PARAMS.h"
18	#include "SHAP_FILT.h"
19
20	C === Routine arguments ===
21	INTEGER myThid
22
23	#ifdef ALLOW_SHAP_FILT
24
25	NAMELIST /SHAP_PARM01/ Shap_funct,
26	& nShapT, nShapTrPhys,
27	& nShapUV, nShapUVPhys,
28	& Shap_Trtau, Shap_TrLength,
29	& Shap_uvtau, Shap_uvLength
30
31	C === Local variables ===
32	C msgBuf - Informational/error meesage buffer
33	C iUnit - Work variable for IO unit number
34	CHARACTER*(MAX_LEN_MBUF) msgBuf
35	INTEGER iUnit
36
37	C-- SHAP_FILT_READPARMS has been called so we know that
38	C the package is active.
39	c SHAPIsOn=.TRUE.
40
41	_BEGIN_MASTER(myThid)
42
43	WRITE(msgBuf,'(A)') ' SHAP_FILT_READPARMS: opening data.shap'
44	CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
45	& SQUEEZE_RIGHT , 1)
46
47	CALL OPEN_COPY_DATA_FILE(
48	I 'data.shap', 'SHAP_FILT_READPARMS',
49	O iUnit,
50	I myThid )
51
52	C-- Default flags and values for Shapiro Filter
53	Shap_funct = 2
54	nShapT = 0
55	nShapUV = 0
56	nShapTrPhys = 0
57	nShapUVPhys = 0
58	Shap_Trtau = deltaTtracer
59	Shap_TrLength = 0.
60	Shap_uvtau = deltaTMom
61	Shap_TrLength = 0.
62
63	C-- Read parameters from open data file
64	READ(UNIT=iUnit,NML=SHAP_PARM01)
65
66	WRITE(msgBuf,'(A)')
67	& ' SHAP_FILT_READPARMS: finished reading data.shap'
68	CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
69	& SQUEEZE_RIGHT , 1)
70
71	write(0,*) 'Shap_funct, nShap_Tr,UV _Phys=',
72	& Shap_funct, nShapT, nShapUV, nShapTrPhys, nShapUVPhys
73	write(0,*) 'Shap_Trtau,Shap_uvtau=',Shap_Trtau,Shap_uvtau
74	C-- Close the open data file
75	CLOSE(iUnit)
76	_END_MASTER(myThid)
77
78	C-- Everyone else must wait for the parameters to be loaded
79	_BARRIER
80
81	#endif /* ALLOW_SHAP_FILT */
82	RETURN
83	END