subroutine radial write(*,*) 'This option is not available.' write(*,*) 'Stopping in radstub' stop return end