From: Stephen Sinclair <radarsat1@gmail.com>
Date: Wed, 14 Sep 2022 13:10:00 +0000
Subject: Patch to remove write statements from Fortran code.

For reasons unknown, it appears that if ANY "write" statements appear in a Fortran source
file, then the resulting object file contains a string in the `.rodata` section with the
path to the source file, which harms reproducibility.  This is despite the correct usage
of `-ffile-prefix-map` during compilation, so it could be a bug in the Fortran compiler.
In any case, removing the write statements fixes the issue, though it may suppress useful
error messages.

diff --git a/externals/hairer/hem5.f b/externals/hairer/hem5.f
index 7f2b6f5..dd490cb 100644
--- a/externals/hairer/hem5.f
+++ b/externals/hairer/hem5.f
@@ -184,8 +184,8 @@ C *** *** *** *** *** *** ***
       MODE=IWK(14)
 C -------- MODE, THE CHOICE OF LINEAR ALGEBRA
       IF ((IWK(14).LE.-1).OR.(IWK(14).GE.6)) THEN
-        WRITE(6,*)' WRONG CHOICE OF IWK(14) (MODE):'
-     &             ,IWK(14)
+c$$$        WRITE(6,*)' WRONG CHOICE OF IWK(14) (MODE):'
+c$$$     &             ,IWK(14)
         ARRET=.TRUE.
       ELSE 
         MODE=IWK(14)
@@ -199,27 +199,27 @@ C -------- MODE, THE CHOICE OF LINEAR ALGEBRA
       IF (MODE.EQ.4) LDF=LDG 
       IF (MODE.GE.4) THEN
         IF (LDG.LE.0) THEN
-          WRITE(6,*)' IWK(1) (LDG) MUST BE POSITIVE'
+c$$$          WRITE(6,*)' IWK(1) (LDG) MUST BE POSITIVE'
           ARRET=.TRUE.
         END IF
         IF (LDF.LE.0) THEN
-          WRITE(6,*)' IWK(2) (LDF) MUST BE POSITIVE'
+c$$$          WRITE(6,*)' IWK(2) (LDF) MUST BE POSITIVE'
           ARRET=.TRUE.
         END IF
         IF ((NMRC.EQ.0).OR.(NBLK.EQ.0)) THEN
-          WRITE(6,*)' IWK(21) (NMRC) AND IWK(22) (NBLK) MUST
-     &        BE POSITIVE'
+c$$$          WRITE(6,*)' IWK(21) (NMRC) AND IWK(22) (NBLK) MUST
+c$$$     &        BE POSITIVE'
           ARRET=.TRUE.
         END IF
         IF (IWK(28).EQ.0) IWK(28)=6
         IF (IWK(25).LT.13*NM) THEN
-           WRITE(6,*)' INTEGER WORK SPACE (IWK(25))
-     &           FOR MA28 TOO SMALL'
+c$$$           WRITE(6,*)' INTEGER WORK SPACE (IWK(25))
+c$$$     &           FOR MA28 TOO SMALL'
            IWK(25)=13*NM
         END IF
         IF (IWK(26).LT.(4*NZA+NM)) THEN
-           WRITE(6,*)' REAL WORK SPACE (IWK(26))
-     &           FOR MA28 TOO SMALL'
+c$$$           WRITE(6,*)' REAL WORK SPACE (IWK(26))
+c$$$     &           FOR MA28 TOO SMALL'
            IWK(26)=4*NZA+NM
         END IF
       END IF
@@ -267,8 +267,8 @@ C -------- UROUND   SMALLEST NUMBER SATISFYING 1.D0+UROUND>1.D0
       ELSE
         UROUND=WK(1)
         IF(UROUND.LE.1.D-35.OR.UROUND.GE.1.D0)THEN
-         WRITE(6,*)' WHICH MACHINE DO YOU HAVE? YOUR UROUND WAS:'
-     &                            ,WK(1)
+c$$$         WRITE(6,*)' WHICH MACHINE DO YOU HAVE? YOUR UROUND WAS:'
+c$$$     &                            ,WK(1)
          ARRET=.TRUE.
         END IF
       END IF
@@ -278,8 +278,8 @@ C -------  SAFETY FACTOR -------------
       ELSE
         SAFE=WK(2)
         IF(SAFE.GE.1.D0.OR.SAFE.LE.1.D-4)THEN
-         WRITE(6,*)' CURIOUS INPUT FOR SAFETY FACTOR WK(2)='
-     &                            ,WK(2)
+c$$$         WRITE(6,*)' CURIOUS INPUT FOR SAFETY FACTOR WK(2)='
+c$$$     &                            ,WK(2)
          ARRET=.TRUE.
         END IF
       END IF
@@ -388,7 +388,7 @@ C ------- PREPARE THE ENTRY-POINTS FOR THE ARRAYS IN WK -----
 C ------ TOTAL STORAGE REQUIREMENT -----------
       IS=IDO+LRDO
       IF(IS.GT.LWK)THEN
-       WRITE(6,*)' INSUFFICIENT STORAGE FOR WK, MIN. LWK=',IS
+c$$$       WRITE(6,*)' INSUFFICIENT STORAGE FOR WK, MIN. LWK=',IS
        ARRET=.TRUE.
       END IF
 C ------- PREPARE THE ENTRY-POINTS FOR THE ARRAYS IN IWK -----
@@ -406,7 +406,7 @@ C ------- PREPARE THE ENTRY-POINTS FOR THE ARRAYS IN IWK -----
 C ------ TOTAL STORAGE REQUIREMENT -----------
       IS=ICONT+LIDO
       IF(IS.GT.LIWK)THEN
-       WRITE(6,*)' INSUFFICIENT STORAGE FOR IWK, MIN. LIWK=',IS
+c$$$       WRITE(6,*)' INSUFFICIENT STORAGE FOR IWK, MIN. LIWK=',IS
        ARRET=.TRUE.
       END IF
 C ------ WHEN A FAIL HAS OCCURED, WE RETURN WITH IDID=-1
@@ -1240,27 +1240,27 @@ C --- STEP IS REJECTED WITHOUT PROJECTION
       GOTO 909
 C --- FAIL EXIT
  176  CONTINUE
-      WRITE(6,979)T   
-      WRITE(6,*) ' MATRIX IN ASET IS  SINGULAR, IER=',IER
+c$$$      WRITE(6,979)T   
+c$$$      WRITE(6,*) ' MATRIX IN ASET IS  SINGULAR, IER=',IER
       IDID=-4
       RETURN
  177  CONTINUE
-      WRITE(6,979)T   
-      WRITE(6,*) ' STEP SIZE T0O SMALL, H=',H
+c$$$      WRITE(6,979)T   
+c$$$      WRITE(6,*) ' STEP SIZE T0O SMALL, H=',H
       IDID=-3
       RETURN
  178  CONTINUE
-      WRITE(6,979)T   
-      WRITE(6,*) ' MORE THAN NMAX =',NMAX,'STEPS ARE NEEDED' 
+c$$$      WRITE(6,979)T   
+c$$$      WRITE(6,*) ' MORE THAN NMAX =',NMAX,'STEPS ARE NEEDED' 
       IDID=-2
       RETURN
 C --- EXIT CAUSED BY SOLOUT
  179  CONTINUE
-      WRITE(6,979)T
-      WRITE(6,*) 'INITIAL PROJECTION: NO CONVERGENCE' 
+c$$$      WRITE(6,979)T
+c$$$      WRITE(6,*) 'INITIAL PROJECTION: NO CONVERGENCE' 
       IDID=-5
       RETURN
- 979  FORMAT(' EXIT OF HEM5 AT X=',E18.4) 
+c$$$ 979  FORMAT(' EXIT OF HEM5 AT X=',E18.4) 
       END
 C*****************************************************************
       SUBROUTINE GIINUM(MODE,NQ,NV,NU,NL,NDIM,MDIM,NMDIM,NM, 
@@ -1493,13 +1493,13 @@ C         WRITE(6,*)'MA28AD'
      &           IUMF(IIK),IUMF(IIW),XUMF,ERROR )
          ISTAT(7)=ISTAT(7)+1
          IF(ERROR.LT.0) THEN
-           WRITE(6,*)'ERROR IN UMDREFAC'
+c$$$           WRITE(6,*)'ERROR IN UMDREFAC'
            IREFAC=1
            GOTO 100
          END IF
          IF(RMIN.LT.THRSH2) THEN
            IREFAC=1
-           WRITE(6,*)'RMIN TOO SMALL'
+c$$$           WRITE(6,*)'RMIN TOO SMALL'
            GOTO 100
          END IF
       END IF
@@ -1772,7 +1772,7 @@ C ---    TEST FOR CONVERGENCE    --------------------------------------
 C
       IF(ERR.LT.EPS) GOTO 200
       IF(ERR.GE.ERROLD*2.D0 .OR. K.GE.ITMAX) THEN
-        WRITE(6,*)'NO CONVERGENCE IN PROJECTION'
+c$$$        WRITE(6,*)'NO CONVERGENCE IN PROJECTION'
 C ---    NO CONVERGENCE STEP IS REJECTED   --------------------------------------------
         ACCEPT = .FALSE.
         RETURN
@@ -1808,7 +1808,7 @@ cc
      &     NZA,LIPAR,LIUMF,LXUMF,LISTA,ISTAT,IPAR,INDG,INDG,
      &     INDFL,INDA,IP,IUMF,XUMF,B,GQ,GQ,FL,AVALUE,IER)
       IF (IER.NE.0) THEN
-        WRITE(6,*)'MATRIX SINGULAR IN APROJ'
+c$$$        WRITE(6,*)'MATRIX SINGULAR IN APROJ'
         ACCEPT=.FALSE.
         RETURN
       END IF
@@ -1836,7 +1836,7 @@ C -- ERROR ESTIMATION
       END DO
       ERR=SQRT(ERR/NV)
       IF(ERR.GT.EPS)ACCEPT=.FALSE.
-      IF(ERR.GT.EPS)write(6,*)'reject after proj',err
+c$$$      IF(ERR.GT.EPS)write(6,*)'reject after proj',err
       RETURN
       END
 C 
@@ -1989,7 +1989,7 @@ c -- qd et vd contiennent q3 et v3
      &     LIPAR,LIUMF,LXUMF,LISTA,ISTAT,IPAR,INDG1,INDGD,
      &     INDFL,INDA,IP,IUMF,XUMF,B,GQ1,GD,FL,AVALUE,IER)
       IF (IER.NE.0) THEN
-         WRITE(6,*)'MATRIX SINGULAR IN CMAT4'
+c$$$         WRITE(6,*)'MATRIX SINGULAR IN CMAT4'
          RETURN
       END IF
       CALL ASOL(MODE,NM,NV,NL,LIPAR,LIUMF,LXUMF,IPAR,IUMF,IP,
diff --git a/externals/hairer/linsp.f b/externals/hairer/linsp.f
index 341b1a7..ca85d1e 100644
--- a/externals/hairer/linsp.f
+++ b/externals/hairer/linsp.f
@@ -190,28 +190,28 @@ c     range  zero to one  and  is thus altered by ma30a/ad.
 c simple data check on input variables and array dimensions.
       IF (N.GT.0) GO TO 10
       IFLAG = -8
-      IF (LP.NE.0) WRITE (LP,99999) N
-      GO TO 210
+c$$$      IF (LP.NE.0) WRITE (LP,99999) N
+      GO TO 220
    10 IF (NZ.GT.0) GO TO 20
       IFLAG = -9
-      IF (LP.NE.0) WRITE (LP,99998) NZ
-      GO TO 210
+c$$$      IF (LP.NE.0) WRITE (LP,99998) NZ
+      GO TO 220
    20 IF (LICN.GE.NZ) GO TO 30
       IFLAG = -10
-      IF (LP.NE.0) WRITE (LP,99997) LICN
-      GO TO 210
+c$$$      IF (LP.NE.0) WRITE (LP,99997) LICN
+      GO TO 220
    30 IF (LIRN.GE.NZ) GO TO 40
       IFLAG = -11
-      IF (LP.NE.0) WRITE (LP,99996) LIRN
-      GO TO 210
+c$$$      IF (LP.NE.0) WRITE (LP,99996) LIRN
+      GO TO 220
 c
 c data check to see if all indices lie between 1 and n.
    40 DO 50 I=1,NZ
         IF (IRN(I).GT.0 .AND. IRN(I).LE.N .AND. ICN(I).GT.0 .AND.
      *   ICN(I).LE.N) GO TO 50
-        IF (IFLAG.EQ.0 .AND. LP.NE.0) WRITE (LP,99995)
+c$$$        IF (IFLAG.EQ.0 .AND. LP.NE.0) WRITE (LP,99995)
         IFLAG = -12
-        IF (LP.NE.0) WRITE (LP,99994) I, A(I), IRN(I), ICN(I)
+c$$$        IF (LP.NE.0) WRITE (LP,99994) I, A(I), IRN(I), ICN(I)
    50 CONTINUE
       IF (IFLAG.LT.0) GO TO 220
 c
@@ -257,7 +257,7 @@ c duplicate element.
   110     MOVE = MOVE + 1
           LENGTH = LENGTH - 1
           JAY = IKEEP(J,3) + NEWJ1
-          IF (MP.NE.0) WRITE (MP,99993) I, J, A(JJ)
+c$$$          IF (MP.NE.0) WRITE (MP,99993) I, J, A(JJ)
           A(JAY) = A(JAY) + A(JJ)
           THEMAX = DMAX1(THEMAX,DABS(A(JAY)))
   120   CONTINUE
@@ -276,8 +276,8 @@ c perform block triangularisation.
       IF (IDISP(1).GT.0) GO TO 170
       IFLAG = -7
       IF (IDISP(1).EQ.-1) IFLAG = -1
-      IF (LP.NE.0) WRITE (LP,99992)
-      GO TO 210
+c$$$      IF (LP.NE.0) WRITE (LP,99992)
+      GO TO 220
 c
 c block triangularization not requested.
 c move structure to end of data arrays in preparation for
@@ -319,8 +319,8 @@ c transfer common block information.
       NDROP = NDROP1
       IF (LBIG) BIG = BIG1
       IF (IFLAG.GE.0) GO TO 200
-      IF (LP.NE.0) WRITE (LP,99991)
-      GO TO 210
+c$$$      IF (LP.NE.0) WRITE (LP,99991)
+      GO TO 220
 c
 c reorder off-diagonal blocks according to pivot permutation.
   200 I1 = IDISP(1) - 1
@@ -337,7 +337,7 @@ c increment growth estimate by original maximum element.
 c set flag if the only error is due to duplicate elements.
       IF (IFLAG.GE.0 .AND. MOVE.NE.0) IFLAG = -14
       GO TO 220
-  210 IF (LP.NE.0) WRITE (LP,99990)
+c$$$  210 IF (LP.NE.0) WRITE (LP,99990)
   220 RETURN
 99999 FORMAT (36X, 17HN OUT OF RANGE = , I10)
 99998 FORMAT (36X, 18HNZ NON POSITIVE = , I10)
@@ -426,23 +426,23 @@ c
 c check to see if elements were dropped in previous ma28a/ad call.
       IF (NDROP.EQ.0) GO TO 10
       IFLAG = -15
-      WRITE (6,99999) IFLAG, NDROP
-      GO TO 70
+c$$$      WRITE (6,99999) IFLAG, NDROP
+      GO TO 60
    10 IFLAG = 0
       MEPS = EPS
       NLP = LP
 c simple data check on variables.
       IF (N.GT.0) GO TO 20
       IFLAG = -11
-      IF (LP.NE.0) WRITE (LP,99998) N
+c$$$      IF (LP.NE.0) WRITE (LP,99998) N
       GO TO 60
    20 IF (NZ.GT.0) GO TO 30
       IFLAG = -10
-      IF (LP.NE.0) WRITE (LP,99997) NZ
+c$$$      IF (LP.NE.0) WRITE (LP,99997) NZ
       GO TO 60
    30 IF (LICN.GE.NZ) GO TO 40
       IFLAG = -9
-      IF (LP.NE.0) WRITE (LP,99996) LICN
+c$$$      IF (LP.NE.0) WRITE (LP,99996) LICN
       GO TO 60
 c
    40 CALL MA28DD(N, A, LICN, IVECT, JVECT, NZ, ICN, IKEEP, IKEEP(1,4),
@@ -464,7 +464,7 @@ c transfer common block information.
       RMIN = MRMIN
       IF (IFLAG.GE.0) GO TO 50
       IFLAG = -2
-      IF (LP.NE.0) WRITE (LP,99995)
+c$$$      IF (LP.NE.0) WRITE (LP,99995)
       GO TO 60
 c
 c optionally calculate the growth parameter.
@@ -476,9 +476,9 @@ c increment estimate by largest element in input matrix.
       IF (GROW .AND. N.GT.1) W(2) = THEMAX
 c set flag if the only error is due to duplicate elements.
       IF (IDUP.EQ.1 .AND. IFLAG.GE.0) IFLAG = -14
-      GO TO 70
-   60 IF (LP.NE.0) WRITE (LP,99994)
-   70 RETURN
+      GO TO 60
+c$$$   60 IF (LP.NE.0) WRITE (LP,99994)
+   60 RETURN
 99999 FORMAT (39H ERROR RETURN FROM MA28B/BD WITH IFLAG=, I4/I7, 4H ENT,
      * 39HRIES DROPPED FROM STRUCTURE BY MA28A/AD)
 99998 FORMAT (36X, 17HN OUT OF RANGE = , I10)
@@ -591,7 +591,7 @@ c   it will be executed nz times in total.
 c perform some validity checks on iold and jold.
           IF (IOLD.LE.N .AND. IOLD.GT.0 .AND. JOLD.LE.N .AND.
      *     JOLD.GT.0) GO TO 30
-          IF (LP.NE.0) WRITE (LP,99999) I, A(I), IOLD, JOLD
+c$$$          IF (LP.NE.0) WRITE (LP,99999) I, A(I), IOLD, JOLD
           IFLAG = -12
           GO TO 180
    30     INEW = IW1(IOLD,1)
@@ -599,7 +599,7 @@ c perform some validity checks on iold and jold.
 c are we in a valid block and is it diagonal or off-diagonal?
           IF (IW1(INEW,3)-IW1(JNEW,3)) 40, 60, 50
    40     IFLAG = -13
-          IF (LP.NE.0) WRITE (LP,99998) IOLD, JOLD
+c$$$          IF (LP.NE.0) WRITE (LP,99998) IOLD, JOLD
           GO TO 180
    50     J1 = IW(INEW,1)
           J2 = J1 + LENOFF(INEW) - 1
@@ -621,14 +621,14 @@ c binary search of ordered list  .. element in l part of row.
    90       J1 = MIDPT
   100     CONTINUE
           IFLAG = -13
-          IF (LP.NE.0) WRITE (LP,99997) IOLD, JOLD
+c$$$          IF (LP.NE.0) WRITE (LP,99997) IOLD, JOLD
           GO TO 180
 c linear search ... element in l part of row or off-diagonal blocks.
   110     DO 120 MIDPT=J1,J2
             IF (IABS(ICN(MIDPT)+0).EQ.JNEW) GO TO 130
   120     CONTINUE
           IFLAG = -13
-          IF (LP.NE.0) WRITE (LP,99997) IOLD, JOLD
+c$$$          IF (LP.NE.0) WRITE (LP,99997) IOLD, JOLD
           GO TO 180
 c equivalent element of icn is in position midpt.
   130     IF (ICN(MIDPT).LT.0) GO TO 160
@@ -926,7 +926,7 @@ c block is structurally singular.
         IF (.NOT.ABORT1) GO TO 80
         IDISP(2) = IACTIV
         IFLAG = -1
-        IF (LP.NE.0) WRITE (LP,99999)
+c$$$        IF (LP.NE.0) WRITE (LP,99999)
 c     return
         GO TO 1120
    50   SCALE = DABS(A(IACTIV))
@@ -940,7 +940,7 @@ c     return
         IF (.NOT.ABORT2) GO TO 70
         IDISP(2) = IACTIV
         IFLAG = -2
-        IF (LP.NE.0) WRITE (LP,99998)
+c$$$        IF (LP.NE.0) WRITE (LP,99998)
         GO TO 1120
    70   A(IBEG) = A(IACTIV)
         ICN(IBEG) = ICN(IACTIV)
@@ -1159,7 +1159,7 @@ c matrix is structurally singular.
           IF (.NOT.ABORT1) GO TO 350
           IDISP(2) = IACTIV
           IFLAG = -1
-          IF (LP.NE.0) WRITE (LP,99999)
+c$$$          IF (LP.NE.0) WRITE (LP,99999)
           GO TO 1120
   350     K = PIVOT - 1
           DO 390 I=ISTART,ILAST
@@ -1175,7 +1175,7 @@ c check now to see if ma30d/dd has created enough available space.
 c create more space by destroying previously created lu factors.
             MOREI = MOREI + IBEG - IDISP(1)
             IBEG = IDISP(1)
-            IF (LP.NE.0) WRITE (LP,99997)
+c$$$            IF (LP.NE.0) WRITE (LP,99997)
             IFLAG = -5
             IF (ABORT3) GO TO 1090
   360       J1 = IPTR(I)
@@ -1210,7 +1210,7 @@ c numerical singularity is recorded here.
           IF (.NOT.ABORT2) GO TO 430
           IDISP(2) = IACTIV
           IFLAG = -2
-          IF (LP.NE.0) WRITE (LP,99998)
+c$$$          IF (LP.NE.0) WRITE (LP,99998)
           GO TO 1120
   430     OLDPIV = IPTR(IPIV) + LENRL(IPIV)
           OLDEND = IPTR(IPIV) + LENR(IPIV) - 1
@@ -1286,7 +1286,7 @@ c check now to see if ma30d/dd has created enough available space.
 c create more space by destroying previously created lu factors.
           MOREI = MOREI + IBEG - IDISP(1)
           IBEG = IDISP(1)
-          IF (LP.NE.0) WRITE (LP,99997)
+c$$$          IF (LP.NE.0) WRITE (LP,99997)
           IFLAG = -5
           IF (ABORT3) GO TO 1090
           IF (IACTIV-IBEG.GE.LENR(IPIV)) GO TO 500
@@ -1498,7 +1498,7 @@ c jump if there is space immediately available for the shifted row.
 c no space available so try to create some by throwing away previous
 c     lu decomposition.
             MOREI = MOREI + IBEG - IDISP(1) - LENPIV - 1
-            IF (LP.NE.0) WRITE (LP,99997)
+c$$$            IF (LP.NE.0) WRITE (LP,99997)
             IFLAG = -5
             IF (ABORT3) GO TO 1090
 c keep record of current pivot row.
@@ -1739,17 +1739,17 @@ c
 c   ***    error returns    ***
  1090 IDISP(2) = IACTIV
       IF (LP.EQ.0) GO TO 1120
-      WRITE (LP,99996)
+c$$$      WRITE (LP,99996)
       GO TO 1110
  1100 IF (IFLAG.EQ.-5) IFLAG = -6
       IF (IFLAG.NE.-6) IFLAG = -3
       IDISP(2) = IACTIV
       IF (LP.EQ.0) GO TO 1120
-      IF (IFLAG.EQ.-3) WRITE (LP,99995)
-      IF (IFLAG.EQ.-6) WRITE (LP,99994)
+c$$$      IF (IFLAG.EQ.-3) WRITE (LP,99995)
+c$$$      IF (IFLAG.EQ.-6) WRITE (LP,99994)
  1110 PIVOT = PIVOT - ISTART + 1
-      WRITE (LP,99993) PIVOT, NBLOCK, ISTART, ILAST
-      IF (PIVOT.EQ.0) WRITE (LP,99992) MINIRN
+c$$$      WRITE (LP,99993) PIVOT, NBLOCK, ISTART, ILAST
+c$$$      IF (PIVOT.EQ.0) WRITE (LP,99992) MINIRN
 c
 c
  1120 RETURN
@@ -1929,8 +1929,8 @@ c   ****    end of main loop    ****
 c
       GO TO 180
 c   ***   error return   ***
-  170 IF (LP.NE.0) WRITE (LP,99999) I
-      IFLAG = -I
+c$$$  170 IF (LP.NE.0) WRITE (LP,99999) I
+  170 IFLAG = -I
 c
   180 RETURN
 99999 FORMAT (54H ERROR RETURN FROM MA30B/BD SINGULARITY DETECTED IN RO,
@@ -2682,7 +2682,7 @@ c find row permutation ip to make diagonal zero-free.
       CALL MC21A(N,ICN,LICN,IW1,LENR,IP,NUMNZ,IW)
 c
 c possible error return for structurally singular matrices.
-      IF (NUMNZ.NE.N.AND.ABORT) GO TO 170
+      IF (NUMNZ.NE.N.AND.ABORT) GO TO 180
 c
 c iw1(.,2) and lenr are permutations of iw1(.,1) and lenr/lenoff
 c     suitable for entry
@@ -2713,14 +2713,14 @@ c idisp(2) is the position in a/icn of the first element in the
 c     diagonal blocks.
       IDISP(2)=LICN-NZ+1
       LARGE=N
-      IF (NZ.EQ.LICN) GO TO 230
+      IF (NZ.EQ.LICN) GO TO 220
       DO 50 K=1,NZ
       J=NZ-K+1
       JJ=LICN-K+1
       A(JJ)=A(J)
    50 ICN(JJ)=ICN(J)
 c 230 = return
-      GO TO 230
+      GO TO 220
 c
 c data structure reordered.
 c
@@ -2769,7 +2769,7 @@ c moves separated off-diagonal elements and untreated rows to
 c     front of storage.
       JNPOS=IBEG
       ILEND=IDISP(1)-1
-      IF (ILEND.LT.IBEG) GO TO 190
+      IF (ILEND.LT.IBEG) GO TO 200
       DO 90 J=IBEG,ILEND
       IF (ICN(J).EQ.0) GO TO 90
       ICN(JNPOS)=ICN(J)
@@ -2777,7 +2777,7 @@ c     front of storage.
       JNPOS=JNPOS+1
    90 CONTINUE
       IDISP(1)=JNPOS
-      IF (IEND-JNPOS.LT.LENOFF(IOLD)) GO TO 190
+      IF (IEND-JNPOS.LT.LENOFF(IOLD)) GO TO 200
       IBEG=LICN+1
 c reset pointers to the beginning of the rows.
       DO 100 I=2,N
@@ -2820,7 +2820,7 @@ c idisp(2) is position of first element in diagonal blocks.
 c
 c this compress is used to move all off-diagonal elements to the
 c     front of the storage.
-      IF (IBEG.GT.LICN) GO TO 230
+      IF (IBEG.GT.LICN) GO TO 220
       JNPOS=IBEG
       ILEND=IDISP(1)-1
       DO 160 J=IBEG,ILEND
@@ -2831,21 +2831,21 @@ c     front of the storage.
   160 CONTINUE
 c idisp(1) is first position after last element of off-diagonal blocks.
       IDISP(1)=JNPOS
-      GO TO 230
+      GO TO 220
 c
 c
 c error return
-  170 IF (LP.NE.0) WRITE(LP,180) NUMNZ
-  180 FORMAT(33X,41H MATRIX IS STRUCTURALLY SINGULAR, RANK = ,I6)
-      IDISP(1)=-1
-      GO TO 210
-  190 IF (LP.NE.0) WRITE(LP,200) N
-  200 FORMAT(33X,33H LICN NOT BIG ENOUGH INCREASE BY ,I6)
-      IDISP(1)=-2
-  210 IF (LP.NE.0) WRITE(LP,220)
-  220 FORMAT(33H+ERROR RETURN FROM MC23AD BECAUSE)
-c
-  230 RETURN
+c$$$  170 IF (LP.NE.0) WRITE(LP,180) NUMNZ
+c$$$  180 FORMAT(33X,41H MATRIX IS STRUCTURALLY SINGULAR, RANK = ,I6)
+  180 IDISP(1)=-1
+      GO TO 220
+c$$$  190 IF (LP.NE.0) WRITE(LP,200) N
+c$$$  200 FORMAT(33X,33H LICN NOT BIG ENOUGH INCREASE BY ,I6)
+  200 IDISP(1)=-2
+c$$$  210 IF (LP.NE.0) WRITE(LP,220)
+c$$$  220 FORMAT(33H+ERROR RETURN FROM MC23AD BECAUSE)
+c
+  220 RETURN
       END
 c######date   01 jan 1984     copyright ukaea, harwell.
 c######alias mc24ad
diff --git a/externals/hairer/radau.f b/externals/hairer/radau.f
index a8cbe8d..f6a40fd 100644
--- a/externals/hairer/radau.f
+++ b/externals/hairer/radau.f
@@ -427,7 +427,7 @@ C -------- NUMBER MAXIMAL AND MINIMAL OF STAGES  NS
       ELSE
          NSUS=IWORK(13)
          IF (NSUS.LE.0.OR.NS.GE.8.OR.NS.EQ.2.OR.NS.EQ.4.OR.NS.EQ.6) THEN
-            WRITE(6,*)' WRONG INPUT IWORK(13)=',IWORK(13)
+c$$$            WRITE(6,*)' WRONG INPUT IWORK(13)=',IWORK(13)
             ARRET=.TRUE.
          END IF
       END IF
@@ -437,7 +437,7 @@ C -------- NMAX , THE MAXIMAL NUMBER OF STEPS -----
       ELSE
          NMAX=IWORK(2)
          IF (NMAX.LE.0) THEN
-            WRITE(6,*)' WRONG INPUT IWORK(2)=',IWORK(2)
+c$$$            WRITE(6,*)' WRONG INPUT IWORK(2)=',IWORK(2)
             ARRET=.TRUE.
          END IF
       END IF
@@ -447,7 +447,7 @@ C -------- NIT    MAXIMAL NUMBER OF NEWTON ITERATIONS
       ELSE
          NIT=IWORK(3)
          IF (NIT.LE.0.OR.NIT.GT.50) THEN
-            WRITE(6,*)' CURIOUS INPUT IWORK(3)=',IWORK(3)
+c$$$            WRITE(6,*)' CURIOUS INPUT IWORK(3)=',IWORK(3)
             ARRET=.TRUE.
          END IF
       END IF
@@ -463,7 +463,7 @@ C -------- PARAMETER FOR DIFFERENTIAL-ALGEBRAIC COMPONENTS
       NIND3=IWORK(7)
       IF (NIND1.EQ.0) NIND1=N
       IF (NIND1+NIND2+NIND3.NE.N) THEN
-       WRITE(6,*)' CURIOUS INPUT FOR IWORK(5,6,7)=',NIND1,NIND2,NIND3
+c$$$       WRITE(6,*)' CURIOUS INPUT FOR IWORK(5,6,7)=',NIND1,NIND2,NIND3
        ARRET=.TRUE.
       END IF
 C -------- PRED   STEP SIZE CONTROL
@@ -479,7 +479,7 @@ C -------- PARAMETER FOR SECOND ORDER EQUATIONS
       IF (M1.EQ.0) M2=N
       IF (M2.EQ.0) M2=M1
       IF (M1.LT.0.OR.M2.LT.0.OR.M1+M2.GT.N) THEN
-       WRITE(6,*)' CURIOUS INPUT FOR IWORK(9,10)=',M1,M2
+c$$$       WRITE(6,*)' CURIOUS INPUT FOR IWORK(9,10)=',M1,M2
        ARRET=.TRUE.
       END IF
 C -------- UROUND   SMALLEST NUMBER SATISFYING 1.0D0+UROUND>1.0D0
@@ -488,20 +488,20 @@ C -------- UROUND   SMALLEST NUMBER SATISFYING 1.0D0+UROUND>1.0D0
       ELSE
          UROUND=WORK(1)
          IF (UROUND.LE.1.0D-19.OR.UROUND.GE.1.0D0) THEN
-            WRITE(6,*)' COEFFICIENTS HAVE 20 DIGITS, UROUND=',WORK(1)
+c$$$            WRITE(6,*)' COEFFICIENTS HAVE 20 DIGITS, UROUND=',WORK(1)
             ARRET=.TRUE.
          END IF
       END IF
 C --------- CHECK IF TOLERANCES ARE O.K.
       IF (ITOL.EQ.0) THEN
           IF (ATOL(1).LE.0.D0.OR.RTOL(1).LE.10.D0*UROUND) THEN
-              WRITE (6,*) ' TOLERANCES ARE TOO SMALL'
+c$$$              WRITE (6,*) ' TOLERANCES ARE TOO SMALL'
               ARRET=.TRUE.
           END IF
       ELSE
           DO I=1,N
           IF (ATOL(I).LE.0.D0.OR.RTOL(I).LE.10.D0*UROUND) THEN
-              WRITE (6,*) ' TOLERANCES(',I,') ARE TOO SMALL'
+c$$$              WRITE (6,*) ' TOLERANCES(',I,') ARE TOO SMALL'
               ARRET=.TRUE.
           END IF
           END DO
@@ -512,7 +512,7 @@ C --------- SAFE     SAFETY FACTOR IN STEP SIZE PREDICTION
       ELSE
          SAFE=WORK(2)
          IF (SAFE.LE.0.001D0.OR.SAFE.GE.1.0D0) THEN
-            WRITE(6,*)' CURIOUS INPUT FOR WORK(2)=',WORK(2)
+c$$$            WRITE(6,*)' CURIOUS INPUT FOR WORK(2)=',WORK(2)
             ARRET=.TRUE.
          END IF
       END IF
@@ -522,7 +522,7 @@ C ------ THET     DECIDES WHETHER THE JACOBIAN SHOULD BE RECOMPUTED;
       ELSE
          THET=WORK(3)
          IF (THET.GE.1.0D0) THEN
-            WRITE(6,*)' CURIOUS INPUT FOR WORK(3)=',WORK(3)
+c$$$            WRITE(6,*)' CURIOUS INPUT FOR WORK(3)=',WORK(3)
             ARRET=.TRUE.
          END IF
       END IF
@@ -538,7 +538,7 @@ C --- QUOT1 AND QUOT2: IF QUOT1 < HNEW/HOLD < QUOT2, STEP SIZE = CONST.
          QUOT2=WORK(6)
       END IF
       IF (QUOT1.GT.1.0D0.OR.QUOT2.LT.1.0D0) THEN
-         WRITE(6,*)' CURIOUS INPUT FOR WORK(5,6)=',QUOT1,QUOT2
+c$$$         WRITE(6,*)' CURIOUS INPUT FOR WORK(5,6)=',QUOT1,QUOT2
          ARRET=.TRUE.
       END IF
 C -------- MAXIMAL STEP SIZE
@@ -559,7 +559,7 @@ C -------  FACL,FACR     PARAMETERS FOR STEP SIZE SELECTION
          FACR=1.D0/WORK(9)
       END IF
       IF (FACL.LT.1.0D0.OR.FACR.GT.1.0D0) THEN
-            WRITE(6,*)' CURIOUS INPUT WORK(8,9)=',WORK(8),WORK(9)
+c$$$            WRITE(6,*)' CURIOUS INPUT WORK(8,9)=',WORK(8),WORK(9)
             ARRET=.TRUE.
          END IF
 C -------- PARAMETERS FOR ORDER SELECTION STRATEGY
@@ -615,8 +615,8 @@ C -- MASS MATRIX
           END IF
 C ------ BANDWITH OF "MAS" NOT SMALLER THAN BANDWITH OF "JAC"
           IF (MLMAS.GT.MLJAC.OR.MUMAS.GT.MUJAC) THEN
-             WRITE (6,*) 'BANDWITH OF "MAS" NOT SMALLER THAN BANDWITH OF
-     & "JAC"'
+c$$$             WRITE (6,*) 'BANDWITH OF "MAS" NOT SMALLER THAN BANDWITH OF
+c$$$     & "JAC"'
             ARRET=.TRUE.
           END IF
       ELSE
@@ -631,8 +631,8 @@ C ------ BANDWITH OF "MAS" NOT SMALLER THAN BANDWITH OF "JAC"
       LDMAS2=MAX(1,LDMAS)
 C ------ HESSENBERG OPTION ONLY FOR EXPLICIT EQU. WITH FULL JACOBIAN
       IF ((IMPLCT.OR.JBAND).AND.IJOB.EQ.7) THEN
-         WRITE(6,*)' HESSENBERG OPTION ONLY FOR EXPLICIT EQUATIONS WITH
-     &FULL JACOBIAN'
+c$$$         WRITE(6,*)' HESSENBERG OPTION ONLY FOR EXPLICIT EQUATIONS WITH
+c$$$     &FULL JACOBIAN'
          ARRET=.TRUE.
       END IF
 C ------- PREPARE THE ENTRY-POINTS FOR THE ARRAYS IN WORK -----
@@ -651,7 +651,7 @@ C ------- PREPARE THE ENTRY-POINTS FOR THE ARRAYS IN WORK -----
 C ------ TOTAL STORAGE REQUIREMENT -----------
       ISTORE=IEE+NMEE*LDE1-1
       IF(ISTORE.GT.LWORK)THEN
-         WRITE(6,*)' INSUFFICIENT STORAGE FOR WORK, MIN. LWORK=',ISTORE
+c$$$         WRITE(6,*)' INSUFFICIENT STORAGE FOR WORK, MIN. LWORK=',ISTORE
          ARRET=.TRUE.
       END IF
 C ------- ENTRY POINTS FOR INTEGER WORKSPACE -----
@@ -661,7 +661,7 @@ C ------- ENTRY POINTS FOR INTEGER WORKSPACE -----
 C --------- TOTAL REQUIREMENT ---------------
       ISTORE=IEIPH+NM1-1
       IF (ISTORE.GT.LIWORK) THEN
-         WRITE(6,*)' INSUFF. STORAGE FOR IWORK, MIN. LIWORK=',ISTORE
+c$$$         WRITE(6,*)' INSUFF. STORAGE FOR IWORK, MIN. LIWORK=',ISTORE
          ARRET=.TRUE.
       END IF
 C ------ WHEN A FAIL HAS OCCURED, WE RETURN WITH IDID=-1
@@ -1668,24 +1668,24 @@ C --- UNEXPECTED STEP-REJECTION
       GOTO 10
 C --- FAIL EXIT
  176  CONTINUE
-      WRITE(6,979)X
-      WRITE(6,*) ' MATRIX IS REPEATEDLY SINGULAR, IER=',IER
+c$$$      WRITE(6,979)X
+c$$$      WRITE(6,*) ' MATRIX IS REPEATEDLY SINGULAR, IER=',IER
       IDID=-4
       RETURN
  177  CONTINUE
-      WRITE(6,979)X
-      WRITE(6,*) ' STEP SIZE T0O SMALL, H=',H
+c$$$      WRITE(6,979)X
+c$$$      WRITE(6,*) ' STEP SIZE T0O SMALL, H=',H
       IDID=-3
       RETURN
  178  CONTINUE
-      WRITE(6,979)X
-      WRITE(6,*) ' MORE THAN NMAX =',NMAX,'STEPS ARE NEEDED'
+c$$$      WRITE(6,979)X
+c$$$      WRITE(6,*) ' MORE THAN NMAX =',NMAX,'STEPS ARE NEEDED'
       IDID=-2
       RETURN
 C --- EXIT CAUSED BY SOLOUT
  179  CONTINUE
-      WRITE(6,979)X
- 979  FORMAT(' EXIT OF RADAU AT X=',E18.4)
+c$$$      WRITE(6,979)X
+c$$$ 979  FORMAT(' EXIT OF RADAU AT X=',E18.4)
       IDID=2
       RETURN
       END
diff --git a/externals/hairer/radau5.f b/externals/hairer/radau5.f
index 9261842..1da979c 100644
--- a/externals/hairer/radau5.f
+++ b/externals/hairer/radau5.f
@@ -392,7 +392,7 @@ C -------- UROUND   SMALLEST NUMBER SATISFYING 1.0D0+UROUND>1.0D0
       ELSE
          UROUND=WORK(1)
          IF (UROUND.LE.1.0D-19.OR.UROUND.GE.1.0D0) THEN
-            WRITE(6,*)' COEFFICIENTS HAVE 20 DIGITS, UROUND=',WORK(1)
+c$$$            WRITE(6,*)' COEFFICIENTS HAVE 20 DIGITS, UROUND=',WORK(1)
             ARRET=.TRUE.
          END IF
       END IF
@@ -400,7 +400,7 @@ C -------- CHECK AND CHANGE THE TOLERANCES
       EXPM=2.0D0/3.0D0
       IF (ITOL.EQ.0) THEN
           IF (ATOL(1).LE.0.D0.OR.RTOL(1).LE.10.D0*UROUND) THEN
-              WRITE (6,*) ' TOLERANCES ARE TOO SMALL'
+c$$$              WRITE (6,*) ' TOLERANCES ARE TOO SMALL'
               ARRET=.TRUE.
           ELSE
               QUOT=ATOL(1)/RTOL(1)
@@ -410,7 +410,7 @@ C -------- CHECK AND CHANGE THE TOLERANCES
       ELSE
           DO I=1,N
           IF (ATOL(I).LE.0.D0.OR.RTOL(I).LE.10.D0*UROUND) THEN
-              WRITE (6,*) ' TOLERANCES(',I,') ARE TOO SMALL'
+c$$$              WRITE (6,*) ' TOLERANCES(',I,') ARE TOO SMALL'
               ARRET=.TRUE.
           ELSE
               QUOT=ATOL(I)/RTOL(I)
@@ -425,7 +425,7 @@ C -------- NMAX , THE MAXIMAL NUMBER OF STEPS -----
       ELSE
          NMAX=IWORK(2)
          IF (NMAX.LE.0) THEN
-            WRITE(6,*)' WRONG INPUT IWORK(2)=',IWORK(2)
+c$$$            WRITE(6,*)' WRONG INPUT IWORK(2)=',IWORK(2)
             ARRET=.TRUE.
          END IF
       END IF
@@ -435,7 +435,7 @@ C -------- NIT    MAXIMAL NUMBER OF NEWTON ITERATIONS
       ELSE
          NIT=IWORK(3)
          IF (NIT.LE.0) THEN
-            WRITE(6,*)' CURIOUS INPUT IWORK(3)=',IWORK(3)
+c$$$            WRITE(6,*)' CURIOUS INPUT IWORK(3)=',IWORK(3)
             ARRET=.TRUE.
          END IF
       END IF
@@ -451,7 +451,7 @@ C -------- PARAMETER FOR DIFFERENTIAL-ALGEBRAIC COMPONENTS
       NIND3=IWORK(7)
       IF (NIND1.EQ.0) NIND1=N
       IF (NIND1+NIND2+NIND3.NE.N) THEN
-       WRITE(6,*)' CURIOUS INPUT FOR IWORK(5,6,7)=',NIND1,NIND2,NIND3
+c$$$       WRITE(6,*)' CURIOUS INPUT FOR IWORK(5,6,7)=',NIND1,NIND2,NIND3
        ARRET=.TRUE.
       END IF
 C -------- PRED   STEP SIZE CONTROL
@@ -467,7 +467,7 @@ C -------- PARAMETER FOR SECOND ORDER EQUATIONS
       IF (M1.EQ.0) M2=N
       IF (M2.EQ.0) M2=M1
       IF (M1.LT.0.OR.M2.LT.0.OR.M1+M2.GT.N) THEN
-       WRITE(6,*)' CURIOUS INPUT FOR IWORK(9,10)=',M1,M2
+c$$$       WRITE(6,*)' CURIOUS INPUT FOR IWORK(9,10)=',M1,M2
        ARRET=.TRUE.
       END IF
 C --------- SAFE     SAFETY FACTOR IN STEP SIZE PREDICTION
@@ -476,7 +476,7 @@ C --------- SAFE     SAFETY FACTOR IN STEP SIZE PREDICTION
       ELSE
          SAFE=WORK(2)
          IF (SAFE.LE.0.001D0.OR.SAFE.GE.1.0D0) THEN
-            WRITE(6,*)' CURIOUS INPUT FOR WORK(2)=',WORK(2)
+c$$$            WRITE(6,*)' CURIOUS INPUT FOR WORK(2)=',WORK(2)
             ARRET=.TRUE.
          END IF
       END IF
@@ -486,7 +486,7 @@ C ------ THET     DECIDES WHETHER THE JACOBIAN SHOULD BE RECOMPUTED;
       ELSE
          THET=WORK(3)
          IF (THET.GE.1.0D0) THEN
-            WRITE(6,*)' CURIOUS INPUT FOR WORK(3)=',WORK(3)
+c$$$            WRITE(6,*)' CURIOUS INPUT FOR WORK(3)=',WORK(3)
             ARRET=.TRUE.
          END IF
       END IF
@@ -497,7 +497,7 @@ C --- FNEWT   STOPPING CRITERION FOR NEWTON'S METHOD, USUALLY CHOSEN <1.
       ELSE
          FNEWT=WORK(4)
          IF (FNEWT.LE.UROUND/TOLST) THEN
-            WRITE(6,*)' CURIOUS INPUT FOR WORK(4)=',WORK(4)
+c$$$            WRITE(6,*)' CURIOUS INPUT FOR WORK(4)=',WORK(4)
             ARRET=.TRUE.
          END IF
       END IF
@@ -513,7 +513,7 @@ C --- QUOT1 AND QUOT2: IF QUOT1 < HNEW/HOLD < QUOT2, STEP SIZE = CONST.
          QUOT2=WORK(6)
       END IF
       IF (QUOT1.GT.1.0D0.OR.QUOT2.LT.1.0D0) THEN
-         WRITE(6,*)' CURIOUS INPUT FOR WORK(5,6)=',QUOT1,QUOT2
+c$$$         WRITE(6,*)' CURIOUS INPUT FOR WORK(5,6)=',QUOT1,QUOT2
          ARRET=.TRUE.
       END IF
 C -------- MAXIMAL STEP SIZE
@@ -534,7 +534,7 @@ C -------  FACL,FACR     PARAMETERS FOR STEP SIZE SELECTION
          FACR=1.D0/WORK(9)
       END IF
       IF (FACL.LT.1.0D0.OR.FACR.GT.1.0D0) THEN
-            WRITE(6,*)' CURIOUS INPUT WORK(8,9)=',WORK(8),WORK(9)
+c$$$            WRITE(6,*)' CURIOUS INPUT WORK(8,9)=',WORK(8),WORK(9)
             ARRET=.TRUE.
          END IF
 C *** *** *** *** *** *** *** *** *** *** *** *** ***
@@ -570,8 +570,8 @@ C -- MASS MATRIX
           END IF
 C ------ BANDWITH OF "MAS" NOT SMALLER THAN BANDWITH OF "JAC"
           IF (MLMAS.GT.MLJAC.OR.MUMAS.GT.MUJAC) THEN
-             WRITE (6,*) 'BANDWITH OF "MAS" NOT SMALLER THAN BANDWITH OF
-     & "JAC"'
+c$$$             WRITE (6,*) 'BANDWITH OF "MAS" NOT SMALLER THAN BANDWITH OF
+c$$$     & "JAC"'
             ARRET=.TRUE.
           END IF
       ELSE
@@ -586,8 +586,8 @@ C ------ BANDWITH OF "MAS" NOT SMALLER THAN BANDWITH OF "JAC"
       LDMAS2=MAX(1,LDMAS)
 C ------ HESSENBERG OPTION ONLY FOR EXPLICIT EQU. WITH FULL JACOBIAN
       IF ((IMPLCT.OR.JBAND).AND.IJOB.EQ.7) THEN
-         WRITE(6,*)' HESSENBERG OPTION ONLY FOR EXPLICIT EQUATIONS WITH
-     &FULL JACOBIAN'
+c$$$         WRITE(6,*)' HESSENBERG OPTION ONLY FOR EXPLICIT EQUATIONS WITH
+c$$$     &FULL JACOBIAN'
          ARRET=.TRUE.
       END IF
 C ------- PREPARE THE ENTRY-POINTS FOR THE ARRAYS IN WORK -----
@@ -608,7 +608,7 @@ C ------- PREPARE THE ENTRY-POINTS FOR THE ARRAYS IN WORK -----
 C ------ TOTAL STORAGE REQUIREMENT -----------
       ISTORE=IEE2I+NM1*LDE1-1
       IF(ISTORE.GT.LWORK)THEN
-         WRITE(6,*)' INSUFFICIENT STORAGE FOR WORK, MIN. LWORK=',ISTORE
+c$$$         WRITE(6,*)' INSUFFICIENT STORAGE FOR WORK, MIN. LWORK=',ISTORE
          ARRET=.TRUE.
       END IF
 C ------- ENTRY POINTS FOR INTEGER WORKSPACE -----
@@ -618,7 +618,7 @@ C ------- ENTRY POINTS FOR INTEGER WORKSPACE -----
 C --------- TOTAL REQUIREMENT ---------------
       ISTORE=IEIPH+NM1-1
       IF (ISTORE.GT.LIWORK) THEN
-         WRITE(6,*)' INSUFF. STORAGE FOR IWORK, MIN. LIWORK=',ISTORE
+c$$$         WRITE(6,*)' INSUFF. STORAGE FOR IWORK, MIN. LIWORK=',ISTORE
          ARRET=.TRUE.
       END IF
 C ------ WHEN A FAIL HAS OCCURED, WE RETURN WITH IDID=-1
@@ -1101,24 +1101,24 @@ C --- UNEXPECTED STEP-REJECTION
       GOTO 10
 C --- FAIL EXIT
  176  CONTINUE
-      WRITE(6,979)X
-      WRITE(6,*) ' MATRIX IS REPEATEDLY SINGULAR, IER=',IER
+c$$$      WRITE(6,979)X
+c$$$      WRITE(6,*) ' MATRIX IS REPEATEDLY SINGULAR, IER=',IER
       IDID=-4
       RETURN
  177  CONTINUE
-      WRITE(6,979)X
-      WRITE(6,*) ' STEP SIZE T0O SMALL, H=',H
+c$$$      WRITE(6,979)X
+c$$$      WRITE(6,*) ' STEP SIZE T0O SMALL, H=',H
       IDID=-3
       RETURN
  178  CONTINUE
-      WRITE(6,979)X
-      WRITE(6,*) ' MORE THAN NMAX =',NMAX,'STEPS ARE NEEDED'
+c$$$      WRITE(6,979)X
+c$$$      WRITE(6,*) ' MORE THAN NMAX =',NMAX,'STEPS ARE NEEDED'
       IDID=-2
       RETURN
 C --- EXIT CAUSED BY SOLOUT
  179  CONTINUE
-      WRITE(6,979)X
- 979  FORMAT(' EXIT OF RADAU5 AT X=',E18.4)
+c$$$      WRITE(6,979)X
+c$$$ 979  FORMAT(' EXIT OF RADAU5 AT X=',E18.4)
       IDID=2
       RETURN
       END
diff --git a/externals/hairer/radaup.f b/externals/hairer/radaup.f
index c802528..7ff5aac 100644
--- a/externals/hairer/radaup.f
+++ b/externals/hairer/radaup.f
@@ -397,7 +397,7 @@ C -------- NUMBER OF STAGES  NS
       ELSE
          NS=IWORK(11)
          IF (NS.LE.2.OR.NS.GE.8.OR.NS.EQ.4.OR.NS.EQ.6) THEN
-            WRITE(6,*)' CURIOUS INPUT IWORK(11)=',IWORK(11)
+c$$$            WRITE(6,*)' CURIOUS INPUT IWORK(11)=',IWORK(11)
             ARRET=.TRUE.
          END IF
       END IF
@@ -407,7 +407,7 @@ C -------- NMAX , THE MAXIMAL NUMBER OF STEPS -----
       ELSE
          NMAX=IWORK(2)
          IF (NMAX.LE.0) THEN
-            WRITE(6,*)' WRONG INPUT IWORK(2)=',IWORK(2)
+c$$$            WRITE(6,*)' WRONG INPUT IWORK(2)=',IWORK(2)
             ARRET=.TRUE.
          END IF
       END IF
@@ -417,7 +417,7 @@ C -------- NIT    MAXIMAL NUMBER OF NEWTON ITERATIONS
       ELSE
          NIT=IWORK(3)
          IF (NIT.LE.0) THEN
-            WRITE(6,*)' CURIOUS INPUT IWORK(3)=',IWORK(3)
+c$$$            WRITE(6,*)' CURIOUS INPUT IWORK(3)=',IWORK(3)
             ARRET=.TRUE.
          END IF
       END IF
@@ -433,7 +433,7 @@ C -------- PARAMETER FOR DIFFERENTIAL-ALGEBRAIC COMPONENTS
       NIND3=IWORK(7)
       IF (NIND1.EQ.0) NIND1=N
       IF (NIND1+NIND2+NIND3.NE.N) THEN
-       WRITE(6,*)' CURIOUS INPUT FOR IWORK(5,6,7)=',NIND1,NIND2,NIND3
+c$$$       WRITE(6,*)' CURIOUS INPUT FOR IWORK(5,6,7)=',NIND1,NIND2,NIND3
        ARRET=.TRUE.
       END IF
 C -------- PRED   STEP SIZE CONTROL
@@ -449,7 +449,7 @@ C -------- PARAMETER FOR SECOND ORDER EQUATIONS
       IF (M1.EQ.0) M2=N
       IF (M2.EQ.0) M2=M1
       IF (M1.LT.0.OR.M2.LT.0.OR.M1+M2.GT.N) THEN
-       WRITE(6,*)' CURIOUS INPUT FOR IWORK(9,10)=',M1,M2
+c$$$       WRITE(6,*)' CURIOUS INPUT FOR IWORK(9,10)=',M1,M2
        ARRET=.TRUE.
       END IF
 C -------- UROUND   SMALLEST NUMBER SATISFYING 1.0D0+UROUND>1.0D0
@@ -458,7 +458,7 @@ C -------- UROUND   SMALLEST NUMBER SATISFYING 1.0D0+UROUND>1.0D0
       ELSE
          UROUND=WORK(1)
          IF (UROUND.LE.1.0D-19.OR.UROUND.GE.1.0D0) THEN
-            WRITE(6,*)' COEFFICIENTS HAVE 20 DIGITS, UROUND=',WORK(1)
+c$$$            WRITE(6,*)' COEFFICIENTS HAVE 20 DIGITS, UROUND=',WORK(1)
             ARRET=.TRUE.
          END IF
       END IF
@@ -466,7 +466,7 @@ C --------- CHECK IF TOLERANCES ARE O.K.
       EXPM=(NS+1.0D0)/(2.0D0*NS)
       IF (ITOL.EQ.0) THEN
           IF (ATOL(1).LE.0.D0.OR.RTOL(1).LE.10.D0*UROUND) THEN
-              WRITE (6,*) ' TOLERANCES ARE TOO SMALL'
+c$$$              WRITE (6,*) ' TOLERANCES ARE TOO SMALL'
               ARRET=.TRUE.
           ELSE
               QUOT=ATOL(1)/RTOL(1)
@@ -476,7 +476,7 @@ C --------- CHECK IF TOLERANCES ARE O.K.
       ELSE
           DO I=1,N
           IF (ATOL(I).LE.0.D0.OR.RTOL(I).LE.10.D0*UROUND) THEN
-              WRITE (6,*) ' TOLERANCES(',I,') ARE TOO SMALL'
+c$$$              WRITE (6,*) ' TOLERANCES(',I,') ARE TOO SMALL'
               ARRET=.TRUE.
           ELSE
               QUOT=ATOL(I)/RTOL(I)
@@ -491,7 +491,7 @@ C --------- SAFE     SAFETY FACTOR IN STEP SIZE PREDICTION
       ELSE
          SAFE=WORK(2)
          IF (SAFE.LE.0.001D0.OR.SAFE.GE.1.0D0) THEN
-            WRITE(6,*)' CURIOUS INPUT FOR WORK(2)=',WORK(2)
+c$$$            WRITE(6,*)' CURIOUS INPUT FOR WORK(2)=',WORK(2)
             ARRET=.TRUE.
          END IF
       END IF
@@ -501,7 +501,7 @@ C ------ THET     DECIDES WHETHER THE JACOBIAN SHOULD BE RECOMPUTED;
       ELSE
          THET=WORK(3)
          IF (THET.GE.1.0D0) THEN
-            WRITE(6,*)' CURIOUS INPUT FOR WORK(3)=',WORK(3)
+c$$$            WRITE(6,*)' CURIOUS INPUT FOR WORK(3)=',WORK(3)
             ARRET=.TRUE.
          END IF
       END IF
@@ -513,7 +513,7 @@ C --- FNEWT   STOPPING CRIERION FOR NEWTON'S METHOD, USUALLY CHOSEN <1.
       ELSE
          FNEWT=WORK(4)
          IF (FNEWT.LE.UROUND/TOLST) THEN
-            WRITE(6,*)' CURIOUS INPUT FOR WORK(4)=',WORK(4)
+c$$$            WRITE(6,*)' CURIOUS INPUT FOR WORK(4)=',WORK(4)
             ARRET=.TRUE.
          END IF
       END IF
@@ -529,7 +529,7 @@ C --- QUOT1 AND QUOT2: IF QUOT1 < HNEW/HOLD < QUOT2, STEP SIZE = CONST.
          QUOT2=WORK(6)
       END IF
       IF (QUOT1.GT.1.0D0.OR.QUOT2.LT.1.0D0) THEN
-         WRITE(6,*)' CURIOUS INPUT FOR WORK(5,6)=',QUOT1,QUOT2
+c$$$         WRITE(6,*)' CURIOUS INPUT FOR WORK(5,6)=',QUOT1,QUOT2
          ARRET=.TRUE.
       END IF
 C -------- MAXIMAL STEP SIZE
@@ -550,7 +550,7 @@ C -------  FACL,FACR     PARAMETERS FOR STEP SIZE SELECTION
          FACR=1.D0/WORK(9)
       END IF
       IF (FACL.LT.1.0D0.OR.FACR.GT.1.0D0) THEN
-            WRITE(6,*)' CURIOUS INPUT WORK(8,9)=',WORK(8),WORK(9)
+c$$$            WRITE(6,*)' CURIOUS INPUT WORK(8,9)=',WORK(8),WORK(9)
             ARRET=.TRUE.
          END IF
 C *** *** *** *** *** *** *** *** *** *** *** *** ***
@@ -585,8 +585,8 @@ C -- MASS MATRIX
           END IF
 C ------ BANDWITH OF "MAS" NOT SMALLER THAN BANDWITH OF "JAC"
           IF (MLMAS.GT.MLJAC.OR.MUMAS.GT.MUJAC) THEN
-             WRITE (6,*) 'BANDWITH OF "MAS" NOT SMALLER THAN BANDWITH OF
-     & "JAC"'
+c$$$             WRITE (6,*) 'BANDWITH OF "MAS" NOT SMALLER THAN BANDWITH OF
+c$$$     & "JAC"'
             ARRET=.TRUE.
           END IF
       ELSE
@@ -601,8 +601,8 @@ C ------ BANDWITH OF "MAS" NOT SMALLER THAN BANDWITH OF "JAC"
       LDMAS2=MAX(1,LDMAS)
 C ------ HESSENBERG OPTION ONLY FOR EXPLICIT EQU. WITH FULL JACOBIAN
       IF ((IMPLCT.OR.JBAND).AND.IJOB.EQ.7) THEN
-         WRITE(6,*)' HESSENBERG OPTION ONLY FOR EXPLICIT EQUATIONS WITH
-     &FULL JACOBIAN'
+c$$$         WRITE(6,*)' HESSENBERG OPTION ONLY FOR EXPLICIT EQUATIONS WITH
+c$$$     &FULL JACOBIAN'
          ARRET=.TRUE.
       END IF
 C ------- PREPARE THE ENTRY-POINTS FOR THE ARRAYS IN WORK -----
@@ -621,7 +621,7 @@ C ------- PREPARE THE ENTRY-POINTS FOR THE ARRAYS IN WORK -----
 C ------ TOTAL STORAGE REQUIREMENT -----------
       ISTORE=IEE+NMEE*LDE1-1
       IF(ISTORE.GT.LWORK)THEN
-         WRITE(6,*)' INSUFFICIENT STORAGE FOR WORK, MIN. LWORK=',ISTORE
+c$$$         WRITE(6,*)' INSUFFICIENT STORAGE FOR WORK, MIN. LWORK=',ISTORE
          ARRET=.TRUE.
       END IF
 C ------- ENTRY POINTS FOR INTEGER WORKSPACE -----
@@ -631,7 +631,7 @@ C ------- ENTRY POINTS FOR INTEGER WORKSPACE -----
 C --------- TOTAL REQUIREMENT ---------------
       ISTORE=IEIPH+NM1-1
       IF (ISTORE.GT.LIWORK) THEN
-         WRITE(6,*)' INSUFF. STORAGE FOR IWORK, MIN. LIWORK=',ISTORE
+c$$$         WRITE(6,*)' INSUFF. STORAGE FOR IWORK, MIN. LIWORK=',ISTORE
          ARRET=.TRUE.
       END IF
 C ------ WHEN A FAIL HAS OCCURED, WE RETURN WITH IDID=-1
@@ -1102,24 +1102,24 @@ C --- UNEXPECTED STEP-REJECTION
       GOTO 10
 C --- FAIL EXIT
  176  CONTINUE
-      WRITE(6,979)X
-      WRITE(6,*) ' MATRIX IS REPEATEDLY SINGULAR, IER=',IER
+c$$$      WRITE(6,979)X
+c$$$      WRITE(6,*) ' MATRIX IS REPEATEDLY SINGULAR, IER=',IER
       IDID=-4
       RETURN
  177  CONTINUE
-      WRITE(6,979)X
-      WRITE(6,*) ' STEP SIZE T0O SMALL, H=',H
+c$$$      WRITE(6,979)X
+c$$$      WRITE(6,*) ' STEP SIZE T0O SMALL, H=',H
       IDID=-3
       RETURN
  178  CONTINUE
-      WRITE(6,979)X
-      WRITE(6,*) ' MORE THAN NMAX =',NMAX,'STEPS ARE NEEDED'
+c$$$      WRITE(6,979)X
+c$$$      WRITE(6,*) ' MORE THAN NMAX =',NMAX,'STEPS ARE NEEDED'
       IDID=-2
       RETURN
 C --- EXIT CAUSED BY SOLOUT
  179  CONTINUE
-      WRITE(6,979)X
- 979  FORMAT(' EXIT OF RADAU5 AT X=',E18.4)
+c$$$      WRITE(6,979)X
+c$$$ 979  FORMAT(' EXIT OF RADAU5 AT X=',E18.4)
       IDID=2
       RETURN
       END
diff --git a/externals/hairer/rodas.f b/externals/hairer/rodas.f
index 88dceff..ea3ffba 100644
--- a/externals/hairer/rodas.f
+++ b/externals/hairer/rodas.f
@@ -341,7 +341,7 @@ C -------- NMAX , THE MAXIMAL NUMBER OF STEPS -----
       ELSE
          NMAX=IWORK(1)
          IF(NMAX.LE.0)THEN
-            WRITE(6,*)' WRONG INPUT IWORK(1)=',IWORK(1)
+c$$$            WRITE(6,*)' WRONG INPUT IWORK(1)=',IWORK(1)
             ARRET=.TRUE.
          END IF
       END IF
@@ -351,7 +351,7 @@ C -------- METH   COEFFICIENTS OF THE METHOD
       ELSE
          METH=IWORK(2)
          IF(METH.LE.0.OR.METH.GE.4)THEN
-            WRITE(6,*)' CURIOUS INPUT IWORK(2)=',IWORK(2)
+c$$$            WRITE(6,*)' CURIOUS INPUT IWORK(2)=',IWORK(2)
             ARRET=.TRUE.
          END IF
       END IF
@@ -368,7 +368,7 @@ C -------- PARAMETER FOR SECOND ORDER EQUATIONS
       IF (M1.EQ.0) M2=N
       IF (M2.EQ.0) M2=M1
       IF (M1.LT.0.OR.M2.LT.0.OR.M1+M2.GT.N) THEN
-       WRITE(6,*)' CURIOUS INPUT FOR IWORK(9,10)=',M1,M2
+c$$$       WRITE(6,*)' CURIOUS INPUT FOR IWORK(9,10)=',M1,M2
        ARRET=.TRUE.
       END IF
 C -------- UROUND   SMALLEST NUMBER SATISFYING 1.D0+UROUND>1.D0
@@ -377,7 +377,7 @@ C -------- UROUND   SMALLEST NUMBER SATISFYING 1.D0+UROUND>1.D0
       ELSE
          UROUND=WORK(1)
          IF(UROUND.LT.1.D-16.OR.UROUND.GE.1.D0)THEN
-            WRITE(6,*)' COEFFICIENTS HAVE 16 DIGITS, UROUND=',WORK(1)
+c$$$            WRITE(6,*)' COEFFICIENTS HAVE 16 DIGITS, UROUND=',WORK(1)
             ARRET=.TRUE.
          END IF
       END IF
@@ -399,7 +399,7 @@ C -------  FAC1,FAC2     PARAMETERS FOR STEP SIZE SELECTION
          FAC2=1.D0/WORK(4)
       END IF
       IF (FAC1.LT.1.0D0.OR.FAC2.GT.1.0D0) THEN
-            WRITE(6,*)' CURIOUS INPUT WORK(3,4)=',WORK(3),WORK(4)
+c$$$            WRITE(6,*)' CURIOUS INPUT WORK(3,4)=',WORK(3),WORK(4)
             ARRET=.TRUE.
          END IF
 C --------- SAFE     SAFETY FACTOR IN STEP SIZE PREDICTION
@@ -408,20 +408,20 @@ C --------- SAFE     SAFETY FACTOR IN STEP SIZE PREDICTION
       ELSE
          SAFE=WORK(5)
          IF (SAFE.LE.0.001D0.OR.SAFE.GE.1.0D0) THEN
-            WRITE(6,*)' CURIOUS INPUT FOR WORK(5)=',WORK(5)
+c$$$            WRITE(6,*)' CURIOUS INPUT FOR WORK(5)=',WORK(5)
             ARRET=.TRUE.
          END IF
       END IF
 C --------- CHECK IF TOLERANCES ARE O.K.
       IF (ITOL.EQ.0) THEN
           IF (ATOL(1).LE.0.D0.OR.RTOL(1).LE.10.D0*UROUND) THEN
-              WRITE (6,*) ' TOLERANCES ARE TOO SMALL'
+c$$$              WRITE (6,*) ' TOLERANCES ARE TOO SMALL'
               ARRET=.TRUE.
           END IF
       ELSE
           DO I=1,N
              IF (ATOL(I).LE.0.D0.OR.RTOL(I).LE.10.D0*UROUND) THEN
-                WRITE (6,*) ' TOLERANCES(',I,') ARE TOO SMALL'
+c$$$                WRITE (6,*) ' TOLERANCES(',I,') ARE TOO SMALL'
                 ARRET=.TRUE.
              END IF
           END DO
@@ -459,8 +459,8 @@ C -- MASS MATRIX
           END IF
 C ------ BANDWITH OF "MAS" NOT LARGER THAN BANDWITH OF "JAC"
           IF (MLMAS.GT.MLJAC.OR.MUMAS.GT.MUJAC) THEN
-              WRITE (6,*) 'BANDWITH OF "MAS" NOT LARGER THAN BANDWITH OF
-     & "JAC"'
+c$$$              WRITE (6,*) 'BANDWITH OF "MAS" NOT LARGER THAN BANDWITH OF
+c$$$     & "JAC"'
             ARRET=.TRUE.
           END IF
       ELSE
@@ -490,14 +490,14 @@ C ------- PREPARE THE ENTRY-POINTS FOR THE ARRAYS IN WORK -----
 C ------ TOTAL STORAGE REQUIREMENT -----------
       ISTORE=IEE+NM1*LDE-1
       IF(ISTORE.GT.LWORK)THEN
-         WRITE(6,*)' INSUFFICIENT STORAGE FOR WORK, MIN. LWORK=',ISTORE
+c$$$         WRITE(6,*)' INSUFFICIENT STORAGE FOR WORK, MIN. LWORK=',ISTORE
          ARRET=.TRUE.
       END IF
 C ------- ENTRY POINTS FOR INTEGER WORKSPACE -----
       IEIP=21
       ISTORE=IEIP+NM1-1
       IF(ISTORE.GT.LIWORK)THEN
-         WRITE(6,*)' INSUFF. STORAGE FOR IWORK, MIN. LIWORK=',ISTORE
+c$$$         WRITE(6,*)' INSUFF. STORAGE FOR IWORK, MIN. LIWORK=',ISTORE
          ARRET=.TRUE.
       END IF
 C ------ WHEN A FAIL HAS OCCURED, WE RETURN WITH IDID=-1
@@ -849,24 +849,24 @@ C --- SINGULAR MATRIX
       GOTO 2
 C --- FAIL EXIT
  176  CONTINUE
-      WRITE(6,979)X
-      WRITE(6,*) ' MATRIX IS REPEATEDLY SINGULAR, IER=',IER
+c$$$      WRITE(6,979)X
+c$$$      WRITE(6,*) ' MATRIX IS REPEATEDLY SINGULAR, IER=',IER
       IDID=-4
       RETURN
  177  CONTINUE
-      WRITE(6,979)X
-      WRITE(6,*) ' STEP SIZE T0O SMALL, H=',H
+c$$$      WRITE(6,979)X
+c$$$      WRITE(6,*) ' STEP SIZE T0O SMALL, H=',H
       IDID=-3
       RETURN
  178  CONTINUE
-      WRITE(6,979)X
-      WRITE(6,*) ' MORE THAN NMAX =',NMAX,'STEPS ARE NEEDED'
+c$$$      WRITE(6,979)X
+c$$$      WRITE(6,*) ' MORE THAN NMAX =',NMAX,'STEPS ARE NEEDED'
       IDID=-2
       RETURN
 C --- EXIT CAUSED BY SOLOUT
  179  CONTINUE
-      WRITE(6,979)X
- 979  FORMAT(' EXIT OF RODAS AT X=',E18.4)
+c$$$      WRITE(6,979)X
+c$$$ 979  FORMAT(' EXIT OF RODAS AT X=',E18.4)
       IDID=2
       RETURN
       END
diff --git a/externals/hairer/seulex.f b/externals/hairer/seulex.f
index 799c11a..a396a63 100644
--- a/externals/hairer/seulex.f
+++ b/externals/hairer/seulex.f
@@ -366,7 +366,7 @@ C -------- NMAX , THE MAXIMAL NUMBER OF STEPS -----
       ELSE
          NMAX=IWORK(2)
          IF(NMAX.LE.0)THEN
-            WRITE(6,*)' WRONG INPUT IWORK(2)=',IWORK(2)
+c$$$            WRITE(6,*)' WRONG INPUT IWORK(2)=',IWORK(2)
             ARRET=.TRUE.
          END IF
       END IF
@@ -376,7 +376,7 @@ C -------- KM     MAXIMUM NUMBER OF COLUMNS IN THE EXTRAPOLATION
       ELSE
          KM=IWORK(3)
          IF(KM.LE.2)THEN
-            WRITE(6,*)' CURIOUS INPUT IWORK(3)=',IWORK(3)
+c$$$            WRITE(6,*)' CURIOUS INPUT IWORK(3)=',IWORK(3)
             ARRET=.TRUE.
          END IF
       END IF
@@ -384,19 +384,19 @@ C -------- NSEQU     CHOICE OF STEP SIZE SEQUENCE
       NSEQU=IWORK(4)
       IF(IWORK(4).EQ.0) NSEQU=2
       IF(NSEQU.LE.0.OR.NSEQU.GE.5)THEN
-         WRITE(6,*)' CURIOUS INPUT IWORK(4)=',IWORK(4)
+c$$$         WRITE(6,*)' CURIOUS INPUT IWORK(4)=',IWORK(4)
          ARRET=.TRUE.
       END IF
 C -------- LAMBDA   PARAMETER FOR DENSE OUTPUT
       LAMBDA=IWORK(5)
       IF(LAMBDA.LT.0.OR.LAMBDA.GE.2)THEN
-         WRITE(6,*)' CURIOUS INPUT IWORK(5)=',IWORK(5)
+c$$$         WRITE(6,*)' CURIOUS INPUT IWORK(5)=',IWORK(5)
          ARRET=.TRUE.
       END IF
 C -------- NRDENS   NUMBER OF DENSE OUTPUT COMPONENTS
       NRDENS=IWORK(6)
       IF(NRDENS.LT.0.OR.NRDENS.GT.N)THEN
-         WRITE(6,*)' CURIOUS INPUT IWORK(6)=',IWORK(6)
+c$$$         WRITE(6,*)' CURIOUS INPUT IWORK(6)=',IWORK(6)
          ARRET=.TRUE.
       END IF
 C -------- PARAMETER FOR SECOND ORDER EQUATIONS
@@ -406,7 +406,7 @@ C -------- PARAMETER FOR SECOND ORDER EQUATIONS
       IF (M1.EQ.0) M2=N
       IF (M2.EQ.0) M2=M1
       IF (M1.LT.0.OR.M2.LT.0.OR.M1+M2.GT.N) THEN
-       WRITE(6,*)' CURIOUS INPUT FOR IWORK(9,10)=',M1,M2
+c$$$       WRITE(6,*)' CURIOUS INPUT FOR IWORK(9,10)=',M1,M2
        ARRET=.TRUE.
       END IF
 C -------- UROUND   SMALLEST NUMBER SATISFYING 1.D0+UROUND>1.D0
@@ -415,7 +415,7 @@ C -------- UROUND   SMALLEST NUMBER SATISFYING 1.D0+UROUND>1.D0
       ELSE
          UROUND=WORK(1)
          IF(UROUND.LE.0.D0.OR.UROUND.GE.1.D0)THEN
-            WRITE(6,*)'  UROUND=',WORK(1)
+c$$$            WRITE(6,*)'  UROUND=',WORK(1)
             ARRET=.TRUE.
          END IF
       END IF
@@ -489,13 +489,13 @@ C ------- WKFCN,WKJAC,WKDEC,WKSOL  ESTIMATED WORK FOR  FCN,JAC,DEC,SOL
 C --------- CHECK IF TOLERANCES ARE O.K.
       IF (ITOL.EQ.0) THEN
           IF (ATOL(1).LE.0.D0.OR.RTOL(1).LE.10.D0*UROUND) THEN
-              WRITE (6,*) ' TOLERANCES ARE TOO SMALL'
+c$$$              WRITE (6,*) ' TOLERANCES ARE TOO SMALL'
               ARRET=.TRUE.
           END IF
       ELSE
         DO I=1,N
           IF (ATOL(I).LE.0.D0.OR.RTOL(I).LE.10.D0*UROUND) THEN
-              WRITE (6,*) ' TOLERANCES(',I,') ARE TOO SMALL'
+c$$$              WRITE (6,*) ' TOLERANCES(',I,') ARE TOO SMALL'
               ARRET=.TRUE.
           END IF
         END DO
@@ -533,8 +533,8 @@ C -- MASS MATRIX
           END IF
 C ------ BANDWITH OF "MAS" NOT LARGER THAN BANDWITH OF "JAC"
           IF (MLMAS.GT.MLJAC.OR.MUMAS.GT.MUJAC) THEN
-              WRITE (6,*) 'BANDWITH OF "MAS" NOT LARGER THAN BANDWITH OF
-     & "JAC"'
+c$$$              WRITE (6,*) 'BANDWITH OF "MAS" NOT LARGER THAN BANDWITH OF
+c$$$     & "JAC"'
             ARRET=.TRUE.
           END IF
       ELSE
@@ -549,8 +549,8 @@ C ------ BANDWITH OF "MAS" NOT LARGER THAN BANDWITH OF "JAC"
       LDMAS2=MAX(1,LDMAS)
 C ------ HESSENBERG OPTION ONLY FOR EXPLICIT EQU. WITH FULL JACOBIAN
       IF ((IMPLCT.OR.JBAND).AND.IJOB.EQ.7) THEN
-         WRITE(6,*)' HESSENBERG OPTION ONLY FOR EXPLICIT EQUATIONS WITH
-     &FULL JACOBIAN'
+c$$$         WRITE(6,*)' HESSENBERG OPTION ONLY FOR EXPLICIT EQUATIONS WITH
+c$$$     &FULL JACOBIAN'
          ARRET=.TRUE.
       END IF
 C ------- PREPARE THE ENTRY-POINTS FOR THE ARRAYS IN WORK -----
@@ -576,7 +576,7 @@ C ------- PREPARE THE ENTRY-POINTS FOR THE ARRAYS IN WORK -----
 C ------ TOTAL STORAGE REQUIREMENT -----------
       ISTORE=IFSAFE+KM2*NRDENS-1
       IF(ISTORE.GT.LWORK)THEN
-         WRITE(6,*)' INSUFFICIENT STORAGE FOR WORK, MIN. LWORK=',ISTORE
+c$$$         WRITE(6,*)' INSUFFICIENT STORAGE FOR WORK, MIN. LWORK=',ISTORE
          ARRET=.TRUE.
       END IF
 C ------- ENTRY POINTS FOR INTEGER WORKSPACE -----
@@ -587,7 +587,7 @@ C ------- ENTRY POINTS FOR INTEGER WORKSPACE -----
 C --------- TOTAL REQUIREMENT ---------------
       ISTORE=IECO+KM-1
       IF(ISTORE.GT.LIWORK)THEN
-         WRITE(6,*)' INSUFF. STORAGE FOR IWORK, MIN. LIWORK=',ISTORE
+c$$$         WRITE(6,*)' INSUFF. STORAGE FOR IWORK, MIN. LIWORK=',ISTORE
          ARRET=.TRUE.
       END IF
 C ------ WHEN A FAIL HAS OCCURED, WE RETURN WITH IDID=-1
@@ -992,9 +992,9 @@ C --- SOLUTION EXIT
       IDID=1
       RETURN
 C --- FAIL EXIT
- 120  WRITE (6,979) X,H
- 979  FORMAT(' EXIT OF SEULEX AT X=',D14.7,'   H=',D14.7)
-      IDID=-1
+c$$$ 120  WRITE (6,979) X,H
+c$$$ 979  FORMAT(' EXIT OF SEULEX AT X=',D14.7,'   H=',D14.7)
+ 120  IDID=-1
       RETURN
       END
 C
@@ -1178,7 +1178,7 @@ C ----- COMPUTE PLACE OF II-TH COMPONENT
         IF (IC(J).EQ.II) I=J
       END DO
       IF (I.EQ.0) THEN
-         WRITE (6,*) ' NO DENSE OUTPUT AVAILABLE FOR COMP.',II
+c$$$         WRITE (6,*) ' NO DENSE OUTPUT AVAILABLE FOR COMP.',II
          RETURN
       END IF
 C ----- COMPUTE THE INTERPOLATED VALUE
diff --git a/externals/netlib/odepack/opkda2.f b/externals/netlib/odepack/opkda2.f
index c4b4ff0..2cfb0e1 100644
--- a/externals/netlib/odepack/opkda2.f
+++ b/externals/netlib/odepack/opkda2.f
@@ -1248,16 +1248,16 @@ C***FIRST EXECUTABLE STATEMENT  XERRWD
 C
 C  Write the message.
 C
-      WRITE (LUNIT,10)  MSG
- 10   FORMAT(1X,A)
-      IF (NI .EQ. 1) WRITE (LUNIT, 20) I1
- 20   FORMAT(6X,'In above message,  I1 =',I10)
-      IF (NI .EQ. 2) WRITE (LUNIT, 30) I1,I2
- 30   FORMAT(6X,'In above message,  I1 =',I10,3X,'I2 =',I10)
-      IF (NR .EQ. 1) WRITE (LUNIT, 40) R1
- 40   FORMAT(6X,'In above message,  R1 =',D21.13)
-      IF (NR .EQ. 2) WRITE (LUNIT, 50) R1,R2
- 50   FORMAT(6X,'In above,  R1 =',D21.13,3X,'R2 =',D21.13)
+c$$$      WRITE (LUNIT,10)  MSG
+c$$$ 10   FORMAT(1X,A)
+c$$$      IF (NI .EQ. 1) WRITE (LUNIT, 20) I1
+c$$$ 20   FORMAT(6X,'In above message,  I1 =',I10)
+c$$$      IF (NI .EQ. 2) WRITE (LUNIT, 30) I1,I2
+c$$$ 30   FORMAT(6X,'In above message,  I1 =',I10,3X,'I2 =',I10)
+c$$$      IF (NR .EQ. 1) WRITE (LUNIT, 40) R1
+c$$$ 40   FORMAT(6X,'In above message,  R1 =',D21.13)
+c$$$      IF (NR .EQ. 2) WRITE (LUNIT, 50) R1,R2
+c$$$ 50   FORMAT(6X,'In above,  R1 =',D21.13,3X,'R2 =',D21.13)
 C
 C  Abort the run if LEVEL = 2.
 C