gridspec_api/regrid/nccf_apply_regrid_type.h

/* $Id: nccf_apply_regrid_type.h 254 2010-11-19 22:10:55Z dkindig $ */

_TYPE_ *ori_data, *tgt_data;
_TYPE_ datum;
status = nccf_get_struct_data_pointer(ori_data_id, (void **)&ori_data);
totErr += abs(status);
status = nccf_get_struct_data_pointer(tgt_data_id, (void **)&tgt_data);
totErr += abs(status);

/* Make sure at least one target point is in the domain */
if (self->n_nonmasked > 0) {
  for (i = 0; i < self->ntargets; ++i) {
    /* Compute sum of weights * field values when inside_domain == 1
       Don't change the value if inside_domain == 0 
       This will set tgt_data[i] to zero if valid and keep 
       tgt_data[i] otherwise
    */
    tgt_data[i] = (1 - self->inside_domain[i])*tgt_data[i];
    /* Linear interpolation */
    for (j = 0; j < nNodes; ++j) {
      datum = ori_data[ self->indices[i][j] ];
      tgt_data[i] += (_TYPE_) (self->weights[i][j]*datum*self->inside_domain[i]);
    }
  }
}