This subroutine initializes the values in the phasemod array here instead of in the main program. It occurs here now as it is easier to modify. The call to phasemod_init is:
call phasemod_init(phasemod,nphase)