standalone/random_fire.py

   1 #!/usr/bin/env python
   2
   3 from math import sqrt
   4 from random import normalvariate as randn
   5 from random import randint,seed
   6 import subprocess as sp
   7
   8 sigma=1.
   9 nx=256
  10 ny=256
  11 dx=6.
  12 dy=6.
  13 igntime=2.
  14 runtime=30*60
  15
  16 slpsigma=sigma*.1/sqrt(2)
  17 windsigma=sigma*5./sqrt(2)
  18
  19 cenx=nx*dx/2.
  20 ceny=ny*dy/2.
  21 dz=sqrt(dx**2+dy**2)
  22 ignr=dz*2
  23 historys=runtime
  24 timestep=dz*1.0/6.
  25
  26 ignxsigma=sigma*dx*nx/10.
  27 ignysigma=sigma*dy*ny/10.
  28
  29 seed()
  30 slpx=randn(0,slpsigma)
  31 slpy=randn(0,slpsigma)
  32 windx=randn(0,windsigma)
  33 windy=randn(0,windsigma)
  34 fuel=randint(1,13)
  35 ignx=randn(cenx,ignxsigma)
  36 igny=randn(ceny,ignysigma)
  37
  38 args=['--nx',nx,'--ny',ny,'--dx',dx,'--dy',dy,
  39       '--windx',windx,'--windy',windy,'--slopex',slpx,'--slopey',slpy,
  40       '--fuelcat',fuel,'--timestep',timestep,'--runtime',runtime,'--history',historys,
  41       '--ignx1',ignx,'--ignx2',ignx,'--igny1',igny,'--igny2',igny,
  42       '--ignr',ignr,'--ignt1',igntime,'--ignt2',igntime]
  43
  44 args=[ str(a) for a in args]
  45
  46 print ' '.join(['fire.py'] +args)
  47 p=sp.Popen(['python','fire.py']+args)
  48 p.communicate()
  49
  50 f=open('params.txt','w')
  51 for i in xrange(0,len(args),2):
  52     f.write('%s=%s\n' % (args[i][2:],args[i+1]))
  53 f.close()