/[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.1 by mlosch, Thu Sep 16 11:27:18 2004 UTC revision 1.5 by jmc, Wed Jun 7 16:34:24 2006 UTC
# Line 33  C     iUnit       - Work variable for IO Line 33  C     iUnit       - Work variable for IO
33  C--   GGL90 vertical mixing parameters  C--   GGL90 vertical mixing parameters
34        NAMELIST /GGL90_PARM01/        NAMELIST /GGL90_PARM01/
35       &     GGL90dumpFreq, GGL90taveFreq,       &     GGL90dumpFreq, GGL90taveFreq,
36         &     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 58  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
70          GGL90diffTKEh        = 0.0 _d 0
71        GGL90mixingLengthMin = 1.0 _d -08        GGL90mixingLengthMin = 1.0 _d -08
72        GGL90TKEFile         = ' '        GGL90TKEFile         = ' '
73    
# Line 97  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 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    
138    #ifdef UNSAFE_INTERNAL_WRITE
139    C--   print TKE vertical mixing parameters to stdout for better debugging
140          WRITE(msgBuf,*) 'GGL90: GGL90dumpFreq        = ', GGL90dumpFreq
141          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
142         &                    SQUEEZE_RIGHT , 1)
143          WRITE(msgBuf,*) 'GGL90: GGL90taveFreq        = ', GGL90taveFreq
144          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
145         &                    SQUEEZE_RIGHT , 1)
146          WRITE(msgBuf,*) 'GGL90: GGL90mixingMaps      = ', GGL90mixingMaps
147          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
148         &                    SQUEEZE_RIGHT , 1)
149          WRITE(msgBuf,*) 'GGL90: GGL90writeState      = ', GGL90writeState
150          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
151         &                    SQUEEZE_RIGHT , 1)
152          WRITE(msgBuf,*) 'GGL90: GGL90ck              = ', GGL90ck
153          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
154         &                    SQUEEZE_RIGHT , 1)
155          WRITE(msgBuf,*) 'GGL90: GGL90ceps            = ', GGL90ceps
156          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
157         &                    SQUEEZE_RIGHT , 1)
158          WRITE(msgBuf,*) 'GGL90: GGL90alpha           = ', GGL90alpha
159          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
160         &                    SQUEEZE_RIGHT , 1)
161          WRITE(msgBuf,*) 'GGL90: GGL90m2              = ', GGL90m2
162          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
163         &                    SQUEEZE_RIGHT , 1)
164          WRITE(msgBuf,*) 'GGL90: GGL90TKEmin          = ', GGL90TKEmin
165          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
166         &                    SQUEEZE_RIGHT , 1)
167          WRITE(msgBuf,*) 'GGL90: GGL90TKEsurfMin      = ', GGL90TKEsurfMin
168          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
169         &                    SQUEEZE_RIGHT , 1)
170          WRITE(msgBuf,*) 'GGL90: GGL90TKEbottom       = ', GGL90TKEbottom
171          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
172         &                    SQUEEZE_RIGHT , 1)
173          WRITE(msgBuf,*) 'GGL90: GGL90viscMax         = ', GGL90viscMax
174          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
175         &                    SQUEEZE_RIGHT , 1)
176          WRITE(msgBuf,*) 'GGL90: GGL90diffMax         = ', GGL90diffMax
177          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
178         &                    SQUEEZE_RIGHT , 1)
179          WRITE(msgBuf,*) 'GGL90: GGL90diffTKEh        = ', GGL90diffTKEh
180          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
181         &                    SQUEEZE_RIGHT , 1)
182          WRITE(msgBuf,*) 'GGL90: GGL90mixingLengthMin = ',
183         &                                           GGL90mixingLengthMin
184          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
185         &                    SQUEEZE_RIGHT , 1)
186          WRITE(msgBuf,*) 'GGL90: GGL90TKEFile         = ', GGL90TKEFile
187          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
188         &                    SQUEEZE_RIGHT , 1)
189    #endif /* UNSAFE_INTERNAL_WRITE */
190    
191          _END_MASTER(myThid)
192    
193    C--   Everyone else must wait for the parameters to be loaded
194          _BARRIER
195    
196  #endif /* ALLOW_GGL90 */  #endif /* ALLOW_GGL90 */
197    
198        return        return
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22