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);
}
