Commits

Takeshi Nakazato authored d9f1960a7f4 Merge
Merge pull request #456 in CASA/casa from bugfix/CAS-11263 to master

* commit '692953af62abec46c9394f90609624709a9c366f': 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

code/msvis/MSVis/MSUtil.cc

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 }

Everything looks good. We'll let you know here if there's anything you should know about.

Add shortcut