/[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.4 by mlosch, Tue Jun 6 19:47:36 2006 UTC revision 1.7 by jmc, Thu Jun 8 14:15:26 2006 UTC
# Line 28  C     errIO       - IO error flag Line 28  C     errIO       - IO error flag
28  C     iUnit       - Work variable for IO unit number  C     iUnit       - Work variable for IO unit number
29    
30        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
31        INTEGER errIO, iUnit        INTEGER errIO, iUnit, iL
32    
33          INTEGER  ILNBLNK
34          EXTERNAL ILNBLNK
35    
36  C--   GGL90 vertical mixing parameters  C--   GGL90 vertical mixing parameters
37        NAMELIST /GGL90_PARM01/        NAMELIST /GGL90_PARM01/
# Line 136  C Now set-up any remaining parameters th Line 139  C Now set-up any remaining parameters th
139        ENDIF        ENDIF
140    
141  C--   print TKE vertical mixing parameters to stdout for better debugging  C--   print TKE vertical mixing parameters to stdout for better debugging
142        WRITE(msgBuf,*) 'GGL90: GGL90dumpFreq        = ', GGL90dumpFreq        WRITE(msgBuf,'(A)')
143        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,       &'// ======================================================='
144       &                    SQUEEZE_RIGHT , 1)        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
       WRITE(msgBuf,*) 'GGL90: GGL90taveFreq        = ', GGL90taveFreq  
       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
      &                    SQUEEZE_RIGHT , 1)  
       WRITE(msgBuf,*) 'GGL90: GGL90mixingMaps      = ', GGL90mixingMaps  
       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
      &                    SQUEEZE_RIGHT , 1)  
       WRITE(msgBuf,*) 'GGL90: GGL90writeState      = ', GGL90writeState  
       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
      &                    SQUEEZE_RIGHT , 1)  
       WRITE(msgBuf,*) 'GGL90: GGL90ck              = ', GGL90ck  
       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
      &                    SQUEEZE_RIGHT , 1)  
       WRITE(msgBuf,*) 'GGL90: GGL90ceps            = ', GGL90ceps  
       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
      &                    SQUEEZE_RIGHT , 1)  
       WRITE(msgBuf,*) 'GGL90: GGL90alpha           = ', GGL90alpha  
       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
      &                    SQUEEZE_RIGHT , 1)  
       WRITE(msgBuf,*) 'GGL90: GGL90m2              = ', GGL90m2  
       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
      &                    SQUEEZE_RIGHT , 1)  
       WRITE(msgBuf,*) 'GGL90: GGL90TKEmin          = ', GGL90TKEmin  
       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
145       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
146        WRITE(msgBuf,*) 'GGL90: GGL90TKEsurfMin      = ', GGL90TKEsurfMin        WRITE(msgBuf,'(A)') '// GGL90 configuration'
147        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
148       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
149        WRITE(msgBuf,*) 'GGL90: GGL90TKEbottom       = ', GGL90TKEbottom        WRITE(msgBuf,'(A)')
150         &'// ======================================================='
151        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
152       &                    SQUEEZE_RIGHT , 1)       &  SQUEEZE_RIGHT , 1)
153        WRITE(msgBuf,*) 'GGL90: GGL90viscMax         = ', GGL90viscMax  
154          CALL WRITE_0D_R8( GGL90dumpFreq, INDEX_NONE,'GGL90dumpFreq =',
155         &'   /* GGL90 state write out interval ( s ). */')
156          CALL WRITE_0D_R8( GGL90taveFreq, INDEX_NONE,'GGL90taveFreq =',
157         &'   /* GGL90 averaging interval ( s ). */')
158          CALL WRITE_0D_L(GGL90mixingMaps,INDEX_NONE,
159         &     'GGL90mixingMAPS =', '   /* GGL90 IO flag. */')
160          CALL WRITE_0D_L(GGL90writeState,INDEX_NONE,
161         &     'GGL90writeState =', '   /* GGL90 IO flag. */')
162          CALL WRITE_0D_R8( GGL90ck, INDEX_NONE,'GGL90ck =',
163         &'   /* GGL90 viscosity parameter. */')
164          CALL WRITE_0D_R8( GGL90ceps, INDEX_NONE,'GGL90ceps =',
165         &'   /* GGL90 dissipation parameter. */')
166          CALL WRITE_0D_R8( GGL90alpha, INDEX_NONE,'GGL90alpha =',
167         &'   /* GGL90 TKE diffusivity parameter. */')
168          CALL WRITE_0D_R8( GGL90m2, INDEX_NONE,'GGL90m2 =',
169         &'   /* GGL90 wind stress to vertical stress ratio. */')
170          CALL WRITE_0D_R8( GGL90TKEmin, INDEX_NONE,'GGL90TKEmin =',
171         &'   /* GGL90 minimum kinetic energy ( m^2/s^2 ). */')
172          CALL WRITE_0D_R8( GGL90TKEsurfMin, INDEX_NONE,
173         &     'GGL90TKEsurfMin =',
174         &'   /* GGL90 minimum surface kinetic energy ( m^2/s^2 ). */')
175          CALL WRITE_0D_R8( GGL90TKEbottom, INDEX_NONE,
176         &     'GGL90TKEbottom =',
177         &     '   /* GGL90 bottom kinetic energy ( m^2/s^2 ). */')
178          CALL WRITE_0D_R8( GGL90viscMax, INDEX_NONE,'GGL90viscMax =',
179         &     '   /* GGL90 upper limit for viscosity ( m^2/s ). */')
180          CALL WRITE_0D_R8( GGL90diffMax, INDEX_NONE,'GGL90diffMax =',
181         &     '   /* GGL90 upper limit for diffusivity ( m^2/s ). */')
182          CALL WRITE_0D_R8( GGL90diffTKEh, INDEX_NONE,'GGL90diffTKEh =',
183         &     '   /* GGL90 horizontal diffusivity for TKE ( m^2/s ). */')
184          CALL WRITE_0D_R8( GGL90mixingLengthMin, INDEX_NONE,
185         &     'GGL90mixingLengthMin =',
186         &     '   /* GGL90 minimum mixing length ( m ). */')
187          iL = MAX_LEN_MBUF - 22
188          iL = MIN(iL,ILNBLNK(GGL90TKEFile))
189          WRITE(msgBuf,'(A,A)')'GGL90: GGL90TKEFile = ',GGL90TKEFile(1:iL)
190        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
191       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
192        WRITE(msgBuf,*) 'GGL90: GGL90diffMax         = ', GGL90diffMax  
193        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        WRITE(msgBuf,'(A)')
194         &'// ======================================================='
195          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
196       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
197        WRITE(msgBuf,*) 'GGL90: GGL90diffTKEh        = ', GGL90diffTKEh        WRITE(msgBuf,'(A)') '// End of GGL90 config. summary'
198        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
199       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
200        WRITE(msgBuf,*) 'GGL90: GGL90mixingLengthMin = ',        WRITE(msgBuf,'(A)')
201       &                                           GGL90mixingLengthMin       &'// ======================================================='
202        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
203       &                    SQUEEZE_RIGHT , 1)       &  SQUEEZE_RIGHT , 1)
204        WRITE(msgBuf,*) 'GGL90: GGL90TKEFile         = ', GGL90TKEFile        WRITE(msgBuf,'(A)') ' '
205        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
206       &                    SQUEEZE_RIGHT , 1)       &  SQUEEZE_RIGHT , 1)
207    
208        _END_MASTER(myThid)        _END_MASTER(myThid)
209    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22