Thanks to everyone who responded to my last message. It has been suggested that the error is in the subroutine: equation_interface_MP_method_of_characteristics this is what the code looks like for this routine: !________________________________________________________________________ !//////////////////////////////////////////////////////////////////////// SUBROUTINE method_of_characteristics(ib20,ib21,xnorm0,xnorm1,q0,q1,q) INTEGER, INTENT(in) :: ib20,ib21 DOUBLE PRECISION, INTENT(in) :: xnorm0(2),xnorm1(2) DOUBLE PRECISION, INTENT(in) :: q0(neq),q1(neq) DOUBLE PRECISION, INTENT(out) :: q(neq) SELECT CASE(eqn%title) CASE('linear euler') CALL method_of_chars_linear_euler(ib20,ib21,xnorm0,xnorm1,q0,q1,q) CASE('euler') CALL method_of_chars_euler(ib20,ib21,xnorm0,xnorm1,q0,q1,q) CASE DEFAULT PRINT *, 'ERROR in SUBROUTINE method_of_characteristics , no equation selected' PRINT *, "equation = ",eqn%title STOP ENDSELECT END SUBROUTINE method_of_characteristics !________________________________________________________________________ !//////////////////////////////////////////////////////////////////////// I suspect that the code is actually not crashing in this subroutine, but in euler_MP_method_of_chars_euler which it calls. Either way it seems that I still am going to have a problem figuring out what the problem might be. _______________________________________________________________________ | Pete Bismuti | | Supercomputer Computations Research Institute | | Florida State University - Department of Mathematics | | [log in to unmask] (904)644-6263 | |_____________________________________________________________________| %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%