/[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.9 by heimbach, Wed Sep 22 20:44:37 2004 UTC revision 1.16 by dimitri, Fri Oct 26 02:00:48 2007 UTC
# Line 19  C     === Global variables === Line 19  C     === Global variables ===
19  #ifdef ALLOW_ORLANSKI  #ifdef ALLOW_ORLANSKI
20  #include "ORLANSKI.h"  #include "ORLANSKI.h"
21  #endif  #endif
22    #ifdef ALLOW_PTRACERS.h
23    #include "PTRACERS_SIZE.h"
24    #include "OBCS_PTRACERS.h"
25    #endif /* ALLOW_PTRACERS */
26    
27  C     === Routine arguments ===  C     === Routine arguments ===
28        INTEGER myThid        INTEGER myThid
29    
30  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
31    
32    C     === Local variables ===
33    C     msgBuf      - Informational/error meesage buffer
34    C     iUnit       - Work variable for IO unit number
35          CHARACTER*(MAX_LEN_MBUF) msgBuf
36          INTEGER iUnit
37          INTEGER I,J,iTracer
38    
39        NAMELIST /OBCS_PARM01/        NAMELIST /OBCS_PARM01/
40       &          OB_Jnorth,OB_Jsouth,OB_Ieast,OB_Iwest,       &          OB_Jnorth,OB_Jsouth,OB_Ieast,OB_Iwest,
41       &          useOrlanskiNorth,useOrlanskiSouth,       &          useOrlanskiNorth,useOrlanskiSouth,
42       &          useOrlanskiEast,useOrlanskiWest,       &          useOrlanskiEast,useOrlanskiWest,
43       &          OBNuFile,OBNvFile,OBNtFile,OBNsFile,       &          OBNuFile,OBNvFile,OBNtFile,OBNsFile,OBNaFile,OBNhFile,
44       &          OBSuFile,OBSvFile,OBStFile,OBSsFile,       &          OBSuFile,OBSvFile,OBStFile,OBSsFile,OBSaFile,OBShFile,
45       &          OBEuFile,OBEvFile,OBEtFile,OBEsFile,       &          OBEuFile,OBEvFile,OBEtFile,OBEsFile,OBEaFile,OBEhFile,
46       &          OBWuFile,OBWvFile,OBWtFile,OBWsFile,       &          OBWuFile,OBWvFile,OBWtFile,OBWsFile,OBWaFile,OBWhFile,
47       &          useOBCSsponge, useOBCSbalance, useOBCSprescribe       &          OBNslFile,OBSslFile,OBEslFile,OBWslFile,
48         &          OBNsnFile,OBSsnFile,OBEsnFile,OBWsnFile,
49         &          OBNuiceFile,OBSuiceFile,OBEuiceFile,OBWuiceFile,
50         &          OBNviceFile,OBSviceFile,OBEviceFile,OBWviceFile,
51         &          useOBCSsponge, useOBCSbalance, useOBCSprescribe,
52         &          OBCSprintDiags
53    #ifdef ALLOW_PTRACERS
54         &        , OBNptrFile,OBSptrFile,OBEptrFile,OBWptrFile
55    #endif
56    
57  #ifdef ALLOW_ORLANSKI  #ifdef ALLOW_ORLANSKI
58        NAMELIST /OBCS_PARM02/        NAMELIST /OBCS_PARM02/
# Line 47  C     === Routine arguments === Line 66  C     === Routine arguments ===
66       &          spongeThickness       &          spongeThickness
67  #endif  #endif
68    
69  C     === Local variables ===        _BEGIN_MASTER(myThid)
 C     msgBuf      - Informational/error meesage buffer  
 C     iUnit       - Work variable for IO unit number  
       CHARACTER*(MAX_LEN_MBUF) msgBuf  
       INTEGER iUnit  
       INTEGER I,J  
70    
71  C--   OBCS_READPARMS has been called so we know that  C--   OBCS_READPARMS has been called so we know that
72  C     the package is active.  C     the package is active.
73        OBCSIsOn=.TRUE.        OBCSIsOn=.TRUE.
74    
       _BEGIN_MASTER(myThid)  
   
75        WRITE(msgBuf,'(A)') ' OBCS_READPARMS: opening data.obcs'        WRITE(msgBuf,'(A)') ' OBCS_READPARMS: opening data.obcs'
76        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
77       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
# Line 84  C--   Default flags and values for OBCS Line 96  C--   Default flags and values for OBCS
96        useOrlanskiWest =.FALSE.        useOrlanskiWest =.FALSE.
97        useOBCSsponge   =.FALSE.        useOBCSsponge   =.FALSE.
98        useOBCSbalance  =.FALSE.        useOBCSbalance  =.FALSE.
99        useOBCSprescribe = .FALSE.        useOBCSprescribe=.FALSE.
100          OBCSprintDiags  =.TRUE.
101    
102        OBNuFile = ' '        OBNuFile = ' '
103        OBNvFile = ' '        OBNvFile = ' '
104        OBNtFile = ' '        OBNtFile = ' '
105        OBNsFile = ' '        OBNsFile = ' '
106          OBNaFile = ' '
107          OBNslFile = ' '
108          OBNsnFile = ' '
109          OBNuiceFile = ' '
110          OBNviceFile = ' '
111          OBNhFile = ' '
112        OBSuFile = ' '        OBSuFile = ' '
113        OBSvFile = ' '        OBSvFile = ' '
114        OBStFile = ' '        OBStFile = ' '
115        OBSsFile = ' '        OBSsFile = ' '
116          OBSaFile = ' '
117          OBShFile = ' '
118          OBSslFile = ' '
119          OBSsnFile = ' '
120          OBSuiceFile = ' '
121          OBSviceFile = ' '
122        OBEuFile = ' '        OBEuFile = ' '
123        OBEvFile = ' '        OBEvFile = ' '
124        OBEtFile = ' '        OBEtFile = ' '
125        OBEsFile = ' '        OBEsFile = ' '
126          OBEaFile = ' '
127          OBEhFile = ' '
128          OBEslFile = ' '
129          OBEsnFile = ' '
130          OBEuiceFile = ' '
131          OBEviceFile = ' '
132        OBWuFile = ' '        OBWuFile = ' '
133        OBWvFile = ' '        OBWvFile = ' '
134        OBWtFile = ' '        OBWtFile = ' '
135        OBWsFile = ' '        OBWsFile = ' '
136          OBWaFile = ' '
137          OBWhFile = ' '
138          OBWslFile = ' '
139          OBWsnFile = ' '
140          OBWuiceFile = ' '
141          OBWviceFile = ' '
142    #ifdef ALLOW_PTRACERS
143          DO iTracer = 1, PTRACERS_num
144           OBNptrFile(iTracer) = ' '
145           OBSptrFile(iTracer) = ' '
146           OBEptrFile(iTracer) = ' '
147           OBWptrFile(iTracer) = ' '
148          ENDDO
149    #endif
150    
151  C--   Read parameters from open data file  C--   Read parameters from open data file
152        READ(UNIT=iUnit,NML=OBCS_PARM01)        READ(UNIT=iUnit,NML=OBCS_PARM01)
# Line 128  C     Default Orlanski radiation paramet Line 173  C     Default Orlanski radiation paramet
173        IF (useOrlanskiNorth.OR.        IF (useOrlanskiNorth.OR.
174       &    useOrlanskiSouth.OR.       &    useOrlanskiSouth.OR.
175       &    useOrlanskiEast.OR.       &    useOrlanskiEast.OR.
176       &    useOrlanskiWest)       &    useOrlanskiWest)
177       & READ(UNIT=iUnit,NML=OBCS_PARM02)       & READ(UNIT=iUnit,NML=OBCS_PARM02)
178  #endif  #endif
179          
180  #ifdef ALLOW_OBCS_SPONGE  #ifdef ALLOW_OBCS_SPONGE
181  C     Default sponge layer parameters  C     Default sponge layer parameters:
182        spongeThickness = 2  C     sponge layer is turned off by default
183        Urelaxobcsinner = 5. _d 0        spongeThickness = 0
184        Urelaxobcsbound = 1. _d 0        Urelaxobcsinner = 0. _d 0
185        Vrelaxobcsinner = 5. _d 0        Urelaxobcsbound = 0. _d 0
186        Vrelaxobcsbound = 1. _d 0        Vrelaxobcsinner = 0. _d 0
187          Vrelaxobcsbound = 0. _d 0
188    CML this was the previous default in units of days
189    CML      spongeThickness = 2
190    CML      Urelaxobcsinner = 5. _d 0
191    CML      Urelaxobcsbound = 1. _d 0
192    CML      Vrelaxobcsinner = 5. _d 0
193    CML      Vrelaxobcsbound = 1. _d 0
194        IF (useOBCSsponge)        IF (useOBCSsponge)
195       & READ(UNIT=iUnit,NML=OBCS_PARM03)       & READ(UNIT=iUnit,NML=OBCS_PARM03)
196  #endif  #endif
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22