//---------------------------------------------------------------------
// first off create an array of points
//---------------------------------------------------------------------
#local point = array[n1][n2];

// this is an array of vectors and is filled by the removed code below:

// ...


//---------------------------------------------------------------------
// draw the final cloth
//---------------------------------------------------------------------
union{
 mesh2{



  vertex_vectors{
   n1*n2
   #local px=0;
   #while(px<n1)
   #local py=0;
   #while(py<n2)

    point[px][py]

   #local py=py+1;
   #end
   #local px=px+1;
   #end
  }
  
  
  
  normal_vectors{
   n1*n2
   #local px=0;
   #while(px<n1)
   #local py=0;
   #while(py<n2)

    // work out the normal here 
    #local vNorm  = <0,0,0>;
    #local vNormn = 0;
    
    #if( px>=1 & py<=n2-2)
     #local vNorm = vNorm + vcross( point[px][py] - point[px][py+1] , point[px][py] - point[px-1][py] );
     #local vNormn = vNormn + 1;
    #end
    
    #if( px<=n1-2 & py<=n2-2)
     #local vNorm = vNorm - vcross( point[px][py] - point[px][py+1] , point[px][py] - point[px+1][py] );
     #local vNormn = vNormn + 1;    
    #end
    
    #if( px<=n1-2 & py>=1)
     #local vNorm = vNorm + vcross( point[px][py] - point[px][py-1] , point[px][py] - point[px+1][py] );
     #local vNormn = vNormn + 1;
    #end
   
    #if( px>=1 & py>=1)
     #local vNorm = vNorm - vcross( point[px][py] - point[px][py-1] , point[px][py] - point[px-1][py] );
     #local vNormn = vNormn + 1;
    #end
   
    (vNorm / vNormn)
   
   #local py=py+1;
   #end
   #local px=px+1;
   #end
  }


 uv_vectors{
   n1*n2
   #local px=0;
   #while(px<n1)
   #local py=0;
   #while(py<n2)

    <px/n1,py/n2>

   #local py=py+1;
   #end
   #local px=px+1;
   #end
   }
 
 
 face_indices{
  (n1-1)*(n2-1)*2
  #local px=0;
  #while(px<n1-1)
   #local py=0;
   #while(py<n2-1)

    <py + px*n2 , py + (px+1)*n2 , py+1 + px*n2>
    <py+1 + (px+1)*n2 , py + (px+1)*n2 , py+1 + px*n2>


   #local py=py+1;
   #end
   #local px=px+1;
   #end
   }
 

 normal_indices{
  (n1-1)*(n2-1)*2
  #local px=0;
  #while(px<n1-1)
   #local py=0;
   #while(py<n2-1)

    <py + px*n2 , py + (px+1)*n2 , py+1 + px*n2>
    <py+1 + (px+1)*n2 , py + (px+1)*n2 , py+1 + px*n2>


   #local py=py+1;
   #end
   #local px=px+1;
   #end
   }
 
 uv_indices{
  (n1-1)*(n2-1)*2
  #local px=0;
  #while(px<n1-1)
   #local py=0;
   #while(py<n2-1)

    <py + px*n2 , py + (px+1)*n2 , py+1 + px*n2>
    <py+1 + (px+1)*n2 , py + (px+1)*n2 , py+1 + px*n2>


   #local py=py+1;
   #end
   #local px=px+1;
   #end
   }
 
 
 }



 }