/[MITgcm]/MITgcm/pkg/ggl90/ggl90_readparms.F
ViewVC logotype

Diff of /MITgcm/pkg/ggl90/ggl90_readparms.F

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

revision 1.2 by mlosch, Mon Sep 27 08:02:04 2004 UTC revision 1.6 by mlosch, Thu Jun 8 00:12:45 2006 UTC
# Line 36  C--   GGL90 vertical mixing parameters Line 36  C--   GGL90 vertical mixing parameters
36       &     GGL90diffTKEh,       &     GGL90diffTKEh,
37       &     GGL90mixingMaps, GGL90writeState,       &     GGL90mixingMaps, GGL90writeState,
38       &     GGL90ck, GGL90ceps, GGL90alpha, GGL90m2,       &     GGL90ck, GGL90ceps, GGL90alpha, GGL90m2,
39       &     GGL90TKEmin, GGL90TKEbottom,       &     GGL90TKEmin, GGL90TKEsurfMin, GGL90TKEbottom,
40       &     GGL90mixingLengthMin, GGL90viscMax, GGL90diffMax,       &     GGL90mixingLengthMin, GGL90viscMax, GGL90diffMax,
41       &     GGL90TKEFile       &     GGL90TKEFile
42        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
# Line 59  C--   set default TKE vertical mixing pa Line 59  C--   set default TKE vertical mixing pa
59        GGL90ck              = 0.1 _d 0        GGL90ck              = 0.1 _d 0
60        GGL90ceps            = 0.7 _d 0        GGL90ceps            = 0.7 _d 0
61        GGL90alpha           = 1.0 _d 0        GGL90alpha           = 1.0 _d 0
62        GGL90m2              = 3.0 _d 0  C     Blanke and Delecluse (1993, JPO) use
63          GGL90m2              = 3.75 _d 0
64        GGL90TKEmin          = 1.0 _d -11        GGL90TKEmin          = 1.0 _d -11
65    C     Blanke and Delecluse (1993, JPO) use
66          GGL90TKEsurfMin      = 1.0 _d -04
67        GGL90TKEbottom       = UNSET_RL        GGL90TKEbottom       = UNSET_RL
68        GGL90viscMax         = 1.0 _d 0        GGL90viscMax         = 1.0 _d 0
69        GGL90diffMax         = 1.0 _d 1        GGL90diffMax         = 1.0 _d 1
# Line 99  C      CALL MODELDATA_EXAMPLE( myThid ) Line 102  C      CALL MODELDATA_EXAMPLE( myThid )
102        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
103       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
104    
       _END_MASTER(myThid)  
   
 C--   Everyone else must wait for the parameters to be loaded  
       _BARRIER  
   
105  C Now set-up any remaining parameters that result from the input parameters  C Now set-up any remaining parameters that result from the input parameters
106        IF ( GGL90TKEbottom .EQ. UNSET_RL ) THEN        IF ( GGL90TKEbottom .EQ. UNSET_RL ) THEN
107         GGL90TKEbottom = GGL90TKEmin         GGL90TKEbottom = GGL90TKEmin
# Line 137  C Now set-up any remaining parameters th Line 135  C Now set-up any remaining parameters th
135         STOP 'ABNORMAL END: S/R GGL90_READPARMS'         STOP 'ABNORMAL END: S/R GGL90_READPARMS'
136        ENDIF        ENDIF
137    
138    C--   print TKE vertical mixing parameters to stdout for better debugging
139          WRITE(msgBuf,'(A)')
140         &'// ======================================================='
141          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
142         &                    SQUEEZE_RIGHT , 1)
143          WRITE(msgBuf,'(A)') '// GGL90 configuration'
144          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
145         &                    SQUEEZE_RIGHT , 1)
146          WRITE(msgBuf,'(A)')
147         &'// ======================================================='
148          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
149         &  SQUEEZE_RIGHT , 1)
150    
151          CALL WRITE_0D_R8( GGL90dumpFreq, INDEX_NONE,'GGL90dumpFreq =',
152         &'   /* GGL90 state write out interval ( s ). */')
153          CALL WRITE_0D_R8( GGL90taveFreq, INDEX_NONE,'GGL90taveFreq =',
154         &'   /* GGL90 averaging interval ( s ). */')
155          CALL WRITE_0D_L(GGL90mixingMaps,INDEX_NONE,
156         &     'GGL90mixingMAPS =', '   /* GGL90 IO flag. */')
157          CALL WRITE_0D_L(GGL90writeState,INDEX_NONE,
158         &     'GGL90writeState =', '   /* GGL90 IO flag. */')
159          CALL WRITE_0D_R8( GGL90ck, INDEX_NONE,'GGL90ck =',
160         &'   /* GGL90 viscosity parameter. */')
161          CALL WRITE_0D_R8( GGL90ceps, INDEX_NONE,'GGL90ceps =',
162         &'   /* GGL90 dissipation parameter. */')
163          CALL WRITE_0D_R8( GGL90alpha, INDEX_NONE,'GGL90alpha =',
164         &'   /* GGL90 TKE diffusivity parameter. */')
165          CALL WRITE_0D_R8( GGL90m2, INDEX_NONE,'GGL90m2 =',
166         &'   /* GGL90 wind stress to vertical stress ratio. */')
167          CALL WRITE_0D_R8( GGL90TKEmin, INDEX_NONE,'GGL90TKEmin =',
168         &'   /* GGL90 minimum kinetic energy ( m^2/s^2 ). */')
169          CALL WRITE_0D_R8( GGL90TKEsurfMin, INDEX_NONE,
170         &     'GGL90TKEsurfMin =',
171         &'   /* GGL90 minimum surface kinetic energy ( m^2/s^2 ). */')
172          CALL WRITE_0D_R8( GGL90TKEbottom, INDEX_NONE,
173         &     'GGL90TKEbottom =',
174         &     '   /* GGL90 bottom kinetic energy ( m^2/s^2 ). */')
175          CALL WRITE_0D_R8( GGL90viscMax, INDEX_NONE,'GGL90viscMax =',
176         &     '   /* GGL90 upper limit for viscosity ( m^2/s ). */')
177          CALL WRITE_0D_R8( GGL90diffMax, INDEX_NONE,'GGL90diffMax =',
178         &     '   /* GGL90 upper limit for diffusivity ( m^2/s ). */')
179          CALL WRITE_0D_R8( GGL90diffTKEh, INDEX_NONE,'GGL90diffTKEh =',
180         &     '   /* GGL90 horizontal diffusivity for TKE ( m^2/s ). */')
181          CALL WRITE_0D_R8( GGL90mixingLengthMin, INDEX_NONE,
182         &     'GGL90mixingLengthMin =',
183         &     '   /* GGL90 minimum mixing length ( m ). */')
184          WRITE(msgBuf,'(A,A)') 'GGL90: GGL90TKEFile = ', GGL90TKEFile
185          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
186         &                    SQUEEZE_RIGHT , 1)
187    
188          WRITE(msgBuf,'(A)')
189         &'// ======================================================='
190          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
191         &                    SQUEEZE_RIGHT , 1)
192          WRITE(msgBuf,'(A)') '// End of GGL90 config. summary'
193          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
194         &                    SQUEEZE_RIGHT , 1)
195          WRITE(msgBuf,'(A)')
196         &'// ======================================================='
197          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
198         &  SQUEEZE_RIGHT , 1)
199          WRITE(msgBuf,'(A)') ' '
200          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
201         &  SQUEEZE_RIGHT , 1)
202    
203          _END_MASTER(myThid)
204    
205    C--   Everyone else must wait for the parameters to be loaded
206          _BARRIER
207    
208  #endif /* ALLOW_GGL90 */  #endif /* ALLOW_GGL90 */
209    
210        return        return
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22