<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7652.24">
<TITLE>Bug fix to WRF-NMM v3.0</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/plain format -->
<BR>
<P><FONT SIZE=2>Dear WRF-NMM community:<BR>
<BR>
I have found a bug in the newly released distribution. In subroutine BUCKETS (dyn_nmm/BUCKETS.F), there are<BR>
calls made to function FPVS0 (phys/module_mp_etanew.F) to calculate saturation vapor pressure. FPVS0 uses<BR>
lookup tables to complete the calculation, and the lookup tables must first be initialized by calling subroutine<BR>
GPVS (also in phys/module_mp_etanew.F). Currently this subroutine is only called if the Eta microphysics are being used; for other microphysics options, FPVS0 returns bogus answers.<BR>
<BR>
This bug can be corrected by modifying dyn_nmm/BUCKETS.F<BR>
<BR>
Near line 77:<BR>
<BR>
! EMK BEGIN BUG FIX<BR>
! USE MODULE_MP_ETANEW,ONLY: C1XPVS,C1XPVS0,C2XPVS,C2XPVS0 &<BR>
! ,FPVS,FPVS0,NX,TBPVS,TBPVS0<BR>
USE MODULE_MP_ETANEW,ONLY: C1XPVS,C1XPVS0,C2XPVS,C2XPVS0 &<BR>
,FPVS,FPVS0,NX,TBPVS,TBPVS0 &<BR>
,GPVS ! Initializes lookup table<BR>
! EMK END BUG FIX<BR>
<BR>
Near line 148:<BR>
<BR>
IF(FIRST_PASS)THEN<BR>
FIRST_PASS=.FALSE.<BR>
!<BR>
EPSILON=R_D/R_V<BR>
ONE_MINUS_EPSILON=1.-EPSILON<BR>
R_FACTOR=1./EPSILON-1.<BR>
CP_FACTOR=CPV/CP-1.<BR>
<BR>
! EMK BEGIN BUG FIX<BR>
call gpvs ! Make sure saturation vapor pressure lookup table<BR>
! is initialized<BR>
! EMK END BUG FIX<BR>
ENDIF<BR>
<BR>
<BR>
Happy computing,<BR>
<BR>
-Eric<BR>
<BR>
<BR>
Eric M. Kemp<BR>
Meteorologist<BR>
Northrop Grumman Information Technology<BR>
Intelligence Group (TASC)<BR>
4801 Stonecroft Boulevard<BR>
Chantilly, VA 20151<BR>
(703) 633-8300 x7078 (lab)<BR>
(703) 633-8300 x8278 (office)<BR>
(703) 449-3400 (fax)<BR>
eric.kemp@ngc.com<BR>
<BR>
</FONT>
</P>
</BODY>
</HTML>