added wind reduction factor from 20ft per Baughman&Rothermel(1980)

openwfm · Sep 16, 2010 · 2195482 · 2195482
1 parent f418d05
commit 2195482
Show file tree

Hide file tree

Showing 3 changed files with 14 additions and 5 deletions.
diff --git a/other/Matlab/vis3d/fire_ros.m b/other/Matlab/vis3d/fire_ros.m
@@ -36,6 +36,7 @@
 rhob     = fuelload/fueldepth;             % ovendry bulk density, lb/ft^3
 c        = 7.47 * exp(-0.133 * savr^0.55); % const in wind coef
 bbb      = 0.02526 * savr^0.54;            % const in wind coef
+c        = c * windrf^bbb;                 % jm: wind reduction from 20ft per Baughman&Rothermel(1980)
 e        = 0.715 * exp( -3.59e-4 * savr);  % const in wind coef
 phiwc    = c * (betafl/betaop)^(-e); 
 rtemp2   = savr^1.5;
@@ -71,4 +72,4 @@
     ros = max(.03333,1.2974 * spdms^1.41);       % spread rate, m/s
 end
 ros=min(ros,6);
-end
+end
diff --git a/wrfv2_fire/phys/module_fr_sfire_phys.F b/wrfv2_fire/phys/module_fr_sfire_phys.F
@@ -115,10 +115,11 @@ module module_fr_sfire_phys
    INTEGER, SAVE      :: no_fuel_cat = 14   ! special category outside of 1:nfuelcats
    CHARACTER (len=80), DIMENSION(mfuelcats ), save :: fuel_name
    INTEGER, DIMENSION( mfuelcats ), save :: ichap
-   REAL   , DIMENSION( mfuelcats ), save :: weight,fgi,fci,fci_d,fct,fcbr, &
+   REAL   , DIMENSION( mfuelcats ), save :: windrf,weight,fgi,fci,fci_d,fct,fcbr, &
                                             fueldepthm,fueldens,fuelmce,   &
                                             savr,st,se
-
+   DATA windrf /0.36, 0.36, 0.44,  0.55,  0.42,  0.44,  0.44, &     
+                0.36, 0.36, 0.36,  0.36,  0.43,  0.46,  1e-7, zf*0 /
    DATA fgi / 0.166, 0.897, 0.675, 2.468, 0.785, 1.345, 1.092, &
               1.121, 0.780, 2.694, 2.582, 7.749, 13.024, 1.e-7, zf*0.  /
    DATA fueldepthm /0.305,  0.305,  0.762, 1.829, 0.61,  0.762,0.762, &
@@ -158,7 +159,7 @@ subroutine init_fuel_cats
 
 ! read 
 namelist /fuel_scalars/ cmbcnst,hfgl,fuelmc_g,fuelmc_c,nfuelcats,no_fuel_cat
-namelist /fuel_categories/ fuel_name,fgi,fueldepthm,savr, &
+namelist /fuel_categories/ fuel_name,windrf,fgi,fueldepthm,savr, &
     fuelmce,fueldens,st,se,weight,fci_d,fct,ichap
 
 !$  if (OMP_GET_THREAD_NUM() .ne. 0)then
@@ -189,6 +190,7 @@ subroutine init_fuel_cats
 call wrf_dm_bcast_real(fuelmc_c,1)
 call wrf_dm_bcast_integer(nfuelcats,1)
 call wrf_dm_bcast_integer(no_fuel_cat,1)
+call wrf_dm_bcast_real(windrf,    nfuelcats)
 call wrf_dm_bcast_real(fgi,       nfuelcats)
 call wrf_dm_bcast_real(fueldepthm,nfuelcats)
 call wrf_dm_bcast_real(savr,      nfuelcats)
@@ -249,6 +251,8 @@ subroutine init_fuel_cats
     call message(msg)
     write(msg,9)'fuel name ',(fuel_name(ii),ii=i,k)
     call message(msg)
+    write(msg,8)'windrf    ',(windrf(ii),ii=i,k)
+    call message(msg)
     write(msg,8)'fgi       ',(fgi(ii),ii=i,k)
     call message(msg)
     write(msg,8)'fueldepthm',(fueldepthm(ii),ii=i,k)
@@ -319,6 +323,7 @@ subroutine write_fuels_m(nsteps,maxwind,maxslope)
 10 format('fuel(',i3,').',a,'=',"'",a,"'",';% ',a)
 do k=1,nfuelcats
     write(iounit,10)k,'fuel_name',trim(fuel_name(k)),'FUEL MODEL NAME'
+    call write_var(k,'windrf',windrf(k),'WIND REDUCTION FACTOR FROM 20ft TO MIDFLAME HEIGHT' )
     call write_var(k,'fgi',fgi(k),'INITIAL TOTAL MASS OF SURFACE FUEL (KG/M**2)' )
     call write_var(k,'fueldepthm',fueldepthm(k),'FUEL DEPTH (M)')
     call write_var(k,'savr',savr(k),'FUEL PARTICLE SURFACE-AREA-TO-VOLUME RATIO, 1/FT')
@@ -496,7 +501,8 @@ subroutine set_fire_params( &
         rhob = fuelload/fueldepth    ! ovendry bulk density, lb/ft^3
 
         c = 7.47 * exp( -0.133 * savr(k)**0.55)    ! const in wind coef
-        fp%bbb(i,j) = 0.02526 * savr(k)**0.54                ! const in wind coef
+        fp%bbb(i,j) = 0.02526 * savr(k)**0.54      ! const in wind coef
+        c = c * windrf(k)**fp%bbb(i,j)             ! jm: wind reduction from 20ft per Baughman&Rothermel(1980)
         e = 0.715 * exp( -3.59e-4 * savr(k))       ! const in wind coef
         fp%phiwc(i,j) = c * (fp%betafl(i,j)/betaop)**(-e)
 

diff --git a/wrfv2_fire/test/em_fire/fireflux_small/namelist.fire b/wrfv2_fire/test/em_fire/fireflux_small/namelist.fire
@@ -24,6 +24,8 @@ no_fuel_cat = 14                   ! extra category for no fuel
 '12: Medium logging slash',
 '13: Heavy logging slash',
 '14: no fuel'
+ windrf= 0.36, 0.36, 0.44,  0.55,  0.42,  0.44,  0.44, 
+         0.36, 0.36, 0.36,  0.36,  0.43,  0.46, 1e-7
  fgi =  0.166, 0.897, 1.08, 2.468, 0.785, 1.345, 1.092, 
 !jc fgi=  0.166, 0.897, 1.076, 2.468, 0.785, 1.345, 1.092
         1.121, 0.780, 2.694, 2.582, 7.749, 13.024, 1.e-7,