real rho,cheat,kcond,hconv,ems,SBconst,Tamb
      parameter (rho=1.05e3) ! kg/m^3, density of polystyrene
      parameter (cheat=1.3e3)! J/kg/K, heat capacity for polystyrene
      parameter (kcond=0.08) ! W/m/K, thermal conductivity coef for polystyrene
      parameter (hconv=10)   ! W/m^2/K, heat transfer coef, free air convection
      parameter (ems=0.10)   ! emissivity of polystyrene
      parameter (SBconst=5.67e-8) ! W/m^2/K^4, Stefan-Boltzmann constant
      parameter (Tamb=293)   ! K, ambient laboratory temperature
      real pi
      parameter (pi=3.14159265)

      real fiberdim
      common /fusiong/fiberdim(3)
      data fiberdim/2*2e-3,0.02/ ! m

      integer ndivx,ndivy,ndivz
      parameter (ndivx=11,ndivy=11,ndivz=31)
      integer nstep,nbuf,obuf,pad
      real time,tstep
      double precision Tfield
      common /fusionf/time,tstep,nstep,nbuf,obuf,pad,
     +                Tfield(ndivx,ndivy,ndivz,2)
      data time,tstep/0,0/
      data nbuf,obuf/1,2/

      real Scent,Ssize,Snorm,Spower
      common /fusions/Scent(3),Ssize(3),Snorm,Spower
      data Scent/0,0,0/                  ! m
      data Ssize/3e-4,3e-4,3e-4/         ! m
      data Snorm/0/                      ! computed in Source()
      data Spower/0/                     ! set in setPower()

      double precision Qcond,Qconv,Qrad
      common /fusionq/Qcond,Qconv,Qrad
      data Qcond,Qconv,Qrad/0,0,0/