anchor.diff:       &     ngauss,                                                      ! eltype
        &     delti,nstep,                                                 ! timdat
        &     istatout,istatoutc,nstatout,nstatestot,nout,npts,            ! ioopts
  +     &     statemin,statemax,ibyteg,                                    ! binucd
        &     idout,idsk,iucd,kw,kp,kucd)                                  ! ioinfo
           else if(matmodel.eq.3) then
             call write_state_cmp(
  
positive.code:      &     ngauss,                                                      ! eltype
       &     delti,nstep,                                                 ! timdat
       &     istatout,istatoutc,nstatout,nstatestot,nout,npts,            ! ioopts
       &     idout,idsk,iucd,kw,kp,kucd)                                  ! ioinfo
          else if(matmodel.eq.13) then
            call write_state_cmp(
  
negative.code:      &     ngauss,                                                      ! eltype
       &     delti,nstep,                                                 ! timdat
       &     istatout,istatoutc,nstatout,nstatestot,nout,npts,            ! ioopts
       &     statemin,statemax,ibyteg,                                    ! binucd
       &     idout,idsk,iucd,kw,kp,kucd)                                  ! ioinfo
          else
            ierr=101
  