double precision aa0thet(-maxthetyp1:maxthetyp1,
&-maxthetyp1:maxthetyp1,-maxthetyp1:maxthetyp1,2),
& aathet(maxtheterm,-maxthetyp1:maxthetyp1,
double precision aa0thet(-maxthetyp1:maxthetyp1,
&-maxthetyp1:maxthetyp1,-maxthetyp1:maxthetyp1,2),
& aathet(maxtheterm,-maxthetyp1:maxthetyp1,