1 |
C $Header$ |
C $Header$ |
2 |
C $Name$ |
C $Name$ |
3 |
|
|
4 |
#include "GGL90_OPTIONS.h" |
#include "GGL90_OPTIONS.h" |
5 |
|
|
6 |
SUBROUTINE GGL90_READPARMS( myThid ) |
SUBROUTINE GGL90_READPARMS( myThid ) |
41 |
& GGL90mixingMaps, GGL90writeState, |
& GGL90mixingMaps, GGL90writeState, |
42 |
& GGL90ck, GGL90ceps, GGL90alpha, GGL90m2, |
& GGL90ck, GGL90ceps, GGL90alpha, GGL90m2, |
43 |
& GGL90TKEmin, GGL90TKEsurfMin, GGL90TKEbottom, |
& GGL90TKEmin, GGL90TKEsurfMin, GGL90TKEbottom, |
44 |
& GGL90mixingLengthMin, GGL90viscMax, GGL90diffMax, |
& GGL90mixingLengthMin, mxlMaxFlag, |
45 |
& GGL90TKEFile |
& GGL90viscMax, GGL90diffMax, GGL90TKEFile |
46 |
_BEGIN_MASTER(myThid) |
_BEGIN_MASTER(myThid) |
47 |
|
|
48 |
WRITE(msgBuf,'(A)') ' GGL90_READPARMS: opening data.ggl90' |
WRITE(msgBuf,'(A)') ' GGL90_READPARMS: opening data.ggl90' |
69 |
C Blanke and Delecluse (1993, JPO) use |
C Blanke and Delecluse (1993, JPO) use |
70 |
GGL90TKEsurfMin = 1.0 _d -04 |
GGL90TKEsurfMin = 1.0 _d -04 |
71 |
GGL90TKEbottom = UNSET_RL |
GGL90TKEbottom = UNSET_RL |
72 |
GGL90viscMax = 1.0 _d 0 |
GGL90viscMax = 1. _d 2 |
73 |
GGL90diffMax = 1.0 _d 1 |
GGL90diffMax = 1. _d 2 |
74 |
GGL90diffTKEh = 0.0 _d 0 |
GGL90diffTKEh = 0.0 _d 0 |
75 |
GGL90mixingLengthMin = 1.0 _d -08 |
GGL90mixingLengthMin = 1.0 _d -08 |
76 |
|
mxlMaxFlag = 0 |
77 |
GGL90TKEFile = ' ' |
GGL90TKEFile = ' ' |
78 |
|
|
79 |
C----------------------------------------------------------------------- |
C----------------------------------------------------------------------- |
186 |
CALL WRITE_0D_R8( GGL90mixingLengthMin, INDEX_NONE, |
CALL WRITE_0D_R8( GGL90mixingLengthMin, INDEX_NONE, |
187 |
& 'GGL90mixingLengthMin =', |
& 'GGL90mixingLengthMin =', |
188 |
& ' /* GGL90 minimum mixing length ( m ). */') |
& ' /* GGL90 minimum mixing length ( m ). */') |
189 |
|
CALL WRITE_0D_I(mxlMaxFlag, INDEX_NONE, 'mxlMaxFlag =', |
190 |
|
& ' /* Flag for limiting mixing-length method */') |
191 |
iL = MAX_LEN_MBUF - 22 |
iL = MAX_LEN_MBUF - 22 |
192 |
iL = MIN( iL, MAX(ILNBLNK(GGL90TKEFile),1) ) |
iL = MIN( iL, MAX(ILNBLNK(GGL90TKEFile),1) ) |
193 |
WRITE(msgBuf,'(A,A)')'GGL90: GGL90TKEFile = ',GGL90TKEFile(1:iL) |
WRITE(msgBuf,'(A,A)')'GGL90: GGL90TKEFile = ',GGL90TKEFile(1:iL) |