324 |
& + surfaceForcingV(I, J+1,bi,bj) ) )**2 |
& + surfaceForcingV(I, J+1,bi,bj) ) )**2 |
325 |
& ) |
& ) |
326 |
C Dirichlet surface boundary condition for TKE |
C Dirichlet surface boundary condition for TKE |
327 |
gTKE(I,J,kSurf) = MAX(GGL90TKEmin,GGL90m2*uStarSquare) |
gTKE(I,J,kSurf) = MAX(GGL90TKEsurfMin,GGL90m2*uStarSquare) |
328 |
& *maskC(I,J,kSurf,bi,bj) |
& *maskC(I,J,kSurf,bi,bj) |
329 |
C Dirichlet bottom boundary condition for TKE = GGL90TKEbottom |
C Dirichlet bottom boundary condition for TKE = GGL90TKEbottom |
330 |
kBottom = MIN(MAX(kLowC(I,J,bi,bj),1),Nr) |
kBottom = MIN(MAX(kLowC(I,J,bi,bj),1),Nr) |