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 |
|_____________________________________________________________________|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|