c*************************************
c	JAN2001/mcnew_sub.inc 9.01.2001
c	include file for mc_photo.for
c	29.11.99 HS
c	8.12.99 inc2 for mc_photo_tab2.for
c	20.9.00 inc3 for mc_photo_tab3.for

c       25.9.2000:
c	common file for mcnew_asy, mcnew_tes
c*************************************
c	mc_sub.inc

c	nn = maximum number of particles (including duplicates)
c	nnp = maximum number of photons
c	nnt =maximum number of test directions
c	mmx,mmy = maximum number of subdivision in tab-routines

	implicit double precision (a-h,o-z)

	parameter(nn=1000000)
	parameter(nnp=10000000)
	parameter(nnt=1000)
	parameter(mmx=100,mmy=100)
	parameter(mmm=mmx*mmy)

	dimension xx(nn),yy(nn),zz(nn),ss(nn),ss2(nn)
	double precision ly,lx

	common /com1/ xx,yy,zz,ss,ss2
	common /com2/ xl,yl,sy,ly,lx,zmin,zmax
	common /com3/ npart,ntest,nbox,icheck,nboxvis,npartori
	common /com4/ ssca,w

	real lattab(nnp),fiitab(nnp),wtab(nnp),scatab(nnp)
	common /com5/ lattab,fiitab,wtab,scatab

	dimension extest(nnt),eytest(nnt),eztest(nnt)
	common /com6/ extest,eytest,eztest

	dimension testw(nnt),testd(nnt),testn(nnt),testm(nnt,10)
	common /com7/ testw,testd,testn
	common /com8/ testm

	dimension indtab(nn),numtab(mmm),mintab(mmm)
	dimension dtab(mmm), cxtab(mmm), cytab(mmm),dtab2(mmm)

	common /com9/indtab,numtab,mintab
	common /com10/cxtab,cytab,dtab,dtab2

	common /com11/mx,my,npot_sum,npot_lkm,nbox_lkm,nbox_sum


	common /com12/ ipart_now,method
	common /com13/ sredu,zredu,tau

	dimension testd_apu0(nnt)
	common /com14/ testd_apu0

c	addition jan2001 9.01.2001
	common /com15/ppar1,ppar2,ppar3
	common /com16/exold,eyold,ezold

	data pii/3.141592654/
        data pii2/6.283185308/
	data ar,ra/0.017453293,57.29577951/


c	for power-law phase-function

	parameter(nntables=1001)

	dimension power_inte(nntables),power_alpha(nntables)
	dimension power_dinte(nntables)
	common /com_power1/ power_inte, power_alpha, power_dinte
	common /com_power2/ power_expo,c_n