64 |
CHARACTER*(MAX_LEN_MBUF) msgBuf |
CHARACTER*(MAX_LEN_MBUF) msgBuf |
65 |
CEOP |
CEOP |
66 |
|
|
|
KbryanLewis79=diffKrBL79surf+(diffKrBL79deep-diffKrBL79surf) |
|
|
& *( atan( -( rF(k)-diffKrBL79Ho )/diffKrBL79scl )/PI+0.5 _d 0) |
|
|
|
|
67 |
IF ( trIdentity.EQ.GAD_TEMPERATURE) THEN |
IF ( trIdentity.EQ.GAD_TEMPERATURE) THEN |
68 |
|
|
69 |
DO k = 1,Nr |
DO k = 1,Nr |
70 |
|
KbryanLewis79=diffKrBL79surf+(diffKrBL79deep-diffKrBL79surf) |
71 |
|
& *( atan( -( rF(k)-diffKrBL79Ho )/diffKrBL79scl )/PI+0.5 _d 0) |
72 |
DO j = 1-Oly,sNy+Oly |
DO j = 1-Oly,sNy+Oly |
73 |
DO i = 1-Olx,sNx+Olx |
DO i = 1-Olx,sNx+Olx |
74 |
KappaRTr(i,j,k) = |
KappaRTr(i,j,k) = |
86 |
ELSEIF ( trIdentity.EQ.GAD_SALINITY) THEN |
ELSEIF ( trIdentity.EQ.GAD_SALINITY) THEN |
87 |
|
|
88 |
DO k = 1,Nr |
DO k = 1,Nr |
89 |
|
KbryanLewis79=diffKrBL79surf+(diffKrBL79deep-diffKrBL79surf) |
90 |
|
& *( atan( -( rF(k)-diffKrBL79Ho )/diffKrBL79scl )/PI+0.5 _d 0) |
91 |
DO j = 1-Oly, sNy+Oly |
DO j = 1-Oly, sNy+Oly |
92 |
DO i = 1-Olx, sNx+Olx |
DO i = 1-Olx, sNx+Olx |
93 |
KappaRTr(i,j,k) = |
KappaRTr(i,j,k) = |
108 |
|
|
109 |
iTr = trIdentity - GAD_TR1 + 1 |
iTr = trIdentity - GAD_TR1 + 1 |
110 |
DO k = 1,Nr |
DO k = 1,Nr |
111 |
|
KbryanLewis79=diffKrBL79surf+(diffKrBL79deep-diffKrBL79surf) |
112 |
|
& *( atan( -( rF(k)-diffKrBL79Ho )/diffKrBL79scl )/PI+0.5 _d 0) |
113 |
DO j = 1-Oly, sNy+Oly |
DO j = 1-Oly, sNy+Oly |
114 |
DO i = 1-Olx, sNx+Olx |
DO i = 1-Olx, sNx+Olx |
115 |
KappaRTr(i,j,k) = |
KappaRTr(i,j,k) = |
147 |
IF (trUseKPP) THEN |
IF (trUseKPP) THEN |
148 |
IF (trIdentity.EQ.GAD_TEMPERATURE) THEN |
IF (trIdentity.EQ.GAD_TEMPERATURE) THEN |
149 |
CALL KPP_CALC_DIFF_T( |
CALL KPP_CALC_DIFF_T( |
150 |
I bi,bj,iMin+1,iMax,jMin+1,jMax,0,Nr, |
I bi,bj,iMin,iMax,jMin,jMax,0,Nr, |
151 |
U KappaRTr, |
U KappaRTr, |
152 |
I myThid) |
I myThid) |
153 |
ELSE |
ELSE |
154 |
CALL KPP_CALC_DIFF_S( |
CALL KPP_CALC_DIFF_S( |
155 |
I bi,bj,iMin+1,iMax,jMin+1,jMax,0,Nr, |
I bi,bj,iMin,iMax,jMin,jMax,0,Nr, |
156 |
U KappaRTr, |
U KappaRTr, |
157 |
I myThid) |
I myThid) |
158 |
ENDIF |
ENDIF |
162 |
#ifdef ALLOW_PP81 |
#ifdef ALLOW_PP81 |
163 |
IF (usePP81) THEN |
IF (usePP81) THEN |
164 |
CALL PP81_CALC_DIFF( |
CALL PP81_CALC_DIFF( |
165 |
I bi,bj,iMin+1,iMax,jMin+1,jMax,0,Nr, |
I bi,bj,iMin,iMax,jMin,jMax,0,Nr, |
166 |
U KappaRTr, |
U KappaRTr, |
167 |
I myThid) |
I myThid) |
168 |
ENDIF |
ENDIF |
171 |
#ifdef ALLOW_MY82 |
#ifdef ALLOW_MY82 |
172 |
IF (useMY82) THEN |
IF (useMY82) THEN |
173 |
CALL MY82_CALC_DIFF( |
CALL MY82_CALC_DIFF( |
174 |
I bi,bj,iMin+1,iMax,jMin+1,jMax,0,Nr, |
I bi,bj,iMin,iMax,jMin,jMax,0,Nr, |
175 |
U KappaRTr, |
U KappaRTr, |
176 |
I myThid) |
I myThid) |
177 |
ENDIF |
ENDIF |
180 |
#ifdef ALLOW_GGL90 |
#ifdef ALLOW_GGL90 |
181 |
IF (useGGL90) THEN |
IF (useGGL90) THEN |
182 |
CALL GGL90_CALC_DIFF( |
CALL GGL90_CALC_DIFF( |
183 |
I bi,bj,iMin+1,iMax,jMin+1,jMax,0,Nr, |
I bi,bj,iMin,iMax,jMin,jMax,0,Nr, |
184 |
O KappaRTr, |
O KappaRTr, |
185 |
I myThid) |
I myThid) |
186 |
ENDIF |
ENDIF |