[MITgcm-support] bug in exf_bulkformulae.F ?

[Patrick Heimbach]

Hi Eli,

thanks for the heads-up.
I don't have the references with me right now,
will check tomorrow.

-Patrick

  > -----Original Message-----
  > From: mitgcm-support-bounces at mitgcm.org 
  > [mailto:mitgcm-support-bounces at mitgcm.org] On Behalf Of 
  > eliyahu biton
  > Sent: Sunday, June 26, 2005 6:46 AM
  > To: mitgcm-support at mitgcm.org
  > Subject: [MITgcm-support] bug in exf_bulkformulae.F ?
  > 
  > 
  > 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
  > 
  > 
  > _______________________________________________
  > MITgcm-support mailing list
  > MITgcm-support at mitgcm.org
  > http://mitgcm.org/mailman/listinfo/mitgcm-support
  >