--- MITgcm/pkg/ggl90/ggl90_calc.F 2005/08/05 23:46:09 1.5 +++ MITgcm/pkg/ggl90/ggl90_calc.F 2006/06/06 16:18:18 1.6 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/ggl90/ggl90_calc.F,v 1.5 2005/08/05 23:46:09 ce107 Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/ggl90/ggl90_calc.F,v 1.6 2006/06/06 16:18:18 mlosch Exp $ C $Name: $ #include "GGL90_OPTIONS.h" @@ -324,7 +324,7 @@ & + surfaceForcingV(I, J+1,bi,bj) ) )**2 & ) C Dirichlet surface boundary condition for TKE - gTKE(I,J,kSurf) = MAX(GGL90TKEmin,GGL90m2*uStarSquare) + gTKE(I,J,kSurf) = MAX(GGL90TKEsurfMin,GGL90m2*uStarSquare) & *maskC(I,J,kSurf,bi,bj) C Dirichlet bottom boundary condition for TKE = GGL90TKEbottom kBottom = MIN(MAX(kLowC(I,J,bi,bj),1),Nr)