/[MITgcm]/MITgcm/pkg/obcs/obcs_readparms.F
ViewVC logotype

Diff of /MITgcm/pkg/obcs/obcs_readparms.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by adcroft, Fri Feb 2 21:36:30 2001 UTC revision 1.3 by adcroft, Mon Feb 5 15:45:47 2001 UTC
# Line 25  C     === Routine arguments === Line 25  C     === Routine arguments ===
25    
26  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
27    
       _RL Orlanski_Cmax,Orlanski_TimeScale  
28        NAMELIST /OBCS_PARM01/        NAMELIST /OBCS_PARM01/
29       &          OB_Jnorth,OB_Jsouth,OB_Ieast,OB_Iwest,       &          OB_Jnorth,OB_Jsouth,OB_Ieast,OB_Iwest,
30       &          useOrlanskiNorth,useOrlanskiSouth,       &          useOrlanskiNorth,useOrlanskiSouth,
31       &          useOrlanskiEast,useOrlanskiWest,       &          useOrlanskiEast,useOrlanskiWest
32       &          Orlanski_Cmax,Orlanski_TimeScale  
33    #ifdef ALLOW_ORLANSKI
34          NAMELIST /OBCS_PARM02/
35         & CMAX, cvelTimeScale
36    #endif
37    
38  C     === Local variables ===  C     === Local variables ===
39  C     msgBuf      - Informational/error meesage buffer  C     msgBuf      - Informational/error meesage buffer
# Line 68  C--   Default flags and values for OBCS Line 71  C--   Default flags and values for OBCS
71        useOrlanskiEast =.FALSE.        useOrlanskiEast =.FALSE.
72        useOrlanskiWest =.FALSE.        useOrlanskiWest =.FALSE.
73    
 C--   Defaults for the Orlanksi package  
       Orlanski_Cmax=0.  
       Orlanski_TimeScale=0.  
   
74  C--   Read parameters from open data file  C--   Read parameters from open data file
75        READ(UNIT=iUnit,NML=OBCS_PARM01)        READ(UNIT=iUnit,NML=OBCS_PARM01)
76    
# Line 87  C     Account for periodicity if negativ Line 86  C     Account for periodicity if negativ
86        write(0,*) 'OB Ie =',OB_Ieast        write(0,*) 'OB Ie =',OB_Ieast
87        write(0,*) 'OB Iw =',OB_Iwest        write(0,*) 'OB Iw =',OB_Iwest
88    
89  #ifndef ALLOW_ORLANSKI  #ifdef ALLOW_ORLANSKI
90        IF (  C     Default Orlanski radiation parameters
91       &    Orlanski_Cmax.NE.0. .OR.        CMAX = 0.45 _d 0 /* maximum allowable phase speed-CFL for AB-II */
92       &    Orlanski_TimeScale.NE.0.) THEN        cvelTimeScale = 2000.0 _d 0 /* Averaging period for phase speed in sec. */
93          WRITE(msgBuf,'(A)')        IF (useOrlanskiNorth.OR.
94       &  'OBCS_READPARMS: Orlanski parameters were set!'       &    useOrlanskiSouth.OR.
95          CALL PRINT_ERROR( msgBuf, 1)       &    useOrlanskiEast.OR.
96          STOP 'OBCS_READPARMS: Inconsistent CPP options and parameters'       &    useOrlanskiWest)
97        ENDIF       & READ(UNIT=iUnit,NML=OBCS_PARM02)
 #else  
       Cmax=Orlanski_Cmax  
       cVelTimeScale=Orlanski_TimeScale  
98  #endif  #endif
99          
100        WRITE(msgBuf,'(A)') ' OBCS_READPARMS: finished reading data.obcs'        WRITE(msgBuf,'(A)') ' OBCS_READPARMS: finished reading data.obcs'
101        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
102       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22