anchor.diff:  
      ! calculate the number of digits in 'y_FAST%NOutSteps' (Maximum number of output steps to be written)
      ! this will be used to pad the write-out step in the VTK filename with zeros in calls to MeshWrVTK_...()
  -   Twidth = int(log10(real(y_FAST%NOutSteps))) + 1
  +   if (y_FAST%NOutSteps==0) then
  +      Twidth = 1
  +   else
  +      Twidth = int(log10(real(y_FAST%NOutSteps))) + 1
  +   endif
      
      
      NumBl = 0
  
positive.code: 
     ! calculate the number of digits in 'y_FAST%NOutSteps' (Maximum number of output steps to be written)
     ! this will be used to pad the write-out step in the VTK filename with zeros in calls to MeshWrVTK()
     Twidth = int(log10(real(y_FAST%NOutSteps))) + 1
     
     NumBl = 0
     if (allocated(ED%Output)) then
  
negative.code: 
     ! calculate the number of digits in 'y_FAST%NOutSteps' (Maximum number of output steps to be written)
     ! this will be used to pad the write-out step in the VTK filename with zeros in calls to MeshWrVTK_...()
     if (y_FAST%NOutSteps==0) then
        Twidth = 1
     else
        Twidth = int(log10(real(y_FAST%NOutSteps))) + 1
     endif
  
     VTK_path = get_vtkroot_path( p_FAST%OutFileRoot )
  
  