97 |
_RL dTdx(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL dTdx(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
98 |
_RL dTdy(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL dTdy(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
99 |
#ifdef ALLOW_KPP |
#ifdef ALLOW_KPP |
100 |
_RL hbl (1-OLx:sNx+OLx,1-OLy:sNy+OLy) ! used by KPP mixing scheme |
_RS hbl (1-OLx:sNx+OLx,1-OLy:sNy+OLy) ! used by KPP mixing scheme |
101 |
_RL frac (1-OLx:sNx+OLx,1-OLy:sNy+OLy) ! used by KPP mixing scheme |
_RS frac (1-OLx:sNx+OLx,1-OLy:sNy+OLy) ! used by KPP mixing scheme |
102 |
|
_RS negone ! used as argument to SWFRAC |
103 |
integer jwtype ! index for Jerlov water type |
integer jwtype ! index for Jerlov water type |
104 |
#endif |
#endif |
105 |
|
|
266 |
ENDDO |
ENDDO |
267 |
j=(sNx+2*OLx)*(sNy+2*OLy) |
j=(sNx+2*OLx)*(sNy+2*OLy) |
268 |
jwtype = 3 |
jwtype = 3 |
269 |
|
negone = -1. |
270 |
CALL SWFRAC( |
CALL SWFRAC( |
271 |
I j, -1., hbl, jwtype, |
I j, negone, hbl, jwtype, |
272 |
O frac ) |
O frac ) |
273 |
|
|
274 |
C Add non local transport coefficient (ghat term) to right-hand-side |
C Add non local transport coefficient (ghat term) to right-hand-side |