anchor.diff:  !
   !     WUENSCHE
   !     -ALLE OUTPUTS VON ZENTRALER ROUTINE, DIE RICHTIG FORMATIERT
  +
  +      USE, INTRINSIC :: ISO_FORTRAN_ENV, ONLY : OUTPUT_UNIT, ERROR_UNIT
  +
         IMPLICIT DOUBLE PRECISION (A-H,O-Z)
  -      PARAMETER(PIENI=1D-17)
  -      PARAMETER(ZERO=0.0D0,HALF=5.0D-1,ONE=1.0D0)
  -      PARAMETER(HALFM=-5.0D-1,ONEM=-1.0D0)
  -      PARAMETER(TWO=2.0D0,THREE=3.0D0,FOUR=4.0D0)
  -      PARAMETER(C1E3=1.0D3,C2E3=2.0D3,C4E3=4.0D3,C1E4=1.0D4)
  -      PARAMETER(C1E12=1.0D12,C1E13=1.0D13,C1E15=1.0D15,C1E16=1.0D16)
  -      PARAMETER(C1M1=1.0D-1,C1M3=1.0D-3,C1M6=1.0D-6,C1M7=1.0D-7)
  -      PARAMETER(C1M8=1.0D-8,C1M9=1.0D-9,C1M10=1.0D-10)
  -      PARAMETER(C1M12=1.0D-12,C1M13=1.0D-13,C1M14=1.0D-14)
  -      PARAMETER(C1M15=1.0D-15,C1M17=1.0D-17,C1M18=1.0D-18)
  -      PARAMETER(C1M21=1.0D-21,C1M24=1.0D-24,C1M38=1.0D-38)
  -      PARAMETER(FIVE=5.0D0,SIX=6.0D0,SEVEN=7.0D0,EIGHT=8.0D0)
  -!      PARAMETER(NINE=9.0D0,TEN=10.0D0)
  -      PARAMETER(C24E0=24.0D0,C120E0=120.0D0,C16E0=16.0D0,C40E0=40.0D0)
  -      PARAMETER(C80E0=80.0D0,C72E0=72.0D0)
  -      PARAMETER(C12E0=12.0D0,C32E0=32.0D0,C48E0=48.0D0,C160E0=160.0D0)
  -      PARAMETER(ISZ=200)
  -!      CHARACTER A*160
  -!      INTEGER NA(50)
  -!
  +
         !For command line argument parsing
         INTEGER :: cmdarg_i, cmdarg_length, cmdarg_status
         CHARACTER(len=100) :: cmdarg_arg
  
positive.code:             continue
           else if (cmdarg_i.eq.1) then
              if (cmdarg_status.ne.0) then
                 write(6,*) "Error: Input file name too long."
                 stop 1
              end if
              fname_in = trim(cmdarg_arg)
           else if (cmdarg_i.eq.2) then
              if (cmdarg_status.ne.0) then
                 write(6,*) "Error: Output file name too long."
                 stop 1
              end if
              fname_out = trim(cmdarg_arg)
           else
              write(6,*) "Error: Expected either no arguments OR "
              write(6,*) "two arguements (input and output file name)"
              stop 1
           end if
           cmdarg_i = cmdarg_i+1
  
negative.code:          ITA = ISCRTY(IA)
           ND = 0
           IF(IA.LT.1.OR.IA.GT.100.OR.IA.GT.LSCR) THEN
              WRITE(ERROR_UNIT,'(1X,A,I4)') '!!! ERROR IN VNAM, IA = ',IA
              WRITE(2,'(1X,A,I4)') '!!! ERROR IN VNAM, IA = ',IA
              STOP 11
           ENDIF
  