Writing output with > 1000 orientations
DDSCAT 7.3 added support for writing output with > 1000 orientations. It works very well, but there is a vestigial check in the REAPAR function that disables much of this functionality.
In version 7.3.1's code (version 160601), in reapar.f90:
1470 ! check that user has not requested more than 1000 orientations and
1471 ! IWRKSC=1
1472
1473 IF(IWRKSC>0.AND.(NBETA*NTHETA*NPHI)>1000)THEN
1474 CALL ERRMSG('FATAL','REAPAR', &
1475 'error: if IWRKSC=1, NBETA*NTHETA*NPHI must be .le. 1000')
1476 ENDIF
Given the changes to the code in 7.3, this can be removed.