int model_fk2z(model* m, int vid, int i, int j, double fk, double* z) { grid* g = m->grids[m->vars[vid].gridid]; float** depth; int ni, nj; grid_getdims(g, &ni, &nj, NULL); if (i < 0 || j < 0 || i >= ni || j >= nj) { *z = NaN; return STATUS_OUTSIDEGRID; } grid_fk2z(g, i, j, fk, z); depth = grid_getdepth(g); if (*z > depth[j][i]) { *z = NaN; return STATUS_OUTSIDEGRID; } return STATUS_OK; }
int model_fk2z(model* m, int vid, int i, int j, double fk, double* z) { return grid_fk2z(m->grids[m->vars[vid].gridid], i, j, fk, z); }