anchor.diff:        OPEN(UNIT=11, FILE="dphiw")
         WRITE(11,1001) DPHIW
         CLOSE(UNIT=11)
  + 1001 FORMAT('',F15.6,2X)
  +END SUBROUTINE
  +
  +SUBROUTINE SAVECL(NTime, CL)
  +      integer, intent(IN) :: NTime
  +      real(kind=8), dimension(NTime), intent(IN) :: CL
  +      OPEN(UNIT=14, FILE="cltime")
  +      WRITE(14,1001) CL
  +      CLOSE(UNIT=14)
    1001 FORMAT('',F15.6,2X)
   END SUBROUTINE
   
  
positive.code:       end do
        call SolvePhiTime(Nelem, B, C, NTime, D, Chit, Phit, DPhiW)
        call CalcCl(Nelem, Xnode, TEat1, NTime, DT, Ut, Phit, Chit, cl)
        write(*,*) cl
  !     Save results
        call SavePhi(Nelem, NTime, PhiT)
        call SaveXWnode(NTime, XWNODE)
        call SaveDPhiW(NTime, DPHIW)
  
negative.code:       CLOSE(UNIT=13)
   1001 FORMAT('',1000(F15.6),2X)
  END SUBROUTINE
  
  
  