pkg/matrix/matrix_readparms.F

5d0c92356b Jean*0001 #include "MATRIX_OPTIONS.h"
a9544d2f65 Sama*0002
8385b3bb5f Ed H*0003 C     !ROUTINE: MATRIX_READPARMS
a9544d2f65 Sama*0004
                0005       SUBROUTINE MATRIX_READPARMS( myThid )
                0006
8385b3bb5f Ed H*0007 C     !DESCRIPTION:
a9544d2f65 Sama*0008 C     Initialize MATRIX parameters, read in data.matrix
                0009
                0010       IMPLICIT NONE
                0011 #include "SIZE.h"
                0012 #include "EEPARAMS.h"
5ecb7b2fa8 Jean*0013 #include "PARAMS.h"
a9544d2f65 Sama*0014 #include "PTRACERS_SIZE.h"
                0015 #include "MATRIX.h"
                0016
8385b3bb5f Ed H*0017 C     !INPUT PARAMETERS:
                0018 C     myThid            :: thread number
a9544d2f65 Sama*0019       INTEGER myThid
                0020
8385b3bb5f Ed H*0021 C     !OUTPUT PARAMETERS:
                0022 C     none
5ecb7b2fa8 Jean*0023
                0024 #ifdef ALLOW_MATRIX
8385b3bb5f Ed H*0025 C     !LOCAL VARIABLES:
8f6dc1243d Jean*0026 C     iUnit  :: unit number for I/O
                0027 C     msgBuf :: message buffer
a9544d2f65 Sama*0028       INTEGER iUnit
                0029       CHARACTER*(MAX_LEN_MBUF) msgBuf
                0030
8385b3bb5f Ed H*0031 C     Open and read the data.matrix file
a9544d2f65 Sama*0032
8f6dc1243d Jean*0033       NAMELIST /MATRIX_PARM01/ expMatrixWriteTime, impMatrixWriteTime
                0034
5ecb7b2fa8 Jean*0035       IF ( .NOT.useMATRIX ) THEN
                0036 C-    pkg MATRIX is not used
                0037         _BEGIN_MASTER(myThid)
                0038 C-    Track pkg activation status:
                0039 C     print a (weak) warning if data.matrix is found
                0040          CALL PACKAGES_UNUSED_MSG( 'useMATRIX', ' ', ' ' )
                0041         _END_MASTER(myThid)
                0042         RETURN
                0043       ENDIF
                0044
5d0c92356b Jean*0045       _BEGIN_MASTER(myThid)
                0046
a9544d2f65 Sama*0047       expMatrixWriteTime = 10. _d 10
                0048       impMatrixWriteTime = 10. _d 10
8f6dc1243d Jean*0049
a9544d2f65 Sama*0050       WRITE(msgBuf,'(A)') ' MATRIX_READPARMS: opening data.matrix'
                0051       CALL PRINT_MESSAGE(msgBuf, standardMessageUnit,
                0052      &                   SQUEEZE_RIGHT , 1)
                0053       CALL OPEN_COPY_DATA_FILE(
                0054      I                   'data.matrix', 'MATRIX_READPARMS',
                0055      O                   iUnit,
                0056      I                   myThid )
                0057       READ(UNIT=iUnit,NML=MATRIX_PARM01)
                0058       WRITE(msgBuf,'(A)')
                0059      &  ' MATRIX_READPARMS: finished reading data.matrix'
                0060       CALL PRINT_MESSAGE(msgBuf, standardMessageUnit,
                0061      &                   SQUEEZE_RIGHT , 1)
                0062
                0063 C Close the open data file
7a77863887 Mart*0064 #ifdef SINGLE_DISK_IO
a9544d2f65 Sama*0065       CLOSE(iUnit)
7a77863887 Mart*0066 #else
                0067       CLOSE(iUnit,STATUS='DELETE')
                0068 #endif /* SINGLE_DISK_IO */
a9544d2f65 Sama*0069       _END_MASTER(myThid)
                0070
                0071 C Everyone else must wait for the parameters to be loaded
                0072       _BARRIER
                0073
                0074 #endif /* ALLOW_MATRIX */
                0075       RETURN
                0076       END