phase shifts at caustics
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| logical, | intent(in) | :: | qleq0 |
SUBROUTINE IncPhaseIfCaustic( qleq0 ) !! phase shifts at caustics LOGICAL, INTENT( IN ) :: qleq0 IF ( IsAtCaustic( qleq0 ) ) THEN phase = phase + pi / 2. END IF END SUBROUTINE IncPhaseIfCaustic