int main(int argc,char *argv[]) {
FILE *fp;
struct RadarParm prm;
struct RawData raw;
fp=fopen(argv[1],"r");
if (fp==NULL) {
fprintf(stderr,"File not found.\n");
exit(-1);
}
while(RawFread(fp,&prm,&raw) !=-1) {
fprintf(stderr,"%.4d-%.2d-%.2d %.2d:%.2d:%.2d\n",
prm.time.yr,prm.time.mo,prm.time.dy,
prm.time.hr,prm.time.mt,prm.time.sc);
prm.cp=1000;
RawFwrite(stdout,&prm,&raw);
}
fclose(fp);
return 0;
}
int main(int argc,char *argv[])
{
/* File format transistion
* ------------------------
*
* When we switch to the new file format remove any reference
* to "new". Change the command line option "new" to "old" and
* remove "old=!new".
*/
char *envstr;
int status;
int arg;
unsigned char help=0;
unsigned char option=0;
unsigned char vb=0;
FILE *fp=NULL;
struct OldRawFp *rawfp=NULL;
int irec=1;
int drec=2;
int dnum=0;
time_t ctime;
int c,n;
char command[128];
char tmstr[40];
prm=RadarParmMake();
raw=RawMake();
OptionAdd(&opt,"-help",'x',&help);
OptionAdd(&opt,"-option",'x',&option);
OptionAdd(&opt,"vb",'x',&vb);
arg=OptionProcess(1,argc,argv,&opt,NULL);
if (help==1)
{
OptionPrintInfo(stdout,hlpstr);
exit(0);
}
if (option==1)
{
OptionDump(stdout,&opt);
exit(0);
}
envstr=getenv("SD_RADAR");
if (envstr==NULL)
{
fprintf(stderr,"Environment variable 'SD_RADAR' must be defined.\n");
exit(-1);
}
fp=fopen(envstr,"r");
if (fp==NULL)
{
fprintf(stderr,"Could not locate radar information file.\n");
exit(-1);
}
network=RadarLoad(fp);
fclose(fp);
if (network==NULL)
{
fprintf(stderr,"Failed to read radar information.\n");
exit(-1);
}
envstr=getenv("SD_HDWPATH");
if (envstr==NULL)
{
fprintf(stderr,"Environment variable 'SD_HDWPATH' must be defined.\n");
exit(-1);
}
RadarLoadHardware(envstr,network);
if (arg==argc) fp=stdin;
else fp=fopen(argv[arg],"r");
if (fp==NULL)
{
fprintf(stderr,"File not found.\n");
exit(-1);
}
status=RawFread(fp,prm,raw);
radar=RadarGetRadar(network,prm->stid);
if (radar==NULL)
{
fprintf(stderr,"Failed to get radar information.\n");
exit(-1);
}
site=RadarYMDHMSGetSite(radar,prm->time.yr,prm->time.mo,
prm->time.dy,prm->time.hr,prm->time.mt,
prm->time.sc);
if (site==NULL)
{
fprintf(stderr,"Failed to get site information.\n");
exit(-1);
}
command[0]=0;
n=0;
for (c=0;c<argc;c++)
{
n+=strlen(argv[c])+1;
if (n>127) break;
if (c !=0) strcat(command," ");
strcat(command,argv[c]);
}
if (vb)
fprintf(stderr,"%d-%d-%d %d:%d:%d beam=%d\n",prm->time.yr,prm->time.mo,
prm->time.dy,prm->time.hr,prm->time.mt,prm->time.sc,prm->bmnum);
do
{
ctime = time((time_t) 0);
RadarParmSetOriginCommand(prm,command);
strcpy(tmstr,asctime(gmtime(&ctime)));
tmstr[24]=0;
RadarParmSetOriginTime(prm,tmstr);
prm->lagfr = 2400;
prm->frang = 360;
status=RawFwrite(stdout,prm,raw);
status=RawFread(fp,prm,raw);
if (vb)
fprintf(stderr,"%d-%d-%d %d:%d:%d beam=%d\n",prm->time.yr,prm->time.mo,
prm->time.dy,prm->time.hr,prm->time.mt,prm->time.sc,prm->bmnum);
} while (status==0);
return 0;
}
int main(int argc,char *argv[]) {
FILE *fp;
char *envstr;
int s;
envstr=getenv("SD_RADAR");
if (envstr==NULL) {
fprintf(stderr,"Environment variable 'SD_RADAR' must be defined.\n");
exit(-1);
}
fp=fopen(envstr,"r");
if (fp==NULL) {
fprintf(stderr,"Could not locate radar information file.\n");
exit(-1);
}
network=RadarLoad(fp);
fclose(fp);
if (network==NULL) {
fprintf(stderr,"Failed to read radar information.\n");
exit(-1);
}
envstr=getenv("SD_HDWPATH");
if (envstr==NULL) {
fprintf(stderr,"Environment variable 'SD_HDWPATH' must be defined.\n");
exit(-1);
}
RadarLoadHardware(envstr,network);
fp=fopen(argv[1],"r");
if (fp==NULL) {
fprintf(stderr,"File not found.\n");
exit(-1);
}
s=RawFread(fp,&prm,&raw);
if (s==-1) {
fprintf(stderr,"Error reading file.\n");
exit(-1);
}
radar=RadarGetRadar(network,prm.stid);
if (radar==NULL) {
fprintf(stderr,"Failed to get radar information.\n");
exit(-1);
}
site=RadarYMDHMSGetSite(radar,prm.time.yr,prm.time.mo,
prm.time.dy,prm.time.hr,prm.time.mt,
prm.time.sc);
if (site==NULL) {
fprintf(stderr,"Failed to get site information.\n");
exit(-1);
}
FitACFStart(site,prm.time.yr,&fblk);
do {
fprintf(stderr,"%.4d-%.2d-%.2d %.2d:%.2d:%.2d\n",
prm.time.yr,prm.time.mo,prm.time.dy,
prm.time.hr,prm.time.mt,prm.time.sc);
FitACF(&prm,&raw,&fblk,&fit);
s=FitFwrite(stdout,&prm,&fit);
s=RawFread(fp,&prm,&raw);
} while (s !=-1);
FitACFEnd(&fblk);
fclose(fp);
return 0;
}
int main(int argc,char *argv[]) {
unsigned char help=0;
unsigned char option=0;
unsigned char vb=0;
int arg=0;
int s;
int thr=-1;
struct RadarParm rprm;
struct RawData rawacf;
FILE *fp;
int rtab[ORIG_MAX_RANGE];
float snr[MAX_RANGE];
int inx,l,step;
float maxval;
int cnt=0;
int recnum=0;
char vstring[256];
OptionAdd(&opt,"-help",'x',&help);
OptionAdd(&opt,"-option",'x',&option);
OptionAdd(&opt,"vb",'x',&vb);
OptionAdd(&opt,"t",'i',&thr);
arg=OptionProcess(1,argc,argv,&opt,NULL);
if (help==1) {
OptionPrintInfo(stdout,hlpstr);
exit(0);
}
if (option==1) {
OptionDump(stdout,&opt);
exit(0);
}
if (arg !=argc) {
fp=fopen(argv[arg],"r");
if (fp==NULL) {
fprintf(stderr,"Error opening file.\n");
exit(-1);
}
} else fp=stdin;
s=RawFread(fp,&rprm,&rawacf);
if (s==-1) {
fprintf(stderr,"Error reading file.\n");
exit(-1);
}
if (thr !=-1) rawacf.thr=thr;
sprintf(vstring,"%d.%.3d",rawacf.revision.major,rawacf.revision.minor);
/*
if (OldRawHeaderFwrite(stdout,"rawwrite",vstring,rawacf.thr,
"example file") !=0) {
fprintf(stderr,"Error writing header.\n");
exit(-1);
}
*/
do {
recnum++;
if (thr !=-1) rawacf.thr=thr;
for (l=0;l<rprm.nrang;l++) {
if (rprm.noise.search>0)
snr[l]=rawacf.pwr0[l]/rprm.noise.search;
else snr[l]=0;
}
/*
if (rprm.nrang>ORIG_MAX_RANGE) {
step=rprm.nrang/ORIG_MAX_RANGE;
for (l=0;l<ORIG_MAX_RANGE;l++) {
maxval=0.0;
inx=l*step;
for (s=0;s<step;s++) {
if (snr[l*step+s]>maxval) {
maxval=snr[l*step+s];
inx=l*step+s;
}
}
rtab[l]=inx;
}
rprm.rsep=rprm.rsep*step;
s=OldRawFwrite(stdout,"rawwrite",&rprm,&rawacf,recnum,rtab);
} else s=OldRawFwrite(stdout,"rawwrite",&rprm,&rawacf,recnum,NULL);
if (s !=0) {
cnt=-1;
break;
}
*/
/*
rprm.stid=(int16)40;
*/
/*
fprintf(stderr,"%d %d\n",rprm.lagfr,rprm.frang);
*/
if (rprm.lagfr==600 && rprm.smsep==300)rprm.lagfr=(int16)1200;
if (rprm.lagfr==400 && rprm.smsep==100)rprm.lagfr=(int16)800;
s=RawFwrite(stdout,&rprm,&rawacf);
if (s==-1) {
cnt=-1;
break;
}
if (vb) fprintf(stderr,"%.4d-%.2d-%.2d %.2d:%.2d:%.2d\n",rprm.time.yr,
rprm.time.mo,rprm.time.dy,rprm.time.hr,rprm.time.mt,
rprm.time.sc);
cnt++;
} while ((s=RawFread(fp,&rprm,&rawacf)) !=-1);
if (cnt==-1) exit(EXIT_FAILURE);
exit(EXIT_SUCCESS);
return 0;
}
int main(int argc,char *argv[]) {
FILE *fp;
char *envstr;
int s,i,j,n;
envstr=getenv("SD_RADAR");
if (envstr==NULL) {
fprintf(stderr,"Environment variable 'SD_RADAR' must be defined.\n");
exit(-1);
}
fp=fopen(envstr,"r");
if (fp==NULL) {
fprintf(stderr,"Could not locate radar information file.\n");
exit(-1);
}
network=RadarLoad(fp);
fclose(fp);
if (network==NULL) {
fprintf(stderr,"Failed to read radar information.\n");
exit(-1);
}
envstr=getenv("SD_HDWPATH");
if (envstr==NULL) {
fprintf(stderr,"Environment variable 'SD_HDWPATH' must be defined.\n");
exit(-1);
}
RadarLoadHardware(envstr,network);
fp=fopen(argv[1],"r");
if (fp==NULL) {
fprintf(stderr,"File not found.\n");
exit(-1);
}
s=RawFread(fp,&prm,&raw);
if (s==-1) {
fprintf(stderr,"Error reading file.\n");
exit(-1);
}
radar=RadarGetRadar(network,prm.stid);
if (radar==NULL) {
fprintf(stderr,"Failed to get radar information.\n");
exit(-1);
}
site=RadarYMDHMSGetSite(radar,prm.time.yr,prm.time.mo,
prm.time.dy,prm.time.hr,prm.time.mt,
prm.time.sc);
if (site==NULL) {
fprintf(stderr,"Failed to get site information.\n");
exit(-1);
}
FitACFStart(site,prm.time.yr,&fblk);
do {
fprintf(stdout,"%.4d-%.2d-%.2d %.2d:%.2d:%.2d\n",
prm.time.yr,prm.time.mo,prm.time.dy,
prm.time.hr,prm.time.mt,prm.time.sc);
fblk.prm.xcf=prm.xcf;
fblk.prm.tfreq=prm.tfreq;
fblk.prm.noise=prm.noise.search;
fblk.prm.nrang=prm.nrang;
fblk.prm.smsep=prm.smsep;
fblk.prm.nave=prm.nave;
fblk.prm.mplgs=prm.mplgs;
fblk.prm.mpinc=prm.mpinc;
fblk.prm.txpl=prm.txpl;
fblk.prm.lagfr=prm.lagfr;
fblk.prm.mppul=prm.mppul;
fblk.prm.bmnum=prm.bmnum;
fblk.prm.cp=prm.cp;
fblk.prm.channel=prm.channel;
fblk.prm.offset=prm.offset; /* stereo offset */
for (i=0;i<fblk.prm.mppul;i++) fblk.prm.pulse[i]=prm.pulse[i];
for (i=0;i<fblk.prm.mplgs;i++) {
fblk.prm.lag[0][i]=prm.lag[i][0];
fblk.prm.lag[1][i]=prm.lag[i][1];
}
for (i=0;i<fblk.prm.nrang;i++) {
fblk.prm.pwr0[i]=raw.pwr0[i];
for (j=0;j<fblk.prm.mplgs;j++) {
fblk.acfd[i][j].x=raw.acfd[i][j][0];
fblk.acfd[i][j].y=raw.acfd[i][j][1];
}
if (fblk.prm.xcf) {
for (j=0;j<fblk.prm.mplgs;j++) {
fblk.xcfd[i][j].x=raw.xcfd[i][j][0];
fblk.xcfd[i][j].y=raw.xcfd[i][j][1];
}
}
}
FitACFBadlags(&fblk.prm,&bsamp);
fprintf(stdout,"%d:",bsamp.nbad);
for (n=0;n<bsamp.nbad;n++) fprintf(stdout,"%d ",bsamp.badsmp[n]);
fprintf(stdout,"\n");
s=RawFread(fp,&prm,&raw);
} while (s !=-1);
fclose(fp);
return 0;
}
