[MITgcm-devel] KPP bug for initialisation of minKPPhbl

[Patrick Heimbach]

There's a bug in the initialisation of
minimum KPP boundary layer depth.
minKPPhbl is set in kpp_readparms.F to be
      minKPPhbl     = drF(1)*0.5

but drF(K) is only set afterwards in S/R ini_vertical_grid
Thus, minKPPhbl is always zero, leading to a too slow
establishment of a boundary layer.
I am going to move it to kpp_init.
(It may be somewhat debatable anyways whether a 
drF(1)-dependent minimal depth is at all meaningful.)

-- 
----------------------------------------------------------
Patrick Heimbach     Massachusetts Institute of Technology
FON: +1/617/253-5259                    EAPS, Room 54-1518
FAX: +1/617/253-4464               77 Massachusetts Avenue
mailto:heimbach at mit.edu                 Cambridge MA 02139
http://www.mit.edu/~heimbach/                          USA