pkg/bbl/bbl_readparms.F

53b68b7823 Dimi*0001 #include "BBL_OPTIONS.h"
                0002
                0003 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0004 CBOP 0
                0005 C !ROUTINE: BBL_READPARMS
                0006
                0007 C !INTERFACE:
                0008       SUBROUTINE BBL_READPARMS( myThid )
                0009
                0010 C     !DESCRIPTION:
                0011 C     Initialize BBL variables and constants.
                0012
                0013 C     !USES:
                0014       IMPLICIT NONE
                0015 #include "SIZE.h"
                0016 #include "EEPARAMS.h"
                0017 #include "PARAMS.h"
                0018 #include "BBL.h"
                0019
                0020 C     !INPUT PARAMETERS:
                0021       INTEGER myThid
                0022 CEOP
                0023
                0024 #ifdef ALLOW_BBL
                0025
                0026 C     !LOCAL VARIABLES:
                0027 C     msgBuf     :: Informational/error message buffer
                0028 C     iUnit      :: Work variable for IO unit number
                0029       CHARACTER*(MAX_LEN_MBUF) msgBuf
                0030       INTEGER iUnit
                0031
                0032       NAMELIST /BBL_PARM01/
5f5cc705f4 Dimi*0033      &     bbl_wvel, bbl_hvel, bbl_initEta,
53b68b7823 Dimi*0034      &     bbl_thetaFile, bbl_saltFile, bbl_etaFile
                0035
ae4c29e0db Jean*0036 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0037
                0038       IF ( .NOT.useBBL ) THEN
                0039 C-    pkg BBL is not used
                0040         _BEGIN_MASTER(myThid)
                0041 C-    Track pkg activation status:
                0042 C     print a (weak) warning if data.bbl is found
                0043          CALL PACKAGES_UNUSED_MSG( 'useBBL', ' ', ' ' )
                0044         _END_MASTER(myThid)
                0045         RETURN
                0046       ENDIF
                0047
53b68b7823 Dimi*0048       _BEGIN_MASTER(myThid)
                0049
                0050 C--   Default values for BBL
5f5cc705f4 Dimi*0051       bbl_wvel      = 1. _d -5
                0052       bbl_hvel      = 1. _d 0
                0053       bbl_initEta   = 20. _d 0
53b68b7823 Dimi*0054
                0055 C-    file names for initial conditions:
                0056       bbl_thetaFile = ' '
                0057       bbl_saltFile  = ' '
                0058       bbl_etaFile   = ' '
                0059
                0060       WRITE(msgBuf,'(A)') 'BBL_READPARMS: opening data.bbl'
                0061       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
                0062      &                    SQUEEZE_RIGHT , 1)
                0063       CALL OPEN_COPY_DATA_FILE(
                0064      I                     'data.bbl', 'BBL_READPARMS',
                0065      O                     iUnit,
                0066      I                     myThid )
                0067
                0068 C     Read parameters from open data file
                0069       READ(UNIT=iUnit,NML=BBL_PARM01)
                0070       WRITE(msgBuf,'(A)')
                0071      &    'BBL_READPARMS: finished reading data.bbl'
                0072       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
                0073      &    SQUEEZE_RIGHT , 1)
                0074 C     Close the open data file
7a77863887 Mart*0075 #ifdef SINGLE_DISK_IO
53b68b7823 Dimi*0076       CLOSE(iUnit)
7a77863887 Mart*0077 #else
                0078       CLOSE(iUnit,STATUS='DELETE')
                0079 #endif /* SINGLE_DISK_IO */
53b68b7823 Dimi*0080
                0081       _END_MASTER(myThid)
                0082
                0083 C--   Everyone else must wait for the parameters to be loaded
                0084       _BARRIER
                0085
                0086 #endif /* ALLOW_BBL */
                0087
                0088       RETURN
                0089       END