void gsdac_getMapVars ( const gsdVars *gsdVars, const char *samMode,
mapVars *mapVars, int *status )
{
const char * cell_sys;
/* Check inherited status */
if ( *status != SAI__OK ) return;
/* Translate Switch mode. */
if ( strncmp ( gsdVars->swMode, "POSITION_SWITCH", 15 ) == 0 ) {
if ( gsdVars->chopping ) {
if ( gsdVars->referenceX == 0.0 && gsdVars->referenceY == 0.0 )
strcpy ( mapVars->swMode, "freqsw" );
else
strcpy ( mapVars->swMode, "pssw" );
} else {
if ( gsdVars->referenceX == 0.0 && gsdVars->referenceY == 0.0 ) {
strcpy ( mapVars->swMode, "freqsw" );
/* Print a message, this was likely intended as a freq. sw. */
msgOutif(MSG__VERB," ", "SWITCH_MODE was POSITION_SWITCH and CHOPPING was 0, this was likely intended to be a frequency switch", status);
} else strcpy ( mapVars->swMode, "pssw" );
}
} else if ( strncmp ( gsdVars->swMode, "BEAMSWITCH", 10 ) == 0 ) {
if ( gsdVars->chopping ) {
strcpy ( mapVars->swMode, "chop" );
} else {
strcpy ( mapVars->swMode, "none" );
msgOutif(MSG__VERB," ", "SWITCH_MODE was BEAMSWITCH and CHOPPING was 0, this may be an error...)", status);
}
} else if ( strncmp ( gsdVars->swMode, "CHOPPING", 8 ) == 0 ) {
strcpy ( mapVars->swMode, "freqsw" );
if ( gsdVars->chopping ) {
msgOutif(MSG__VERB," ", "SWITCH_MODE was CHOPPING and CHOPPING was 1, this appears to be a misconfigured frequency switch", status);
}
} else if ( strncmp ( gsdVars->swMode, "NO_SWITCH", 7 ) == 0 ) {
strcpy ( mapVars->swMode, "none" );
if ( gsdVars->chopping ) {
msgOutif(MSG__VERB," ", "SWITCH_MODE was NO_SWITCH and CHOPPING was 1, this may be an error...", status);
}
} else {
*status = SAI__ERROR;
msgSetc ( "SWITCHMODE", gsdVars->swMode );
errRep ( "gsdac_getMapVars", "Couldn't identify switch mode ^SWITCHMODE", status );
return;
}
/* Get the chopping parameters for grid beamswitches
and samples. */
if ( strcmp ( mapVars->swMode, "chop" ) == 0 ) {
if ( strncmp ( gsdVars->chopCoords, "AZ", 2 ) == 0 )
strcpy ( mapVars->chopCrd, "AZEL" );
else if ( strncmp ( gsdVars->chopCoords, "EQ", 2 ) == 0 )
strcpy ( mapVars->chopCrd, "HADEC" );
else if ( strncmp ( gsdVars->chopCoords, "RB", 2 ) == 0 )
strcpy ( mapVars->chopCrd, "B1950" );
else if ( strncmp ( gsdVars->chopCoords, "RJ", 2 ) == 0 )
strcpy ( mapVars->chopCrd, "J2000" );
else if ( strncmp ( gsdVars->chopCoords, "RD", 2 ) == 0 )
strcpy ( mapVars->chopCrd, "APP" );
else if ( strncmp ( gsdVars->chopCoords, "GA", 2 ) == 0 )
strcpy ( mapVars->chopCrd, "GAL" );
else {
strcpy ( mapVars->chopCrd, "" );
msgOutif(MSG__VERB," ",
"Couldn't identify chop coordinates, continuing anyway", status);
}
}
/* Get the local offset coordinate system. */
cell_sys = gsdac_code2tcssys( gsdVars->cellCode, status );
if (*status == SAI__OK && strlen(cell_sys) == 0) {
msgOutif(MSG__VERB," ",
"Couldn't identify cell coordinates, continuing anyway", status);
}
one_strlcpy( mapVars->loclCrd, cell_sys,
sizeof(mapVars->loclCrd), status );
/* Convert to ACSIS formatted string. */
sprintf ( mapVars->skyRefX, "[OFFSET] %.0f [%s]",
gsdVars->referenceX, mapVars->loclCrd );
sprintf ( mapVars->skyRefY, "[OFFSET] %.0f [%s]",
gsdVars->referenceY, mapVars->loclCrd );
/* Get the scanning coordinates. */
strcpy ( mapVars->scanCrd, mapVars->loclCrd );
/* Get the map and scan parameters for rasters. */
if ( strcmp ( samMode, "scan" ) == 0
&& strcmp ( mapVars->swMode, "pssw" ) == 0 ) {
/* Get the map height and width. */
mapVars->mapHght = gsdVars->nMapPtsX * gsdVars->cellX;
mapVars->mapWdth = gsdVars->nMapPtsY * gsdVars->cellY;
/* Get the scan velocity and spacing. */
if ( strncmp ( gsdVars->obsDirection, "HORIZONTAL", 10 ) == 0 ) {
mapVars->scanVel = gsdVars->cellX * gsdVars->nMapPtsX /
gsdVars->scanTime;
mapVars->scanDy = gsdVars->cellY;
} else if ( strncmp ( gsdVars->obsDirection, "VERTICAL", 8 ) == 0 ) {
mapVars->scanVel = gsdVars->cellY * gsdVars->nMapPtsY /
gsdVars->scanTime;
mapVars->scanDy = gsdVars->cellX;
} else {
*status = SAI__ERROR;
errRep ( "gsdac_getMapVars", "Error getting scan velocity",
status );
return;
}
if ( gsdVars->scanRev ) strcpy ( mapVars->scanPat, "BOUSTROPHEDON" );
else strcpy ( mapVars->scanPat, "RASTER" );
}
mapVars->scanPA = gsdVars->cellV2X;
mapVars->mapPA = mapVars->scanPA - 90.0;
}
void gsdac_putFits ( const gsdVars *gsdVars, const int subBandNum,
const int nSubsys, const int obsNum,
const int utDate, const int nSteps,
const char *backend, const int recepsUsed,
char *recepNames[], const char *samMode,
const char *obsType, const dateVars *dateVars,
const mapVars *mapVars, const double *lineFreqs,
const double *IFFreqs, const gsdWCS *wcs,
AstFitsChan *fitschan, int *status )
{
/* Local variables */
char bwMode[SZFITSTR]; /* ACSIS total bandwidth setup */
char curChar; /* character pointer */
int day; /* days for time conversion. */
char doppler[SZFITSTR]; /* doppler velocity definition */
double etal; /* telescope efficiency */
int hour; /* hours for time conversion. */
int i; /* loop counter */
float IFchanSp; /* TOPO IF channel spacing (Hz) */
char instrume[SZFITSTR]; /* front-end receiver */
double intTime; /* total time spent integrating (s) */
int josMult; /* times around the jiggle (i.e. 1) */
int josMin; /* minimum step time psw */
int min; /* minutes for time conversion. */
const char *molecule; /* target molecular species */
int month; /* months for time conversion. */
int nMix; /* number of mixers */
double nRefStep; /* number of nod sets repeated */
char object[SZFITSTR]; /* object of interest */
char object2[SZFITSTR]; /* object of interest (second half of name) */
double obsgeo[3]; /* cartesian coordinates of telescope */
char obsIDSS[SZFITSTR]; /* unique observation number + subsystem
number in format
INSTR_NNNNN_YYYYMMDDTHHMMSS_N */
char obsSB[SZFITSTR]; /* observed sideband */
int parse; /* flag for incorrect date string. */
char recipe[25]; /* ORAC-DR recipe name */
char recptors[SZFITSTR]; /* active FE receptor IDs for this obs */
double refChan; /* reference IF channel no. */
int seeingok; /* True if the seeing measurement is okay */
char seeDatSt[SZFITSTR]; /* time of seeingSt in format
YYYY-MM-DDTHH:MM:SS */
char sSysObs[SZFITSTR]; /* spectral ref. frame during observation */
int standard; /* true for spectral line standards */
int startIdx; /* index in pattern at start of observation */
int stBetRef; /* max number of steps between refs */
int stBetRef_defined; /* Use the stBetRef value? */
double stepTime; /* RTS step time */
char subBands[SZFITSTR]; /* ACSIS sub-band set-up */
int tauok; /* True if the tau measurement is okay */
char tauDatSt[SZFITSTR]; /* time of tau225St observation in
format YYYY-MM-DDTHH:MM:SS */
const char *transiti; /* target transition for molecule */
int year; /* year for time conversion. */
/* Check inherited status */
if ( *status != SAI__OK ) return;
/* Get the cartesian coordinates of the telescope's location. */
smf_terr( gsdVars->telLatitude * DD2R, gsdVars->telHeight * 1000.0,
-gsdVars->telLongitude * DD2R, obsgeo );
/* Get the telescope efficiency and convert from percentage to decimal */
etal = gsdVars->etal / 100.0;
/* Obs Id, Date, Pointing Info */
/* Truncate the object names and concatenate. */
cnfImprt ( gsdVars->object1, 16, object );
if ( strncmp ( gsdVars->object2, " ", 1 ) != 0 ) {
cnfImprt ( gsdVars->object2, 16, object2 );
strcat ( object, ", " );
strcat ( object, object2 );
}
/* Determine if this is a spectral line standard. */
/* Loop through the array of standards and see if the object matches
one of them. */
standard = 0;
i = 0;
while ( strcmp ( standards[i], "" ) != 0 ) {
if ( strcmp ( object, standards[i] ) == 0 ) {
standard = 1;
}
i++;
}
/* Copy the obsID into the obsIDSS and add the subsystem number. */
sprintf ( obsIDSS, "%s_%i", dateVars->obsID,
subBandNum % nSubsys + 1 );
/* Integration time related. */
/* Get the sum of the integration times. This is either the
sum of the elements of the intTimes table for grids, or
the scan time times the number of scans for rasters. */
if ( gsdVars->obsContinuous ) {
intTime = gsdVars->scanTime * gsdVars->nScan;
} else {
intTime = gsdVars->nCycle * gsdVars->cycleTime;
}
/* ACSIS Specific. */
/* Get the molecule and transition. */
smf_get_moltrans ( lineFreqs[subBandNum] * 1000.0, &molecule, &transiti,
status );
/* Get the bandwidth setup. */
/***** NOTE: may be different for rxb widebands *****/
sprintf ( bwMode, "%iMHzx%i", (int)(gsdVars->bandwidths[subBandNum]),
gsdVars->BEChans[subBandNum] );
/***** NOTE: Possibly undef? *****/
strcpy ( subBands, bwMode );
/* Get the reference channel (offset from channel which contains
IFFreq). */
refChan = ( (double)( gsdVars->BEChans[subBandNum] ) / 2.0 ) +
( ( IFFreqs[subBandNum] - gsdVars->totIFs[subBandNum] ) /
( fabs(gsdVars->freqRes[subBandNum]) / 1000.0 ) );
IFchanSp = gsdVars->freqRes[subBandNum] * 1000000.0;
strcpy( recipe, "REDUCE_SCIENCE" );
/* Is this a pointing? */
if ( strncasecmp( gsdVars->backend, "DAS", 3 ) == 0 &&
strncasecmp( gsdVars->obsType, "FIVEPOINT", 9 ) == 0 )
strcpy( recipe, "REDUCE_POINTING" );
else if ( strncasecmp( gsdVars->obsType, "focus", 5 ) == 0 )
strcpy( recipe, "REDUCE_FOCUS" );
/* Solar-system objects have a different default recipe. */
else if ( strncasecmp( gsdVars->object1, "SUN", 3 ) == 0 ||
strncasecmp( gsdVars->object1, "MERCURY", 7 ) == 0 ||
strncasecmp( gsdVars->object1, "VENUS", 5 ) == 0 ||
strncasecmp( gsdVars->object1, "MOON", 4 ) == 0 ||
strncasecmp( gsdVars->object1, "MARS", 4 ) == 0 ||
strncasecmp( gsdVars->object1, "JUPITER", 7 ) == 0 ||
strncasecmp( gsdVars->object1, "SATURN", 6 ) == 0 ||
strncasecmp( gsdVars->object1, "TITAN", 5 ) == 0 ||
strncasecmp( gsdVars->object1, "URANUS", 6 ) == 0 ||
strncasecmp( gsdVars->object1, "NEPTUNE", 7 ) == 0 ||
strncasecmp( gsdVars->object1, "PLUTO", 5 ) == 0 )
strcpy( recipe, "REDUCE_SCIENCE_CONTINUUM" );
/* Try to catch most galaxies but exclude regions of the Galaxy. */
else if ( gsdVars->velocity > 120 )
strcpy( recipe, "REDUCE_SCIENCE_BROADLINE" );
else if ( strncmp( mapVars->swMode, "freq", 4 ) == 0 )
strcpy( recipe, "REDUCE_SCIENCE_FSW" );
/* FE Specific. */
/* We need to take the opportunity to correct the frontend name.
In many cases, especially in B-band, the instrument names are
incorrect. */
gsdac_getRealInstrumentName( gsdVars, instrume, sizeof(instrume), status );
/* Get the number of mixers, 2 for RXB3 & RXW, 1 for others.
but note that RXB3I has 1. We know that instrume is nul-terminated.
Both RXWC and RXWD have two mixers so just check first 3 characters.
*/
if ( strcmp ( "RXB3", instrume ) == 0 ||
strncmp ( "RXW", instrume, 3 ) == 0 ) {
nMix = 2;
} else {
nMix = 1;
}
/* Get the observed sideband (-ve value = LSB, +ve value = USB ). */
if ( gsdVars->sbSigns[subBandNum] > 0 ) strcpy ( obsSB, "USB" );
else strcpy ( obsSB, "LSB" );
/* Get the names of the receptors. */
strcpy ( recptors, recepNames[0] );
for ( i = 1; i < recepsUsed; i++ ) {
strcat ( recptors, " " );
strcat ( recptors, recepNames[i] );
}
/* Truncate the doppler velocity definition and set
to lowercase. */
cnfImprt ( gsdVars->velDefn, 16, doppler );
i = 0;
curChar = doppler[i];
while ( curChar != '\0' ) {
doppler[i] = tolower(curChar);
i++;
curChar = doppler[i];
}
/***** NOTE : Possibly comes from VEL_REF (c12vref). */
strcpy ( sSysObs, "TOPOCENT" );
/* Environmental data. */
/* The GSD magic value for bad seeing and tau seems to be VAL__BADR */
tauDatSt[0] = '\0'; /* make sure it is initialised */
if (gsdVars->tau225 == VAL__BADR) {
tauok = 0;
} else {
tauok = 1;
/* Convert dates from YYMMDDHHMMSS to
YYYY-MM-DDTHH:MM:SS. */
parse = sscanf ( gsdVars->tauTime, "%02d%02d%02d%02d%02d", &year,
&month, &day, &hour, &min );
if ( parse == 0 || parse == EOF ) {
msgOutf ( FUNC_NAME, "Couldn't convert CSO tau time of '%s', continuing anyway.",
status, gsdVars->tauTime );
} else {
/* Inelegant method to get YYYY from YY. */
if ( year > 70 ) year = year + 1900;
else year = year + 2000;
sprintf ( tauDatSt, "%04d-%02d-%02dT%02d:%02d:00",
year, month, day, hour, min );
}
}
/* The GSD magic value for bad seeing and tau seems to be VAL__BADR */
seeDatSt[0] = '\0'; /* Make sure it is initialised */
if (gsdVars->seeing == VAL__BADR) {
seeingok = 0;
} else {
seeingok = 1;
/* Convert dates from YYMMDDHHMMSS to
YYYY-MM-DDTHH:MM:SS. */
parse = sscanf ( gsdVars->seeTime, "%02d%02d%02d%02d%02d", &year,
&month, &day, &hour, &min );
if ( parse == 0 || parse == EOF ) {
msgOutiff(MSG__VERB," ",
"Couldn't convert seeing time of '%s', continuing anyway.", status,
gsdVars->seeTime );
} else {
/* Kludge to get YYYY from YY. */
if ( year > 70 ) year = year + 1900;
else year = year + 2000;
sprintf ( seeDatSt, "%04d-%02d-%02dT%02d:%02d:00",
year, month, day, hour, min );
}
}
/* JOS parameters */
/* STEPTIME is SCAN_TIME divided by the number of points in a scan
for rasters. For PSW spectra it is SCAN_TIME/2 and for BSW
it is SCAN_TIME/4 -- because acsis expects a minimum unit
of ABBA for the nod --*/
if ( strcmp ( samMode, "scan" ) == 0 ) {
if ( strncmp ( gsdVars->obsDirection, "HORIZONTAL", 10 ) == 0 ) {
stepTime = gsdVars->scanTime / gsdVars->nMapPtsX;
} else {
stepTime = gsdVars->scanTime / gsdVars->nMapPtsY;
}
josMin = 1;
josMult = 0;
} else if ( strcmp ( mapVars->swMode, "chop" ) == 0 ) {
stepTime = gsdVars->scanTime / 4.0;
josMin = 0;
josMult = 1;
} else {
stepTime = gsdVars->scanTime / 2.0;
josMin = stepTime;
josMult = 0;
}
/* Get the length of time in the reference, and the number
of steps between references. set up defaults first to avoid compiler
warnings. */
nRefStep = 1;
stBetRef = 1;
if ( gsdVars->obsContinuous ) {
/* For rasters, determine the number of time scanning
each point in one row from the total time for
the row / number of points in the row. The length
of time in the reference is then sqrt (number of
points in the row) * (time per point). */
nRefStep = sqrt ( (double)(gsdVars->nScanPts) ) *
( (double)(gsdVars->scanTime) / (double)(gsdVars->nScanPts) );
stBetRef = gsdVars->nScanPts;
}
/* Set a flag indicating if the stBetRef value is defined. */
stBetRef_defined = ( strcmp ( mapVars->swMode, "chop" ) != 0 );
/* Get the starting index into the pattern. */
gsdac_getStartIdx ( gsdVars, samMode, &startIdx, status );
/************************************/
/* WRITE OUT FITS HEADERS */
/************************************/
astSetFitsS ( fitschan, "TELESCOP", gsdVars->telName,
"Name of Telescope", 0 );
astSetFitsS ( fitschan, "ORIGIN", "Joint Astronomy Centre, Hilo",
"Origin of file", 0 );
astSetFitsF ( fitschan, "ALT-OBS", gsdVars->telHeight * 1000.0,
"[m] Height of observation above sea level", 0 );
astSetFitsF ( fitschan, "LAT-OBS", gsdVars->telLatitude,
"[deg] Latitude of Observatory", 0 );
astSetFitsF ( fitschan, "LONG-OBS", -(gsdVars->telLongitude),
"[deg] East longitude of Observatory", 0 );
astSetFitsF ( fitschan, "OBSGEO-X", obsgeo[0],
"[m]", 0 );
astSetFitsF ( fitschan, "OBSGEO-Y", obsgeo[1],
"[m]", 0 );
astSetFitsF ( fitschan, "OBSGEO-Z", obsgeo[2],
"[m]", 0 );
astSetFitsF ( fitschan, "ETAL", etal,
"Telescope efficiency", 0 );
/* OMP and ORAC-DR Specific */
astSetFitsCM ( fitschan,
"---- OMP and ORAC-DR Specific ----", 0 );
astSetFitsS ( fitschan, "PROJECT", gsdVars->project,
"PATT number", 0 );
/***** NOTE: possibly REDUCE_POINTING for spectral 5 points *****/
astSetFitsS ( fitschan, "RECIPE", recipe, "ORAC-DR recipe", 0 );
astSetFitsU ( fitschan, "DRGROUP", "Data Reduction group ID", 0 );
astSetFitsS ( fitschan, "MSBID", "GSD2ACSIS_V1",
"ID of minimum schedulable block", 0 );
astSetFitsU ( fitschan, "MSBTID", "Transaction ID of MSB", 0 );
astSetFitsU ( fitschan, "SURVEY", "Survey Name", 0 );
astSetFitsU ( fitschan, "RMTAGENT", "name of Remote Agent", 0 );
astSetFitsU ( fitschan, "AGENTID", "Unique identifier for remote agent", 0 );
/* Obs Id, Date, Pointing Info */
astSetFitsCM ( fitschan,
"---- Obs Id, Date, pointing Info ----", 0 );
astSetFitsS ( fitschan, "OBJECT", object,
"Object of interest", 0 );
astSetFitsL ( fitschan, "STANDARD", standard,
"True if the spectral line is a standard", 0 );
astSetFitsI ( fitschan, "OBSNUM", obsNum,
"Observation number", 0 );
astSetFitsI ( fitschan, "NSUBSCAN", 1,
"Sub-scan number", 0 );
astSetFitsL ( fitschan, "OBSEND", 1,
"True if the file is the last in current observation", 0 );
astSetFitsI ( fitschan, "UTDATE", utDate,
"UT Date as integer in yyyymmdd format", 0 );
astSetFitsS ( fitschan, "DATE-OBS", dateVars->dateObs,
"UTC Datetime of start of observation", 0 );
astSetFitsS ( fitschan, "DATE-END", dateVars->dateEnd,
"UTC Datetime of end of observation", 0 );
astSetFitsF ( fitschan, "DUT1", gsdVars->obsUT1C,
"[d] UT1-UTC correction", 0 );
astSetFitsS ( fitschan, "OBSID", dateVars->obsID,
"Unique observation identifier", 0 );
astSetFitsS ( fitschan, "OBSIDSS", obsIDSS,
"Unique observation + subsystem ID", 0 );
/***** NOTE: possibly same as REFRECEP *****/
astSetFitsU ( fitschan, "INSTAP", "Receptor at tracking centre (if any)", 0 );
astSetFitsF ( fitschan, "INSTAP_X", 0.0,
"[arcsec] Aperture X off. rel. to instr centre", 0 );
astSetFitsF ( fitschan, "INSTAP_Y", 0.0,
"[arcsec] Aperture Y off. rel. to instr centre", 0 );
/* The following 6 cards are just placeholders and will be
updated later. */
astSetFitsU ( fitschan, "AMSTART", "Airmass at start of observation", 0 );
astSetFitsU ( fitschan, "AMEND", "Airmass at end of observation", 0 );
astSetFitsU ( fitschan, "AZSTART", "[deg] Azimuth at start of observation", 0 );
astSetFitsU ( fitschan, "AZEND", "[deg] Azimuth at end of observation", 0 );
astSetFitsU ( fitschan, "ELSTART", "[deg] Elevation at start of observation", 0 );
astSetFitsU ( fitschan, "ELEND", "[deg] Elevation at end of observation", 0 );
/* Base position in tracking coordinates */
astSetFitsS( fitschan, "TRACKSYS", gsdac_code2tcssys(gsdVars->centreCode, status),
"TCS tracking coordinate system", 0 );
astSetFitsF ( fitschan, "BASEC1", wcs->baseTr1 * AST__DR2D,
"[deg] TCS BASE position (longitude) in TRACKSYS", 0 );
astSetFitsF ( fitschan, "BASEC2", wcs->baseTr2 * AST__DR2D,
"[deg] TCS BASE position (latitude) in TRACKSYS", 0 );
astSetFitsS ( fitschan, "HSTSTART", dateVars->HSTstart,
"HST at start of observation", 0 );
astSetFitsS ( fitschan, "HSTEND", dateVars->HSTend,
"HST at end of observation", 0 );
astSetFitsS ( fitschan, "LSTSTART", dateVars->LSTstart,
"LST at start of observation", 0 );
astSetFitsS ( fitschan, "LSTEND", dateVars->LSTend,
"LST at end of observation", 0 );
/* Integration time related. */
astSetFitsCM ( fitschan,
"---- Integration time related ----", 0 );
astSetFitsF ( fitschan, "INT_TIME", intTime,
"Time spent integrating, entire", 0 );
/* DAS Specific */
astSetFitsCM ( fitschan,
"---- DAS Specific ----", 0 );
astSetFitsI ( fitschan, "SPECID", gsdVars->BESubsys[subBandNum],
"Spectrum ID for this subband", 0 );
/* ACSIS Specific. */
astSetFitsCM ( fitschan,
"---- ACSIS Specific ----", 0 );
astSetFitsS ( fitschan, "BACKEND", backend,
"Name of the backend", 0 );
astSetFitsS ( fitschan, "MOLECULE", molecule,
"Target molecular species", 0 );
astSetFitsS ( fitschan, "TRANSITI", transiti,
"Target transition for MOLECULE", 0 );
astSetFitsU ( fitschan, "DRRECIPE", "ACSIS-DR recipe name", 0 );
astSetFitsS ( fitschan, "BWMODE", bwMode,
"Bandwidth setup", 0 );
astSetFitsI ( fitschan, "SUBSYSNR", ( subBandNum % nSubsys ) + 1,
"Sub-system number", 0 );
astSetFitsS ( fitschan, "SUBBANDS", bwMode,
"Sub-band setup", 0 );
astSetFitsI ( fitschan, "NSUBBAND", 1,
"Number of subbands", 0 );
astSetFitsF ( fitschan, "SUBREFP1", refChan,
"Reference channel for subband1", 0 );
astSetFitsU ( fitschan, "SUBREFP2", "Reference channel for subband2", 0 );
astSetFitsI ( fitschan, "NCHNSUBS", gsdVars->BEChans[subBandNum],
"Number of subbands", 0 );
astSetFitsF ( fitschan, "REFCHAN", refChan,
"Reference IF channel No.", 0 );
astSetFitsF ( fitschan, "IFCHANSP", IFchanSp,
"[Hz] TOPO IF channel spacing (signed)", 0 );
astSetFitsS ( fitschan, "FFT_WIN", "truncate",
"Type of window used for FFT", 0 );
if ( strncmp ( gsdVars->backend, "DAS", 3 ) == 0 ) {
astSetFitsF ( fitschan, "BEDEGFAC", 1.15,
"Backend degradation factor", 0 );
} else if ( strncmp ( gsdVars->backend, "AOSC", 4 ) == 0 ) {
astSetFitsF ( fitschan, "BEDEGFAC", 1.0,
"Backend degradation factor", 0 );
} else {
astSetFitsU ( fitschan, "BEDEGFAC", "Backend degradation factor", 0 );
}
astSetFitsU ( fitschan, "MSROOT", "Root name of raw measurement sets", 0 );
/* FE Specific. */
astSetFitsCM ( fitschan,
"---- FE Specific ----", 0 );
astSetFitsS ( fitschan, "INSTRUME", instrume,
"Front-end receiver", 0 );
astSetFitsS ( fitschan, "SB_MODE", gsdVars->sbMode,
"Sideband mode", 0 );
astSetFitsF ( fitschan, "IFFREQ", IFFreqs[subBandNum],
"[GHz] IF Frequency", 0 );
astSetFitsI ( fitschan, "N_MIX", nMix,
"No. of mixers", 0 );
astSetFitsS ( fitschan, "OBS_SB", obsSB,
"The observed sideband", 0 );
astSetFitsF ( fitschan, "LOFREQS", gsdVars->LOFreqs[subBandNum],
"[GHz] LO Frequency at start of obs", 0 );
astSetFitsF ( fitschan, "LOFREQE", gsdVars->LOFreqs[subBandNum],
"[GHz] LO Frequency at end of obs", 0 );
astSetFitsS ( fitschan, "RECPTORS", recptors,
"Active FE receptor IDs for this obs", 0 );
astSetFitsS ( fitschan, "REFRECEP", recepNames[0],
"Receptor with unit sensitivity", 0 );
if ( strcmp ( samMode, "grid" ) == 0 ) {
astSetFitsU ( fitschan, "MEDTSYS", "[K] Median of the T-sys across all receptors", 0 );
} else {
astSetFitsF ( fitschan, "MEDTSYS",
gsdVars->sourceSysTemps[subBandNum],
"[K] Median of the T-sys across all receptors", 0 );
}
astSetFitsS ( fitschan, "TEMPSCAL", "TA*",
"Temperature scale in use", 0 );
astSetFitsS ( fitschan, "DOPPLER", doppler,
"Doppler velocity definition", 0 );
astSetFitsS ( fitschan, "SSYSOBS", sSysObs,
"Spectral ref. frame during observation", 0 );
/* Environmental data. */
astSetFitsCM ( fitschan,
"---- Environmental Data ----", 0 );
astSetFitsF ( fitschan, "ATSTART", gsdVars->tamb,
"[degC] Air temp at start of observation", 0 );
astSetFitsF ( fitschan, "ATEND", gsdVars->tamb,
"[degC] Air temp at end of observation", 0 );
astSetFitsF ( fitschan, "HUMSTART",gsdVars->hamb ,
"Rel Humidity at observation start", 0 );
astSetFitsF ( fitschan, "HUMEND", gsdVars->hamb,
"Rel Humidity observation end", 0 );
astSetFitsF ( fitschan, "BPSTART", gsdVars->pamb,
"[mbar] Pressure at observation start", 0 );
astSetFitsF ( fitschan, "BPEND", gsdVars->pamb,
"[mbar] Pressure at observation end", 0 );
astSetFitsU ( fitschan, "WNDSPDST", "[km/h] Wind Speed at obs start", 0 );
astSetFitsU ( fitschan, "WNDSPDEN", "[km/h] Wind Speed at obs end", 0 );
astSetFitsU ( fitschan, "WNDDIRST", "[deg] Wind direction, azimuth at obs start", 0 );
astSetFitsU ( fitschan, "WNDDIREN", "[deg] Wind direction, azimuth at obs end", 0 );
if (tauok) {
astSetFitsF ( fitschan, "TAU225ST", gsdVars->tau225,
"Tau at 225 GHz from CSO at start", 0 );
astSetFitsF ( fitschan, "TAU225EN", gsdVars->tau225,
"Tau at 225 GHz from CSO at end", 0 );
} else {
astSetFitsU ( fitschan, "TAU225ST",
"Tau at 225 GHz from CSO at start", 0 );
astSetFitsU ( fitschan, "TAU225EN",
"Tau at 225 GHz from CSO at end", 0 );
}
if (tauok && tauDatSt[0]) {
astSetFitsS ( fitschan, "TAUDATST", tauDatSt,
"Time of TAU225ST observation", 0 );
astSetFitsS ( fitschan, "TAUDATEN", tauDatSt,
"Time of TAU225EN observation", 0 );
} else {
astSetFitsU ( fitschan, "TAUDATST",
"Time of TAU225ST observation", 0 );
astSetFitsU ( fitschan, "TAUDATEN",
"Time of TAU225EN observation", 0 );
}
astSetFitsS ( fitschan, "TAUSRC", "CSO225GHZ",
"Source of the TAU225 value", 0 );
astSetFitsU ( fitschan, "WVMTAUST", "186GHz Tau from JCMT WVM at start", 0 );
astSetFitsU ( fitschan, "WVMTAUEN", "185GHz Tau from JCMT WVM at end", 0 );
astSetFitsU ( fitschan, "WVMDATST", "Time of WVMTAUST", 0 );
astSetFitsU ( fitschan, "WVMDATEN", "Time of WVMTAUEN", 0 );
if ( seeingok ) {
astSetFitsF ( fitschan, "SEEINGST", gsdVars->seeing,
"[arcsec] SAO atmospheric seeing (start)", 0 );
astSetFitsF ( fitschan, "SEEINGEN", gsdVars->seeing,
"[arcsec] SAO atmospheric seeing (end)", 0 );
} else {
/* Use undefined header value if there was no measurement */
astSetFitsU ( fitschan, "SEEINGST",
"[arcsec] SAO atmospheric seeing (start)", 0 );
astSetFitsU ( fitschan, "SEEINGEN",
"[arcsec] SAO atmospheric seeing (end)", 0 );
}
if (seeingok && seeDatSt[0] ) {
astSetFitsS ( fitschan, "SEEDATST", seeDatSt,
"Date/Time of SEEINGST", 0 );
astSetFitsS ( fitschan, "SEEDATEN", seeDatSt,
"Date/Time of SEEINGEN", 0 );
} else {
astSetFitsU ( fitschan, "SEEDATST",
"Date/Time of SEEINGST", 0 );
astSetFitsU ( fitschan, "SEEDATEN",
"Date/Time of SEEINGEN", 0 );
}
astSetFitsU ( fitschan, "FRLEGTST", "[degC] Mean Front leg temperature - Start", 0 );
astSetFitsU ( fitschan, "FRLEGTEN", "[degC] Mean Front leg temperature - End", 0 );
astSetFitsU ( fitschan, "BKLEGTST", "[degC] Mean Back leg temperature - Start", 0 );
astSetFitsU ( fitschan, "BKLEGTEN", "[degC] Mean Back leg temperature - End", 0 );
/* Switching and Map setup for the observation. */
astSetFitsCM ( fitschan,
"---- Switching and Map setup for the observationi ----", 0 );
astSetFitsS ( fitschan, "SAM_MODE", samMode,
"Sampling Mode", 0 );
astSetFitsS ( fitschan, "SW_MODE", mapVars->swMode,
"Switch Mode", 0 );
astSetFitsS ( fitschan, "SKYREFX", mapVars->skyRefX,
"X co-ord of Reference Position", 0 );
astSetFitsS ( fitschan, "SKYREFY", mapVars->skyRefY,
"Y co-ord of Reference Position", 0 );
astSetFitsS ( fitschan, "OBS_TYPE", obsType,
"Type of observation", 0 );
if ( strcmp ( samMode, "grid" ) == 0
&& strcmp ( mapVars->swMode, "chop" ) == 0 ) {
astSetFitsS ( fitschan, "CHOP_CRD", mapVars->chopCrd,
"Chopping co-ordinate system", 0 );
astSetFitsF ( fitschan, "CHOP_FRQ", gsdVars->chopFrequency,
"[Hz] Chop frequency", 0 );
astSetFitsF ( fitschan, "CHOP_PA", gsdVars->chopPA,
"[deg] Chop PA; 0=in lat, 90=in long", 0 );
astSetFitsF ( fitschan, "CHOP_THR", gsdVars->chopThrow,
"[arcsec] Chop throw", 0 );
} else {
astSetFitsU ( fitschan, "CHOP_CRD", "Chopping co-ordinate system", 0 );
astSetFitsU ( fitschan, "CHOP_FRQ", "[Hz] Chop frequency", 0 );
astSetFitsU ( fitschan, "CHOP_PA", "[deg] Chop PA; 0=in lat, 90=in long", 0 );
astSetFitsU ( fitschan, "CHOP_THR", "[arcsec] Chop throw", 0 );
}
astSetFitsU ( fitschan, "ROT_CRD", "Coordinate frame of image rotator", 0 );
astSetFitsU ( fitschan, "ROT_PA", "[[deg] Angle of image rotator", 0 );
astSetFitsU ( fitschan, "JIGL_CNT", "Number of offsets in jiggle pattern", 0 );
astSetFitsU ( fitschan, "JIGL_NAM", "File containing the jiggle offsets", 0 );
astSetFitsU ( fitschan, "JIG_PA", "[deg] Jiggle PA; 0=in lat, 90=in long", 0 );
astSetFitsU ( fitschan, "JIG_CRD", "Jiggling co-ordinate system", 0 );
astSetFitsU ( fitschan, "JIG_SCAL", "Scale size of jiggle pattern", 0 );
if ( strcmp ( samMode, "scan" ) == 0
&& strcmp ( mapVars->swMode, "pssw" ) == 0 ) {
astSetFitsF ( fitschan, "MAP_HGHT", mapVars->mapHght,
"[arcsec] Requested height of map", 0 );
astSetFitsF ( fitschan, "MAP_PA", mapVars->mapPA,
"[deg] Requested PA of map", 0 );
astSetFitsF ( fitschan, "MAP_WDTH", mapVars->mapWdth,
"[arcsec] Requested width of map", 0 );
astSetFitsS ( fitschan, "LOCL_CRD", mapVars->loclCrd,
"Local offset/map PA co-ordinate system", 0 );
astSetFitsF ( fitschan, "MAP_X", gsdVars->centreOffsetX,
"[arcsec] Requested map offset from telescope centre", 0 );
astSetFitsF ( fitschan, "MAP_Y", gsdVars->centreOffsetY,
"[arcsec] Requested map offset from telescope centre", 0 );
astSetFitsS ( fitschan, "SCAN_CRD", mapVars->scanCrd,
"Co-ordinate system for scan", 0 );
astSetFitsF ( fitschan, "SCAN_VEL", mapVars->scanVel,
"[arcsec/sec] Scan velocity along scan direction", 0 );
astSetFitsF ( fitschan, "SCAN_DY", mapVars->scanDy,
"[arcsec] Scan spacing perp. to scan", 0 );
astSetFitsS ( fitschan, "SCAN_PAT", mapVars->scanPat,
"Scan pattern name", 0 );
} else {
astSetFitsU ( fitschan, "MAP_HGHT", "[arcsec] Requested height of map", 0 );
astSetFitsU ( fitschan, "MAP_PA", "[deg] Requested PA of map", 0 );
astSetFitsU ( fitschan, "MAP_WDTH", "[arcsec] Requested width of map", 0 );
astSetFitsU ( fitschan, "LOCL_CRD", "Local offset/map PA co-ordinate system", 0 );
astSetFitsU ( fitschan, "MAP_X", "[arcsec] Requested map offset from telescope centre", 0 );
astSetFitsU ( fitschan, "MAP_Y", "[arcsec] Requested map offset from telescope centre", 0 );
astSetFitsU ( fitschan, "SCAN_CRD", "Co-ordinate system for scan", 0 );
astSetFitsU ( fitschan, "SCAN_VEL", "[arcsec/sec] Scan velocity along scan direction", 0 );
astSetFitsU ( fitschan, "SCAN_DY", "[arcsec] Scan spacing perp. to scan", 0 );
astSetFitsU ( fitschan, "SCAN_PAT", "Scan pattern name", 0 );
}
/* SMU */
astSetFitsCM ( fitschan,
"---- SMU ----", 0 );
astSetFitsF ( fitschan, "ALIGN_DX", gsdVars->smuDX,
"SMU tables X axis focus offset", 0 );
astSetFitsF ( fitschan, "ALIGN_DY", gsdVars->smuDY,
"SMU tables Y axis focus offset", 0 );
astSetFitsF ( fitschan, "FOCUS_DZ", gsdVars->smuDZ,
"SMU tables Z axis focus offset", 0 );
astSetFitsF ( fitschan, "DAZ", gsdVars->smuOffsEW,
"SMU azimuth pointing offset", 0 );
astSetFitsF ( fitschan, "DEL", gsdVars->smuOffsNS,
"SMU elevation pointing offset", 0 );
astSetFitsF ( fitschan, "UAZ", gsdVars->errAz,
"User azimuth pointing offset", 0 );
astSetFitsF ( fitschan, "UEL", gsdVars->errEl,
"User elevation pointing offset", 0 );
/* JOS parameters */
astSetFitsCM ( fitschan,
"---- JOS parameters ----", 0 );
astSetFitsF ( fitschan, "STEPTIME", stepTime,
"RTS step time during an RTS sequence", 0 );
/* Since coadding is happening in the correlator before writing out
** the spectrum, in ACSIS land the number of cycles is 1 instead
** of gsdVars->nCycle
*/
astSetFitsI ( fitschan, "NUM_CYC", 1,
"Number of times to repeat entire recipe", 0 );
astSetFitsI ( fitschan, "NUM_NODS", 1,
"Number of times to repeat nod set", 0 );
if (josMult == 0) {
astSetFitsU ( fitschan, "JOS_MULT", "Nr of steps between bsw", 0 );
} else {
astSetFitsI ( fitschan, "JOS_MULT", josMult, "Nr of steps between bsw", 0 );
}
if (josMin == 0) {
astSetFitsU ( fitschan, "JOS_MIN", "Nr of steps between psw", 0 );
} else {
astSetFitsI ( fitschan, "JOS_MIN", josMin, "Nr of steps between psw", 0 );
}
astSetFitsU ( fitschan, "NCALSTEP", "Number of RTS steps for each CAL", 0 );
astSetFitsF ( fitschan, "NREFSTEP", nRefStep,
"Mean no. of RTS steps for each REF", 0 );
if( stBetRef_defined ) {
astSetFitsI ( fitschan, "STBETREF", stBetRef,
"Target number of RTS steps between REFs", 0 );
} else {
astSetFitsU ( fitschan, "STBETREF", "Target number of RTS steps between REFs", 0 );
}
astSetFitsU ( fitschan, "STBETCAL", "Target number of RTS steps between CALs",0 );
astSetFitsI ( fitschan, "STARTIDX", startIdx,
"Index into pattern at start of obs", 0 );
astSetFitsU ( fitschan, "FOCAXIS", "Focus Axis to move (X, Y, Z)", 0 );
astSetFitsU ( fitschan, "NFOCSTEP", "Number of focal position steps", 0 );
astSetFitsU ( fitschan, "FOCSTEP", "Distance between focal steps", 0 );
/* Miscellaneous parameters */
astSetFitsCM ( fitschan,
"---- Miscellaneous ----", 0 );
astSetFitsU ( fitschan, "OCSCFG", "OCS config filename", 0 );
astSetFitsL ( fitschan, "SIMULATE", 0,
"True if any data are simulated", 0 );
astSetFitsL ( fitschan, "SIM_CORR", 0,
"True if any CORRTASK is simualted", 0 );
astSetFitsL ( fitschan, "SIM_SMU", 0,
"True if SMU data is simulated", 0 );
astSetFitsL ( fitschan, "SIM_TCS", 0,
"True if TCS data is simulated", 0 );
astSetFitsL ( fitschan, "SIM_RTS", 0,
"True if RTS data is simulated", 0 );
astSetFitsL ( fitschan, "SIM_IF", 0,
"True if IF data is simulated", 0 );
astSetFitsS ( fitschan, "STATUS", "NORMAL",
"Status at end of observation", 0 );
/* Rover=specific parameters */
astSetFitsCM ( fitschan,
"---- ROVER Specific ----", 0 );
astSetFitsL ( fitschan, "POL_CONN", 0,
"True if ROVER is connected", 0 );
astSetFitsU ( fitschan, "POL_MODE", "Step-and-integrate (STEPINT) or spinning (SPIN)", 0 );
astSetFitsU ( fitschan, "ROTAFREQ", "[Hz] Spin frequency (if spinning)", 0 );
astSetFitsU ( fitschan, "POL_CRD", "Coordinate frame of polarimeter angles", 0 );
astSetFitsU ( fitschan, "POL_PA", "[[deg] Angle of pol fast axis", 0 );
astSetFitsS ( fitschan, "BUNIT", "K", "", 0 );
}
