EQUIL - UPDATE - Officially deprecates old Hansen and Marklin eq_types to avoid convusion with correct modern interfaces

Author: logan-nc

equil/equil.f

@@ -174,10 +174,13 @@ SUBROUTINE equil_read(unit, op_psihigh, op_psilow)
       CASE("t7")
          CALL read_eq_t7
       CASE("marklin_direct")
+         CALL program_stop("Error: Deprecated eq_type "//TRIM(eq_type))
          CALL read_eq_marklin_direct
       CASE("marklin_inverse")
+         CALL program_stop("Error: Deprecated eq_type "//TRIM(eq_type))
          CALL read_eq_marklin_inverse
       CASE("hansen_inverse")
+         CALL program_stop("Error: Deprecated eq_type "//TRIM(eq_type))
          CALL read_eq_hansen_inverse
       CASE("pfrc")
          CALL read_eq_pfrc