Skip to content

Commit

Permalink
Mods to the fire routines before the release. Files
Browse files Browse the repository at this point in the history 
from Janice Coen and Jonathan Beezley.

M    phys/module_fr_sfire_phys.F
M    phys/module_fr_sfire_atm.F
M    phys/module_fr_sfire_core.F
M    phys/module_fr_sfire_util.F
M    phys/module_fr_sfire_driver.F
M    phys/module_fr_sfire_model.F



git-svn-id: https://svn-wrf-model.cgd.ucar.edu/trunk@4752 b0b5d27b-6f0f-0410-a2a3-cb1e977edc3d
  • Loading branch information
@davegill
davegill committed Mar 4, 2011
1 parent 21f723f commit 54a1084
Show file tree
Hide file tree
Showing 6 changed files with 95 additions and 198 deletions.
20 changes: 9 additions & 11 deletions phys/module_fr_sfire_atm.F
View file
@@ -1,8 +1,4 @@
!WRF:MEDIATION_LAYER:FIRE_MODEL

!*** Jan Mandel August 2007 - February 2008
!*** email: jmandel@ucar.edu or Jan.Mandel@gmail.com or Jan.Mandel@cudenver.edu

! Routines dealing with the atmosphere

module module_fr_sfire_atm
Expand All @@ -22,9 +18,8 @@ SUBROUTINE fire_tendency( &
rthfrten,rqvfrten) ! theta and Qv tendencies

! This routine is atmospheric physics
! it does NOT go into module_fr_sfire_phys because it is not fire physics
! it does NOT go into module_fr_sfire_phys because it is not related to fire physical processes

! taken from the code by Ned Patton, only order of arguments change to the convention here
! --- this routine takes fire generated heat and moisture fluxes and
! calculates their influence on the theta and water vapor
! --- note that these tendencies are valid at the Arakawa-A location
Expand All @@ -46,7 +41,7 @@ SUBROUTINE fire_tendency( &
REAL, INTENT(in), DIMENSION( ims:ime,kms:kme,jms:jme ) :: dz8w ! dz across w-lvl
REAL, INTENT(in), DIMENSION( ims:ime,kms:kme,jms:jme ) :: rho ! density

REAL, INTENT(in) :: alfg ! extinction depth ground fire heat (m)
REAL, INTENT(in) :: alfg ! extinction depth surface fire heat (m)
REAL, INTENT(in) :: alfc ! extinction depth crown fire heat (m)
REAL, INTENT(in) :: z1can ! height of crown fire heat release (m)

Expand All @@ -65,6 +60,7 @@ SUBROUTINE fire_tendency( &
REAL :: xlv_i
REAL :: z_w
REAL :: fact_g, fact_c
REAL :: alfg_i, alfc_i

REAL, DIMENSION( its:ite,kts:kte,jts:jte ) :: hfx,qfx

Expand All @@ -85,6 +81,8 @@ SUBROUTINE fire_tendency( &

cp_i = 1./cp ! inverse of specific heat
xlv_i = 1./xlv ! inverse of latent heat
alfg_i = 1./alfg
alfc_i = 1./alfc

!!write(msg,'(8e11.3)')cp,cp_i,xlv,xlv_i,alfg,alfc,z1can
!!call message(msg)
Expand Down Expand Up @@ -113,11 +111,11 @@ SUBROUTINE fire_tendency( &

! --- heat flux

fact_g = cp_i * EXP( - alfg * z_w )
fact_g = cp_i * EXP( - alfg_i * z_w )
IF ( z_w < z1can ) THEN
fact_c = cp_i
ELSE
fact_c = cp_i * EXP( - alfc * (z_w - z1can) )
fact_c = cp_i * EXP( - alfc_i * (z_w - z1can) )
END IF
hfx(i,k,j) = fact_g * grnhfx(i,j) + fact_c * canhfx(i,j)

Expand All @@ -127,11 +125,11 @@ SUBROUTINE fire_tendency( &

! --- vapor flux

fact_g = xlv_i * EXP( - alfg * z_w )
fact_g = xlv_i * EXP( - alfg_i * z_w )
IF (z_w < z1can) THEN
fact_c = xlv_i
ELSE
fact_c = xlv_i * EXP( - alfc * (z_w - z1can) )
fact_c = xlv_i * EXP( - alfc_i * (z_w - z1can) )
END IF
qfx(i,k,j) = fact_g * grnqfx(i,j) + fact_c * canqfx(i,j)

Expand Down
22 changes: 8 additions & 14 deletions phys/module_fr_sfire_core.F
View file
@@ -1,7 +1,4 @@
!
!*** Jan Mandel August-October 2007 email: jmandel@ucar.edu or Jan.Mandel@gmail.com
!
! With contributions by Minjeong Kim.
#define DEBUG_OUT
#define DEBUG_PRINT
!#define FUEL_LEFT
Expand Down Expand Up @@ -355,8 +352,6 @@ subroutine fuel_left(&
!*** purpose: determine fraction of fuel remaining
!*** NOTE: because variables are cell centered, need halo/sync width 1 before

!*** Jan Mandel August 2007 email: jmandel@ucar.edu or Jan.Mandel@gmail.com

!*** arguments

integer, intent(in) :: its,ite,jts,jte,ims,ime,jms,jme,ifs,ife,jfs,jfe
Expand Down Expand Up @@ -1292,8 +1287,6 @@ subroutine prop_ls( id, & ! for debug

!*** purpose: advance level function in time

! Jan Mandel August 2007 - February 2008

!*** description
!
! Propagation of closed curve by a level function method. The level function
Expand Down Expand Up @@ -1657,7 +1650,8 @@ subroutine tend_ls( id, &
!*** local
real:: te,diffLx,diffLy,diffRx,diffRy, &
diffCx,diffCy,diff2x,diff2y,grad,rr, &
ros_back,ros_wind,ros_slope,advx,advy,scale,nvx,nvy,speed,tanphi
ros_base,ros_wind,ros_slope,ros_back,advx,advy,scale,nvx,nvy, &
speed,tanphi
integer::i,j,itso,iteo,jtso,jteo
character(len=128)msg

Expand Down Expand Up @@ -1743,10 +1737,10 @@ subroutine tend_ls( id, &

! get rate of spread from wind speed and slope

call fire_ros(ros_back,ros_wind,ros_slope, &
call fire_ros(ros_base,ros_wind,ros_slope, &
nvx,nvy,i,j,fp)

rr=ros_back + ros_wind + ros_slope
rr=ros_base + ros_wind + ros_slope
if(fire_grows_only.gt.0)rr=max(rr,0.)

! set ros for output
Expand All @@ -1760,7 +1754,7 @@ subroutine tend_ls( id, &
else

! normal direction backing rate only
te = - ros_back*grad
te = - ros_base*grad

! advection in wind direction
if (abs(speed)> eps) then
Expand Down Expand Up @@ -1923,7 +1917,7 @@ real function speed_func(diffCx,diffCy,dx,dy,i,j,fp)
type(fire_params),intent(in)::fp
!*** local
real::scale,nvx,nvy,r
real::ros_back , ros_wind , ros_slope
real::ros_base , ros_wind , ros_slope
real, parameter:: eps=epsilon(0.0)
!*** executable
! normal direction, from central differences
Expand All @@ -1933,10 +1927,10 @@ real function speed_func(diffCx,diffCy,dx,dy,i,j,fp)

! get rate of spread from wind speed and slope

call fire_ros(ros_back,ros_wind,ros_slope, &
call fire_ros(ros_base,ros_wind,ros_slope, &
nvx,nvy,i,j,fp)

r=ros_back + ros_wind + ros_slope
r=ros_base + ros_wind + ros_slope
if(fire_grows_only.gt.0)r=max(r,0.)
speed_func=r

Expand Down
96 changes: 22 additions & 74 deletions phys/module_fr_sfire_driver.F
View file
@@ -1,67 +1,17 @@
! SFIRE - Spread fire model in WRF-Fire
!
!*** Jan Mandel August 2007 - July 2010
!*** email: Jan.Mandel@gmail.com

! For support please subscribe to the wrf-fire mailing list at NCAR at
! http://mailman.ucar.edu/mailman/listinfo/wrf-fire
! or go to http://www.openwfm.org/wiki/WRF-Fire_user_support

! This module is the only entry point from WRF-ARW to the wildland
! fire model. The call to sfire_driver advances the fire model by
! one timestep. The fire model inputs the wind and outputs
! temperature and humidity tendencies. The fire model also inputs a
! This module is the entry point from WRF ARW to the wildland
! fire module. The call to sfire_driver advances the fire module by
! one timestep. The fire module inputs the wind and outputs
! temperature and humidity tendencies. The fire module also inputs a
! number of constant arrays (fuel data, topography). Additional
! arguments are model state (for data assimilation) and constant arrays
! the model gives to WRF for safekeeping because it is not allowed
! to save anything.

! This code is described in [1]. The fire model is coupled with WRF
! but the fire code itself is not dependent on WRF in any way other
! than calls to few WRF utilities from module_fr_sfire_util. This
! model uses a level set function method for advancing the fireline.
! It is a reimplementation of an earlier model, which used fireline
! propagation by tracers and was coupled with the Clark-Hall
! atmospheric code, described in [2]. For WRF documentation see [3].

! If you use this code, please acknowledge our work by citing [1].
! Thank you.

! Acknowledgements
!
! Contributions to development of the level set
! method by Mijeong Kim. Contribution to fuel calculation by Volodymyr
! Kondratenko.
! Contributions to this wildland fire module have come from individuals at
! NCAR, the U.S.D.A. Forest Service, the Australian Bureau of Meteorology,
! and the University of Colorado at Denver.
!
! The fire physics is adapted from an earlier code by Terry
! L. Clark, Janice L. Coen, and Don Latham [2]. The interface with
! WRF is adapted from a code by Ned Patton for coupling of the earlier
! fire model with WRF, based on the foundation outlined in [4].
! The changes in WRF infrastructure and support of
! refined fire grids were provided by John Michalakes.

! Jonathan D. Beezley has set up and maintained the WRF build and
! execution environment, provided software engineering infrastructure
! including synchronization with the WRF repository, and was responsibe
! for all aspects of WRF modification.

! Refefences
!
! [1] Jan Mandel, Jonathan D. Beezley, Janice L. Coen, and Minjeong Kim,
! Data Asimilation for Wildland Fires: Ensemble Kalman filters in
! coupled atmosphere-surface models, IEEE Control Systems Magazine 29,
! Issue 3, June 2009, 47-65.

! [2] T. L. Clark, J. Coen, and D. Latham, Description of a coupled
! atmosphere-fire model, Intl. J. Wildland Fire, vol. 13, pp. 49-64,
! 2004
!
! [3] http://www.mmm.ucar.edu/wrf/OnLineTutorial/Introduction/index.html
!
! [4] Edward G. Patton and Janice L. Coen, WRF-Fire: A Coupled
! Atmosphere-Fire Module for WRF, Preprints of Joint MM5/Weather
! Research and Forecasting Model Users' Workshop, Boulder, CO,
! June 22-25, 2004, pp. 221-223, NCAR

module module_fr_sfire_driver
! use this module for standalone call, you only need to provide some mock-up wrf modules
Expand Down Expand Up @@ -135,22 +85,20 @@ subroutine sfire_driver_em ( grid , config_flags &

!*** executable

if(fire_ifun_start.le.1)call print_id ! print id only once, during initialization
! populate our structures from wrf

! pointers to be passed to fire spread formulas
fp%vx => grid%uf ! fire winds
fp%vy => grid%vf ! fire winds
! pointers to be passed to fire rate of spread formulas
fp%vx => grid%uf ! W-E winds used in fire module
fp%vy => grid%vf ! S-N winds used in fire module
fp%zsf => grid%zsf ! terrain height
fp%dzdxf => grid%dzdxf ! terrain grad
fp%dzdyf => grid%dzdyf ! terrain grad
fp%bbb => grid%bbb ! spread formula coeff
fp%betafl => grid%betafl ! spread formula coeff
fp%phiwc => grid%phiwc ! spread formula coeff
fp%r_0 => grid%r_0 ! spread formula coeff
fp%fgip => grid%fgip ! spread formula coeff
fp%ischap => grid%ischap ! spread formula coeff
fp%bbb => grid%bbb ! a rate of spread formula coeff
fp%betafl => grid%betafl ! a rate of spread formula variable
fp%phiwc => grid%phiwc ! a rate of spread formula coeff
fp%r_0 => grid%r_0 ! a rate of spread formula variable
fp%fgip => grid%fgip ! a rate of spread formula coeff
fp%ischap => grid%ischap ! a rate of spread formula switch

! get ignition data
call fire_ignition_convert (config_flags,fire_max_ignitions,fire_ignition_longlat, &
Expand Down Expand Up @@ -352,7 +300,7 @@ subroutine sfire_driver_phys (ifun,need_lfn_update, &
u_frame,v_frame, & ! velocity offset
unit_fxlong,unit_fxlat, & ! fxlong, fxlat units in m
fire_crwn_hgt, & ! lowest height crown fire heat is released (m)
fire_ext_grnd, & ! extinction depth of ground fire heat (m)
fire_ext_grnd, & ! extinction depth of surface fire heat flux (m)
fire_ext_crwn, & ! wind height for vertical interploation to fire spread
fire_wind_height

Expand Down Expand Up @@ -383,14 +331,14 @@ subroutine sfire_driver_phys (ifun,need_lfn_update, &

real, intent(out), dimension(ims:ime, jms:jme):: & ! redundant arrays, for display purposes only (atm grid)
avg_fuel_frac, & ! average fuel fraction
grnhfx, & ! heat flux from ground fire (W/m^2)
grnqfx, & ! moisture flux from ground fire (W/m^2)
grnhfx, & ! heat flux from surface fire (W/m^2)
grnqfx, & ! moisture flux from surface fire (W/m^2)
canhfx, & ! heat flux from crown fire (W/m^2)
canqfx ! moisture flux from crown fire (W/m^2)

real, intent(out), dimension(ifms:ifme, jfms:jfme):: & ! redundant arrays, for display only, fire grid
fgrnhfx, & ! heat flux from ground fire (W/m^2)
fgrnqfx, & ! moisture flux from ground fire (W/m^2)
fgrnhfx, & ! heat flux from surface fire (W/m^2)
fgrnqfx, & ! moisture flux from surface fire (W/m^2)
fcanhfx, & ! heat flux from crown fire (W/m^2)
fcanqfx, & ! moisture flux from crown fire (W/m^2)
ros ! fire rate of spread (m/s)
Expand Down Expand Up @@ -698,7 +646,7 @@ subroutine sfire_driver_phys (ifun,need_lfn_update, &
s = 1./(ir*jr)
do j=jts,jte
do i=its,ite
! scale ground fluxes to get the averages
! scale surface fluxes to get the averages
avg_fuel_frac(i,j)=avg_fuel_frac(i,j)*s
grnhfx(i,j)=fire_atm_feedback*grnhfx(i,j)*s
grnqfx(i,j)=fire_atm_feedback*grnqfx(i,j)*s
Expand Down
3 changes: 0 additions & 3 deletions phys/module_fr_sfire_model.F
View file
@@ -1,7 +1,4 @@
!
!*** Jan Mandel October 2007 email: jmandel@ucar.edu or Jan.Mandel@gmail.com
!

#define DEBUG_OUT

module module_fr_sfire_model
Expand Down

0 comments on commit 54a1084

Please sign in to comment.