wrfv2_fire/phys/fuel_burnt_test.F

   1 program fuel_burnt_test
   2 use module_fr_sfire_core
   3 use module_fr_sfire_phys
   4 implicit none
   5
   6 ! ims,ime,jms,jme   in   memory dimensions
   7 ! its,ite,jts,jte   in   tile dimensions (cells where fuel_frac computed)
   8 ! ifs,ife,jfs,jfe   in   fuel_frac memory dimensions
   9 ! ir,jr,            in   refinement - quadrature cells per fire cell, even, at least 2
  10 ! lfn               in   level function, at nodes at midpoints of cells
  11 ! tign              in   ignition time, at nodes at nodes at midpoints of cells
  12 ! fuel_time         in   time constant of fuel, per cell
  13 ! time_now          in   time now
  14 ! fuel_frac         out  fraction of fuel remaining, per cell
  15
  16 !*** local
  17 integer:: ids,ide,jds,jde,ims,ime,jms,jme,ir,jr,k
  18 parameter(ims=1,ime=2,jms=1,jme=2,ir=2,jr=2)
  19 real:: time_now,off,fuel_left_f
  20 real,dimension(ims:ime,jms:jme)::lff,tif,lfn0,fuel_time
  21 integer::i,iu
  22
  23 lfn0(1,1)=0.419279
  24 lfn0(1,2)=-1.9766893E-02
  25 lfn0(2,1)=6.310914
  26 lfn0(2,2)=5.983462
  27 time_now=2.0
  28 fuel_time(1,1)= 8.23529
  29 iu=1
  30 open(iu,file='tmp.txt')
  31 do k=1,23
  32
  33   off = (k-1)*0.05
  34   lff(1,1)=lfn0(1,1)-off
  35   lff(1,2)=lfn0(1,2)-off
  36   lff(2,1)=lfn0(2,1)-off
  37   lff(2,2)=lfn0(2,2)-off
  38
  39   tif(1,1)=2.000000
  40   tif(1,2)=1.565244
  41   tif(2,1)=2.000000
  42   tif(2,2)=2.000000
  43
  44   fuel_left_f=fuel_left_cell_1( &
  45                  lff(1,1),lff(1,2),lff(2,1),lff(2,2),&
  46                  tif(1,1),tif(1,2),tif(2,1),tif(2,2),&
  47                  time_now, fuel_time(1,1))
  48                                               ! premultiply by weight
  49   fuel_left_f=1-fuel_left_f
  50   print *,'fuel_left_f=',fuel_left_f
  51   write(iu,*) fuel_left_f
  52 end do
  53 close(iu)
  54 end program
  55