grid.c

	#include <stdio.h>
	#define N1 1690
	#define N2 285
	int main(void) {
	double gridsp_x[4] = {
	0.32,
	0.16,
	0.08,
	0.04
	};
	double gridsp_y[4] = {
	0.32,
	0.16,
	0.08,
	0.04
	};
	double gridset_x[7];
	double gridset_y[7];
	double iii[N1][N2], jjj[N1][N2], dx[N1][N2], dy[N1][N2];
	int sp_x[10] = {
	-1000,
	-1000,
	30,
	40,
	50,
	1540,
	1550,
	1560,
	1000,
	1000
	};
	int sp_y[10] = {
	-1000,
	-1000,
	60,
	70,
	80,
	205,
	215,
	225,
	1000,
	1000
	};
	int i, j, im, jm;
	jm = N2;
	im = N1;
	gridset_x[0] = -46;
	gridset_x[1] = gridset_x[0] + gridsp_x[0] * sp_x[2];
	gridset_y[0] = -24.1;
	gridset_y[1] = gridset_y[0] + gridsp_y[0] * sp_y[2];
	for (i = 2; i < 7; i++) {
	gridset_x[i] = gridset_x[i - 1] + gridsp_x[i - 1] * (sp_x[i + 1] - sp_x[i]);
	gridset_y[i] = gridset_y[i - 1] + gridsp_y[i - 1] * (sp_y[i + 1] - sp_y[i]);
	}
	for (j = 0; j < jm; j++) {
	for (i = 0; i < im; i++) {
	if (i < sp_x[2] || i > sp_x[7]) {
	dx[i][j] = gridsp_x[0];
	} else if (i < sp_x[3] || i > sp_x[6]) {
	dx[i][j] = gridsp_x[1];
	} else if (i < sp_x[4] || i > sp_x[5]) {
	dx[i][j] = gridsp_x[2];
	} else {
	dx[i][j] = gridsp_x[3];
	}
	}
	}
	for (j = 0; j < jm; j++) {
	for (i = 0; i < im; i++) {
	if (j < sp_y[2] || j > sp_y[7]) {
	dy[i][j] = gridsp_y[0];
	} else if (j < sp_y[3] || j > sp_y[6]) {
	dy[i][j] = gridsp_y[1];
	} else if (j < sp_y[4] || j > sp_y[5]) {
	dy[i][j] = gridsp_y[2];
	} else {
	dy[i][j] = gridsp_y[3];
	}
	}
	}
	for (j = 0; j < jm; j++) {
	for (i = 0; i < im; i++) {
	if (i <= sp_x[2]) {
	iii[i][j] = gridset_x[0] + gridsp_x[0] * i;
	}
	if (i <= sp_x[3] && i > sp_x[2]) {
	iii[i][j] = gridset_x[1] + gridsp_x[1] * (i - sp_x[2]);
	}
	if (i <= sp_x[4] && i > sp_x[3]) {
	iii[i][j] = gridset_x[2] + gridsp_x[2] * (i - sp_x[3]);
	}
	if (i <= sp_x[5] && i > sp_x[4]) {
	iii[i][j] = gridset_x[3] + gridsp_x[3] * (i - sp_x[4]);
	}
	if (i <= sp_x[6] && i > sp_x[5]) {
	iii[i][j] = gridset_x[4] + gridsp_x[2] * (i - sp_x[5]);
	}
	if (i <= sp_x[7] && i > sp_x[6]) {
	iii[i][j] = gridset_x[5] + gridsp_x[1] * (i - sp_x[6]);
	}
	if (i > sp_x[7]) {
	iii[i][j] = gridset_x[6] + gridsp_x[0] * (i - sp_x[7]);
	}
	}
	}
	for (j = 0; j < jm; j++) {
	for (i = 0; i < im; i++) {
	if (j <= sp_y[2]) {
	jjj[i][j] = gridset_y[0] + gridsp_y[0] * j;
	}
	if (j <= sp_y[3] && j > sp_y[2]) {
	jjj[i][j] = gridset_y[1] + gridsp_y[1] * (j - sp_y[2]);
	}
	if (j <= sp_y[4] && j > sp_y[3]) {
	jjj[i][j] = gridset_y[2] + gridsp_y[2] * (j - sp_y[3]);
	}
	if (j <= sp_y[5] && j > sp_y[4]) {
	jjj[i][j] = gridset_y[3] + gridsp_y[3] * (j - sp_y[4]);
	}
	if (j <= sp_y[6] && j > sp_y[5]) {
	jjj[i][j] = gridset_y[4] + gridsp_y[2] * (j - sp_y[5]);
	}
	if (j <= sp_y[7] && j > sp_y[6]) {
	jjj[i][j] = gridset_y[5] + gridsp_y[1] * (j - sp_y[6]);
	}
	if (j > sp_y[7]) {
	jjj[i][j] = gridset_y[6] + gridsp_y[0] * (j - sp_y[7]);
	}
	}
	}
	return 0;
	}

view raw grid.c hosted with ❤ by GitHub