sourcePos

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  
  function [sourceRA sourceDEC] = sourcePos(sourceName, mjd)

   this function is just a simple correspondance between a source name and
   anumber, for posterity, just to make saving things as a text file
   easier.

  sjcm
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

0001 function [sourceRA sourceDEC] = sourcePos(sourceName)
0002 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0003 %
0004 %  function [sourceRA sourceDEC] = sourcePos(sourceName, mjd)
0005 %
0006 %   this function is just a simple correspondance between a source name and
0007 %   anumber, for posterity, just to make saving things as a text file
0008 %   easier.
0009 %
0010 %  sjcm
0011 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0012 
0013 if(nargin<2)
0014   year = 2012;
0015 else
0016   year = mjd2fracyear(mjd);
0017 end
0018 [casS casSerr cygS cygSerr tauS tauSerr] = CalibFluxDen(5.0, year);
0019 
0020 
0021 switch(lower(sourceName))
0022   case 'cyga'
0023     sourceNum = 1;
0024     sourceFlux = cygS;
0025     fluxErr    = cygSerr; 
0026     sourceRA = 299.8667 ;
0027     sourceDEC = 40.6667 ;    
0028 %    sourceFlux = 365;
0029 %    fluxErr    = 0.02; %
0030     
0031   case 'cygaoff'
0032     sourceNum = 1;
0033     sourceFlux = cygS;
0034     fluxErr    = cygSerr; %
0035     
0036   case 'taua'
0037     sourceNum = 2;
0038     sourceFlux = tauS;
0039     fluxErr    = tauSerr;
0040     sourceRA = 83.6333 ;
0041     sourceDEC= 22.0144 ;
0042 %    sourceFlux = 615;
0043 %    fluxErr    = 0.02; %
0044     
0045   case 'tauaoff'
0046     sourceNum = 2;
0047     sourceFlux = tauS;
0048     fluxErr    = tauSerr;
0049     
0050   case 'j0437+296' %3c123
0051     sourceNum = 3;
0052     sourceFlux = 16.56;
0053     fluxErr    = 0.03; %
0054     sourceRA =  0;
0055     sourceDEC = 0;
0056   case 'j0915-119' %3C218
0057     sourceNum = 4;
0058     sourceFlux = 13.65;
0059     fluxErr    = 0.03; %
0060     sourceRA =  0;
0061     sourceDEC = 0;
0062   case 'j1230+123' %3c274
0063     sourceNum = 5;
0064     sourceFlux = 68;
0065     fluxErr    = 2; %
0066     sourceRA =  0;
0067     sourceDEC = 0;
0068   case '3c274' %3c274; virgoA
0069     sourceNum = 5;
0070     sourceFlux = 68;
0071     fluxErr    = 2; %
0072     sourceRA =  187.7042;
0073     sourceDEC =  12.3911 ;
0074      case 'virgoa' %3c274; virgoA
0075     sourceNum = 5;
0076     sourceFlux = 68;
0077     fluxErr    = 2; %
0078     sourceRA =  187.7042;
0079     sourceDEC =  12.3911 ;
0080   case 'j1331+305' %3C286'
0081     sourceNum = 6;
0082     sourceFlux = 7.6;
0083     fluxErr    = 0.03; %
0084     sourceRA =  0;
0085     sourceDEC = 0;
0086   case 'j1720-009' %3C353
0087     sourceNum = 7;
0088     sourceFlux = 23.49;
0089     fluxErr    = 0.03; %
0090     sourceRA =  0;
0091     sourceDEC = 0;
0092   case '3c353' %3C353
0093     sourceNum = 7;
0094     sourceFlux = 23.49;
0095     fluxErr    = 0.03; %
0096     sourceRA =  0;
0097     sourceDEC = 0;
0098   case 'casa'
0099     sourceNum = 8;
0100     sourceFlux = casS;
0101     fluxErr    = casSerr;
0102     sourceRA =  350.8667;
0103     sourceDEC = 58.8117;
0104 %    sourceFlux = 670;
0105 %    fluxErr    = 0.02; %
0106     
0107   case 'casaoff'
0108     sourceNum = 8;
0109     sourceFlux = casS;
0110     fluxErr    = casSerr;
0111       
0112   case 'j0225+621'
0113     sourceNum = 9;
0114     sourceFlux = 135.5;  % my best estimate
0115     fluxErr    = 0.05; %
0116     sourceRA = 0 ;
0117     sourceDEC =0;
0118     
0119     
0120   case 'j0927+390'
0121     sourceNum = 10;
0122     sourceFlux = 8.5;
0123     fluxErr    = 0.50; %
0124     sourceRA =  0;
0125     sourceDEC = 0;
0126     case 'm42' %3C145
0127     sourceNum = 11;
0128     sourceFlux = 0;  %
0129     fluxErr    = 0; %
0130     sourceRA = 83.8208 ;
0131     sourceDEC = -5.3911 ;
0132   otherwise
0133     sourceNum = 99;
0134     sourceFlux = 0;
0135     fluxErr    = 0;
0136 end
0137 
0138 return;