anchor.diff:        f1=half*deltp*(one-alfap)/vis
         f2=one/((one+pr)/e+half*alfap*deltp/vis)
         emean=third*(ee(1)+ee(2)+ee(3))
  +      smean=e*emean/(one-two*pr)
  +      smeantp=third*(state(1,1)+state(2,1)+state(3,1))
  +cdebug      write(6,*) "e,pr,vis,rmu,tmax,f1,f2,emean,smean,smeantp:",
  +cdebug     & e,pr,vis,rmu,tmax,f1,f2,emean,smean,smeantp
         do i=1,nstr
           eet=ee(i)-diag(i)*emean-state(i,3)
  -        dstate(i,1)=f2*(eet-f1*state(i,1))
  -        stau=(one-alfap)*state(i,1)+alfap*dstate(i,1)
  +        sdevtp=state(i,1)-diag(i)*smeantp
  +        sdev=f2*(eet-f1*sdevtp)
  +        dstate(i,1)=sdev+diag(i)*smean
  +        stau=(one-alfap)*sdevtp+alfap*sdev
           dstate(i,3)=half*deltp*stau/vis
           dstate(i,2)=ee(i)
  +cdebug        write(6,*) "i,eet,sdevtp,sdev,ds1,ds2,ds3,s1,s2,s3,stau:",
  +cdebug     &   i,eet,sdevtp,sdev,dstate(i,1),dstate(i,2),dstate(i,3),
  +cdebug     &   state(i,1),state(i,2),state(i,3),stau
         end do
   c
         return
  
positive.code:         errstrng="mat_prt_5"
          return
  c
   700  format("Material number:       ",i7,/,
       &       "Material type:         ",a37,/,
       &       "Number of properties:  ",i7,/)
   710  format(15x,a15,3x,1pe15.8)
  c
        end
  
negative.code:           dmat(iddmat(i,j))=od
          end do
        end do
  cdebug      write(6,*) "dmat:",(dmat(idb),idb=1,nddmat)
        return
        end
  c
  