Using FORTRAN write statement with parallel NONMEM
Dear all,
I'm using the BAYES method and attempting to write out MCMC samples for
individual parameters per the method shown in example 8 of the NONMEM 7.2
documentation. This involves using a FORTRAN write statement in verbatim
code at the end of the PK block:
"IF(BAYES_EXTRA==1 .AND. ITER_REPORT>=0 .AND. TIME==0.0) THEN
"WRITE(50,98) ITER_REPORT,ID,CL,V1,Q,V2
"98 FORMAT(I12,1X,F14.0,4(1X,1PG12.5))
"ENDIF
This works as expected when running on a single processor. The desired
samples are written to a file called fort.50. When I try to run this in
parallel per "Introduction to NONMEM 7.2.0" p. 70, the run completes
successfully, but only a subset of the MCMC samples are written to fort.50.
I'm guessing this is the subset that are computed on the master node.
Comparable files for the worker nodes are nowhere to be found. In
particular they do not appear in the directories named worker1, worker2,
etc.
Any idea where the samples are being written or alternative approaches to
access those values?
Thanks,
Bill