Changeset 2550 for branches/hpc33/src/MSFiller.cpp
- Timestamp:
- 05/29/12 11:24:15 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/hpc33/src/MSFiller.cpp
r2549 r2550 304 304 // set dummy epoch 305 305 mf.set( currentTime ) ; 306 307 // initialize dirType308 //dirType = MDirection::N_Types ;309 306 310 307 // … … 957 954 pointingDirection.xyPlane( idx-1 ), pointingDirection.xyPlane( idx ) ) ; 958 955 } 959 //mf.resetEpoch( currentTime ) ;960 956 mf.set( currentTime ) ; 961 957 Quantum< Vector<Double> > tmp( d.column( 0 ), Unit( "rad" ) ) ; 962 958 if ( dirType != MDirection::J2000 ) { 963 //MDirection::Convert toj2000( dirType, MDirection::Ref( MDirection::J2000, mf ) ) ;964 959 dir = toj2000( tmp ).getAngle( "rad" ).getValue() ; 965 960 } … … 968 963 } 969 964 if ( dirType != MDirection::AZELGEO ) { 970 //MDirection::Convert toazel( dirType, MDirection::Ref( MDirection::AZELGEO, mf ) ) ;971 //MDirection::Convert toazel( dirType, MDirection::Ref( MDirection::AZEL, mf ) ) ;972 965 azel = toazel( tmp ).getAngle( "rad" ).getValue() ; 973 966 } … … 981 974 { 982 975 dir = sourceDir.getAngle( "rad" ).getValue() ; 983 mf.resetEpoch( currentTime ) ; 984 //MDirection::Convert toazel( dirType, MDirection::Ref( MDirection::AZELGEO, mf ) ) ; 976 mf.set( currentTime ) ; 985 977 azel = toazel( Quantum< Vector<Double> >( dir, Unit("rad") ) ).getAngle( "rad" ).getValue() ; 986 978 if ( dirType != MDirection::J2000 ) { 987 //MDirection::Convert toj2000( dirType, MDirection::Ref( MDirection::J2000, mf ) ) ;988 979 dir = toj2000( Quantum< Vector<Double> >( dir, Unit("rad") ) ).getAngle( "rad" ).getValue() ; 989 980 }
Note: See TracChangeset
for help on using the changeset viewer.