/[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.5 by jmc, Wed Jun 7 16:34:24 2006 UTC revision 1.6 by mlosch, Thu Jun 8 00:12:45 2006 UTC
# Line 135  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    
 #ifdef UNSAFE_INTERNAL_WRITE  
138  C--   print TKE vertical mixing parameters to stdout for better debugging  C--   print TKE vertical mixing parameters to stdout for better debugging
139        WRITE(msgBuf,*) 'GGL90: GGL90dumpFreq        = ', GGL90dumpFreq        WRITE(msgBuf,'(A)')
140        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,       &'// ======================================================='
141       &                    SQUEEZE_RIGHT , 1)        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
       WRITE(msgBuf,*) 'GGL90: GGL90taveFreq        = ', GGL90taveFreq  
       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
142       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
143        WRITE(msgBuf,*) 'GGL90: GGL90mixingMaps      = ', GGL90mixingMaps        WRITE(msgBuf,'(A)') '// GGL90 configuration'
144        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
145       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
146        WRITE(msgBuf,*) 'GGL90: GGL90writeState      = ', GGL90writeState        WRITE(msgBuf,'(A)')
147         &'// ======================================================='
148        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
149       &                    SQUEEZE_RIGHT , 1)       &  SQUEEZE_RIGHT , 1)
150        WRITE(msgBuf,*) 'GGL90: GGL90ck              = ', GGL90ck  
151        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL WRITE_0D_R8( GGL90dumpFreq, INDEX_NONE,'GGL90dumpFreq =',
152       &                    SQUEEZE_RIGHT , 1)       &'   /* GGL90 state write out interval ( s ). */')
153        WRITE(msgBuf,*) 'GGL90: GGL90ceps            = ', GGL90ceps        CALL WRITE_0D_R8( GGL90taveFreq, INDEX_NONE,'GGL90taveFreq =',
154        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,       &'   /* GGL90 averaging interval ( s ). */')
155       &                    SQUEEZE_RIGHT , 1)        CALL WRITE_0D_L(GGL90mixingMaps,INDEX_NONE,
156        WRITE(msgBuf,*) 'GGL90: GGL90alpha           = ', GGL90alpha       &     'GGL90mixingMAPS =', '   /* GGL90 IO flag. */')
157        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL WRITE_0D_L(GGL90writeState,INDEX_NONE,
158       &                    SQUEEZE_RIGHT , 1)       &     'GGL90writeState =', '   /* GGL90 IO flag. */')
159        WRITE(msgBuf,*) 'GGL90: GGL90m2              = ', GGL90m2        CALL WRITE_0D_R8( GGL90ck, INDEX_NONE,'GGL90ck =',
160        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,       &'   /* GGL90 viscosity parameter. */')
161       &                    SQUEEZE_RIGHT , 1)        CALL WRITE_0D_R8( GGL90ceps, INDEX_NONE,'GGL90ceps =',
162        WRITE(msgBuf,*) 'GGL90: GGL90TKEmin          = ', GGL90TKEmin       &'   /* GGL90 dissipation parameter. */')
163        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL WRITE_0D_R8( GGL90alpha, INDEX_NONE,'GGL90alpha =',
164       &                    SQUEEZE_RIGHT , 1)       &'   /* GGL90 TKE diffusivity parameter. */')
165        WRITE(msgBuf,*) 'GGL90: GGL90TKEsurfMin      = ', GGL90TKEsurfMin        CALL WRITE_0D_R8( GGL90m2, INDEX_NONE,'GGL90m2 =',
166        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,       &'   /* GGL90 wind stress to vertical stress ratio. */')
167       &                    SQUEEZE_RIGHT , 1)        CALL WRITE_0D_R8( GGL90TKEmin, INDEX_NONE,'GGL90TKEmin =',
168        WRITE(msgBuf,*) 'GGL90: GGL90TKEbottom       = ', GGL90TKEbottom       &'   /* GGL90 minimum kinetic energy ( m^2/s^2 ). */')
169        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL WRITE_0D_R8( GGL90TKEsurfMin, INDEX_NONE,
170       &                    SQUEEZE_RIGHT , 1)       &     'GGL90TKEsurfMin =',
171        WRITE(msgBuf,*) 'GGL90: GGL90viscMax         = ', GGL90viscMax       &'   /* 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,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
186       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
187        WRITE(msgBuf,*) 'GGL90: GGL90diffMax         = ', GGL90diffMax  
188        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        WRITE(msgBuf,'(A)')
189         &'// ======================================================='
190          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
191       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
192        WRITE(msgBuf,*) 'GGL90: GGL90diffTKEh        = ', GGL90diffTKEh        WRITE(msgBuf,'(A)') '// End of GGL90 config. summary'
193        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
194       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
195        WRITE(msgBuf,*) 'GGL90: GGL90mixingLengthMin = ',        WRITE(msgBuf,'(A)')
196       &                                           GGL90mixingLengthMin       &'// ======================================================='
197        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
198       &                    SQUEEZE_RIGHT , 1)       &  SQUEEZE_RIGHT , 1)
199        WRITE(msgBuf,*) 'GGL90: GGL90TKEFile         = ', GGL90TKEFile        WRITE(msgBuf,'(A)') ' '
200        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
201       &                    SQUEEZE_RIGHT , 1)       &  SQUEEZE_RIGHT , 1)
 #endif /* UNSAFE_INTERNAL_WRITE */  
202    
203        _END_MASTER(myThid)        _END_MASTER(myThid)
204    
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22