// And compute indices:
              for (i = 0; i < nlat-1; i++) {
                 for (j = 0; j < nl; j++) {
                    idx = 5*(i*nl+j);
                    ci[idx+0] = i*nl+j;
                    ci[idx+1] = i*nl+(j+1)\%nl;
                    ci[idx+2] = (i+1)*nl+(j+1)\%nl;
                    ci[idx+3] = (i+1)*nl+j;
                    ci[idx+4] = -1;  // end-of-polygon
                 }
              }
           }"
        ]
     }