static void ps_draw_bond(t_psdata ps,
atom_id ai,atom_id aj,iv2 vec2[],
rvec x[],char **atomnm[],int size[],bool bBalls)
{
char *ic,*jc;
int xi,yi,xj,yj;
int xm,ym;
if (bBalls) {
ps_draw_atom(ps,ai,vec2,atomnm);
ps_draw_atom(ps,aj,vec2,atomnm);
}
else {
if (local_pbc_dx(x[ai],x[aj])) {
ic=*atomnm[ai];
jc=*atomnm[aj];
xi=vec2[ai][XX];
yi=vec2[ai][YY];
xj=vec2[aj][XX];
yj=vec2[aj][YY];
if (ic != jc) {
xm=(xi+xj) >> 1;
ym=(yi+yj) >> 1;
ps_rgb(ps,Type2RGB(ic));
ps_line(ps,xi,yi,xm,ym);
ps_rgb(ps,Type2RGB(jc));
ps_line(ps,xm,ym,xj,yj);
}
else {
ps_rgb(ps,Type2RGB(ic));
ps_line(ps,xi,yi,xj,yj);
}
}
static void ps_draw_atom(t_psdata ps,atom_id ai,iv2 vec2[],char **atomnm[])
{
int xi,yi;
xi=vec2[ai][XX];
yi=vec2[ai][YY];
ps_rgb(ps,Type2RGB(*atomnm[ai]));
ps_line(ps,xi-MSIZE,yi,xi+MSIZE+1,yi);
ps_line(ps,xi,yi-MSIZE,xi,yi+MSIZE+1);
}
void paint_ellipse (double x0, double y0, double angle, double major, double minor, double scale, double t11,double t12,double t21,double t22, GMT_LONG polygon, int rgb[3], GMT_LONG outline)
/* Make an ellipse at center x0,y0 */
{
#define NPOINTS 362
GMT_LONG npoints = NPOINTS;
GMT_LONG i;
/* relative to center of ellipse */
double dxe[NPOINTS],dye[NPOINTS];
/* absolute paper coordinates */
double axe[NPOINTS],aye[NPOINTS];
trace_ellipse(angle, major, minor, npoints, dxe, dye);
for (i = 0; i < npoints - 2; i++) {
transform_local (x0,y0,dxe[i],dye[i],scale,t11,t12,t21,t22,&axe[i],&aye[i]);
}
if(polygon)
ps_polygon(axe, aye, npoints - 2, rgb, outline);
else
ps_line(axe, aye, npoints - 2, 3, TRUE);
}
static void draw_zerolines(t_psdata out, real x0, real y0, real w,
int nmat, t_matrix mat[], t_psrec *psr)
{
real xx, yy, dy, xx00, yy00;
int i, x, y;
xx00 = x0-1.5;
yy00 = y0-1.5;
ps_linewidth(out, psr->zerolinewidth);
for (i = 0; (i < nmat); i++)
{
dy = box_height(&(mat[i]), psr);
/* mat[i].axis_x and _y were already set by draw_boxes */
if (psr->X.lineatzero)
{
ps_rgb(out, linecolors[psr->X.lineatzero]);
for (x = 0; (x < mat[i].nx); x++)
{
xx = xx00+(x+0.7)*psr->xboxsize;
/* draw lines whenever tick label almost zero (e.g. next trajectory) */
if (x != 0 && x < mat[i].nx-1 &&
abs(mat[i].axis_x[x]) <
0.1*abs(mat[i].axis_x[x+1]-mat[i].axis_x[x]) )
{
ps_line (out, xx, yy00, xx, yy00+dy+2);
}
}
}
if (psr->Y.lineatzero)
{
ps_rgb(out, linecolors[psr->Y.lineatzero]);
for (y = 0; (y < mat[i].ny); y++)
{
yy = yy00+(y+0.7)*psr->yboxsize;
/* draw lines whenever tick label almost zero (e.g. next trajectory) */
if (y != 0 && y < mat[i].ny-1 &&
abs(mat[i].axis_y[y]) <
0.1*abs(mat[i].axis_y[y+1]-mat[i].axis_y[y]) )
{
ps_line (out, xx00, yy, xx00+w+2, yy);
}
}
}
yy00 += box_height(&(mat[i]), psr)+box_dh(psr)+box_dh_top(IS_ONCE, psr);
}
}
int ps_main(int argc, char **argv)
{
DIR *d;
struct dirent *de;
char *namefilter = 0;
int pidfilter = 0;
int threads = 0;
d = opendir("/proc");
if(d == 0) return -1;
while(argc > 1){
if(!strcmp(argv[1],"-t")) {
threads = 1;
} else if(!strcmp(argv[1],"-x")) {
display_flags |= SHOW_TIME;
} else if(!strcmp(argv[1], "-Z")) {
display_flags |= SHOW_MACLABEL;
} else if(!strcmp(argv[1],"-P")) {
display_flags |= SHOW_POLICY;
} else if(!strcmp(argv[1],"-p")) {
display_flags |= SHOW_PRIO;
} else if(!strcmp(argv[1],"-c")) {
display_flags |= SHOW_CPU;
} else if(!strcmp(argv[1],"--abi")) {
display_flags |= SHOW_ABI;
} else if(isdigit(argv[1][0])){
pidfilter = atoi(argv[1]);
} else {
namefilter = argv[1];
}
argc--;
argv++;
}
if (display_flags & SHOW_MACLABEL) {
printf("LABEL USER PID PPID NAME\n");
} else {
printf("USER PID PPID VSIZE RSS %s%s %s WCHAN PC %sNAME\n",
(display_flags&SHOW_CPU)?"CPU ":"",
(display_flags&SHOW_PRIO)?"PRIO NICE RTPRI SCHED ":"",
(display_flags&SHOW_POLICY)?"PCY " : "",
(display_flags&SHOW_ABI)?"ABI " : "");
}
while((de = readdir(d)) != 0){
if(isdigit(de->d_name[0])){
int pid = atoi(de->d_name);
if(!pidfilter || (pidfilter == pid)) {
ps_line(pid, 0, namefilter);
if(threads) ps_threads(pid, namefilter);
}
}
}
closedir(d);
return 0;
}
static int
ps_main(OBJ *obj, int argc, char **argv)
{
DIR *d;
char buf[256];
struct dirent *de;
char *namefilter = 0;
int pidfilter = 0;
int threads = 0;
d = opendir("/proc");
if(d == 0) return -1;
while(argc > 1){
if(!strcmp(argv[1],"-t")) {
threads = 1;
} else if(!strcmp(argv[1],"-x")) {
display_flags |= SHOW_TIME;
} else if(!strcmp(argv[1],"-p")) {
display_flags |= SHOW_PRIO;
} else if(isdigit(argv[1][0])){
pidfilter = atoi(argv[1]);
} else {
namefilter = argv[1];
}
argc--;
argv++;
}
PRINTF("USER PID PPID VSIZE RSS WCHAN PC NAME\n");
while((de = readdir(d)) != 0){
if(isdigit(de->d_name[0])){
int pid = atoi(de->d_name);
if(!pidfilter || (pidfilter == pid)) {
ps_line(obj, pid, 0, namefilter);
if(threads) ps_threads(obj, pid, namefilter);
}
}
}
closedir(d);
return 0;
}
void ps_threads(int pid, char *namefilter)
{
char tmp[128];
DIR *d;
struct dirent *de;
sprintf(tmp,"/proc/%d/task",pid);
d = opendir(tmp);
if(d == 0) return;
while((de = readdir(d)) != 0){
if(isdigit(de->d_name[0])){
int tid = atoi(de->d_name);
if(tid == pid) continue;
ps_line(pid, tid, namefilter);
}
}
closedir(d);
}
static void draw_boxes(t_psdata ps,real x0,real y0,real w,
int nmat,t_matrix mat[],t_psrec *psr)
{
char buf[12];
char *mylab;
real xxx;
char **xtick,**ytick;
real xx,yy,dy,xx00,yy00,offset_x,offset_y;
int i,j,x,y,ntx,nty,strlength;
/* Only necessary when there will be no y-labels */
strlength = 0;
/* Draw the box */
ps_rgb(ps,BLACK);
ps_linewidth(ps,psr->boxlinewidth);
yy00=y0;
for(i=0; (i<nmat); i++) {
dy=box_height(&(mat[i]),psr);
ps_box(ps,x0-1,yy00-1,x0+w+1,yy00+dy+1);
yy00+=dy+box_dh(psr)+box_dh_top(IS_ONCE,psr);
}
/* Draw the ticks on the axes */
ps_linewidth(ps,psr->ticklinewidth);
xx00=x0-1;
yy00=y0-1;
for (i=0; (i<nmat); i++) {
if (mat[i].flags & MAT_SPATIAL_X) {
ntx = mat[i].nx + 1;
offset_x = 0.1;
} else {
ntx = mat[i].nx;
offset_x = 0.6;
}
if (mat[i].flags & MAT_SPATIAL_Y) {
nty = mat[i].ny + 1;
offset_y = 0.1;
} else {
nty = mat[i].ny;
offset_y = 0.6;
}
snew(xtick,ntx);
for(j=0; (j<ntx); j++) {
sprintf(buf,"%g",mat[i].axis_x[j]);
xtick[j]=strdup(buf);
}
ps_strfont(ps,psr->X.tickfont,psr->X.tickfontsize);
for(x=0; (x<ntx); x++) {
xx = xx00 + (x + offset_x)*psr->xboxsize;
if ( ( bRmod(mat[i].axis_x[x], psr->X.offset, psr->X.major) ||
(psr->X.first && (x==0))) &&
( (i == 0) || box_do_all_x_maj_ticks(psr) ) ) {
/* Longer tick marks */
ps_line (ps,xx,yy00,xx,yy00-psr->X.majorticklen);
/* Plot label on lowest graph only */
if (i == 0)
ps_ctext(ps,xx,
yy00-DDD-psr->X.majorticklen-psr->X.tickfontsize*0.8,
xtick[x],eXCenter);
} else if ( bRmod(mat[i].axis_x[x], psr->X.offset, psr->X.minor) &&
( (i == 0) || box_do_all_x_min_ticks(psr) ) ){
/* Shorter tick marks */
ps_line(ps,xx,yy00,xx,yy00-psr->X.minorticklen);
} else if ( bRmod(mat[i].axis_x[x], psr->X.offset, psr->X.major) ) {
/* Even shorter marks, only each X.major */
ps_line(ps,xx,yy00,xx,yy00-(psr->boxspacing/2));
}
}
ps_strfont(ps,psr->Y.tickfont,psr->Y.tickfontsize);
snew(ytick,nty);
for(j=0; (j<nty); j++) {
sprintf(buf,"%g",mat[i].axis_y[j]);
ytick[j]=strdup(buf);
}
for(y=0; (y<nty); y++) {
yy = yy00 + (y + offset_y)*psr->yboxsize;
if ( bRmod(mat[i].axis_y[y], psr->Y.offset, psr->Y.major) ||
(psr->Y.first && (y==0))) {
/* Major ticks */
strlength=max(strlength,(int)strlen(ytick[y]));
ps_line (ps,xx00,yy,xx00-psr->Y.majorticklen,yy);
ps_ctext(ps,xx00-psr->Y.majorticklen-DDD,
yy-psr->Y.tickfontsize/3.0,ytick[y],eXRight);
}
else if ( bRmod(mat[i].axis_y[y], psr->Y.offset, psr->Y.minor) ) {
/* Minor ticks */
ps_line(ps,xx00,yy,xx00-psr->Y.minorticklen,yy);
}
}
sfree(xtick);
sfree(ytick);
/* Label on Y-axis */
if (!psr->bYonce || i==nmat/2) {
if (strlen(psr->Y.label) > 0)
mylab = psr->Y.label;
else
mylab = mat[i].label_y;
if (strlen(mylab) > 0) {
ps_strfont(ps,psr->Y.font,psr->Y.fontsize);
ps_flip(ps,TRUE);
xxx=x0-psr->X.majorticklen-psr->X.tickfontsize*strlength-DDD;
ps_ctext(ps,yy00+box_height(&mat[i],psr)/2.0,612.5-xxx,
mylab,eXCenter);
ps_flip(ps,FALSE);
}
}
yy00+=box_height(&(mat[i]),psr)+box_dh(psr)+box_dh_top(IS_ONCE,psr);
}
/* Label on X-axis */
if (strlen(psr->X.label) > 0)
mylab = psr->X.label;
else
mylab = mat[0].label_x;
if (strlen(mylab) > 0) {
ps_strfont(ps,psr->X.font,psr->X.fontsize);
ps_ctext(ps,x0+w/2,y0-DDD-psr->X.majorticklen-psr->X.tickfontsize*FUDGE-
psr->X.fontsize,mylab,eXCenter);
}
}
void zzpsli_( int *ix1 , int *iy1 , int *ix2 , int *iy2 )
{ ps_line( *ix1,*iy1 , *ix2,*iy2 ) ; }
int main (int argc, char **argv)
{
char program_name[] = "MBPS";
char help_message[] = "MBPS reads a swath bathymetry data file and creates a postscript 3-d mesh plot";
char usage_message[] = "mbps [-Iinfile -Fformat -Nnpings -Ppings\n\t-Byr/mo/da/hr/mn/sc -Eyr/mo/da/hr/mn/sc \n\t-Aalpha -Keta -Dviewdir -Xvertexag \n\t-T\"title\" -Wmetersperinch \n\t-Sspeedmin -Ggap -Ydisplay_stats \n\t-Zdisplay_scales -V -H]";
extern char *optarg;
int errflg = 0;
int c;
int help = 0;
int flag = 0;
/*ALBERTO definitions */
int gap=1;
double *xp, *yp;
double xl[4], yl[4];
double alpha = ALPHA_DEF;
double eta = ETA_DEF;
double ve = VE_DEF;
char viewdir = VIEWDIR_DEF;
int display_stats = MB_YES;
int display_scales = MB_YES;
double sin_eta, cos_eta;
double sin_alpha, cos_alpha;
double track_length, xscale, zscale, zscale_inch;
double mean_xp=0.0, mean_yp=0.0, min_xp, max_xp, min_yp, max_yp;
double scaling, x_off, y_off;
double min_z, max_z, range_z, meters_per_inch=(-1.0);
double mean_lat=0.0;
double mean_lon=0.0;
double mean_latmin;
double mean_lonmin;
double mean_hdg=0.0;
int done, mean_knt=0;
int orient;
char label[100];
int a, b, rotate;
double x, y, z;
/* MBIO status variables */
int status = MB_SUCCESS;
int verbose = 0;
int error = MB_ERROR_NO_ERROR;
char *message;
/* MBIO read control parameters */
int format;
int lonflip;
double bounds[4];
int btime_i[7];
int etime_i[7];
double btime_d;
double etime_d;
double speedmin;
double timegap;
char file[MB_PATH_MAXLINE];
int pings = 1;
int beams_bath;
int beams_amp;
int pixels_ss;
int num_pings_max = MBPS_MAXPINGS;
/* MBIO read values */
void *mbio_ptr = NULL;
int kind;
struct ping data[MBPS_MAXPINGS+3];
int time_i[7];
double time_d;
double navlon;
double navlat;
double speed;
double heading;
double distance;
double altitude;
double sonardepth;
char *beamflag;
double *bath = NULL;
double *bathacrosstrack = NULL;
double *bathalongtrack = NULL;
double *amp = NULL;
double *ss = NULL;
double *ssacrosstrack = NULL;
double *ssalongtrack = NULL;
char comment[MB_COMMENT_MAXLINE];
int timbeg_i[7];
int timend_i[7];
double distot = 0.0;
int nread;
char title[MB_COMMENT_MAXLINE];
int forward;
double xx, yy, zz;
double heading_start, dheading, dheadingx, dheadingy;
int i, j, jj, k;
void Polygon_Fill();
void Good_Polygon();
/* initialize some time variables */
for (i=0;i<7;i++)
{
timbeg_i[i] = 0;
timend_i[i] = 0;
}
/* get current default values */
status = mb_defaults(verbose,&format,&pings,&lonflip,bounds,
btime_i,etime_i,&speedmin,&timegap);
/* set default input to stdin */
strcpy (file, "stdin");
/* process argument list */
while ((c = getopt(argc, argv, "VvHhF:f:B:b:E:e:S:s:T:t:I:i:A:a:X:x:K:k:D:d:N:n:P:p:W:w:G:g:YyZz")) != -1)
switch (c)
{
case 'H':
case 'h':
help++;
break;
case 'V':
case 'v':
verbose++;
break;
case 'A':
case 'a':
sscanf (optarg, "%lf", &alpha);
flag++;
break;
case 'B':
case 'b':
sscanf (optarg,"%d/%d/%d/%d/%d/%d",
&btime_i[0],&btime_i[1],&btime_i[2],
&btime_i[3],&btime_i[4],&btime_i[5]);
btime_i[6] = 0;
flag++;
break;
case 'D':
case 'd':
sscanf (optarg, "%c", &viewdir);
flag++;
break;
case 'E':
case 'e':
sscanf (optarg,"%d/%d/%d/%d/%d/%d",
&etime_i[0],&etime_i[1],&etime_i[2],
&etime_i[3],&etime_i[4],&etime_i[5]);
etime_i[6] = 0;
flag++;
break;
case 'F':
case 'f':
sscanf (optarg,"%d", &format);
flag++;
break;
case 'G':
case 'g':
sscanf (optarg, "%d", &gap);
flag++;
break;
case 'I':
case 'i':
sscanf (optarg,"%s", file);
flag++;
break;
case 'K':
case 'k':
sscanf (optarg, "%lf", &eta);
flag++;
break;
case 'N':
case 'n':
sscanf (optarg, "%d", &num_pings_max);
if (num_pings_max < 2
|| num_pings_max > MBPS_MAXPINGS)
num_pings_max = MBPS_MAXPINGS;
flag++;
break;
case 'P':
case 'p':
sscanf (optarg, "%d", &pings);
flag++;
break;
case 'S':
case 's':
sscanf (optarg,"%lf", &speedmin);
flag++;
break;
case 'T':
case 't':
sscanf (optarg,"%s", title);
flag++;
break;
case 'X':
case 'x':
sscanf (optarg, "%lf", &ve);
flag++;
break;
case 'W':
case 'w':
sscanf (optarg, "%lf", &meters_per_inch);
flag++;
break;
case 'Y':
case 'y':
display_stats = MB_NO;
flag++;
break;
case 'Z':
case 'z':
display_scales = MB_NO;
flag++;
break;
case '?':
errflg++;
break;
} /* switch */
/* Process the title of the plot */
for (i = 1; i < argc; i++)
{
if (argv[i][0] == '-'&& ((argv[i][1]=='T')||(argv[i][1]=='t')) ) {
strcpy(title,argv[i]);
title[0]=' ';
title[1]=' ';
}
}
/* check that otions are allowed */
if ((viewdir!='P') && (viewdir!='S') && (viewdir!='B') &&
(viewdir!='p') && (viewdir!='s') && (viewdir!='b'))
{
fprintf(stderr,"viewdir must be either P/p (port) S/s (stbd) or B/b (back)\n");
errflg++;
}
/* if error flagged then print it and exit */
if (errflg)
{
fprintf(stderr,"usage: %s\n", usage_message);
fprintf(stderr,"\nProgram <%s> Terminated\n",
program_name);
error = MB_ERROR_BAD_USAGE;
exit(error);
}
/* print starting message */
if (verbose == 1 || help)
{
fprintf(stderr,"\nProgram %s\n",program_name);
fprintf(stderr,"MB-system Version %s\n",MB_VERSION);
}
/* print starting debug statements */
if (verbose >= 2)
{
fprintf(stderr,"\ndbg2 Program <%s>\n",program_name);
fprintf(stderr,"dbg2 MB-system Version %s\n",MB_VERSION);
fprintf(stderr,"dbg2 Control Parameters:\n");
fprintf(stderr,"dbg2 verbose: %d\n",verbose);
fprintf(stderr,"dbg2 help: %d\n",help);
fprintf(stderr,"dbg2 format: %d\n",format);
fprintf(stderr,"dbg2 btime_i[0]: %d\n",btime_i[0]);
fprintf(stderr,"dbg2 btime_i[1]: %d\n",btime_i[1]);
fprintf(stderr,"dbg2 btime_i[2]: %d\n",btime_i[2]);
fprintf(stderr,"dbg2 btime_i[3]: %d\n",btime_i[3]);
fprintf(stderr,"dbg2 btime_i[4]: %d\n",btime_i[4]);
fprintf(stderr,"dbg2 btime_i[5]: %d\n",btime_i[5]);
fprintf(stderr,"dbg2 btime_i[6]: %d\n",btime_i[6]);
fprintf(stderr,"dbg2 etime_i[0]: %d\n",etime_i[0]);
fprintf(stderr,"dbg2 etime_i[1]: %d\n",etime_i[1]);
fprintf(stderr,"dbg2 etime_i[2]: %d\n",etime_i[2]);
fprintf(stderr,"dbg2 etime_i[3]: %d\n",etime_i[3]);
fprintf(stderr,"dbg2 etime_i[4]: %d\n",etime_i[4]);
fprintf(stderr,"dbg2 etime_i[5]: %d\n",etime_i[5]);
fprintf(stderr,"dbg2 etime_i[6]: %d\n",etime_i[6]);
fprintf(stderr,"dbg2 speedmin: %f\n",speedmin);
fprintf(stderr,"dbg2 file: %s\n",file);
}
/* if help desired then print it and exit */
if (help)
{
fprintf(stderr,"\n%s\n",help_message);
fprintf(stderr,"\nusage: %s\n", usage_message);
exit(error);
}
/* get format if required */
if (format == 0)
mb_get_format(verbose,file,NULL,&format,&error);
/* initialize reading the swath file */
if ((status = mb_read_init(
verbose,file,format,pings,lonflip,bounds,
btime_i,etime_i,speedmin,timegap,
&mbio_ptr,&btime_d,&etime_d,
&beams_bath,&beams_amp,&pixels_ss,&error)) != MB_SUCCESS)
{
mb_error(verbose,error,&message);
fprintf(stderr,"\nMBIO Error returned from function <mb_read_init>:\n%s\n",message);
fprintf(stderr,"\nMultibeam File <%s> not initialized for reading\n",file);
fprintf(stderr,"\nProgram <%s> Terminated\n",
program_name);
exit(error);
}
/* initialize values */
sin_alpha = sin(alpha*DTR);
cos_alpha = cos(alpha*DTR);
sin_eta = sin(eta*DTR);
cos_eta = cos(eta*DTR);
min_z = 0.0;
max_z = -9999.0;
/* allocate memory for data arrays */
beamflag = NULL;
bath = NULL;
bathacrosstrack = NULL;
bathalongtrack = NULL;
amp = NULL;
ss = NULL;
ssacrosstrack = NULL;
ssalongtrack = NULL;
if (error == MB_ERROR_NO_ERROR)
status = mb_register_array(verbose, mbio_ptr, MB_MEM_TYPE_BATHYMETRY,
sizeof(char), (void **)&beamflag, &error);
if (error == MB_ERROR_NO_ERROR)
status = mb_register_array(verbose, mbio_ptr, MB_MEM_TYPE_BATHYMETRY,
sizeof(double), (void **)&bath, &error);
if (error == MB_ERROR_NO_ERROR)
status = mb_register_array(verbose, mbio_ptr, MB_MEM_TYPE_BATHYMETRY,
sizeof(double), (void **)&bathacrosstrack, &error);
if (error == MB_ERROR_NO_ERROR)
status = mb_register_array(verbose, mbio_ptr, MB_MEM_TYPE_BATHYMETRY,
sizeof(double), (void **)&bathalongtrack, &error);
if (error == MB_ERROR_NO_ERROR)
status = mb_register_array(verbose, mbio_ptr, MB_MEM_TYPE_AMPLITUDE,
sizeof(double), (void **)&, &error);
if (error == MB_ERROR_NO_ERROR)
status = mb_register_array(verbose, mbio_ptr, MB_MEM_TYPE_SIDESCAN,
sizeof(double), (void **)&ss, &error);
if (error == MB_ERROR_NO_ERROR)
status = mb_register_array(verbose, mbio_ptr, MB_MEM_TYPE_SIDESCAN,
sizeof(double), (void **)&ssacrosstrack, &error);
if (error == MB_ERROR_NO_ERROR)
status = mb_register_array(verbose, mbio_ptr, MB_MEM_TYPE_SIDESCAN,
sizeof(double), (void **)&ssalongtrack, &error);
for (i=0;i<num_pings_max+3;i++)
{
data[i].beams_bath = 0;
data[i].beamflag = NULL;
data[i].bath = NULL;
data[i].bathacrosstrack = NULL;
data[i].bathalongtrack = NULL;
data[i].xp = NULL;
data[i].yp = NULL;
}
/* if error initializing memory then quit */
if (error != MB_ERROR_NO_ERROR)
{
mb_error(verbose,error,&message);
fprintf(stderr,"\nMBIO Error allocating data arrays:\n%s\n",message);
fprintf(stderr,"\nProgram <%s> Terminated\n",
program_name);
exit(error);
}
/* read and process data */
nread = 0;
done = MB_NO;
error = MB_ERROR_NO_ERROR;
while (done == MB_NO && error <= MB_ERROR_NO_ERROR)
{
/* read a ping of data */
status = mb_get(verbose,mbio_ptr,&kind,&pings,
time_i,&time_d,
&navlon,&navlat,
&speed,&heading,
&distance,&altitude,&sonardepth,
&beams_bath,&beams_amp,&pixels_ss,
beamflag,bath,amp,
bathacrosstrack,
bathalongtrack,
ss,
ssacrosstrack,
ssalongtrack,
comment,&error);
/* only work with survey data */
if (error == MB_ERROR_NO_ERROR && kind == MB_DATA_DATA)
{
/* allocate arrays */
data[nread].beams_bath = beams_bath;
status = mb_mallocd(verbose, __FILE__, __LINE__, beams_bath*sizeof(char), (void **)&(data[nread].beamflag), &error);
status = mb_mallocd(verbose, __FILE__, __LINE__, beams_bath*sizeof(double), (void **)&(data[nread].bath), &error);
status = mb_mallocd(verbose, __FILE__, __LINE__, beams_bath*sizeof(double), (void **)&(data[nread].bathacrosstrack), &error);
status = mb_mallocd(verbose, __FILE__, __LINE__, beams_bath*sizeof(double), (void **)&(data[nread].bathalongtrack), &error);
status = mb_mallocd(verbose, __FILE__, __LINE__, beams_bath*sizeof(double), (void **)&(data[nread].xp), &error);
status = mb_mallocd(verbose, __FILE__, __LINE__, beams_bath*sizeof(double), (void **)&(data[nread].yp), &error);
/* copy data to storage arrays */
for (i=0;i<beams_bath;i++)
{
data[nread].beamflag[i] = beamflag[i];
data[nread].bath[i] = bath[i];
data[nread].bathacrosstrack[i] = bathacrosstrack[i];
data[nread].bathalongtrack[i] = bathalongtrack[i];
data[nread].xp[i] = BAD;
data[nread].yp[i] = BAD;
}
/* ignore time gaps */
if (error == MB_ERROR_TIME_GAP)
{
error = MB_ERROR_NO_ERROR;
status = MB_SUCCESS;
}
/* output error messages */
if (error == MB_ERROR_COMMENT)
{
/* do nothing */
}
else if (verbose >= 1 && error < MB_ERROR_NO_ERROR
&& error >= MB_ERROR_OTHER)
{
mb_error(verbose,error,&message);
fprintf(stderr,"\nNonfatal MBIO Error:\n%s\n",
message);
fprintf(stderr,"Time: %d %d %d %d %d %d %d\n",
time_i[0],time_i[1],time_i[2],
time_i[3],time_i[4],time_i[5],
time_i[6]);
}
else if (verbose >= 1 && error < MB_ERROR_NO_ERROR)
{
mb_error(verbose,error,&message);
fprintf(stderr,"\nNonfatal MBIO Error:\n%s\n",
message);
fprintf(stderr,"Number of good records so far: %d\n",nread);
}
else if (verbose >= 1 && error > MB_ERROR_NO_ERROR
&& error != MB_ERROR_EOF)
{
mb_error(verbose,error,&message);
fprintf(stderr,"\nFatal MBIO Error:\n%s\n",
message);
fprintf(stderr,"Last Good Time: %d %d %d %d %d %d %d\n",
time_i[0],time_i[1],time_i[2],
time_i[3],time_i[4],time_i[5],
time_i[6]);
}
/* calculate raw x,y locations for each beam */
if (status == MB_SUCCESS)
{
/* set initial heading */
if (nread == 0)
heading_start = heading;
/* get heading x and y components */
dheading = heading - heading_start;
if (dheading > 360.0)
dheading -= 360.0;
else if (dheading < 0.0)
dheading += 360.0;
dheadingx = sin(DTR * dheading);
dheadingy = cos(DTR * dheading);
/* get alongtrack distance in nav */
distot += distance * 1000.0; /* distance in meters */
/* loop over the beams */
for (j=0; j<beams_bath; j++)
{
if (j >= data[nread].beams_bath)
{
data[nread].beamflag[j] = MB_FLAG_NULL;
data[nread].xp[j] = BAD;
data[nread].yp[j] = BAD;
}
else if (mb_beam_ok(beamflag[j]))
{
xx = dheadingy * bathacrosstrack[j]
+ dheadingx * bathalongtrack[j];
yy = distot
- dheadingx * bathacrosstrack[j]
+ dheadingy * bathalongtrack[j];
zz = -bath[j];
if (viewdir=='S' || viewdir=='s')
{
data[nread].xp[j] = yy
+ xx * sin_eta * cos_alpha;
data[nread].yp[j] = zz * cos_eta * ve
- xx * sin_eta * sin_alpha;
}
else if (viewdir=='P' || viewdir=='p')
{
data[nread].xp[j]= -yy
- xx * sin_eta * cos_alpha;
data[nread].yp[j]= zz * cos_eta * ve
+ xx * sin_eta * sin_alpha;
}
else if (viewdir=='B' || viewdir=='b')
{
data[nread].xp[j] = xx
+ yy * sin_eta * cos_alpha;
data[nread].yp[j]= zz * cos_eta * ve
+ yy * sin_eta * sin_alpha;
}
mean_lat += navlat;
mean_lon += navlon;
mean_hdg += heading;
mean_xp += data[nread].xp[j];
mean_yp += data[nread].yp[j];
mean_knt++;
if (-data[nread].bath[j] < min_z)
min_z= -data[nread].bath[j];
if (-data[nread].bath[j] > max_z)
max_z= -data[nread].bath[j];
}
else
{
data[nread].xp[j] = BAD;
data[nread].yp[j] = BAD;
}
} /* for j=0 ... */
if (nread == 0)
{
for (k=0; k<7; k++)
timbeg_i[k] = time_i[k];
}
else
{
for (k=0; k<7; k++)
timend_i[k]=time_i[k];
}
} /* if status==MB_SUCCESS */
/* increment counters */
if (error == MB_ERROR_NO_ERROR)
{
nread++;
}
}
/* print debug statements */
if (verbose >= 2)
{
fprintf(stderr,"\ndbg2 Reading loop finished in program <%s>\n",
program_name);
fprintf(stderr,"dbg2 status: %d\n",status);
fprintf(stderr,"dbg2 error: %d\n",error);
fprintf(stderr,"dbg2 nread: %d\n",nread);
fprintf(stderr,"dbg2 pings: %d\n",pings);
}
/* test if done */
if (nread >= num_pings_max
&& verbose >= 1)
{
fprintf(stderr, "%s: Maximum number of pings [%d] read before end of file reached...\n",
program_name, num_pings_max);
done = MB_YES;
}
if (nread >= num_pings_max || error > MB_ERROR_NO_ERROR)
{
done = MB_YES;
}
} /* end of processing data, 1'st while under read/process data */
/* close the swath file */
status = mb_close(verbose,&mbio_ptr,&error);
/* print debug statements */
if (verbose >= 2)
{
fprintf(stderr,"\ndbg2 Reading loop finished in program <%s>\n",
program_name);
fprintf(stderr,"dbg2 status: %d\n",status);
fprintf(stderr,"dbg2 error: %d\n",error);
fprintf(stderr,"dbg2 nread: %d\n",nread);
fprintf(stderr,"dbg2 pings: %d\n",pings);
}
/* total track length in m */
track_length = distot;
mean_lat /= mean_knt;
mean_latmin = fabs(mean_lat - (int) mean_lat) * 60.0;
mean_lon /= mean_knt;
mean_lonmin = fabs(mean_lon - (int) mean_lon) * 60.0;
mean_hdg /= mean_knt;
mean_xp /= mean_knt;
mean_yp /= mean_knt;
/* rescale xp[],yp[] to zero mean; get min and max */
max_yp = min_yp = max_xp = min_xp = 0.0;
for (i=0; i<nread; i++)
{
beamflag = data[i].beamflag;
xp = data[i].xp;
yp = data[i].yp;
for (j=0; j<data[i].beams_bath; j++)
{
if (mb_beam_ok(beamflag[j]))
{
yp[j] -= mean_yp;
xp[j] -= mean_xp;
min_xp = MIN(min_xp, xp[j]);
max_xp = MAX(max_xp, xp[j]);
min_yp = MIN(min_yp, yp[j]);
max_yp = MAX(max_yp, yp[j]);
} /* if yp[][] */
} /* for j */
} /* for i */
/* get page orientation, scaling(in/m) factor and startup plot */
if ((viewdir=='P') || (viewdir=='S') || (viewdir=='p') || (viewdir=='s'))
{
/* Landscape */
orient = 0;
if (meters_per_inch > 0.0)
{
scaling = 1.0 / meters_per_inch;
x_off = 11. / 2;
y_off = 8.5 / 2.;
}
else
{
if ( (5.2 / (max_yp - min_yp)) < (8.5 / (max_xp - min_xp)) )
scaling = (5.2 / (max_yp - min_yp));
else
scaling = (8.5 / (max_xp - min_xp));
x_off=(-(max_xp + min_xp) * scaling / 2.0) + (11. / 2);
y_off=(-(max_yp + min_yp) * scaling / 2.0) + (8.5 / 2) - .2;
}
}
else
{
/* Portrait */
orient = 1;
if (meters_per_inch > 0.0)
{
scaling = 1.0 / meters_per_inch;
x_off = 8.5 / 2.0;
y_off = 11. / 2.0;
}
else
{
if ( (8./(max_yp-min_yp))<(6.5/(max_xp-min_xp)) )
scaling = (8./(max_yp-min_yp));
else
scaling = (6.5/(max_xp-min_xp));
x_off=(-(max_xp+min_xp)*scaling/2.0)+(8.5/2);
y_off=(-(max_yp+min_yp)*scaling/2.0)+(11./2)-.2;
}
}
/* initialize the Postscript plotting */
#ifdef GMT_MINOR_VERSION
ps_plotinit_hires(NULL,0,orient,x_off,y_off,1.0,1.0,1,300,1,
gmtdefs.paper_width, gmtdefs.page_rgb,
gmtdefs.encoding.name,
GMT_epsinfo (argv[0]));
#else
ps_plotinit(NULL,0,orient,x_off,y_off,1.0,1.0,1,300,1,
gmtdefs.paper_width, gmtdefs.page_rgb,
gmtdefs.encoding.name,
GMT_epsinfo (argv[0]));
#endif
GMT_echo_command (argc, argv);
/* now loop over the data in the appropriate order
laying down white filled boxes with black outlines
wherever the data is good */
if ((viewdir=='S') || (viewdir=='s'))
forward = MB_YES;
else if ((viewdir=='P') || (viewdir=='p'))
forward = MB_NO;
else if ((viewdir=='B') || (viewdir=='b'))
{
if (alpha < 90.0)
forward = MB_YES;
else
forward = MB_NO;
}
for (j=0;j<beams_bath-1;j++)
{
for (i=0;i<nread-1;i++)
{
if (forward == MB_YES)
jj = j;
else
jj = beams_bath - 2 - j;
/* make box */
if (mb_beam_ok(data[i].beamflag[jj])
&& mb_beam_ok(data[i+1].beamflag[jj])
&& mb_beam_ok(data[i].beamflag[jj+1])
&& mb_beam_ok(data[i+1].beamflag[jj+1]))
{
xl[0] = scaling * data[i].xp[jj];
yl[0] = scaling * data[i].yp[jj];
xl[1] = scaling * data[i+1].xp[jj];
yl[1] = scaling * data[i+1].yp[jj];
xl[2] = scaling * data[i+1].xp[jj+1];
yl[2] = scaling * data[i+1].yp[jj+1];
xl[3] = scaling * data[i].xp[jj+1];
yl[3] = scaling * data[i].yp[jj+1];
ps_polygon(xl,yl,4,rgb_white,1);
}
}
}
/* titles and such */
ps_setline(2); /* set line width */
if (display_stats == MB_NO)
{
/* plot a title */
xl[0]=0;
yl[0]=max_yp*scaling+.6;
sprintf(label,"%s",title);
ps_text(xl[0],yl[0],20.,label,0.,6,0);
}
else
{
/* plot a title */
xl[0]=0;
yl[0]=max_yp*scaling+1.3;
sprintf(label,"%s",title);
ps_text(xl[0],yl[0],20.,label,0.,6,0);
/*xl[0]-=3.25;*/
yl[0]-=0.3;
sprintf(label,"Mean Lat.: %3d@+o@+ %4.1f' Mean Lon.: %4d@+o @+%4.1f' Heading: %.1lf@+o @+",(int)mean_lat, mean_latmin, (int)mean_lon, mean_lonmin, mean_hdg);
ps_text(xl[0],yl[0],15.,label,0.,6,0);
yl[0]-=0.3;
sprintf(label,"View Angle: %.1lf@+o @+ V.E.: %.1lfX Scale: %.0lf m/inch Track Length: %.1lf km",eta,ve,1.0/scaling,track_length/1000.0);
ps_text(xl[0],yl[0],15.,label,0.,6,0);
yl[0]-=0.3;
sprintf(label,
"From %.4d/%.2d/%.2d %.2d:%.2d:%.2d to %.4d/%.2d/%.2d %.2d:%.2d:%.2d",
timbeg_i[0],timbeg_i[1],timbeg_i[2],timbeg_i[3],
timbeg_i[4],timbeg_i[5],timend_i[0],timend_i[1],
timend_i[2],timend_i[3],timend_i[4],timend_i[5]);
ps_text(xl[0],yl[0],15.,label,0.,6,0);
} /* else after if display_stats */
if (display_scales == MB_YES)
{
/* plot the x-scale */
xscale=10000; /* x scale in m */
if (track_length < 50000) xscale=5000;
if (track_length < 20000) xscale=2000;
if (track_length < 10000) xscale=1000;
xl[0]=xl[1]= (-xscale*scaling/2.0);
xl[2]=xl[3]= (-xl[0]);
xl[0]+=2.;xl[1]+=2.;xl[2]+=2.;xl[3]+=2.;
yl[1]=yl[2]= min_yp*scaling-1.;
yl[0]=yl[3]= yl[1]+0.1;
#ifdef GMT_MINOR_VERSION
ps_line(xl,yl,4,3,0);
#else
ps_line(xl,yl,4,3,0,0);
#endif
sprintf(label,"%.0f km",xscale/1000.0);
ps_text(xl[0]+.5,yl[0]+.05,15.,label,0.,6,0);
/* plot the z-scale */
range_z=(max_z-min_z);
zscale=2000; /* z scale in m */
if (range_z < 3000) zscale=1000;
if (range_z < 1000) zscale=500;
if (range_z < 500) zscale=200;
if (range_z < 250) zscale=100;
zscale_inch= zscale*scaling*cos_eta*ve;
xl[1]=xl[2]+0.5;
xl[2]=xl[1];
xl[0]=xl[3]= xl[1]+.1;
yl[0]=yl[1]= min_yp*scaling-1.;
yl[2]=yl[3]= yl[0]+zscale_inch;
#ifdef GMT_MINOR_VERSION
ps_line(xl,yl,4,3,0);
#else
ps_line(xl,yl,4,3,0,0);
#endif
sprintf(label,"%.0f m",zscale);
ps_text(xl[0]+0.3,yl[0]+zscale_inch/2.0,15.,label,0.,6,0);
/* plot an arrow in the ship's direction */
a=0;
b=beams_bath/2;
while (!mb_beam_ok(data[a++].beamflag[b])) {}
xl[0] = data[--a].xp[b];
yl[0] = data[a].yp[b];
a = nread - 1;
while (!mb_beam_ok(data[a--].beamflag[b])) {}
xl[1] = data[++a].xp[b];
yl[1] = data[a].yp[b];
xl[1] = ((xl[1]-xl[0])/distot/2)+.6;
yl[1] = ((yl[1]-yl[0])/distot/2) + min_yp*scaling-1.;
xl[0] = 0.+.6;
yl[0] = 0.+min_yp*scaling-0.85;
ps_vector(xl[0],yl[0],xl[1],yl[1],
0.01,0.25,0.1,1.0,rgb_black,0);
ps_text(xl[0]-1.7,yl[0]+.2,15.,"ship heading",0.,1,0);
ps_text(xl[0]-1.7,yl[0],15.,"direction",0.,1,0);
/* plot the three axes */
for (i=0;i<3;i++)
{
xl[0]=0.; /* point in center of page */
yl[0]=0.;
rotate=0; /* set to 1 if arrow is rotated below */
if (i==0)
{
/* x-axis */
x=1.;
y=0;
z=0;
}
else if (i==1)
{
/* y-axis */
x=0;
y=1.;
z=0;
}
else if (i==2)
{
/* z-axis */
x=0;
y=0;
z=-1.;
}
if (viewdir=='P' || viewdir=='p')
{
xl[1]=-y-x*sin_eta*cos_alpha+xl[0];
yl[1]= -z*cos_eta+x*sin_eta*sin_alpha+yl[0];
}
else if (viewdir=='B' || viewdir=='b')
{
xl[1]=(x+y*sin_eta*cos_alpha)+xl[0];
yl[1]=-z*cos_eta+y*sin_eta*sin_alpha+yl[0];
}
else if (viewdir=='S' || viewdir=='s')
{
xl[1]=y+x*sin_eta*cos_alpha+xl[0];
yl[1]=z*cos_eta-x*sin_eta*sin_alpha+yl[0];
}
if (yl[1]<yl[0])
{
/* rotate arrows 180 if facing downward */
xl[1]=-xl[1];
yl[1]=-yl[1];
rotate=1;
}
xl[0]=(-3.); /* move arrows from center to lower left corner */
yl[0]=(min_yp*scaling-1.);
xl[1]=xl[0]+xl[1];
yl[1]=yl[0]+yl[1];
ps_vector(xl[0],yl[0],xl[1],yl[1],
0.01,0.25,0.1,1.0,rgb_black,0);
if (i==0&&rotate==0)
ps_text(xl[1],yl[1]+.15,15.,"x",0.,6,0);
else if (i==1&&rotate==0)
ps_text(xl[1],yl[1]+.15,15.,"y",0.,6,0);
else if (i==2&&rotate==0)
ps_text(xl[1],yl[1]+.15,15.,"z",0.,6,0);
else if (i==0&&rotate==1)
ps_text(xl[1],yl[1]+.15,15.,"-x",0.,6,0);
else if (i==1&&rotate==1)
ps_text(xl[1],yl[1]+.15,15.,"-y",0.,6,0);
else if (i==2&&rotate==1)
ps_text(xl[1],yl[1]+.15,15.,"z",0.,6,0);
} /* (i=0;i<3;i++) */
} /* if display_scales */
/* end the postscript file */
ps_plotend(1);
/* deallocate arrays */
for (i=0;i<nread;i++)
{
mb_freed(verbose,__FILE__, __LINE__, (void **)&(data[i].beams_bath), &error);
mb_freed(verbose,__FILE__, __LINE__, (void **)&(data[i].bath), &error);
mb_freed(verbose,__FILE__, __LINE__, (void **)&(data[i].bathacrosstrack), &error);
mb_freed(verbose,__FILE__, __LINE__, (void **)&(data[i].bathalongtrack), &error);
mb_freed(verbose,__FILE__, __LINE__, (void **)&(data[i].xp), &error);
mb_freed(verbose,__FILE__, __LINE__, (void **)&(data[i].yp), &error);
}
/* check memory */
if (verbose >= 4)
status = mb_memory_list(verbose,&error);
/* print output debug statements */
if (verbose >= 2)
{
fprintf(stderr,"\ndbg2 Program <%s> completed\n",
program_name);
fprintf(stderr,"dbg2 Ending status:\n");
fprintf(stderr,"dbg2 status: %d\n",status);
}
/* end it all */
exit(error);
} /* main */
int ps_main(int argc, char **argv)
{
DIR *d;
struct dirent *de;
int pidfilter = 0;
int threads = 0;
while(argc > 1){
if(!strcmp(argv[1],"-t")) {
threads = 1;
} else if(!strcmp(argv[1],"-n")) {
display_flags |= SHOW_NUMERIC_UID;
} else if(!strcmp(argv[1],"-x")) {
display_flags |= SHOW_TIME;
} else if(!strcmp(argv[1], "-Z")) {
display_flags |= SHOW_MACLABEL;
} else if(!strcmp(argv[1],"-P")) {
display_flags |= SHOW_POLICY;
} else if(!strcmp(argv[1],"-p")) {
display_flags |= SHOW_PRIO;
} else if(!strcmp(argv[1],"-c")) {
display_flags |= SHOW_CPU;
} else if(!strcmp(argv[1],"--abi")) {
display_flags |= SHOW_ABI;
} else if(!strcmp(argv[1],"--ppid")) {
ppid_filter = atoi(argv[2]);
if (ppid_filter == 0) {
/* Bug 26554285: Use printf because some apps require at least
* one line of output to stdout even for errors.
*/
printf("bad ppid '%s'\n", argv[2]);
return 1;
}
argc--;
argv++;
} else {
pidfilter = atoi(argv[1]);
if (pidfilter == 0) {
/* Bug 26554285: Use printf because some apps require at least
* one line of output to stdout even for errors.
*/
printf("bad pid '%s'\n", argv[1]);
return 1;
}
}
argc--;
argv++;
}
if (display_flags & SHOW_MACLABEL) {
printf("LABEL ");
}
printf("USER PID PPID VSIZE RSS %s%s %sWCHAN %*s %sNAME\n",
(display_flags&SHOW_CPU)?"CPU ":"",
(display_flags&SHOW_PRIO)?"PRIO NICE RTPRI SCHED ":"",
(display_flags&SHOW_POLICY)?"PCY " : "",
(int) PC_WIDTH, "PC",
(display_flags&SHOW_ABI)?"ABI " : "");
d = opendir("/proc");
if(d == 0) return -1;
while((de = readdir(d)) != 0){
if(isdigit(de->d_name[0])){
int pid = atoi(de->d_name);
if(!pidfilter || (pidfilter == pid)) {
ps_line(pid, 0);
if(threads) ps_threads(pid);
}
}
}
closedir(d);
return 0;
}