Changeset 2940 for trunk/external-alma/atnf/PKSIO/GBTFITSreader.cc
- Timestamp:
- 05/13/14 14:39:48 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/external-alma/atnf/PKSIO/GBTFITSreader.cc
r2907 r2940 170 170 // Move to the SDFITS extension. 171 171 cALFA = cALFA_BD = cALFA_CIMA = 0; 172 char *extName = "SINGLE DISH" ;173 if (fits_movnam_hdu(cSDptr, BINARY_TBL, "SINGLE DISH", 0, &cStatus)) {172 char *extName = (char *)"SINGLE DISH" ; 173 if (fits_movnam_hdu(cSDptr, BINARY_TBL, extName, 0, &cStatus)) { 174 174 // No SDFITS table, look for BDFITS or CIMAFITS. 175 175 cStatus = 0; 176 if (fits_movnam_hdu(cSDptr, BINARY_TBL, "BDFITS", 0, &cStatus) == 0) {176 if (fits_movnam_hdu(cSDptr, BINARY_TBL, (char *)"BDFITS", 0, &cStatus) == 0) { 177 177 cALFA_BD = 1; 178 extName = "BDFITS" ;178 extName = (char *)"BDFITS" ; 179 179 180 180 } else { 181 181 cStatus = 0; 182 if (fits_movnam_hdu(cSDptr, BINARY_TBL, "CIMAFITS", 0, &cStatus) == 0) {182 if (fits_movnam_hdu(cSDptr, BINARY_TBL, (char *)"CIMAFITS", 0, &cStatus) == 0) { 183 183 cALFA_CIMA = 1; 184 extName = "CIMAFITS" ;184 extName = (char *)"CIMAFITS" ; 185 185 186 186 // Check for later versions of CIMAFITS. 187 187 float version; 188 readParm( "VERSION", TFLOAT, &version);188 readParm((char *)(char *)"VERSION", TFLOAT, &version); 189 189 if (version >= 2.0f) cALFA_CIMA = int(version); 190 190 … … 282 282 // GBT data. 283 283 char telescope[32]; 284 readParm( "TELESCOP", TSTRING, telescope); // Core.284 readParm((char *)(char *)"TELESCOP", TSTRING, telescope); // Core. 285 285 cGBT = strncmp(telescope, "GBT", 3) == 0 || 286 286 strncmp(telescope, "NRAO_GBT", 8) == 0; … … 288 288 289 289 // Check that the DATA array column is present. 290 findData(DATA, "DATA", TFLOAT);290 findData(DATA, (char *)(char *)"DATA", TFLOAT); 291 291 haveSpectra = cHaveSpectra = cData[DATA].colnum > 0; 292 292 … … 305 305 // no TDIM (or MAXISn) card; use the LAGS_IN value. 306 306 cNAxes = 5; 307 readParm( "LAGS_IN", TLONG, cNAxis);307 readParm((char *)"LAGS_IN", TLONG, cNAxis); 308 308 cNAxis[1] = 1; 309 309 cNAxis[2] = 1; … … 325 325 } 326 326 327 findData(FLAGGED, "FLAGGED", TBYTE);327 findData(FLAGGED, (char *)(char *)"FLAGGED", TBYTE); 328 328 329 329 } else { 330 330 // DATA column not present, check for a DATAXED keyword. 331 findData(DATAXED, "DATAXED", TSTRING);331 findData(DATAXED, (char *)(char *)"DATAXED", TSTRING); 332 332 if (cData[DATAXED].colnum < 0) { 333 333 log(LogOrigin( className, methodName, WHERE ), LogIO::SEVERE, "DATA array column absent from binary table."); … … 338 338 // Determine the number of axes and their length. 339 339 char dataxed[32]; 340 readParm( "DATAXED", TSTRING, dataxed);340 readParm((char *)(char *)"DATAXED", TSTRING, dataxed); 341 341 342 342 for (int iaxis = 0; iaxis < 5; iaxis++) cNAxis[iaxis] = 0; … … 348 348 } 349 349 350 char *CTYPE[5] = { "CTYPE1", "CTYPE2", "CTYPE3", "CTYPE4","CTYPE5"};351 char *CRPIX[5] = { "CRPIX1", "CRPIX2", "CRPIX3", "CRPIX4","CRPIX5"};352 char *CRVAL[5] = { "CRVAL1", "CRVAL2", "CRVAL3", "CRVAL4","CRVAL5"};353 char *CDELT[5] = { "CDELT1", "CDELT2", "CDELT3", "CDELT4","CDELT5"};350 char *CTYPE[5] = {(char *)"CTYPE1", (char *)"CTYPE2", (char *)"CTYPE3", (char *)"CTYPE4", (char *)"CTYPE5"}; 351 char *CRPIX[5] = {(char *)"CRPIX1", (char *)"CRPIX2", (char *)"CRPIX3", (char *)"CRPIX4", (char *)"CRPIX5"}; 352 char *CRVAL[5] = {(char *)"CRVAL1", (char *)"CRVAL2", (char *)"CRVAL3", (char *)"CRVAL4", (char *)"CRVAL5"}; 353 char *CDELT[5] = {(char *)"CDELT1", (char *)"CDELT2", (char *)"CDELT3", (char *)"CDELT4", (char *)"CDELT5"}; 354 354 355 355 // Find required DATA array axes. … … 429 429 // Fixed in ALFA CIMAFITS. 430 430 cRaAxis = 2; 431 raCRVAL = "CRVAL2A";431 raCRVAL = (char *)"CRVAL2A"; 432 432 433 433 cDecAxis = 3; 434 decCRVAL = "CRVAL3A";434 decCRVAL = (char *)"CRVAL3A"; 435 435 } 436 436 … … 444 444 445 445 // Set up machinery for data retrieval. 446 findData(SCAN, "SCAN", TINT); // Shared.447 findData(CYCLE, "CYCLE", TINT); // Additional.448 findData(DATE_OBS, "DATE-OBS", TSTRING); // Core.446 findData(SCAN, (char *)(char *)"SCAN", TINT); // Shared. 447 findData(CYCLE, (char *)(char *)"CYCLE", TINT); // Additional. 448 findData(DATE_OBS, (char *)"DATE-OBS", TSTRING); // Core. 449 449 450 450 if (cTimeAxis >= 0) { … … 461 461 462 462 } else { 463 findData(TIME, "TIME", TDOUBLE); // Core.464 } 465 466 findData(EXPOSURE, "EXPOSURE", TFLOAT); // Core.467 findData(OBJECT, "OBJECT", TSTRING); // Core.468 findData(OBJ_RA, "OBJ-RA", TDOUBLE); // Additional.469 findData(OBJ_DEC, "OBJ-DEC", TDOUBLE); // Additional.470 findData(RESTFRQ, "RESTFRQ", TDOUBLE); // Additional.471 findData(OBSMODE, "OBSMODE", TSTRING); // Shared.472 473 findData(BEAM, "BEAM", TSHORT); // Additional.474 findData(IF, "IF", TSHORT); // Additional.463 findData(TIME, (char *)(char *)"TIME", TDOUBLE); // Core. 464 } 465 466 findData(EXPOSURE, (char *)(char *)"EXPOSURE", TFLOAT); // Core. 467 findData(OBJECT, (char *)(char *)"OBJECT", TSTRING); // Core. 468 findData(OBJ_RA, (char *)"OBJ-RA", TDOUBLE); // Additional. 469 findData(OBJ_DEC, (char *)"OBJ-DEC", TDOUBLE); // Additional. 470 findData(RESTFRQ, (char *)(char *)"RESTFRQ", TDOUBLE); // Additional. 471 findData(OBSMODE, (char *)(char *)"OBSMODE", TSTRING); // Shared. 472 473 findData(BEAM, (char *)(char *)"BEAM", TSHORT); // Additional. 474 findData(IF, (char *)(char *)"IF", TSHORT); // Additional. 475 475 findData(FqRefVal, fqCRVAL, TDOUBLE); // Frequency reference value. 476 476 findData(FqDelt, fqCDELT, TDOUBLE); // Frequency increment. … … 478 478 findData(RA, raCRVAL, TDOUBLE); // Right ascension. 479 479 findData(DEC, decCRVAL, TDOUBLE); // Declination. 480 findData(SCANRATE, "SCANRATE", TFLOAT); // Additional.481 482 findData(TSYS, "TSYS", TFLOAT); // Core.483 findData(CALFCTR, "CALFCTR", TFLOAT); // Additional.484 findData(XCALFCTR, "XCALFCTR", TFLOAT); // Additional.485 findData(BASELIN, "BASELIN", TFLOAT); // Additional.486 findData(BASESUB, "BASESUB", TFLOAT); // Additional.487 findData(XPOLDATA, "XPOLDATA", TFLOAT); // Additional.488 489 findData(REFBEAM, "REFBEAM", TSHORT); // Additional.490 findData(TCAL, "TCAL", TFLOAT); // Shared.491 findData(TCALTIME, "TCALTIME", TSTRING); // Additional.492 findData(AZIMUTH, "AZIMUTH", TFLOAT); // Shared.493 findData(ELEVATIO, "ELEVATIO", TFLOAT); // Shared.494 findData(PARANGLE, "PARANGLE", TFLOAT); // Additional.495 findData(FOCUSAXI, "FOCUSAXI", TFLOAT); // Additional.496 findData(FOCUSTAN, "FOCUSTAN", TFLOAT); // Additional.497 findData(FOCUSROT, "FOCUSROT", TFLOAT); // Additional.498 findData(TAMBIENT, "TAMBIENT", TFLOAT); // Shared.499 findData(PRESSURE, "PRESSURE", TFLOAT); // Shared.500 findData(HUMIDITY, "HUMIDITY", TFLOAT); // Shared.501 findData(WINDSPEE, "WINDSPEE", TFLOAT); // Shared.502 findData(WINDDIRE, "WINDDIRE", TFLOAT); // Shared.480 findData(SCANRATE, (char *)(char *)"SCANRATE", TFLOAT); // Additional. 481 482 findData(TSYS, (char *)(char *)"TSYS", TFLOAT); // Core. 483 findData(CALFCTR, (char *)(char *)"CALFCTR", TFLOAT); // Additional. 484 findData(XCALFCTR, (char *)(char *)"XCALFCTR", TFLOAT); // Additional. 485 findData(BASELIN, (char *)(char *)"BASELIN", TFLOAT); // Additional. 486 findData(BASESUB, (char *)(char *)"BASESUB", TFLOAT); // Additional. 487 findData(XPOLDATA, (char *)(char *)"XPOLDATA", TFLOAT); // Additional. 488 489 findData(REFBEAM, (char *)(char *)"REFBEAM", TSHORT); // Additional. 490 findData(TCAL, (char *)(char *)"TCAL", TFLOAT); // Shared. 491 findData(TCALTIME, (char *)(char *)"TCALTIME", TSTRING); // Additional. 492 findData(AZIMUTH, (char *)(char *)"AZIMUTH", TFLOAT); // Shared. 493 findData(ELEVATIO, (char *)(char *)"ELEVATIO", TFLOAT); // Shared. 494 findData(PARANGLE, (char *)(char *)"PARANGLE", TFLOAT); // Additional. 495 findData(FOCUSAXI, (char *)(char *)"FOCUSAXI", TFLOAT); // Additional. 496 findData(FOCUSTAN, (char *)(char *)"FOCUSTAN", TFLOAT); // Additional. 497 findData(FOCUSROT, (char *)(char *)"FOCUSROT", TFLOAT); // Additional. 498 findData(TAMBIENT, (char *)(char *)"TAMBIENT", TFLOAT); // Shared. 499 findData(PRESSURE, (char *)(char *)"PRESSURE", TFLOAT); // Shared. 500 findData(HUMIDITY, (char *)(char *)"HUMIDITY", TFLOAT); // Shared. 501 findData(WINDSPEE, (char *)(char *)"WINDSPEE", TFLOAT); // Shared. 502 findData(WINDDIRE, (char *)(char *)"WINDDIRE", TFLOAT); // Shared. 503 503 504 504 findData(STOKES, polCRVAL, TINT); 505 findData(SIG, "SIG", TSTRING);506 findData(CAL, "CAL", TSTRING);507 508 findData(RVSYS, "RVSYS", TDOUBLE);509 findData(VFRAME, "VFRAME", TDOUBLE);510 findData(VELDEF, "VELDEF", TSTRING);511 512 findData(TIMESTAMP, "TIMESTAMP", TSTRING);513 findData(DURATION, "DURATION", TDOUBLE);514 findData(SAMPLER, "SAMPLER", TSTRING);505 findData(SIG, (char *)(char *)"SIG", TSTRING); 506 findData(CAL, (char *)(char *)"CAL", TSTRING); 507 508 findData(RVSYS, (char *)(char *)"RVSYS", TDOUBLE); 509 findData(VFRAME, (char *)(char *)"VFRAME", TDOUBLE); 510 findData(VELDEF, (char *)(char *)"VELDEF", TSTRING); 511 512 findData(TIMESTAMP, (char *)(char *)"TIMESTAMP", TSTRING); 513 findData(DURATION, (char *)(char *)"DURATION", TDOUBLE); 514 findData(SAMPLER, (char *)(char *)"SAMPLER", TSTRING); 515 515 516 516 if (cStatus) { … … 527 527 cScanNo = 0; 528 528 if (cALFA_CIMA) { 529 findData(SCAN, "SCAN_ID", TINT);529 findData(SCAN, (char *)"SCAN_ID", TINT); 530 530 if (cALFA_CIMA > 1) { 531 531 // Note that RECNUM increases by cNAxisTime per row. 532 findData(CYCLE, "RECNUM", TINT);532 findData(CYCLE, (char *)(char *)"RECNUM", TINT); 533 533 } else { 534 findData(CYCLE, "SUBSCAN", TINT);534 findData(CYCLE, (char *)(char *)"SUBSCAN", TINT); 535 535 } 536 536 } else if (cALFA_BD) { 537 findData(SCAN, "SCAN_NUMBER", TINT);538 findData(CYCLE, "PATTERN_NUMBER", TINT);537 findData(SCAN, (char *)"SCAN_NUMBER", TINT); 538 findData(CYCLE, (char *)"PATTERN_NUMBER", TINT); 539 539 } 540 540 } else { … … 554 554 if (cALFA) { 555 555 // ALFA INPUT_ID, 0-relative (overrides BEAM column if present). 556 findData(BEAM, "INPUT_ID", TSHORT);556 findData(BEAM, (char *)"INPUT_ID", TSHORT); 557 557 cBeam_1rel = 0; 558 558 … … 563 563 } else { 564 564 // ms2sdfits output, 0-relative "feed" number. 565 findData(BEAM, "MAIN_FEED1", TSHORT);565 findData(BEAM, (char *)"MAIN_FEED1", TSHORT); 566 566 cBeam_1rel = 0; 567 567 } … … 573 573 // ALFA data, 0-relative. 574 574 if (cALFA_CIMA > 1) { 575 findData(IF, "IFN", TSHORT);575 findData(IF, (char *)(char *)"IFN", TSHORT); 576 576 } else { 577 findData(IF, "IFVAL", TSHORT);577 findData(IF, (char *)(char *)"IFVAL", TSHORT); 578 578 } 579 579 cIF_1rel = 0; … … 582 582 // ms2sdfits writes a scalar "TSYS" column that averages the polarizations. 583 583 int colnum; 584 findCol( "SYSCAL_TSYS", &colnum);584 findCol((char *)"SYSCAL_TSYS", &colnum); 585 585 if (colnum > 0) { 586 586 // This contains the vector Tsys. 587 findData(TSYS, "SYSCAL_TSYS", TFLOAT);587 findData(TSYS, (char *)"SYSCAL_TSYS", TFLOAT); 588 588 } 589 589 … … 591 591 592 592 if (cData[SCANRATE].colnum < 0) { 593 findData(SCANRATE, "FIELD_POINTING_DIR_RATE", TFLOAT);593 findData(SCANRATE, (char *)"FIELD_POINTING_DIR_RATE", TFLOAT); 594 594 } 595 595 596 596 if (cData[RESTFRQ].colnum < 0) { 597 findData(RESTFRQ, "RESTFREQ", TDOUBLE);597 findData(RESTFRQ, (char *)(char *)"RESTFREQ", TDOUBLE); 598 598 if (cData[RESTFRQ].colnum < 0) { 599 findData(RESTFRQ, "SPECTRAL_WINDOW_REST_FREQUENCY", TDOUBLE);599 findData(RESTFRQ, (char *)"SPECTRAL_WINDOW_REST_FREQUENCY", TDOUBLE); 600 600 } 601 601 } 602 602 603 603 if (cData[OBJ_RA].colnum < 0) { 604 findData(OBJ_RA, "SOURCE_DIRECTION", TDOUBLE);604 findData(OBJ_RA, (char *)"SOURCE_DIRECTION", TDOUBLE); 605 605 } 606 606 if (cData[OBJ_DEC].colnum < 0) { 607 findData(OBJ_DEC, "SOURCE_DIRECTION", TDOUBLE);607 findData(OBJ_DEC, (char *)"SOURCE_DIRECTION", TDOUBLE); 608 608 } 609 609 … … 611 611 612 612 if (cData[TCAL].colnum < 0) { 613 findData(TCAL, "SYSCAL_TCAL", TFLOAT);613 findData(TCAL, (char *)"SYSCAL_TCAL", TFLOAT); 614 614 } else if (cALFA_BD) { 615 615 // ALFA BDFITS has a different TCAL with 64 elements - kill it! 616 findData(TCAL, "NO NO NO", TFLOAT);616 findData(TCAL, (char *)"NO NO NO", TFLOAT); 617 617 } 618 618 619 619 if (cALFA_BD) { 620 620 // ALFA BDFITS. 621 findData(AZIMUTH, "CRVAL2B", TFLOAT);622 findData(ELEVATIO, "CRVAL3B", TFLOAT);621 findData(AZIMUTH, (char *)(char *)"CRVAL2B", TFLOAT); 622 findData(ELEVATIO, (char *)(char *)"CRVAL3B", TFLOAT); 623 623 } 624 624 625 625 if (cALFA) { 626 626 // ALFA data. 627 findData(PARANGLE, "PARA_ANG", TFLOAT);627 findData(PARANGLE, (char *)"PARA_ANG", TFLOAT); 628 628 } 629 629 630 630 if (cData[TAMBIENT].colnum < 0) { 631 findData(TAMBIENT, "WEATHER_TEMPERATURE", TFLOAT);631 findData(TAMBIENT, (char *)"WEATHER_TEMPERATURE", TFLOAT); 632 632 } 633 633 634 634 if (cData[PRESSURE].colnum < 0) { 635 findData(PRESSURE, "WEATHER_PRESSURE", TFLOAT);635 findData(PRESSURE, (char *)"WEATHER_PRESSURE", TFLOAT); 636 636 } 637 637 638 638 if (cData[HUMIDITY].colnum < 0) { 639 findData(HUMIDITY, "WEATHER_REL_HUMIDITY", TFLOAT);639 findData(HUMIDITY, (char *)"WEATHER_REL_HUMIDITY", TFLOAT); 640 640 } 641 641 642 642 if (cData[WINDSPEE].colnum < 0) { 643 findData(WINDSPEE, "WEATHER_WIND_SPEED", TFLOAT);643 findData(WINDSPEE, (char *)"WEATHER_WIND_SPEED", TFLOAT); 644 644 } 645 645 646 646 if (cData[WINDDIRE].colnum < 0) { 647 findData(WINDDIRE, "WEATHER_WIND_DIRECTION", TFLOAT);647 findData(WINDDIRE, (char *)"WEATHER_WIND_DIRECTION", TFLOAT); 648 648 } 649 649 … … 768 768 if (cData[DATAXED].colnum > 0) { 769 769 char dataxed[32]; 770 readParm( "DATAXED", TSTRING, dataxed);770 readParm((char *)(char *)"DATAXED", TSTRING, dataxed); 771 771 772 772 sscanf(dataxed, "(%ld,%ld,%ld,%ld,%ld)", cNAxis, cNAxis+1, … … 1014 1014 if (cALFA_CIMA > 1) { 1015 1015 // FFTs per second when the Mock correlator operates in RFI blanking mode. 1016 readData( "PHFFTACC", TFLOAT, 0, &cALFAacc);1016 readData((char *)"PHFFTACC", TFLOAT, 0, &cALFAacc); 1017 1017 } 1018 1018 … … 1051 1051 1052 1052 // Read parameter values. 1053 readParm( "OBSERVER", TSTRING, observer); // Shared.1054 readParm( "PROJID", TSTRING, project); // Shared.1055 readParm( "TELESCOP", TSTRING, telescope); // Core.1053 readParm((char *)(char *)"OBSERVER", TSTRING, observer); // Shared. 1054 readParm((char *)(char *)"PROJID", TSTRING, project); // Shared. 1055 readParm((char *)(char *)"TELESCOP", TSTRING, telescope); // Core. 1056 1056 1057 1057 antPos[0] = 0.0; 1058 1058 antPos[1] = 0.0; 1059 1059 antPos[2] = 0.0; 1060 if (readParm( "ANTENNA_POSITION", TDOUBLE, antPos)) {1061 readParm( "OBSGEO-X", TDOUBLE, antPos); // Additional.1062 readParm( "OBSGEO-Y", TDOUBLE, antPos + 1); // Additional.1063 readParm( "OBSGEO-Z", TDOUBLE, antPos + 2); // Additional.1060 if (readParm((char *)"ANTENNA_POSITION", TDOUBLE, antPos)) { 1061 readParm((char *)"OBSGEO-X", TDOUBLE, antPos); // Additional. 1062 readParm((char *)"OBSGEO-Y", TDOUBLE, antPos + 1); // Additional. 1063 readParm((char *)"OBSGEO-Z", TDOUBLE, antPos + 2); // Additional. 1064 1064 } 1065 1065 1066 1066 if (antPos[0] == 0.0) { 1067 1067 if (cGBT) { 1068 readParm( "SITELONG", TDOUBLE, antPos ) ;1069 readParm( "SITELAT", TDOUBLE, antPos+1 ) ;1070 readParm( "SITEELEV", TDOUBLE, antPos+2 ) ;1068 readParm( (char *)(char *)"SITELONG", TDOUBLE, antPos ) ; 1069 readParm( (char *)(char *)"SITELAT", TDOUBLE, antPos+1 ) ; 1070 readParm( (char *)(char *)"SITEELEV", TDOUBLE, antPos+2 ) ; 1071 1071 Vector<Double> pos( 2 ) ; 1072 1072 pos[0] = Double( antPos[0] ) ; … … 1116 1116 } 1117 1117 1118 readParm( "EQUINOX", TFLOAT, &equinox); // Shared.1118 readParm((char *)(char *)"EQUINOX", TFLOAT, &equinox); // Shared. 1119 1119 if (cStatus == 405) { 1120 1120 // EQUINOX was written as string value in early versions. 1121 1121 cStatus = 0; 1122 1122 char strtmp[32]; 1123 readParm( "EQUINOX", TSTRING, strtmp);1123 readParm((char *)(char *)"EQUINOX", TSTRING, strtmp); 1124 1124 sscanf(strtmp, "%f", &equinox); 1125 1125 } 1126 1126 1127 if (readParm( "RADESYS", TSTRING, radecsys)) { // Additional.1128 if (readParm( "RADECSYS", TSTRING, radecsys)) { // Additional.1127 if (readParm((char *)(char *)"RADESYS", TSTRING, radecsys)) { // Additional. 1128 if (readParm((char *)(char *)"RADECSYS", TSTRING, radecsys)) { // Additional. 1129 1129 strcpy(radecsys, ""); 1130 1130 } 1131 1131 } 1132 1132 1133 if (readParm( "SPECSYS", TSTRING, dopplerFrame)) { // Additional.1133 if (readParm((char *)(char *)"SPECSYS", TSTRING, dopplerFrame)) { // Additional. 1134 1134 // Fallback value. 1135 1135 strcpy(dopplerFrame, "TOPOCENT"); … … 1139 1139 // Added few more codes currently (as of 2009 Oct) used in the GBT 1140 1140 // SDFITS (based io_sdfits_define.pro of GBTIDL). - TT 1141 if (readParm( "VELFRAME", TSTRING, dopplerFrame)) { // Additional.1141 if (readParm((char *)(char *)"VELFRAME", TSTRING, dopplerFrame)) { // Additional. 1142 1142 // No, try digging it out of the CTYPE card (AIPS convention). 1143 1143 char keyw[9], ctype[9]; … … 1211 1211 readTime(1, 1, datobs, utc); 1212 1212 readData(FqRefVal, 1, &refFreq); 1213 readParm( "BANDWID", TDOUBLE, &bandwidth); // Core.1213 readParm((char *)(char *)"BANDWID", TDOUBLE, &bandwidth); // Core. 1214 1214 1215 1215 if (cStatus) { … … 1263 1263 if (cALFA_BD) { 1264 1264 unsigned char invert; 1265 readData( "UPPERSB", TBYTE, irow, &invert);1265 readData((char *)"UPPERSB", TBYTE, irow, &invert); 1266 1266 1267 1267 if (invert) { … … 1867 1867 unsigned char invert; 1868 1868 int anynul, colnum; 1869 findCol( "UPPERSB", &colnum);1869 findCol((char *)"UPPERSB", &colnum); 1870 1870 fits_read_col(cSDptr, TBYTE, colnum, cRow, 1, 1, 0, &invert, &anynul, 1871 1871 &cStatus); … … 2011 2011 // Rescale according to the number of unblanked accumulations. 2012 2012 int colnum, naccum; 2013 findCol( "STAT", &colnum);2013 findCol((char *)"STAT", &colnum); 2014 2014 fits_read_col(cSDptr, TINT, colnum, cRowR, 10*(cTimeIdx-1)+2, 1, 0, 2015 2015 &naccum, &anynul, &cStatus); … … 2631 2631 char chars[32]; 2632 2632 if (cALFA_BD) { 2633 readData( "OBS_NAME", TSTRING, cRow, chars);2633 readData((char *)"OBS_NAME", TSTRING, cRow, chars); 2634 2634 } else { 2635 readData( "SCANTYPE", TSTRING, cRow, chars);2635 readData((char *)"SCANTYPE", TSTRING, cRow, chars); 2636 2636 } 2637 2637 … … 2691 2691 if (cALFA_CIMA > 1) { 2692 2692 int colnum, naccum; 2693 findCol( "STAT", &colnum);2693 findCol((char *)"STAT", &colnum); 2694 2694 fits_read_col(cSDptr, TINT, colnum, cRow, 2, 1, 0, &naccum, &anynul, 2695 2695 &cStatus); … … 2916 2916 MEpoch me( Quantity( dd, "d" ), MEpoch::UTC ) ; 2917 2917 double antPos[3] ; 2918 readParm( "SITELONG", TDOUBLE, antPos ) ;2919 readParm( "SITELAT", TDOUBLE, antPos+1 ) ;2920 readParm( "SITEELEV", TDOUBLE, antPos+2 ) ;2918 readParm( (char *)(char *)"SITELONG", TDOUBLE, antPos ) ; 2919 readParm( (char *)(char *)"SITELAT", TDOUBLE, antPos+1 ) ; 2920 readParm( (char *)(char *)"SITEELEV", TDOUBLE, antPos+2 ) ; 2921 2921 Vector<Double> pos( 2 ) ; 2922 2922 pos[0] = Double( antPos[0] ) ;
Note: See TracChangeset
for help on using the changeset viewer.