Skip to content
Permalink
Browse files

fix when all is burning

  • Loading branch information...
Jan Mandel authored and janmandel committed Mar 22, 2011
1 parent 0778c78 commit 5f6709b1638308600b1ba78a75ed72a893182a9f
Showing with 6 additions and 5 deletions.
  1. +5 −4 wrfv2_fire/phys/module_fr_sfire_core.F
  2. +1 −1 wrfv2_fire/test/em_fire/hill/namelist.input
@@ -1056,9 +1056,7 @@ real, intent(in)::lfn00,lfn01,lfn10,lfn11 ! level set function at 4 corners o
DREAL::ps,aps,area,ta,out
DREAL::t00,t01,t10,t11
DREAL,parameter::safe=tiny(aps)
character(len=128)::msg
DREAL::dx,dy ! mesh sizes
!*** local
integer::i,j,k

DREAL,dimension(3)::u
@@ -1078,6 +1076,7 @@ real, intent(in)::lfn00,lfn01,lfn10,lfn11 ! level set function at 4 corners o
DREAL::upper,lower,ah,ch,aa,cc,aupp,cupp,alow,clow
DREAL,dimension(2,2)::mQ
DREAL,dimension(2)::ut
character(len=128)::msg

!calls
intrinsic epsilon
@@ -1120,6 +1119,7 @@ real, intent(in)::lfn00,lfn01,lfn10,lfn11 ! level set function at 4 corners o
!*** case0 Do nothing
if ( lfn00>=0 .and. lfn10>=0 .and. lfn01>=0 .and. lfn11>=0 ) then
out = 1.0 ! fuel_left, no burning
area= 0.
!*** case4 all four coners are burning
else if (lfn00<=0 .and. lfn10<=0 .and. lfn01<=0 .and. lfn11<=0) then
! All burning
@@ -1153,9 +1153,10 @@ real, intent(in)::lfn00,lfn01,lfn10,lfn11 ! level set function at 4 corners o
u(3)=-u(3)/fuel_time_cell
s1=u(1)
s2=u(2)
out=1-exp(u(3))*intexp(s1)*intexp(s2)
out=exp(u(3))*intexp(s1)*intexp(s2)
area=1
if ( out<0 .or. out>1.0 ) then
print *,'case4, out should be between 0 and 1'
call message('WARNING: fuel_left_cell: case all burning: out should be between 0 and 1')
end if
!*** case 1,2,3- other cases
!*** part of cell is burning
@@ -185,7 +185,7 @@
fire_grows_only=1, ! if >0 level set function cannot increase = fire can only grow
fire_viscosity=0.4, ! artificial viscosity in level set method (max 1, needed with fire_upwinding=0)
fire_upwinding=3, ! 0=none, 1=standard, 2=godunov, 3=eno, 4=sethian
fire_fuel_left_method=1, ! for now, use 1 only
fire_fuel_left_method=2, ! for now, use 1 only
fire_lfn_ext_up=1.0, ! 0.=extend level set function at boundary by reflection, 1.=always up
fire_advection=0, ! 0 = cawfe, 1 = use abs speed/slope in spread rate, then project on normal to fireline
/

0 comments on commit 5f6709b

Please sign in to comment.
You can’t perform that action at this time.