If you want to do this yourself, then you could follow these steps to add the functionality that you are looking for.
In the raman_scattering subroutine of the raman.F file, add
#include "mafdecls.fh"
to the list of other included header files. Then comment out the line
parameter (numpts = 1000)
Finally, add the lines
if (.not.rtdb_get(rtdb,'raman:numpts',mt_int,1,numpts))
& numpts=1000
After recompiling/linking the code, you should be able to control the number of points in the plot by adding
set raman:numpts <integer>
to your input file. If you don't have the set directive in your input, then it will default to 1000, which is what it is currently hardcoded as.
|