using fwh from fuel categories not fire_wind_height when fire_wind_lo…

…g_interp=2 or 3 fwh = fire_wind_height is now used only for fire_wind_log_interp=4
openwfm · Apr 5, 2011 · d71244e · d71244e
1 parent 03bb9bb
commit d71244e
Showing 1 changed file with 8 additions and 6 deletions.
diff --git 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)