From d71244e7f8ddaa6b8c4417a4d6b05072f1ba013c Mon Sep 17 00:00:00 2001 From: Jan Mandel Date: Tue, 5 Apr 2011 17:53:19 -0600 Subject: [PATCH] using fwh from fuel categories not fire_wind_height when fire_wind_log_interp=2 or 3 fwh = fire_wind_height is now used only for fire_wind_log_interp=4 --- wrfv2_fire/phys/module_fr_sfire_atm.F | 14 ++++++++------ 1 file changed, 8 insertions(+), 6 deletions(-) diff --git a/wrfv2_fire/phys/module_fr_sfire_atm.F b/wrfv2_fire/phys/module_fr_sfire_atm.F index 43beda35..e1d5b975 100644 --- a/wrfv2_fire/phys/module_fr_sfire_atm.F +++ b/wrfv2_fire/phys/module_fr_sfire_atm.F @@ -690,9 +690,10 @@ subroutine setup_wind_log_interpolation( & enddo enddo do j=jfts,jfte - do i=ifts,ifte - fwh(i,j)=fire_wind_height - enddo + do i=ifts,ifte + k = int(nfuel_cat(i,j)) + fwh(i,j)=fcwh(k) + enddo enddo case(3) @@ -709,9 +710,10 @@ subroutine setup_wind_log_interpolation( & z0, & ! atm grid arrays in fz0) ! fire grid arrays out do j=jfts,jfte - do i=ifts,ifte - fwh(i,j)=fire_wind_height - enddo + do i=ifts,ifte + k = int(nfuel_cat(i,j)) + fwh(i,j)=fcwh(k) + enddo enddo case(4)