casa / casa6 / cc478d62719

Commits

casa5/gcwrap/python/scripts/simutil.py

Modified

-986
+                eta_b geometrical blockage efficiency,
-987
+                eta_t taper efficiency,
-988
+                eta_q correlator efficiency including quantization,
-989
+                t_rx  reciever temperature.
-990
+        antenna efficiency = eta_p * eta_s * eta_b * eta_t
-991
+        Notes: VLA correlator efficiency includes waveguide loss
-992
+               EVLA correlator efficiency is probably optimistic at 0.88
-993
+        """
-994
-995
+        if telescope==None: telescope=self.telescopename
-+
+        # Noisetemp only knows about 6 observatories,
-+
+        # none of which have measure.observatory dicts that
-+
+        # contain lowercase characters so str.upper is harmless here
-999
+        telescope=str.upper(telescope)
-1000
-1001
+        obs =['ALMASD','ALMA','ACA','EVLA','VLA','SMA']
-1002
+        d   =[ 12.    ,12.   ,7.   ,25.   ,25.  , 6. ]
-1003
+        ds  =[ 0.75   ,0.75  ,0.75 ,0.364 ,0.364,0.35] # subreflector size for ACA?
-1004
+        eps =[ 25.    ,25.   ,20.  ,300   ,300  ,15. ] # antenna surface accuracy
-1005
-1006
+        cq  =[ 0.845, 0.845,  0.88,  0.79, 0.86, 0.88] # correlator eff
-1007
+        # SMA is from Wright http://astro.berkeley.edu/~wright/obsrms.py
-1008
+        # ALMA includes quantization eff of 0.96
-1172
+        if antennalist==None:
-1173
+            if telescope==None:
-1174
+                self.msg("Telescope name has not been set.",priority="error")
-1175
+                return False
-1176
+            if diam==None:
-1177
+                self.msg("Antenna diameter has not been set.",priority="error")
-1178
+                return False
-1179
+            if nant==None:
-1180
+                self.msg("Number of antennas has not been set.",priority="error")
-1181
+                return False
--
--
+            found=False
--
+            t=telescope.upper()
--
+            for l in pl.arange(len(t)-1)+2:
--
+                if t[0:l] in me.obslist(): found=True
--
+            if found:
--
+                posobs=me.measure(me.observatory(telescope),'WGS84')
-+
-+
+            known=False
-+
+            # check if telescope is known to measures tool
-+
+            # ensure case insensitivity - CAS-12753
-+
+            obslist_lower = [obs.lower() for obs in me.obslist()]
-+
+            if self.telescope.lower() in obslist_lower:
-+
+                t = self.telescope
-+
+                known = True
-+
-+
+            if known == True:
-+
+                posobs = me.measure(me.observatory(t), 'WGS84')
-1193
+            else:
--
+                self.msg("Unknown telescope and no antenna list.",priority="error")
-+
+                self.msg("Unknown telescope and no antenna list.",
-+
+                         priority="error")
-1196
+                return False
-1197
-1198
+            obslat=qa.convert(posobs['m1'],'deg')['value']
-1199
+            obslon=qa.convert(posobs['m0'],'deg')['value']
-1200
+            obsalt=qa.convert(posobs['m2'],'m')['value']
-1201
+            stnx,stny,stnz = self.locxyz2itrf(obslat,obslon,obsalt,0,0,0)
-1202
+            antnames="A00"
-1203
+            # use a single dish of diameter diam - must be set if antennalist
-1204
+            # is not set.
-1205
+            stnd=[diam]
-1664
+            self.coordsys=params["coordsys"]
-1665
-1666
+        if "observatory" in params:
-1667
+            self.telescopename=params["observatory"]
-1668
+        else:
-1669
+            self.telescopename="SIMULATED"
-1670
+        if self.verbose:
-1671
+            self.msg("Using observatory= %s" % self.telescopename,
-1672
+                     origin="readantenna")
-1673
--
+        # me.observatory has partial matching implemented
--
+        found=False
--
+        t=self.telescopename.upper()
--
+        for l in pl.arange(len(t)-1)+2:
--
+            if t[0:l] in me.obslist(): found=True
--
+        if found:
--
+            posobs=me.measure(me.observatory(self.telescopename),'WGS84')
-+
+        known = False
-+
+        # case insensitive check if telescopename is known
-+
+        obslist_lower = [obs.lower() for obs in me.obslist()]
-+
+        if self.telescopename.lower() in obslist_lower:
-+
+            t = self.telescopename
-+
+            known = True
-+
+            posobs=me.measure(me.observatory(t),'WGS84')
-1681
-1682
+        if "COFA" in params:
-+
+            if known:
-+
+                self.msg("antenna config file specifies COFA for a known observatory "+
-+
+                         self.telescopename+", overriding with specified COFA.",priority="warn")
-1686
+            obs_latlon=params["COFA"].split(",")
-1687
+            cofa_lon=float(obs_latlon[0])
-1688
+            cofa_lat=float(obs_latlon[1])
-1689
+            cofa_alt=0.
-1690
+            posobs=me.position("WGS84",qa.quantity(cofa_lon,"deg"),qa.quantity(cofa_lat,"deg"),qa.quantity(cofa_alt,"m"))
--
+            if found:
--
+                self.msg("antenna config file specifies COFA but a known observatory "+self.telescopename+", so ignoring specified COFA.",priority="warn")
--
+        elif not found:
-+
+        elif not known:
-1692
+            if params["coordsys"].upper()[0:3]=="LOC":
-1693
+                self.msg("To use local coordinates in the antenna position file, you must either use a known observatory name, or provide the COFA explicitly",priority="error")
-1694
+                return -1
-1695
+            else:
-1696
+                # we have absolute coords, so can create the posobs from their
-1697
+                # average at the end
-1698
+                posobs={}
-1699
-1700
-1701

casa6

Commits

Jira Issues

Add shortcut