C $Header: /u/gcmpack/MITgcm/pkg/obcs/obcs_readparms.F,v 1.26 2010/01/31 20:25:49 jmc Exp $
C $Name: $
#include "OBCS_OPTIONS.h"
SUBROUTINE OBCS_READPARMS( myThid )
C *==========================================================*
C | SUBROUTINE OBCS_READPARMS
C | o Routine to initialize OBCS variables and constants.
C *==========================================================*
C *==========================================================*
IMPLICIT NONE
C === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "OBCS.h"
#ifdef ALLOW_ORLANSKI
#include "ORLANSKI.h"
#endif
#ifdef ALLOW_PTRACERS
#include "PTRACERS_SIZE.h"
#include "OBCS_PTRACERS.h"
#endif /* ALLOW_PTRACERS */
C === Routine arguments ===
INTEGER myThid
#ifdef ALLOW_OBCS
C === Local variables ===
C msgBuf - Informational/error message buffer
C iUnit - Work variable for IO unit number
CHARACTER*(MAX_LEN_MBUF) msgBuf
INTEGER iUnit
INTEGER I, J
INTEGER bi, bj, iG, jG, iGm, jGm
#ifdef ALLOW_PTRACERS
INTEGER iTracer
#endif
NAMELIST //OBCS_PARM01
& OB_Jnorth,OB_Jsouth,OB_Ieast,OB_Iwest,
& useOrlanskiNorth,useOrlanskiSouth,
& useOrlanskiEast,useOrlanskiWest,
& OBNuFile,OBNvFile,OBNtFile,OBNsFile,OBNaFile,OBNhFile,
& OBSuFile,OBSvFile,OBStFile,OBSsFile,OBSaFile,OBShFile,
& OBEuFile,OBEvFile,OBEtFile,OBEsFile,OBEaFile,OBEhFile,
& OBWuFile,OBWvFile,OBWtFile,OBWsFile,OBWaFile,OBWhFile,
& OBNslFile,OBSslFile,OBEslFile,OBWslFile,
& OBNsnFile,OBSsnFile,OBEsnFile,OBWsnFile,
& OBNuiceFile,OBSuiceFile,OBEuiceFile,OBWuiceFile,
& OBNviceFile,OBSviceFile,OBEviceFile,OBWviceFile,
& OBNetaFile, OBSetaFile, OBEetaFile, OBWetaFile,
& OBNwFile, OBSwFile, OBEwFile, OBWwFile,
#ifdef ALLOW_PTRACERS
& OBNptrFile,OBSptrFile,OBEptrFile,OBWptrFile,
#endif
& useOBCSsponge, useOBCSbalance, useOBCSprescribe,
& OBCSprintDiags, useOBCSYearlyFields, OBCSfixTopo
#ifdef ALLOW_ORLANSKI
NAMELIST //OBCS_PARM02
& CMAX, cvelTimeScale, CFIX, useFixedCEast, useFixedCWest
#endif
#ifdef ALLOW_OBCS_SPONGE
NAMELIST //OBCS_PARM03
& Urelaxobcsinner,Urelaxobcsbound,
& Vrelaxobcsinner,Vrelaxobcsbound,
& spongeThickness
#endif
_BEGIN_MASTER(myThid)
C-- OBCS_READPARMS has been called so we know that
C the package is active.
c OBCSIsOn=.TRUE.
IF ( debugLevel .GE. debLevB ) THEN
WRITE(msgBuf,'(A)') ' OBCS_READPARMS: opening data.obcs'
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
& SQUEEZE_RIGHT , 1)
ENDIF
CALL OPEN_COPY_DATA_FILE(
I 'data.obcs', 'OBCS_READPARMS',
O iUnit,
I myThid )
C-- Default flags and values for OBCS
DO I=1,Nx
OB_Jnorth(I)=0
OB_Jsouth(I)=0
ENDDO
DO J=1,Ny
OB_Ieast(J)=0
OB_Iwest(J)=0
ENDDO
useOrlanskiNorth =.FALSE.
useOrlanskiSouth =.FALSE.
useOrlanskiEast =.FALSE.
useOrlanskiWest =.FALSE.
useOBCSsponge =.FALSE.
useOBCSbalance =.FALSE.
useOBCSprescribe =.FALSE.
OBCSprintDiags = debugLevel.GE.debLevB
useOBCSYearlyFields=.FALSE.
OBCSfixTopo =.TRUE.
OBNuFile = ' '
OBNvFile = ' '
OBNtFile = ' '
OBNsFile = ' '
OBNaFile = ' '
OBNslFile = ' '
OBNsnFile = ' '
OBNuiceFile = ' '
OBNviceFile = ' '
OBNhFile = ' '
OBSuFile = ' '
OBSvFile = ' '
OBStFile = ' '
OBSsFile = ' '
OBSaFile = ' '
OBShFile = ' '
OBSslFile = ' '
OBSsnFile = ' '
OBSuiceFile = ' '
OBSviceFile = ' '
OBEuFile = ' '
OBEvFile = ' '
OBEtFile = ' '
OBEsFile = ' '
OBEaFile = ' '
OBEhFile = ' '
OBEslFile = ' '
OBEsnFile = ' '
OBEuiceFile = ' '
OBEviceFile = ' '
OBWuFile = ' '
OBWvFile = ' '
OBWtFile = ' '
OBWsFile = ' '
OBWaFile = ' '
OBWhFile = ' '
OBWslFile = ' '
OBWsnFile = ' '
OBWuiceFile = ' '
OBWviceFile = ' '
OBNetaFile = ' '
OBSetaFile = ' '
OBEetaFile = ' '
OBWetaFile = ' '
OBNwFile = ' '
OBSwFile = ' '
OBEwFile = ' '
OBWwFile = ' '
#ifdef ALLOW_PTRACERS
DO iTracer = 1, PTRACERS_num
OBNptrFile(iTracer) = ' '
OBSptrFile(iTracer) = ' '
OBEptrFile(iTracer) = ' '
OBWptrFile(iTracer) = ' '
ENDDO
#endif
C-- Read parameters from open data file
READ(UNIT=iUnit,NML=OBCS_PARM01)
C Account for periodicity if negative indices were supplied
DO J=1,Ny
IF (OB_Ieast(J).lt.0) OB_Ieast(J)=OB_Ieast(J)+Nx+1
ENDDO
DO I=1,Nx
IF (OB_Jnorth(I).lt.0) OB_Jnorth(I)=OB_Jnorth(I)+Ny+1
ENDDO
IF ( debugLevel .GE. debLevB ) THEN
write(*,*) 'OB Jn =',OB_Jnorth
write(*,*) 'OB Js =',OB_Jsouth
write(*,*) 'OB Ie =',OB_Ieast
write(*,*) 'OB Iw =',OB_Iwest
ENDIF
#ifdef ALLOW_ORLANSKI
C Default Orlanski radiation parameters
CMAX = 0.45 _d 0 /* maximum allowable phase speed-CFL for AB-II */
cvelTimeScale = 2000.0 _d 0 /* Averaging period for phase speed in sec. */
CFIX = 0.8 _d 0 /* Fixed boundary phase speed in m/s */
useFixedCEast=.FALSE.
useFixedCWest=.FALSE.
IF (useOrlanskiNorth.OR.
& useOrlanskiSouth.OR.
& useOrlanskiEast.OR.
& useOrlanskiWest)
& READ(UNIT=iUnit,NML=OBCS_PARM02)
#endif
#ifdef ALLOW_OBCS_SPONGE
C Default sponge layer parameters:
C sponge layer is turned off by default
spongeThickness = 0
Urelaxobcsinner = 0. _d 0
Urelaxobcsbound = 0. _d 0
Vrelaxobcsinner = 0. _d 0
Vrelaxobcsbound = 0. _d 0
CML this was the previous default in units of days
CML spongeThickness = 2
CML Urelaxobcsinner = 5. _d 0
CML Urelaxobcsbound = 1. _d 0
CML Vrelaxobcsinner = 5. _d 0
CML Vrelaxobcsbound = 1. _d 0
IF (useOBCSsponge)
& READ(UNIT=iUnit,NML=OBCS_PARM03)
#endif
IF ( debugLevel .GE. debLevB ) THEN
WRITE(msgBuf,'(A)') ' OBCS_READPARMS: finished reading data.obcs'
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
& SQUEEZE_RIGHT , 1)
ENDIF
C-- Close the open data file
CLOSE(iUnit)
_END_MASTER(myThid)
C-- Everyone else must wait for the parameters to be loaded
_BARRIER
C-- Calculate the tiled index arrays OB_Jn/Js/Ie/Iw here from the
C global arrays OB_Jnorth/Jsouth/Ieast/Iwest.
C Note: This part of the code has been moved from obcs_init_fixed to
C routine routine because the OB_Jn/Js/Ie/Iw index arrays are
C required by ini_depth which is called befoer obcs_init_fixed
DO bj = myByLo(myThid), myByHi(myThid)
DO bi = myBxLo(myThid), myBxHi(myThid)
DO I=1-Olx,sNx+Olx
OB_Jn(I,bi,bj)=0
OB_Js(I,bi,bj)=0
ENDDO
DO J=1-Oly,sNy+Oly
OB_Ie(J,bi,bj)=0
OB_Iw(J,bi,bj)=0
ENDDO
DO J=1-Oly,sNy+Oly
C convert from local y index J to global y index jG
jG = myYGlobalLo-1+(bj-1)*sNy+J
C use periodicity to deal with out of range points caused by the overlaps.
C they will be excluded by the mask in any case, but this saves array
C out-of-bounds errors here.
jGm = 1+mod( jG-1+Ny , Ny )
C loop over local x index I
DO I=1,sNx
iG = myXGlobalLo-1+(bi-1)*sNx+I
iGm = 1+mod( iG-1+Nx , Nx )
C OB_Ieast(jGm) allows for the eastern boundary to be at variable x locations
IF (iG.EQ.OB_Ieast(jGm)) OB_Ie(J,bi,bj)=I
IF (iG.EQ.OB_Iwest(jGm)) OB_Iw(J,bi,bj)=I
ENDDO
ENDDO
DO J=1,sNy
jG = myYGlobalLo-1+(bj-1)*sNy+J
jGm = 1+mod( jG-1+Ny , Ny )
DO I=1-Olx,sNx+Olx
iG = myXGlobalLo-1+(bi-1)*sNx+I
iGm = 1+mod( iG-1+Nx , Nx )
C OB_Jnorth(iGm) allows for the northern boundary to be at variable y locations
IF (jG.EQ.OB_Jnorth(iGm)) OB_Jn(I,bi,bj)=J
IF (jG.EQ.OB_Jsouth(iGm)) OB_Js(I,bi,bj)=J
ENDDO
ENDDO
C bi,bj-loops
ENDDO
ENDDO
#endif /* ALLOW_OBCS */
RETURN
END