22 fprintf(
ioQQQ,
"PROBLEM the illuminate command has come before any "
23 "radiation field shape commands.\nThis must come after the field"
24 " is specified.\nSorry.\n");
32 else if( p.
nMatch(
"REVE" ) )
34 else if( p.
nMatch(
"SYMM" ) )
46 double AngleIllumRadian = 0.;
52 p.
NoNumb(
"illumination angle");
58 AngleIllumRadian = a/
RADIAN;
59 if( AngleIllumRadian < 0. || AngleIllumRadian >=
PI/2. )
61 fprintf(
ioQQQ,
" Angle of illumination must be between 0 and 90 degrees "
62 "or 0 and pi/2 radians.\n" );
71 (
realnum)(1./cos(AngleIllumRadian));