pkg/kl10/kl10_readparms.F

d8d1486ca1 Jean*0001 #include "KL10_OPTIONS.h"
                0002
                0003 CBOP
                0004 C !ROUTINE: KL10_READPARMS
                0005
                0006 C !INTERFACE: ==========================================================
                0007       SUBROUTINE KL10_READPARMS( myThid )
                0008
                0009 C !DESCRIPTION:
                0010 C     Initialize KL10 parameters, read in data.kl10
                0011
                0012 C !USES: ===============================================================
                0013       IMPLICIT NONE
                0014 #include "SIZE.h"
                0015 #include "EEPARAMS.h"
                0016 #include "KL10.h"
                0017 #include "PARAMS.h"
                0018
                0019 C !INPUT PARAMETERS: ===================================================
                0020 C  myThid               :: thread number
                0021       INTEGER myThid
                0022
                0023 C !OUTPUT PARAMETERS: ==================================================
                0024 C  none
                0025
                0026 #ifdef ALLOW_KL10
                0027
                0028 C !LOCAL VARIABLES: ====================================================
                0029 C  iUnit                :: unit number for I/O
                0030 C  msgBuf               :: message buffer
                0031       INTEGER iUnit
                0032       CHARACTER*(MAX_LEN_MBUF) msgBuf
                0033 CEOP
                0034
                0035       NAMELIST /KL10_PARM01/
                0036      &     KLviscMax,
                0037      &     KLdumpFreq,
                0038      &     KLwriteState,
                0039      &     KLtaveFreq
                0040
                0041       IF ( .NOT.useKL10 ) THEN
                0042 C-    pkg KL10 is not used
                0043         _BEGIN_MASTER(myThid)
                0044 C-    Track pkg activation status:
                0045          KL10isON = .FALSE.
                0046 C     print a (weak) warning if data.kl10 is found
                0047          CALL PACKAGES_UNUSED_MSG( 'useKL10', ' ', ' ' )
                0048         _END_MASTER(myThid)
                0049         RETURN
                0050       ENDIF
                0051
                0052       _BEGIN_MASTER(myThid)
                0053
                0054 C This routine has been called by the main model so we set our
                0055 C internal flag to indicate we are in business
                0056       KL10isON = .TRUE.
                0057
                0058 C Set defaults values for parameters in KL10.h
                0059       KLviscMax    = 300. _d 0
                0060       KLdumpFreq   = dumpFreq
                0061       KLtaveFreq   = taveFreq
                0062       KLwriteState = .FALSE.
                0063
                0064 C Open and read the data.kl10 file
                0065       WRITE(msgBuf,'(A)') ' KL10_READPARMS: opening data.kl10'
                0066       CALL PRINT_MESSAGE(msgBuf, standardMessageUnit,
                0067      &                   SQUEEZE_RIGHT , 1)
                0068       CALL OPEN_COPY_DATA_FILE(
                0069      I                   'data.kl10', 'KL10_READPARMS',
                0070      O                   iUnit,
                0071      I                   myThid )
                0072       READ(UNIT=iUnit,NML=KL10_PARM01)
                0073       WRITE(msgBuf,'(A)')
                0074      &  ' KL10_READPARMS: finished reading data.kl10'
                0075       CALL PRINT_MESSAGE(msgBuf, standardMessageUnit,
                0076      &                   SQUEEZE_RIGHT , 1)
                0077
                0078 C Close the open data file
7a77863887 Mart*0079 #ifdef SINGLE_DISK_IO
d8d1486ca1 Jean*0080       CLOSE(iUnit)
7a77863887 Mart*0081 #else
                0082       CLOSE(iUnit,STATUS='DELETE')
                0083 #endif /* SINGLE_DISK_IO */
d8d1486ca1 Jean*0084       _END_MASTER(myThid)
                0085
                0086 C Everyone else must wait for the parameters to be loaded
                0087       _BARRIER
                0088
                0089 C Now set-up any remaining parameters that result from the input parameters
                0090       IF ( KLviscMax .LE. 0. ) THEN
                0091        WRITE(msgBuf,'(A)') 'KLviscMax must be greater than zero'
                0092        CALL PRINT_ERROR( msgBuf , 1)
                0093        STOP 'ABNORMAL END: S/R KL10_READPARMS'
                0094       ENDIF
                0095 #endif /* ALLOW_KL10 */
                0096
                0097       RETURN
                0098       END