[MITgcm-support] bug in exf_bulkformulae.F ?

[eliyahu biton]

Hello everyone!
I tried to follow the calculation of the turbulent heat fluxes, as they 
appear within the subroutine 'exf_bulkformulae.F' at the 'exf' pkg and I 
think that there is a problem with the calculation of the term 'psimh' 
(line 394).
            psimh  = -psim_fac*huol*stable +
    &                     (exf_one - stable)*
    &                     log((exf_one + x*(exf_two + x))*
    &                     (exf_one + xsq)/8.) - exf_two*atan(x) +
    &                     pi*exf_half
According to Large & Pond paper from 1981 (eq 7) it should be (the last 
two arguments should be included within the unstable case) ;
            psimh  = -psim_fac*huol*stable +
    &                     (exf_one - stable)*
    &                     (log((exf_one + x*(exf_two + x))*
    &                     (exf_one + xsq)/8.) - exf_two*atan(x) +
    &                     pi*exf_half)
I checked how this term calculated at subroutine bulkf_formula_lanl.F 
(line  159 ) in 'bulk_force' pkg,  and there it  looks o.k. ( the last 
two argument are included within the unstable case).
             psimh = -5. _d 0*huol*stable + (1. _d 0-stable)*
    &                    (2. _d 0*log(5. _d -1*(1. _d 0+x)) +
    &                     2. _d 0*log(5. _d -1*(1. _d 0+xsq)) -
    &                     2. _d 0*atan(x) + pi*.5 _d 0)


Eli

-- 
Eli Biton
Dept. of Environmental Sciences,
Weizmann Institute of Sciences,
Rehovot, 76100
Israel
Tel: +972-8-9344924
Fax: +972-8-9344124