version 1.8, 2012/08/21 23:29:34
|
version 1.9, 2012/09/12 18:26:58
|
|
|
} | } |
} | } |
| |
|
if (verbflag) |
printf("gapfile read, gapsize = %d\n",gapsize); | printf("gapfile read, gapsize = %d\n",gapsize); |
| |
data_dim=gapsize; | data_dim=gapsize; |
|
|
} | } |
} | } |
| |
|
if (verbflag) |
printf("num_sections = %d\n",num_sections); | printf("num_sections = %d\n",num_sections); |
| |
DRMS_RecordSet_t *inrecset = drms_open_records(drms_env, inrecquery, &status); | DRMS_RecordSet_t *inrecset = drms_open_records(drms_env, inrecquery, &status); |
|
|
tsstart[1]=m; | tsstart[1]=m; |
tsend[1]=m; | tsend[1]=m; |
outarr = drms_array_create(usetype, 2, tslength, tmpptr, &status); | outarr = drms_array_create(usetype, 2, tslength, tmpptr, &status); |
|
outarr->bzero=segoutarr[TSOUT]->bzero; |
|
outarr->bscale=segoutarr[TSOUT]->bscale; |
status=drms_segment_writeslice(segoutarr[TSOUT], outarr, tsstart, tsend, 0); | status=drms_segment_writeslice(segoutarr[TSOUT], outarr, tsstart, tsend, 0); |
free(outarr); | free(outarr); |
if (status != DRMS_SUCCESS) | if (status != DRMS_SUCCESS) |
|
|
tsstart[1]=m; | tsstart[1]=m; |
tsend[1]=m; | tsend[1]=m; |
outarr = drms_array_create(usetype, 2, tslength, tmpptr, &status); | outarr = drms_array_create(usetype, 2, tslength, tmpptr, &status); |
|
outarr->bzero=segoutarr[FFTOUT]->bzero; |
|
outarr->bscale=segoutarr[FFTOUT]->bscale; |
status=drms_segment_writeslice(segoutarr[FFTOUT], outarr, tsstart, tsend, 0); | status=drms_segment_writeslice(segoutarr[FFTOUT], outarr, tsstart, tsend, 0); |
free(outarr); | free(outarr); |
if (status != DRMS_SUCCESS) | if (status != DRMS_SUCCESS) |
|
|
fft1start[1]=lmode+m*msignarr[FFT1OUT]; | fft1start[1]=lmode+m*msignarr[FFT1OUT]; |
fft1end[1]=lmode+m*msignarr[FFT1OUT]; | fft1end[1]=lmode+m*msignarr[FFT1OUT]; |
outarr = drms_array_create(usetype, 2, fft1length, data+2*imin, &status); | outarr = drms_array_create(usetype, 2, fft1length, data+2*imin, &status); |
|
outarr->bzero=segoutarr[FFT1OUT]->bzero; |
|
outarr->bscale=segoutarr[FFT1OUT]->bscale; |
status=drms_segment_writeslice(segoutarr[FFT1OUT], outarr, fft1start, fft1end, 0); | status=drms_segment_writeslice(segoutarr[FFT1OUT], outarr, fft1start, fft1end, 0); |
free(outarr); | free(outarr); |
if (status != DRMS_SUCCESS) | if (status != DRMS_SUCCESS) |
|
|
fft1start[1]=lmode-m*msignarr[FFT1OUT]; | fft1start[1]=lmode-m*msignarr[FFT1OUT]; |
fft1end[1]=lmode-m*msignarr[FFT1OUT]; | fft1end[1]=lmode-m*msignarr[FFT1OUT]; |
outarr = drms_array_create(usetype, 2, fft1length, datahold, &status); | outarr = drms_array_create(usetype, 2, fft1length, datahold, &status); |
|
outarr->bzero=segoutarr[FFT1OUT]->bzero; |
|
outarr->bscale=segoutarr[FFT1OUT]->bscale; |
status=drms_segment_writeslice(segoutarr[FFT1OUT], outarr, fft1start, fft1end, 0); | status=drms_segment_writeslice(segoutarr[FFT1OUT], outarr, fft1start, fft1end, 0); |
free(outarr); | free(outarr); |
if (status != DRMS_SUCCESS) | if (status != DRMS_SUCCESS) |
|
|
powstart[1]=lmode+m*msignarr[POWOUT]; | powstart[1]=lmode+m*msignarr[POWOUT]; |
powend[1]=lmode+m*msignarr[POWOUT]; | powend[1]=lmode+m*msignarr[POWOUT]; |
outarr = drms_array_create(usetype, 2, powlength, pow+imin, &status); | outarr = drms_array_create(usetype, 2, powlength, pow+imin, &status); |
|
outarr->bzero=segoutarr[POWOUT]->bzero; |
|
outarr->bscale=segoutarr[POWOUT]->bscale; |
status=drms_segment_writeslice(segoutarr[POWOUT], outarr, powstart, powend, 0); | status=drms_segment_writeslice(segoutarr[POWOUT], outarr, powstart, powend, 0); |
free(outarr); | free(outarr); |
if (status != DRMS_SUCCESS) | if (status != DRMS_SUCCESS) |
|
|
powstart[1]=lmode-m*msignarr[POWOUT]; | powstart[1]=lmode-m*msignarr[POWOUT]; |
powend[1]=lmode-m*msignarr[POWOUT]; | powend[1]=lmode-m*msignarr[POWOUT]; |
outarr = drms_array_create(usetype, 2, powlength, datahold, &status); | outarr = drms_array_create(usetype, 2, powlength, datahold, &status); |
|
outarr->bzero=segoutarr[POWOUT]->bzero; |
|
outarr->bscale=segoutarr[POWOUT]->bscale; |
status=drms_segment_writeslice(segoutarr[POWOUT], outarr, powstart, powend, 0); | status=drms_segment_writeslice(segoutarr[POWOUT], outarr, powstart, powend, 0); |
free(outarr); | free(outarr); |
if (status != DRMS_SUCCESS) | if (status != DRMS_SUCCESS) |
|
|
powstart[1]=lmode+m*msignarr[ARPOWOUT]; | powstart[1]=lmode+m*msignarr[ARPOWOUT]; |
powend[1]=lmode+m*msignarr[ARPOWOUT]; | powend[1]=lmode+m*msignarr[ARPOWOUT]; |
outarr = drms_array_create(usetype, 2, powlength, pow+imin, &status); | outarr = drms_array_create(usetype, 2, powlength, pow+imin, &status); |
|
outarr->bzero=segoutarr[ARPOWOUT]->bzero; |
|
outarr->bscale=segoutarr[ARPOWOUT]->bscale; |
status=drms_segment_writeslice(segoutarr[ARPOWOUT], outarr, powstart, powend, 0); | status=drms_segment_writeslice(segoutarr[ARPOWOUT], outarr, powstart, powend, 0); |
free(outarr); | free(outarr); |
if (status != DRMS_SUCCESS) | if (status != DRMS_SUCCESS) |
|
|
powstart[1]=lmode-m*msignarr[ARPOWOUT]; | powstart[1]=lmode-m*msignarr[ARPOWOUT]; |
powend[1]=lmode-m*msignarr[ARPOWOUT]; | powend[1]=lmode-m*msignarr[ARPOWOUT]; |
outarr = drms_array_create(usetype, 2, powlength, datahold, &status); | outarr = drms_array_create(usetype, 2, powlength, datahold, &status); |
|
outarr->bzero=segoutarr[ARPOWOUT]->bzero; |
|
outarr->bscale=segoutarr[ARPOWOUT]->bscale; |
status=drms_segment_writeslice(segoutarr[ARPOWOUT], outarr, powstart, powend, 0); | status=drms_segment_writeslice(segoutarr[ARPOWOUT], outarr, powstart, powend, 0); |
free(outarr); | free(outarr); |
if (status != DRMS_SUCCESS) | if (status != DRMS_SUCCESS) |
|
|
for (i=0; i<npow; i++) | for (i=0; i<npow; i++) |
datahold[i] = msum[imin+i] * c; | datahold[i] = msum[imin+i] * c; |
outarr = drms_array_create(usetype, 2, powlength, datahold, &status); | outarr = drms_array_create(usetype, 2, powlength, datahold, &status); |
|
outarr->bzero=segoutarr[MAVGOUT]->bzero; |
|
outarr->bscale=segoutarr[MAVGOUT]->bscale; |
status=drms_segment_write(segoutarr[MAVGOUT], outarr, 0); | status=drms_segment_write(segoutarr[MAVGOUT], outarr, 0); |
free(outarr); | free(outarr); |
if (status != DRMS_SUCCESS) | if (status != DRMS_SUCCESS) |