static void grid_setcoords(grid* g, int htype, int hnodetype, int nx, int ny, int nz, void* x, void* y, double* z)
{
g->htype = htype;
if (htype == GRIDHTYPE_LATLON)
g->gridnodes_xy = gnxy_simple_create(nx, ny, x, y);
#if !defined(NO_GRIDUTILS)
else if (htype == GRIDHTYPE_CURVILINEAR) {
g->gridnodes_xy = gnxy_curv_create(hnodetype, nx, ny, x, y, g->maptype);
#if defined(GRIDNODES_WRITE)
{
char fname[MAXSTRLEN];
snprintf(fname, MAXSTRLEN, "gridnodes-%d.txt", grid_getid(g));
if (!file_exists(fname))
gridnodes_write(((gnxy_curv*) g->gridnodes_xy)->gn, fname, CT_XY);
}
#endif
}
#endif
else
enkf_quit("programming error");
grid_setlonbase(g);
if (g->vtype == GRIDVTYPE_Z || g->vtype == GRIDVTYPE_SIGMA)
g->gridnodes_z = gnz_simple_create(nz, z);
else
enkf_quit("not implemented");
}
int main(int argc, char* argv[])
{
char* gridfname = NULL;
char* maskfname = NULL;
NODETYPE ntin = NT_DD;
NODETYPE ntout = NT_DD;
COORDTYPE ct = CT_XY;
int tweaknpolar = 0;
gridnodes* gn = NULL;
parse_commandline(argc, argv, &gridfname, &maskfname, &ntin, &ntout, &ct, &tweaknpolar);
gn = gridnodes_read(gridfname, ntin);
gridnodes_validate(gn);
if (ntin != ntout) {
gridnodes* gnnew = gridnodes_transform(gn, ntout);
gridnodes_destroy(gn);
gn = gnnew;
}
if (maskfname != NULL) {
int nx = gridnodes_getnce1(gn);
int ny = gridnodes_getnce2(gn);
int** mask = gu_readmask(maskfname, nx, ny);
gridnodes_applymask(gn, mask);
gu_free2d(mask);
}
if (tweaknpolar)
gridnodes_tweaknpolar(gn);
if (gu_verbose)
gridnodes_calcstats(gn);
gridnodes_write(gn, "stdout", ct);
gridnodes_destroy(gn);
return 0;
}
int main(int argc, char* argv[])
{
char* fname1 = NULL;
char* fname2 = NULL;
gridnodes* grid1 = NULL;
gridnodes* grid2 = NULL;
gridnodes* grid3 = NULL;
NODETYPE type = NT_NONE;
int i0 = INT_MAX;
int j0 = INT_MAX;
int nx1 = 0;
int ny1 = 0;
int nx2 = 0;
int ny2 = 0;
int nx3 = 0;
int ny3 = 0;
double** gx1 = NULL;
double** gy1 = NULL;
double** gx2 = NULL;
double** gy2 = NULL;
double** gx3 = NULL;
double** gy3 = NULL;
int i1start, i2start, j1start, j2start;
int i1, i2, j1, j2, i3, j3;
parse_commandline(argc, argv, &fname1, &fname2, &i0, &j0, &type);
grid1 = gridnodes_read(fname1, NT_NONE);
grid2 = gridnodes_read(fname2, NT_NONE);
nx1 = gridnodes_getnx(grid1);
ny1 = gridnodes_getny(grid1);
nx2 = gridnodes_getnx(grid2);
ny2 = gridnodes_getny(grid2);
gx1 = gridnodes_getx(grid1);
gy1 = gridnodes_gety(grid1);
gx2 = gridnodes_getx(grid2);
gy2 = gridnodes_gety(grid2);
nx3 = (i0 + nx2 > nx1) ? i0 + nx2 : nx1;
ny3 = (j0 + ny2 > ny1) ? j0 + ny2 : ny1;
grid3 = gridnodes_create(nx3, ny3, type);
gx3 = gridnodes_getx(grid3);
gy3 = gridnodes_gety(grid3);
if (i0 < 0) {
i1start = i0;
i2start = 0;
} else {
i1start = 0;
i2start = -i0;
}
if (j0 < 0) {
j1start = j0;
j2start = 0;
} else {
j1start = 0;
j2start = -j0;
}
if (!merge) {
for (j3 = 0, j1 = j1start, j2 = j2start; j3 < ny3; ++j3, ++j1, ++j2) {
for (i3 = 0, i1 = i1start, i2 = i2start; i3 < nx3; ++i3, ++i1, ++i2) {
int ok2 = (j2 >= 0 && j2 < ny2 && i2 >= 0 && i2 <= nx2);
if (ok2) {
gx3[j3][i3] = gx2[j2][i2];
gy3[j3][i3] = gy2[j2][i2];
} else {
int ok1 = (j1 >= 0 && j1 < ny1 && i1 >= 0 && i1 <= nx1);
if (ok1) {
gx3[j3][i3] = gx1[j1][i1];
gy3[j3][i3] = gy1[j1][i1];
} else {
gx3[j3][i3] = NaN;
gy3[j3][i3] = NaN;
}
}
}
}
} else {
for (j3 = 0, j1 = j1start, j2 = j2start; j3 < ny3; ++j3, ++j1, ++j2) {
for (i3 = 0, i1 = i1start, i2 = i2start; i3 < nx3; ++i3, ++i1, ++i2) {
int ok2 = (j2 >= 0 && j2 < ny2 && i2 >= 0 && i2 <= nx2);
int ok1 = (j1 >= 0 && j1 < ny1 && i1 >= 0 && i1 <= nx1);
if (!ok2 && !ok1) {
gx3[j3][i3] = NaN;
gy3[j3][i3] = NaN;
} else if (ok2 && !ok1) {
gx3[j3][i3] = gx2[j2][i2];
gy3[j3][i3] = gy2[j2][i2];
} else if (!ok2 && ok1) {
gx3[j3][i3] = gx1[j1][i1];
gy3[j3][i3] = gy1[j1][i1];
} else {
if (!isnan(gx2[j2][i2])) {
gx3[j3][i3] = gx2[j2][i2];
gy3[j3][i3] = gy2[j2][i2];
} else {
gx3[j3][i3] = gx1[j1][i1];
gy3[j3][i3] = gy1[j1][i1];
}
}
}
}
}
gridnodes_validate(grid3);
gridnodes_write(grid3, "stdout", CT_XY);
gridnodes_destroy(grid1);
gridnodes_destroy(grid2);
gridnodes_destroy(grid3);
return 0;
}
