real pi_ parameter (pi_=3.1415926) real K_,me_c2,c_ parameter (K_=0.307e6,me_c2=511.e+3,c_=3.0e10) real Z_,A_,I_ parameter (Z_=6,A_=12.01,I_=80.) real density,refindex,heatcap parameter (density=3.52,refindex=2.42,heatcap=0.47) real rhoCp,kappa,ems,SBc,Tambient parameter (rhoCp=density*heatcap) c parameter (kappa=20.) ! for type IIa diamonds parameter (kappa=10.) ! for type Ib diamonds (typical range 5-15) parameter (ems=0.02) parameter (SBc=5.67e-12) parameter (Tambient=300.) real thickness,width,height real kappaWire,lengthWire,diameterWire parameter (kappaWire=1.3) parameter (lengthWire=0.8) parameter (diameterWire=18e-3) ! drain all heat at boundaries real beamsizeX,beamsizeY,hpower common /changeables/thickness,width,height, + beamsizeX,beamsizeY,hpower integer NX,NY parameter (NX=30,NY=30) double precision Tgrid_(30,30) double precision time_ common /dprec/Tgrid_,time_