コード例 #1
0
ファイル: holonDualTP.c プロジェクト: nimeshpatel/smaSoftware
int main(int argc, char **argv)
{
    FILE       *fp_raw_data, *fp_raw_data2;
    union 	ad_union amp_vvm[10], pha_vvm[10];
    union 	ad_union amp_vvm2[10], pha_vvm2[10];
    char	command[35],file_name[80],file_name2[80],*cfreq,c;
    int		j, k, l, ant_scn, ant_scn2, ref_ant, n_size, n_start, offsetunit=0;
    int antennas[11] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
    long	i, i_delay, n_delay;
    int 	adc0_fd, adc1_fd, adc0_rd, adc1_rd, count;
    int 	adc2_fd, adc3_fd, adc2_rd, adc3_rd;
    int 	OFLAG;
    int		istatus,icount;
    int		amp_read[2000],phase_read[2000],truetime2ms[20000];
    int		amp_read2[2000],phase_read2[2000];
    int         ant;
    time_t      start_time,cur_time,Start_time,time1[2000],time2[2000],time3[2000],time4[2000];
    float       az_step, el_step, azoff, eloff, freq,az_bore,el_bore, az_bore2, el_bore2; 
    float       az_off,el_off, az_ref, el_ref;
    float       contDet_a1[2], contDet_a2[2], contDet_ar[2];
    float       contDet1_a1, contDet1_a2, contDet1_ar;
    double      az_enc,el_enc,az_enc2, el_enc2;
    double      az_read[2000],el_read[2000], az_read2[2000], el_read2[2000];
    int rm_status, antlist[RM_ARRAY_SIZE];/*variables for encoderClient*/
    short	currentAzoff, currentAzoff2;
    short	handshake;
    struct tm	ts;	
    double      currentAzoffD;
    unsigned short m1009,m1010;
    smapoptContext optCon;
    struct  smapoptOption optionsTable[] = {
	    {"help",'h',SMAPOPT_ARG_NONE,0,'h'},
	    {"primary antenna",'a',SMAPOPT_ARG_INT,&ant_scn,'a'}, 
	    {"secondary antenna",'d',SMAPOPT_ARG_INT,&ant_scn2,'d'}, 
/*	    {"antennas", 'a', SMAPOPT_ARG_ANTENNAS, &antennas[0], 0, "Specify a list of antennas to move"},*/
	    {"refant",'r',SMAPOPT_ARG_INT,&ref_ant,'r'},
	    {"frequency",'f',SMAPOPT_ARG_STRING,&cfreq,'f'},
	    {"map_size",'m',SMAPOPT_ARG_INT,&n_size,'m'},
	    {"scan_speed",'s',SMAPOPT_ARG_INT,&offsetunit,'s'},
	    {"line_start",'l',SMAPOPT_ARG_INT,&n_start,'l'},
	    {NULL,0,0,NULL,0}
    };
/*     ant_scn=0;
     ant_scn2=0;
     for (ant = 0; ant < N_ANTENNAS; ant++){
	printf ("antennas[ant] %d\n", antennas[ant]);
     }
     for (ant = 0; ant < N_ANTENNAS; ant++){
        if (antennas[ant]){
          printf("%d ", ant);
        if (ant_scn == 0){
        ant_scn = ant;
	}
	else {
	ant_scn2 = ant;
	}
	}
	}
*/
	printf ("antennas: %d %d \n", ant_scn, ant_scn2);

/**********get command line arguments and initializations**************/

	count=2; /* set how many counts we read from A/D board device */
	n_start=0; /* start from the beginning as default setting */

	if(argc<7) usage(-1,"Insufficient number of arguments","");
	optCon = smapoptGetContext("holon", argc, argv, optionsTable,0);

	while ((c = smapoptGetNextOpt(optCon)) >= 0) {
	switch(c) {
		case 'h':
		usage(0,NULL,NULL);
		break;
	}
	}
	if(c<-1) {
	fprintf(stderr, "%s: %s\n",
		smapoptBadOption(optCon, SMAPOPT_BADOPTION_NOALIAS),
		smapoptStrerror(c));
	}

	if(ant_scn<1||ant_scn>8)
	{
	printf("ant should be 1-8\n");
	exit(0);
	}

	if(ref_ant<1||ref_ant>8)
	{
	printf("reference ant should be 1-8\n");
	exit(0);
	}

	freq=atof(cfreq);
	if(freq<50.0||freq>700.0)
	{
	printf("200 < freq < 700\n");
	exit(0);
	}

	if(n_size<n_start)
	{
	printf("n_start is larger than n_size.\n");
	exit(0);
	}

	if(offsetunit<1)
	{
	printf("scan_speed should be larger than 1.\n");
	exit(0);
	}

	if(freq < 500.0) n_delay=(long) 8e7/offsetunit;
	else n_delay=(long) 4e7/offsetunit;
	az_step=7734.0/freq;el_step=az_step; /*23.3 for 332 GHz, and 33.3 for 232.4 GHz*/
printf("ant_scn=%d ref_ant=%d freq=%f size=%d scanspeed=%d n_delay=%ld step=%f\n",ant_scn,ref_ant,freq,n_size,offsetunit,n_delay,az_step);
/* Dual */
printf("ant_scn=%d ref_ant=%d freq=%f size=%d scanspeed=%d n_delay=%ld step=%f\n",ant_scn2,ref_ant,freq,n_size,offsetunit,n_delay,az_step);
	azoff=0.; eloff=0.;
    	start_time = time(NULL);
	ts = *localtime(&start_time);
	sprintf(file_name, "/data/engineering/holo/ant%1d/%02d%02d_%02d%02d.holo%1d",ant_scn,ts.tm_mon+1,ts.tm_mday,ts.tm_hour,ts.tm_min,ant_scn);
	printf("freq=%4.1fGHz grid=%4.1f:%4.1f filename %s\n",freq,az_step,el_step,file_name);
/* Dual */
	sprintf(file_name2, "/data/engineering/holo/ant%1d/%02d%02d_%02d%02d.holo%1d",ant_scn2,ts.tm_mon+1,ts.tm_mday,ts.tm_hour,ts.tm_min,ant_scn2);
	printf("freq=%4.1fGHz grid=%4.1f:%4.1f filename %s\n",freq,az_step,el_step,file_name2);

	if(ant_scn==1){
		if(freq < 500.0){
			az_bore=AZ_BORE1;
			el_bore=EL_BORE1;
		}
		else{
			az_bore=AZ_BORE1H;
			el_bore=EL_BORE1H;
		}
	}

	if(ant_scn==2){
		if(freq < 500.0){
			az_bore=AZ_BORE2;
			el_bore=EL_BORE2;
		}
		else{
			az_bore=AZ_BORE2H;
			el_bore=EL_BORE2H;
		}
	}
	else if(ant_scn==3){
		if(freq < 500.0){
			az_bore=AZ_BORE3;
			el_bore=EL_BORE3;
		}
		else{
			az_bore=AZ_BORE3H;
			el_bore=EL_BORE3H;
		}
	}
	else if(ant_scn==4){
		if(freq < 500.0){
			az_bore=AZ_BORE4;
			el_bore=EL_BORE4;
		}
		else{
			az_bore=AZ_BORE4H;
			el_bore=EL_BORE4H;
		}
	}
	else if(ant_scn==5){
		if(freq < 500.0){
			az_bore=AZ_BORE5;
			el_bore=EL_BORE5;
		}
		else{
			az_bore=AZ_BORE5H;
			el_bore=EL_BORE5H;
		}
	}
	else if(ant_scn==6){
		if(freq < 500.0){
			az_bore=AZ_BORE6;
			el_bore=EL_BORE6;
		}
		else{
			az_bore=AZ_BORE6H;
			el_bore=EL_BORE6H;
		}
	}
	else if(ant_scn==7){
		if(freq < 500.0){
			az_bore=AZ_BORE7;
			el_bore=EL_BORE7;
		}
		else{
			az_bore=AZ_BORE7H;
			el_bore=EL_BORE7H;
		}
	}
	else if(ant_scn==8){
		if(freq < 500.0){
			az_bore=AZ_BORE8;
			el_bore=EL_BORE8;
		}
		else{
			az_bore=AZ_BORE8H;
			el_bore=EL_BORE8H;
		}
	}
	else if(ant_scn==0){ /* for debug */
		az_bore=AZ_BORE0;
		el_bore=EL_BORE0;
		ant_scn=4;
	}
/* Dual */
	if(ant_scn2==1){
		if(freq < 500.0){
			az_bore2=AZ_BORE1;
			el_bore2=EL_BORE1;
		}
		else{
			az_bore2=AZ_BORE1H;
			el_bore2=EL_BORE1H;
		}
	}

	if(ant_scn2==2){
		if(freq < 500.0){
			az_bore2=AZ_BORE2;
			el_bore2=EL_BORE2;
		}
		else{
			az_bore2=AZ_BORE2H;
			el_bore2=EL_BORE2H;
		}
	}
	else if(ant_scn2==3){
		if(freq < 500.0){
			az_bore2=AZ_BORE3;
			el_bore2=EL_BORE3;
		}
		else{
			az_bore2=AZ_BORE3H;
			el_bore2=EL_BORE3H;
		}
	}
	else if(ant_scn2==4){
		if(freq < 500.0){
			az_bore2=AZ_BORE4;
			el_bore2=EL_BORE4;
		}
		else{
			az_bore2=AZ_BORE4H;
			el_bore2=EL_BORE4H;
		}
	}
	else if(ant_scn2==5){
		if(freq < 500.0){
			az_bore2=AZ_BORE5;
			el_bore2=EL_BORE5;
		}
		else{
			az_bore2=AZ_BORE5H;
			el_bore2=EL_BORE5H;
		}
	}
	else if(ant_scn2==6){
		if(freq < 500.0){
			az_bore2=AZ_BORE6;
			el_bore2=EL_BORE6;
		}
		else{
			az_bore2=AZ_BORE6H;
			el_bore2=EL_BORE6H;
		}
	}
	else if(ant_scn2==7){
		if(freq < 500.0){
			az_bore2=AZ_BORE7;
			el_bore2=EL_BORE7;
		}
		else{
			az_bore2=AZ_BORE7H;
			el_bore2=EL_BORE7H;
		}
	}
	else if(ant_scn2==8){
		if(freq < 500.0){
			az_bore2=AZ_BORE8;
			el_bore2=EL_BORE8;
		}
		else{
			az_bore2=AZ_BORE8H;
			el_bore2=EL_BORE8H;
		}
	}
	else if(ant_scn2==0){ /* for debug */
		az_bore2=AZ_BORE0;
		el_bore2=EL_BORE0;
		ant_scn2=2;
	}
/* Dual End */

	if(ref_ant==1){
		if(freq < 500.0){
			az_ref=AZ_BORE1;
			el_ref=EL_BORE1;
		}
		else{
			az_ref=AZ_BORE1H;
			el_ref=EL_BORE1H;
		}
	}

	if(ref_ant==2){
		if(freq < 500.0){
			az_ref=AZ_BORE2;
			el_ref=EL_BORE2;
		}
		else{
			az_ref=AZ_BORE2H;
			el_ref=EL_BORE2H;
		}
	}
	else if(ref_ant==3){
		if(freq < 500.0){
			az_ref=AZ_BORE3;
			el_ref=EL_BORE3;
		}
		else{
			az_ref=AZ_BORE3H;
			el_ref=EL_BORE3H;
		}
	}
	else if(ref_ant==4){
		if(freq < 500.0){
			az_ref=AZ_BORE4;
			el_ref=EL_BORE4;
		}
		else{
			az_ref=AZ_BORE4H;
			el_ref=EL_BORE4H;
		}
	}
	else if(ref_ant==5){
		if(freq < 500.0){
			az_ref=AZ_BORE5;
			el_ref=EL_BORE5;
		}
		else{
			az_ref=AZ_BORE5H;
			el_ref=EL_BORE5H;
		}
	}
	else if(ref_ant==6){
		if(freq < 500.0){
			az_ref=AZ_BORE6;
			el_ref=EL_BORE6;
		}
		else{
			az_ref=AZ_BORE6H;
			el_ref=EL_BORE6H;
		}
	}
	else if(ref_ant==7){
		if(freq < 500.0){
			az_ref=AZ_BORE7;
			el_ref=EL_BORE7;
		}
		else{
			az_ref=AZ_BORE7H;
			el_ref=EL_BORE7H;
		}
	}
	else if(ref_ant==8){
		if(freq < 500.0){
			az_ref=AZ_BORE8;
			el_ref=EL_BORE8;
		}
		else{
			az_ref=AZ_BORE8H;
			el_ref=EL_BORE8H;
		}
	}
	else if(ref_ant==0){ /* for debug */
		az_ref=AZ_BORE0;
		el_ref=EL_BORE0;
		ref_ant=5;
	}

	fp_raw_data = fopen(file_name,"w");
	if(fp_raw_data==NULL){
		printf("cannot open the data file\n");
		exit(1);
	}
/* Dual */
	fp_raw_data2 = fopen(file_name2,"w");
	if(fp_raw_data2==NULL){
		printf("cannot open the data file 2\n");
		exit(1);
	}
/* End Dual */

/* initializing ref. mem. */
	rm_status=rm_open(antlist);
	if(rm_status != RM_SUCCESS) {
		rm_error_message(rm_status,"rm_open()");
		exit(1);
	}

/* configure for interrupt through RM */
 	rm_status=rm_monitor(ant_scn,"RM_HOLO_ENC_HNDSHK_S");
	if(rm_status != RM_SUCCESS) {
		rm_error_message(rm_status,"rm_monitor()");
		exit(1);
	}
/* Dual */
 	rm_status=rm_monitor(ant_scn2,"RM_HOLO_ENC_HNDSHK_S");
	if(rm_status != RM_SUCCESS) {
		rm_error_message(rm_status,"rm_monitor()");
		exit(1);
	}
/* */
	printf("checking the encoderServer status\n");

	handshake = 0;

	rm_status=rm_write_notify(ant_scn,"RM_HOLO_ENC_HNDSHK_S",&handshake);
		if(rm_status != RM_SUCCESS) {
		rm_error_message(rm_status,"rm_write_notify()");
		exit(1);
		} 

	sleep(1);

	rm_status=rm_read(ant_scn,"RM_HOLO_ENC_HNDSHK_S",&handshake);
		if(rm_status != RM_SUCCESS) {
		rm_error_message(rm_status,"rm_read()");
		exit(1);
		}  

/*	if (handshake==0) 
	{
	printf("encoder Server is not running on acc%1d!\n",ant_scn);
	exit(-1);
	}
*/

/* Dual */
/* Should check encoderServer for ant_scn2 here */

/* open Device tables */
	OFLAG=O_RDWR|O_NDELAY;
	adc0_fd = open("/dev/xVME564-0", OFLAG, 0); 
	if(adc0_fd<0){printf("cannot open device 0\n");exit(1);}
	adc1_fd = open("/dev/xVME564-1", OFLAG, 0); 
	if(adc1_fd<0){printf("cannot open device 1\n");exit(1);}
/* Dual */
	adc2_fd = open("/dev/xVME564-2", OFLAG, 0); 
	if(adc2_fd<0){printf("cannot open device 0\n");exit(1);}
	adc3_fd = open("/dev/xVME564-3", OFLAG, 0);
	if(adc3_fd<0){printf("cannot open device 1\n");exit(1);}
/* */
/*********************Get the antenna ready for*************************/
/* go to starting position  */
/* numbers for ant6 frozen in for now 30 Jul 2002 
	printf("frozen az_bore & el_bore for ant6\n");
	az_bore=48.25;
	el_bore=18.60; */
	printf ("Going to bore-sight position \n");
	sprintf(command, "azoff -a %d -s 0", ant_scn);
	system (command);
	sprintf(command, "eloff -a %d -s 0", ant_scn);
	system (command);
	sprintf(command, "azel -a %d -z %f -e %f", ref_ant, az_ref, el_ref);
	printf ("command: %s \n", command);
	system (command);
	sprintf(command, "azel -a %d -z %f -e %f", ant_scn, az_bore, el_bore);
	printf ("command: %s \n", command);
	system (command);
	sleep(1);
	sprintf(command, "offsetUnit -a %d -s %d", ant_scn, -offsetunit); 
	printf ("command: %s \n", command);
	system (command);
	sleep(1);
	sprintf(command, "resume -a %d", ant_scn); 
	system (command);

/* Dual */
	printf ("Going to bore-sight position \n");
	sprintf(command, "azoff -a %d -s 0", ant_scn2);
	system (command);
	sprintf(command, "eloff -a %d -s 0", ant_scn2);
	system (command);
	sprintf(command, "azel -a %d -z %f -e %f", ref_ant, az_ref, el_ref);
	printf ("command: %s \n", command);
	system (command);
	sprintf(command, "azel -a %d -z %f -e %f", ant_scn2, az_bore2, el_bore2);
	printf ("command: %s \n", command);
	system (command);
	sleep(1);
	sprintf(command, "offsetUnit -a %d -s %d", ant_scn2, -offsetunit); 
	printf ("command: %s \n", command);
	system (command);
	sleep(1);
	sprintf(command, "resume -a %d", ant_scn2); 
	system (command);
/* */
/*	printf ("Make sure encoderServer is running on acc%1d.\n",ant_scn); */
/* Dual */
	printf ("Make sure encoderServer is running on acc%1d and acc%1d\n",ant_scn,ant_scn2);
	printf ("Hit Return to continue:");
	fgetc(stdin);
	printf ("Starting raster on antennas %d and %d\n", ant_scn, ant_scn2);
	azoff=((float)n_size/2.+1.)*az_step;
	eloff=((float)n_size/2.)*el_step - n_start*el_step;
/*	sprintf(command, "azoff -a %d -s %f", ant_scn, azoff+offsetunit*3); */
/* Dual */
	sprintf(command, "azoff -a %d,%d -s %f", ant_scn, ant_scn2,azoff+offsetunit*3); 
/* */
	printf ("command: %s ", command);
	system (command);
	sleep(1);
/*	sprintf(command, "eloff -a %d -s %f", ant_scn, eloff);  */
/* Dual */
	sprintf(command, "eloff -a %d,%d -s %f", ant_scn, ant_scn2, eloff); 
/* */
	printf ("%s \n", command);
	system (command);
	sleep(5);

	
	/* The following is from Bob's ttsubs.c main */
        ttime.timeout_ticks = 2;        /* This margin should be safe */
        ttfd = open("/dev/vme_sg_simple", O_RDWR, 0);
        if(ttfd <= 0) {
            fprintf(stderr, "Error opening TrueTime - /dev/vme_sg_simple\n");
            exit(SYSERR_RTN);
        }


/*********************Start taking data*************************/
	printf ("Start raster! \n");
    	start_time = time(NULL);
	Start_time = start_time;

/* Elevation Loop */
	for (j=n_start; j<n_size; j++){  /*do one row*/ 
		sprintf(command, "azel -a %d -z %f -e %f", ref_ant, az_ref, el_ref);
		printf ("%s \n", command);
		system (command);
		i=0; 
	currentAzoff=0;
    		start_time = time(NULL);
/*		sprintf(command, "azscan -a %d", ant_scn);*/
/* Dual */
		sprintf(command, "azscan -a %d,%d", ant_scn, ant_scn2);
/* */
		printf ("%s \n", command);
		sleep(1);
		istatus=system (command);

		printf("%6d %6d  %6f\n",istatus,currentAzoff,azoff);
		printf("%6d  %6f \n",currentAzoff,azoff); fflush(stdout);
/*		while((1.0*currentAzoff)>-azoff){ */
/*		while((1.0*currentAzoff)>(-(azoff+offsetunit))){ */
/* Dual */
		while(((1.0*currentAzoff)>(-(azoff+offsetunit)))||((1.0*currentAzoff2)>(-(azoff+offsetunit)))){
/* */

/*			 for (i_delay=1; i_delay<n_delay; i_delay++){} */ /*  delay using GetPosTime() from ttsubs.c- using interrupts from TrueTime */
                        cur_time=time(NULL);
                        time1[i]=cur_time-Start_time;

			GetPosTime();
/*
            printf("Day %d, msec %d\n", day, msec);
*/
			truetime2ms[i]=msec;
                        cur_time=time(NULL);
                        time2[i]=cur_time-Start_time;
			adc0_rd=read(adc0_fd, amp_vvm[1].byte, count);
			if(adc0_rd<=0){
				printf("cannot read device 0\n");
				exit(1);
				close(adc0_fd);
				close(adc1_fd);
			} 
/* Read total power values 10 May 2010 */
			rm_status=rm_read(ant_scn,"RM_CONT_DET_MUWATT_V2_F",&contDet_a1);
			if(rm_status != RM_SUCCESS) {
			rm_error_message(rm_status,"rm_read()");
			exit(1);
       			}/*if end*/
			rm_status=rm_read(ant_scn2,"RM_CONT_DET_MUWATT_V2_F",&contDet_a2);
			if(rm_status != RM_SUCCESS) {
			rm_error_message(rm_status,"rm_read()");
			exit(1);
       			}/*if end*/
			rm_status=rm_read(ref_ant,"RM_CONT_DET_MUWATT_V2_F",&contDet_ar);
			if(rm_status != RM_SUCCESS) {
			rm_error_message(rm_status,"rm_read()");
			exit(1);
       			}/*if end*/
 /* usleep(1); */
			adc1_rd=read(adc1_fd, pha_vvm[1].byte, count);
			if(adc1_rd<=0){
				printf("cannot read device 1\n");
				exit(1);
				close(adc0_fd);
				close(adc1_fd);
			}
/* Dual */
			adc2_rd=read(adc2_fd, amp_vvm2[1].byte, count);
			if(adc2_rd<=0){
				printf("cannot read device 2\n");
				exit(1);
				close(adc2_fd);
				close(adc3_fd);
			} 
 /* usleep(1); */
			adc3_rd=read(adc3_fd, pha_vvm2[1].byte, count);
			if(adc3_rd<=0){
				printf("cannot read device 3\n");
				exit(1);
				close(adc2_fd);
				close(adc3_fd);
			}
/* End Dual */
                        cur_time=time(NULL);
                        time3[i]=cur_time-Start_time;
/* read encoders using encoderClient */
			encoderClient(ant_scn,&az_enc,&el_enc,&m1009,&m1010); 
/*		        printf("First read done\n %f %f %f\n", az_enc, el_enc, currentAzoffD); */
/* Dual */
			encoderClient(ant_scn2,&az_enc2,&el_enc2,&m1009,&m1010); 
/*		        printf("Second read done\n %f %f %f\n", az_enc2, el_enc2, currentAzoffD); */
/* */
                        cur_time=time(NULL);
                        time4[i]=cur_time-Start_time;

			az_read[i]=az_enc;el_read[i]=el_enc;
			amp_read[i]=amp_vvm[1].word;
			phase_read[i]=pha_vvm[1].word;
/* Dual */
			az_read2[i]=az_enc2;el_read2[i]=el_enc2;
			amp_read2[i]=amp_vvm2[1].word;
			phase_read2[i]=pha_vvm2[1].word;
/* End Dual */
/* adding time stamping of samples: 01 Nov 04 */

			rm_status=rm_read(ant_scn,"RM_AZOFF_D",&currentAzoffD);
			if(rm_status != RM_SUCCESS) {
			rm_error_message(rm_status,"rm_read()");
			exit(1);
       			}/*if end*/
			currentAzoff=(short)currentAzoffD;
/* Dual */
			rm_status=rm_read(ant_scn2,"RM_AZOFF_D",&currentAzoffD);
			if(rm_status != RM_SUCCESS) {
			rm_error_message(rm_status,"rm_read()");
			exit(1);
       			}/*if end*/
			currentAzoff2=(short)currentAzoffD;
/* End Dual */
/* printf("line# %3d %3d %4d %6.2f %6.2f %6.2f %5d %d\n",j,i,currentAzoff,azoff,az_enc,el_enc,amp_read[i],icount); fflush(stdout); */

			i++;
		}/*while end*/

/*		sprintf(command, "stopScan -a %d", ant_scn); */
/* Dual */
		sprintf(command, "stopScan -a %d,%d", ant_scn, ant_scn2);
/* */
	      	system (command);
		sleep(1); 
		printf ("command: %s \n", command); 
	/*	system (command);
		usleep(10); */
		eloff = eloff - el_step;
/*		sprintf(command, "eloff -a %d -s %f", ant_scn, eloff); */
/* Dual */
		sprintf(command, "eloff -a %d,%d -s %f", ant_scn, ant_scn2, eloff); 
/* */
		system (command);
		sleep(1);
		printf ("command: %s ", command);fflush(stdout);
	/*	system (command);
		usleep(10); */
/*		sprintf(command, "azoff -a %d -s %f", ant_scn, azoff+offsetunit*3); */
/* Dual */
		sprintf(command, "azoff -a %d,%d -s %f", ant_scn, ant_scn2, azoff+offsetunit*3); 
/* */
		system (command);
		sleep(1);
		printf ("command: %s ", command);fflush(stdout);
	/*	system (command); */
/* write out also the total powers 10 May 2010 */
		for(k=1;k<i;k++){
			fprintf(fp_raw_data,"%10.6f %10.6f %6d %6d %6.4f %3d %d %d %d %d %d\n" ,\
			el_read[k],az_read[k],amp_read[k],phase_read[k],contDet_a1[0],j,time1[k],time2[k],time3[k],time4[k],truetime2ms[k]);
			fflush(fp_raw_data);
    		}/*for end*/
/* Dual */
		for(k=1;k<i;k++){
			fprintf(fp_raw_data2,"%10.6f %10.6f %6d %6d %6.4f %3d %d %d %d %d %d\n" ,\
			el_read2[k],az_read2[k],amp_read2[k],phase_read2[k],contDet_a2[0],j,time1[k],time2[k],time3[k],time4[k],truetime2ms[k]);
			fflush(fp_raw_data2);
    		}/*for end*/
/* End Dual */
		/* Place for including online dispaly command */
		printf ("Flush the data\n"); 
    		cur_time=time(NULL);
    		printf ("%3dth line: Stop time: %d time taken: %d secs samples: %d\n", \
		j, cur_time, cur_time-start_time, i);
		sleep(4);
/*		sprintf(command, "azoff -a %d -s %f", ant_scn, azoff);
		printf ("command: %s ", command);fflush(stdout);
		system (command);
		usleep(10);
		sprintf(command, "eloff -a %d -s %f", ant_scn, eloff);
		printf ("command: %s ", command);fflush(stdout);
		system (command);
		usleep(10);*/
    	} 
/* for end; El loop done */

	printf ("Start: %d Stop: %d Time taken: %d secs\n", Start_time, cur_time, cur_time - Start_time);
/*	sprintf(command, "azoff -a %d -s 0", ant_scn); */
/* Dual */
	sprintf(command, "azoff -a %d,%d -s 0", ant_scn,ant_scn2); 
/* */
	system (command);
	sleep(10);
	printf ("command: %s ", command);
/*	sprintf(command, "eloff -a %d -s 0", ant_scn); */
/* Dual */
	sprintf(command, "eloff -a %d,%d -s 0", ant_scn,ant_scn2);
/* */
	system (command);
	sleep(10);
	printf ("command: %s\n", command);

/*********************clear up everything*************************/
	fclose(fp_raw_data);
	close(adc0_fd);
	close(adc1_fd);
/* Dual */
	fclose(fp_raw_data2);
	close(adc2_fd);
	close(adc3_fd);
/* End Dual */
/* release the RM interrupt waiting */
	rm_status=rm_clear_monitor();
	if(rm_status != RM_SUCCESS) {
	rm_error_message(rm_status,"rm_clear()");
	}
	printf ("DONE!\n");
	exit(0);
}				/* end main */
コード例 #2
0
int main(int argc, char *argv[])  {

	char c,command_n[30];

	short pmac_command_flag=0;
	int dsm_status;
	smapoptContext optCon;
	int i ;

        int trackStatus=0;
        int tracktimestamp,timestamp;
	time_t dsmtimestamp;

	 struct  smapoptOption optionsTable[] = {
                {"help",'h',SMAPOPT_ARG_NONE,0,'h'},
                {NULL,0,0,NULL,0}
        };


        optCon = smapoptGetContext("stop", argc, argv, optionsTable,0);
 
        while ((c = smapoptGetNextOpt(optCon)) >= 0) {
 
            switch(c) {
                    case 'h':
                    usage(0,NULL,NULL);
                    break;
            }
        }


	for(i=0;i<30;i++) {
        command_n[i]=0x0;
        };
	command_n[0]='@'; /* send the command */
        pmac_command_flag=0;

 
        if(c<-1) {
        fprintf(stderr, "%s: %s\n",
                smapoptBadOption(optCon, SMAPOPT_BADOPTION_NOALIAS),
                smapoptStrerror(c));
        }
         smapoptFreeContext(optCon);
 
              /* initializing DSM */
        dsm_status=dsm_open();
        if(dsm_status != DSM_SUCCESS) {
                dsm_error_message(dsm_status,"dsm_open()");
                exit(1);
        }


#if 0
          /* check if track is running on this antenna */

        dsm_status=dsm_read(DSM_HOST,"DSM_UNIX_TIME_L",&timestamp,&dsmtimestamp);
	dsm_status=dsm_read(DSM_HOST,"DSM_TRACK_TIMESTAMP_L",&tracktimestamp,&dsmtimestamp);
        if(abs(tracktimestamp-timestamp)>3L) {
        trackStatus=0;
        printf("Track is not running.\n");
        }
        if(abs(tracktimestamp-timestamp)<=3L) trackStatus=1;

        if(trackStatus==1) {
        dsm_status=dsm_write(DSM_HOST,"DSM_COMMANDED_TRACK_COMMAND_C30",
                                        &command_n);
        if(dsm_status != DSM_SUCCESS) {
                dsm_error_message(dsm_status,"dsm_write()");
                exit(1);
        }
#endif
 
        dsm_status=dsm_write(DSM_HOST,"DSM_COMMANDED_TRACK_COMMAND_C30",
                                        &command_n);
        if(dsm_status != DSM_SUCCESS) {
                dsm_error_message(dsm_status,"dsm_write()");
                exit(1);
        }
 
/* the following dsm_write should actually be dsm_write_notify- but due 
to a bug in gltTrack, dsm_write_notify does not work */
     dsm_status=dsm_write(DSM_HOST,"DSM_COMMAND_FLAG_S",
                                        &pmac_command_flag);
 
     if(dsm_status != DSM_SUCCESS) {
                dsm_error_message(dsm_status,"dsm_write()");
                exit(1);
        }
 
#if 0
 
	} /* if track is running */
#endif

	return 0;
 
}