pkg/flt/flt_readparms.F

C $Header: /u/gcmpack/MITgcm/pkg/flt/flt_readparms.F,v 1.1 2009/01/31 22:21:31 jmc Exp $
C $Name:  $

#include "FLT_OPTIONS.h"

      SUBROUTINE FLT_READPARMS( myThid  )

C     ==================================================================
C     SUBROUTINE FLT_READPARMS
C     ==================================================================
C     o read float-pkg parameters from file "data.flt"
C     ==================================================================

C     !USES:
      IMPLICIT NONE

#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "FLT.h"

C     == routine arguments ==
C     myThid - thread number for this instance of the routine.
      INTEGER myThid

C     == local variables ==
      CHARACTER*(MAX_LEN_MBUF) msgBuf
      INTEGER iUnit

      NAMELIST /FLT_NML/ flt_int_traj, flt_int_prof,
     &                   flt_noise,
     &                   flt_file

C     == end of interface ==

      _BEGIN_MASTER(myThid)

C     Set default values.
      flt_int_traj =  3600.
      flt_int_prof = 43200.
      flt_noise    = 0.0
      flt_file     = 'float_pos'

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

C     Read parameters from open data file
      READ(UNIT=iUnit,NML=FLT_NML)
      WRITE(msgBuf,'(A)') ' FLT_READPARMS: finished reading data.flt'
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &                    SQUEEZE_RIGHT, myThid )
C     Close the open data file
      CLOSE(iUnit)

C DO some checks
      IF ( useFLT .AND. useOBCS ) THEN
        WRITE(msgBuf,'(A,A)')
     &   'S/R FLT_READPARMS: floats currently not implemented',
     &   ' with open boundaries'
        CALL PRINT_ERROR( msgBuf , myThid)
        STOP 'ABNORMAL END: S/R FLT_READPARMS'
      ENDIF

      _END_MASTER(myThid)
      _BARRIER

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/flt/flt_readparms.F,v 1.1 2009/01/31 22:21:31 jmc Exp $
2	C $Name: $
3
4	#include "FLT_OPTIONS.h"
5
6	SUBROUTINE FLT_READPARMS( myThid )
7
8	C ==================================================================
9	C SUBROUTINE FLT_READPARMS
10	C ==================================================================
11	C o read float-pkg parameters from file "data.flt"
12	C ==================================================================
13
14	C !USES:
15	IMPLICIT NONE
16
17	#include "SIZE.h"
18	#include "EEPARAMS.h"
19	#include "PARAMS.h"
20	#include "FLT.h"
21
22	C == routine arguments ==
23	C myThid - thread number for this instance of the routine.
24	INTEGER myThid
25
26	C == local variables ==
27	CHARACTER*(MAX_LEN_MBUF) msgBuf
28	INTEGER iUnit
29
30	NAMELIST /FLT_NML/ flt_int_traj, flt_int_prof,
31	& flt_noise,
32	& flt_file
33
34	C == end of interface ==
35
36	_BEGIN_MASTER(myThid)
37
38	C Set default values.
39	flt_int_traj = 3600.
40	flt_int_prof = 43200.
41	flt_noise = 0.0
42	flt_file = 'float_pos'
43
44	WRITE(msgBuf,'(A)') ' FLT_READPARMS: opening data.flt'
45	CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
46	& SQUEEZE_RIGHT, myThid )
47	CALL OPEN_COPY_DATA_FILE(
48	I 'data.flt', 'FLT_READPARMS',
49	O iUnit,
50	I myThid )
51
52	C Read parameters from open data file
53	READ(UNIT=iUnit,NML=FLT_NML)
54	WRITE(msgBuf,'(A)') ' FLT_READPARMS: finished reading data.flt'
55	CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
56	& SQUEEZE_RIGHT, myThid )
57	C Close the open data file
58	CLOSE(iUnit)
59
60	C DO some checks
61	IF ( useFLT .AND. useOBCS ) THEN
62	WRITE(msgBuf,'(A,A)')
63	& 'S/R FLT_READPARMS: floats currently not implemented',
64	& ' with open boundaries'
65	CALL PRINT_ERROR( msgBuf , myThid)
66	STOP 'ABNORMAL END: S/R FLT_READPARMS'
67	ENDIF
68
69	_END_MASTER(myThid)
70	_BARRIER
71
72	RETURN
73	END