casa / casa5 / 692953af62a

Commits

Kumar Golap authored 692953af62a 04 Apr 2018

Fixed an issue when there are more fields than what is used in the ms..and picking the right one for frequency frame conversion calculations

Modified


401 401	if(freqStart > freqEndObs) freqStart=freqEndObs;
402 402	if(freqEnd < freqEndObs) freqEnd=freqEndObs;
403 403	retval=True;
404 404	}
405 405	else{
406 406	MFrequency::Convert toframe(obsMFreqType,
407 407	MFrequency::Ref(freqframe, frame));
408 408	for (uInt j=0; j< nTimes; ++j){
409 409	if((useFieldsInMS \|\| anyEQ(fieldIds, fldId[elindx[uniqIndx[j]]])) && anyEQ(ddOfSpw, ddId[elindx[uniqIndx[j]]])){
410 410	timeCol.get(elindx[uniqIndx[j]], ep);
411 +	dir=fieldCol.phaseDirMeas(fldId[elindx[uniqIndx[j]]]);
411 412	frame.resetEpoch(ep);
413 +	frame.resetDirection(dir);
412 414	Double freqTmp=toframe(Quantity(freqStartObs, "Hz")).get("Hz").getValue();
413 415	if(freqStart > freqTmp) freqStart=freqTmp;
414 416	if(freqEnd < freqTmp) freqEnd=freqTmp;
415 417	freqTmp=toframe(Quantity(freqEndObs, "Hz")).get("Hz").getValue();
416 418	if(freqStart > freqTmp) freqStart=freqTmp;
417 419	if(freqEnd < freqTmp) freqEnd=freqTmp;
418 420	retval=True;
419 421	}
420 422	}
421 423	}