V3FIT
errray1.f
1  SUBROUTINE errray1(errcode,psi,thet,x,y,err,tol,tola)
2  USE precision
3  IMPLICIT NONE
4  INTEGER :: errcode
5  REAL(rprec) :: psi,thet,x,y,err,tol,tola
6 c.......................................................................
7  IF (errcode .eq. 1)then
8  WRITE (6,'(a)')'rayatt: sliding window search out of range'
9  WRITE (6,100)' psi = ',psi,' thet = ',thet
10  WRITE (6,100)' x = ',x, ' y = ',y
11  RETURN
12  ENDIF
13 c
14  IF (errcode .eq. 2)then
15  WRITE (6,'(a)')'rayatt: sliding window search fails (10)'
16  WRITE (6,100)' psi = ',psi,' thet = ',thet
17  WRITE (6,100)' x = ',x, ' y = ',y
18  RETURN
19  ENDIF
20 c
21  IF (errcode .eq. 3)then
22  WRITE (6,'(a)')'rayatt: Newton-Raphson INTerations fails (20)'
23  WRITE (6,100)' psi = ',psi,' thet = ',thet
24  WRITE (6,100)' err = ',err,' tol = ',tol
25  IF (abs(err) .gt. abs(tola))then
26  WRITE (6,'(a)')'rayatt: err > max.tol'
27  RETURN
28  ELSE
29  WRITE (6,'(a)')'rayatt: tol < err < max. tol (warning)'
30  RETURN
31  ENDIF
32  ENDIF
33 c
34  WRITE(6,*) 'error in errray1--errray999'
35  stop
36  100 FORMAT(a,1pe12.4,a,1pe12.4)
37  END