subroutine gfssflx::frh2o | ( | ) |
New version (June 2001): much faster and more accurate Newton iteration achieved by first taking log of eqn cited above – less than 4 (typically 1 or 2) iterations achieves convergence. Also, explicit 1-step solution option for special case of paramter ck=0, which reduces the orginal implicit equation to a simpler explicit form, known as the "flerchinger eqn". Improved handling of solution in the limit of freezing point temperature t0.
Referenced by snksrc().