C $Header: /u/gcmpack/MITgcm/pkg/atm_ocn_coupler/cpl_read_params.F,v 1.6 2016/01/06 00:32:44 jmc Exp $ C $Name: $ #include "CPP_OPTIONS.h" CBOP 0 C !ROUTINE: CPL_READ_PARAMS C !INTERFACE: SUBROUTINE CPL_READ_PARAMS( msgUnit ) C !DESCRIPTION: C *==========================================================* C | SUBROUTINE CPL_READ_PARAMS C | o Read coupler and mapping parameters C *==========================================================* C !USES: IMPLICIT NONE C == Global variables == #include "CPL_PARAMS.h" C !INPUT/OUTPUT PARAMETERS: C msgUnit :: log-file I/O unit INTEGER msgUnit C !LOCAL VARIABLES: INTEGER ioUnit LOGICAL fileExist CEOP NAMELIST //COUPLER_PARAMS & cpl_sequential, & cpl_exchange_RunOff, & cpl_exchange1W_sIce, cpl_exchange2W_sIce, & cpl_exchange_SaltPl, cpl_exchange_DIC, & runOffMapSize, runOffMapFile C-- Default values for Coupler parameters cpl_sequential = 0 C default coupling-exch selector value: fit a simple warm C aqua-planet (no land, no seaice, no CO2/DIC coupling) cpl_exchange_RunOff = 0 cpl_exchange1W_sIce = 0 cpl_exchange2W_sIce = 0 cpl_exchange_SaltPl = 0 cpl_exchange_DIC = 0 runOffMapSize = 0 runOffMapFile = ' ' C- other parameters: cplErrorCount = 0 fileExist = .FALSE. C-- Read-in parameter file: INQUIRE( FILE='data.cpl', EXIST=fileExist ) IF ( fileExist ) THEN WRITE(msgUnit,'(2A)') 'CPL_READ_PARAMS: ', & 'Reading parameter file "data.cpl"' ioUnit = 88 OPEN( ioUnit, FILE='data.cpl',STATUS='old') READ( ioUnit, COUPLER_PARAMS ) CLOSE(ioUnit ) ELSE cplErrorCount = 1 WRITE(msgUnit,'(2A)') ' *** ERROR *** CPL_READ_PARAMS: ', & 'parameter file "data.cpl" not found' WRITE(msgUnit,'(2A,I4,A)') ' *** ERROR *** CPL_READ_PARAMS: ', & '==> set ErrorCount=', cplErrorCount, ' (Fatal)' ENDIF C-- Check parameter consistency: IF ( cpl_exchange2W_sIce.EQ.3 .AND. cpl_sequential.EQ.0 ) THEN cplErrorCount = cplErrorCount + 1 WRITE(msgUnit,'(2A)') ' *** ERROR *** CPL_READ_PARAMS: ', & 'needs "cpl_sequential=1" to use: cpl_exchange2W_sIce=3' WRITE(msgUnit,'(2A,I4,A)') ' *** ERROR *** CPL_READ_PARAMS: ', & '==> set ErrorCount=', cplErrorCount, ' (Fatal)' ENDIF IF ( cpl_exchange2W_sIce.EQ.1 .AND. cpl_sequential.EQ.1 ) THEN cplErrorCount = cplErrorCount + 1 WRITE(msgUnit,'(2A)') ' *** ERROR *** CPL_READ_PARAMS: ', & 'needs "cpl_sequential=0" to use: cpl_exchange2W_sIce=1' WRITE(msgUnit,'(2A,I4,A)') ' *** ERROR *** CPL_READ_PARAMS: ', & '==> set ErrorCount=', cplErrorCount, ' (Fatal)' ENDIF C-- Print summary of Coupler parameters value IF ( cplErrorCount.EQ.0 ) THEN WRITE(msgUnit,'(2A)') 'CPL_READ_PARAMS: ', & '------ Coupler parameter Summary: -----------------' C- main coupling time-stepping selector: WRITE(msgUnit,'(2A,I4)') 'CPL_READ_PARAMS: ', & ' cpl_sequential = ', cpl_sequential C- cpl_exchange field selector: WRITE(msgUnit,'(2A,I4)') 'CPL_READ_PARAMS: ', & ' cpl_exchange_RunOff =', cpl_exchange_RunOff WRITE(msgUnit,'(2A,I4)') 'CPL_READ_PARAMS: ', & ' cpl_exchange1W_sIce =', cpl_exchange1W_sIce WRITE(msgUnit,'(2A,I4)') 'CPL_READ_PARAMS: ', & ' cpl_exchange2W_sIce =', cpl_exchange2W_sIce WRITE(msgUnit,'(2A,I4)') 'CPL_READ_PARAMS: ', & ' cpl_exchange_SaltPl =', cpl_exchange_SaltPl WRITE(msgUnit,'(2A,I4)') 'CPL_READ_PARAMS: ', & ' cpl_exchange_DIC =', cpl_exchange_DIC C- run-off mapping parameters: WRITE(msgUnit,'(2A,I8)') 'CPL_READ_PARAMS: ', & ' runOffMapSize=', runOffMapSize WRITE(msgUnit,'(3A)') 'CPL_READ_PARAMS: ', & ' runOffMapFile=', runOffMapFile WRITE(msgUnit,'(2A)') 'CPL_READ_PARAMS: ', & '------ End of Coupler parameter Summary ------------' ENDIF RETURN END