Line data Source code
1 : //# PBMath1DEVLA.cc: Implementation for PBMath1DEVLA
2 : //# Copyright (C) 2016
3 : //# Associated Universities, Inc. Washington DC, USA.
4 : //#
5 : //# This library is free software; you can redistribute it and/or modify it
6 : //# under the terms of the GNU General Public License as published by
7 : //# the Free Software Foundation; either version 2 of the License, or (at your
8 : //# option) any later version.
9 : //#
10 : //# This library is distributed in the hope that it will be useful, but WITHOUT
11 : //# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
12 : //# FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public
13 : //# License for more details.
14 : //#
15 : //# You should have received a copy of the GNU General Public License
16 : //# along with this library; if not, write to the Free Software Foundation,
17 : //# Inc., 675 Massachusetts Ave, Cambridge, MA 02139, USA.
18 : //#
19 : //# Correspondence concerning AIPS++ should be adressed as follows:
20 : //# Internet email: aips2-request@nrao.edu.
21 : //# Postal address: AIPS++ Project Office
22 : //# National Radio Astronomy Observatory
23 : //# 520 Edgemont Road
24 : //# Charlottesville, VA 22903-2475 USA
25 : //#
26 : //#
27 : //# $Id$
28 :
29 : #include <synthesis/TransformMachines/PBMath1DEVLA.h>
30 : #include <synthesis/TransformMachines2/Utils.h>
31 : #include <casacore/measures/Measures.h>
32 : #include <casacore/measures/Measures/MeasConvert.h>
33 : #include <casacore/measures/Measures/MDirection.h>
34 :
35 : using namespace casacore;
36 :
37 : namespace casa { //# NAMESPACE CASA - BEGIN
38 :
39 0 : PBMath1DEVLA::PBMath1DEVLA(): pbMathPoly_p(0) {
40 0 : wideFit_p = true;
41 0 : init();
42 0 : maxRad_p=casacore::Quantity(1.0, "deg");
43 0 : }
44 :
45 271 : PBMath1DEVLA::PBMath1DEVLA(Quantity maxRad, Bool useSymmetricBeam, Double freqToUse): pbMathPoly_p(0), maxRad_p(maxRad), useSymmetric_p(useSymmetricBeam){
46 271 : wideFit_p=True;
47 271 : init();
48 271 : nearestVPArray(freqToUse, False);
49 271 : }
50 :
51 200 : PBMathInterface::PBClass PBMath1DEVLA::whichPBClass(){
52 200 : if(pbMathPoly_p)
53 200 : return pbMathPoly_p->whichPBClass();
54 :
55 0 : return PBMathInterface::COMMONPB;
56 : }
57 :
58 0 : void PBMath1DEVLA::fillPBArray(){
59 : ///This is a dummy to make compiler happy
60 : ///as in this class we recalculate vp_p when needed
61 0 : nearestVPArray(1.425e9, False);
62 :
63 0 : }
64 :
65 1553 : void PBMath1DEVLA::nearestVPArray(Double freq, bool printINFO){
66 4659 : LogIO os(LogOrigin("PBMATH1DEVLA", "nearestVPArray"));
67 3106 : String band="";
68 :
69 1553 : if(refim::SynthesisUtils::parseBandName(bandOrFeedName_p).nelements() > 0){
70 1553 : band=refim::SynthesisUtils::parseBandName(bandOrFeedName_p)[0];
71 1553 : band=band.after("EVLA_");
72 : }
73 1553 : if(band==""){
74 1553 : band=feed(freq);
75 : }
76 1553 : Double origFreq=freq;
77 1553 : limitFreqForBand(band, freq);
78 :
79 1553 : auto confiter=feedConf_p.find(band);
80 1553 : Float mag = 1.21;
81 1553 : if(confiter != feedConf_p.end()){
82 3106 : squint_p=BeamSquint(MDirection(Quantity((mag*cos(confiter->second)), "'"),
83 3106 : Quantity((mag*sin(confiter->second)), "'"),
84 3106 : MDirection::Ref(MDirection::AZEL)),
85 4659 : Quantity(1.0, "GHz"));
86 : }
87 : else{
88 : //No squint
89 0 : squint_p=BeamSquint(MDirection(Quantity(0.0, "'"),
90 0 : Quantity(0.0, "'"),
91 0 : MDirection::Ref(MDirection::AZEL)),
92 0 : Quantity(1.0, "GHz"));
93 : }
94 3106 : Vector<Double> coeff;
95 1553 : Double freqMHz=freq*1e-6;
96 1553 : Double freqUsed=freq*1e-6;
97 1553 : std::map<Double, std::vector<Double> >::iterator low, prev;
98 1553 : low = coeffmap_p.lower_bound(freqMHz);
99 1553 : if (low == coeffmap_p.end()) {
100 0 : --low;
101 0 : coeff=Vector<Double>(low->second);
102 0 : freqUsed=low->first;
103 : }
104 1553 : else if (low == coeffmap_p.begin()) {
105 0 : coeff=Vector<Double>(low->second);
106 0 : freqUsed=low->first;
107 : }
108 : else{
109 1553 : prev = low;
110 1553 : --prev;
111 : //cerr << "freqMHz " << freqMHz << " prev " << prev->first << " low " <<low->first << endl;
112 1553 : if (fabs(freqMHz - prev->first) < fabs(low->first - freqMHz)){
113 686 : coeff=Vector<Double>(prev->second);
114 686 : freqUsed=prev->first;
115 : }
116 : else{
117 867 : coeff=Vector<Double>(low->second);
118 867 : freqUsed=low->first;
119 : }
120 :
121 : }
122 1553 : if(printINFO)
123 1282 : os << LogIO::NORMAL1 << "Using EVLA "+band+" band beam model of frequency " << freqUsed << " MHz scaled to original frequency "<< origFreq*1e-6 << LogIO::POST;
124 1553 : pbMathPoly_p= new PBMath1DPoly(coeff, maxRad_p, Quantity(1.0, "GHz"), False, squint_p, useSymmetric_p);
125 1553 : (this->vp_p).resize();
126 1553 : (this->vp_p)=pbMathPoly_p->vp_p;
127 : //cerr << "vp_p " << vp_p << endl;
128 1553 : maximumRadius_p=pbMathPoly_p->maximumRadius_p;
129 1553 : refFreq_p=pbMathPoly_p->refFreq_p;
130 1553 : fScale_p=pbMathPoly_p->fScale_p;
131 1553 : inverseIncrementRadius_p=pbMathPoly_p->inverseIncrementRadius_p;
132 1553 : scale_p=pbMathPoly_p->scale_p;
133 1553 : composite_p=pbMathPoly_p->composite_p;
134 : //cerr << "maximumRad " << maximumRadius_p << " inverse " << inverseIncrementRadius_p << endl;
135 1553 : }
136 1553 : String PBMath1DEVLA::feed(const Double freq, const Double freqCenter){
137 1553 : if((freqCenter > 224e6 && freqCenter<480e6) || (freq > 224e6 && freq<480e6))
138 0 : return "P";
139 1553 : if((freqCenter > 900e6 && freqCenter<2003.0e6) || (freq > 900e6 && freq<2003e6))
140 1461 : return "L";
141 92 : if((freqCenter > 1990e6 && freqCenter<4001.0e6) || (freq > 1990e6 && freq<4001e6))
142 2 : return "S";
143 90 : if((freqCenter > 3990e6 && freqCenter<8001.0e6) || (freq > 3990e6 && freq<8001e6))
144 90 : return "C";
145 0 : if((freqCenter > 7990e6 && freqCenter<12001.0e6) || (freq > 7990e6 && freq<12001e6))
146 0 : return "X";
147 0 : if((freqCenter > 12000e6 && freqCenter<18000.0e6) || (freq > 12000e6 && freq<18000e6))
148 0 : return "U";
149 0 : if((freqCenter > 19000e6 && freqCenter<26000.0e6) || (freq > 19000e6 && freq<26000e6))
150 0 : return "K";
151 0 : if((freqCenter > 28000e6 && freqCenter<38000.0e6) || (freq > 28000e6 && freq<38000e6))
152 0 : return "A";
153 0 : if((freqCenter > 41000e6 && freqCenter<50000.0e6) || (freq > 41000e6 && freq<50000e6))
154 0 : return "Q";
155 :
156 0 : return "";
157 : }
158 1553 : void PBMath1DEVLA::limitFreqForBand(const casacore::String& band, casacore::Double& freq){
159 :
160 1553 : if(band=="P"){
161 0 : if(freq <= 224e6) freq=232e6;
162 0 : if(freq >=480e6) freq=470e6;
163 : }
164 1553 : else if(band=="L"){
165 1461 : if(freq <= 900e6) freq=1040e6;
166 1461 : if(freq >=2000e6) freq=2000e6;
167 : }
168 92 : else if(band=="S"){
169 2 : if(freq < 2052e6) freq=2052e6;
170 2 : if(freq >=3948e6) freq=3948e6;
171 : }
172 90 : else if(band=="C"){
173 90 : if(freq < 4052e6) freq=4052e6;
174 90 : if(freq >=7948e6) freq=7948e6;
175 : }
176 0 : else if(band=="X"){
177 0 : if(freq < 8052e6) freq=8052e6;
178 0 : if(freq >=11948e6) freq=11948e6;
179 : }
180 0 : else if(band=="U"){
181 0 : if(freq < 12052e6) freq=12052e6;
182 0 : if(freq >=17948e6) freq=17948e6;
183 : }
184 0 : else if(band=="K"){
185 0 : if(freq < 19052e6) freq=19052e6;
186 0 : if(freq >=25948e6) freq=25948e6;
187 : }
188 0 : else if(band=="A"){
189 0 : if(freq < 28052e6) freq=28052e6;
190 0 : if(freq >=38048e6) freq=38048e6;
191 : }
192 0 : else if(band=="Q"){
193 0 : if(freq < 41052e6) freq=41052e6;
194 0 : if(freq >=43948e6) freq=43948e6;
195 : }
196 1553 : }
197 :
198 271 : void PBMath1DEVLA::init(){
199 542 : wFreqs_p=Vector<Double>(
200 542 : std::vector<Double>({232., 246., 281., 296., 312., 328., 344., 357., 382., 392., 403., 421., 458., 470., 1040, 1104, 1168, 1232, 1296, 1360, 1424, 1488, 1552, 1680, 1744, 1808, 1872, 1936, 2000})
201 271 : );
202 271 : wFreqs_p *=1e6;
203 : //Float mag = 1.21; // half-squint magnitude in arcmin at 1 GHz)
204 271 : feedConf_p["L"] = (-185.9)*C::pi/180.0; // squint orientation, rads, North of +AZ axis
205 271 : feedConf_p["S"] = (-11.61)*C::pi/180.0;
206 271 : feedConf_p["C"]=(-104.8)*C::pi/180.0;
207 271 : feedConf_p["X"]=(-113.7)*C::pi/180.0;
208 271 : feedConf_p["U"]=(42.4)*C::pi/180.0;
209 271 : feedConf_p["K"]=(64.4)*C::pi/180.0;
210 271 : feedConf_p["A"]=(106.9)*C::pi/180.0;
211 271 : feedConf_p["Q"]=(85.5)*C::pi/180.0;
212 :
213 542 : squint_p=BeamSquint( MDirection( Quantity(1.25, "arcmin"), Quantity(90.0, "deg"),
214 813 : MDirection::Ref(MDirection::AZEL)), Quantity(1.0, "GHz"));
215 : ////P
216 271 : coeffmap_p[232]={1.0, -1.137e-3, 5.19e-7, -1.04e-10, 0.71e-14};
217 271 : coeffmap_p[246]={1.0, -1.130e-3, 5.04e-7, -1.02e-10, 0.77e-14};
218 271 : coeffmap_p[281]={1.0, -1.106e-3, 5.11e-7, -1.10e-10, 0.91e-14};
219 271 : coeffmap_p[296]={1.0, -1.125e-3, 5.27e-7, -1.14e-10, 0.96e-14};
220 271 : coeffmap_p[312]={1.0, -1.030e-3, 4.44e-7, -0.89e-10, 0.68e-14};
221 271 : coeffmap_p[328]={1.0, -0.980e-3, 4.25e-7, -0.87e-10, 0.69e-14};
222 271 : coeffmap_p[344]={1.0, -0.974e-3, 4.09e-7, -0.76e-10, 0.53e-14};
223 271 : coeffmap_p[357]={1.0, -0.996e-3, 4.23e-7, -0.79e-10, 0.51e-14};
224 271 : coeffmap_p[382]={1.0, -1.002e-3, 4.39e-7, -0.88e-10, 0.64e-14};
225 271 : coeffmap_p[392]={1.0, -1.067e-3, 5.13e-7, -1.12e-10, 0.90e-14};
226 271 : coeffmap_p[403]={1.0, -1.057e-3, 4.90e-7, -1.06e-10, 0.87e-14};
227 271 : coeffmap_p[421]={1.0, -1.154e-3, 5.85e-7, -1.33e-10, 1.08e-14};
228 271 : coeffmap_p[458]={1.0, -0.993e-3, 4.67e-7, -1.04e-10, 0.88e-14};
229 271 : coeffmap_p[470]={1.0, -1.010e-3, 4.85e-7, -1.07e-10, 0.86e-14};
230 : /////////L
231 271 : coeffmap_p[1040]={ 1.000, -1.529e-3, 8.69e-7, -1.88e-10};
232 271 : coeffmap_p[1104]={ 1.000, -1.486e-3, 8.15e-7, -1.68e-10};
233 271 : coeffmap_p[1168]={ 1.000, -1.439e-3, 7.53e-7, -1.45e-10};
234 271 : coeffmap_p[1232]={ 1.000, -1.450e-3, 7.87e-7, -1.63e-10};
235 271 : coeffmap_p[1296]={ 1.000, -1.428e-3, 7.62e-7, -1.54e-10};
236 271 : coeffmap_p[1360]={ 1.000, -1.449e-3, 8.02e-7, -1.74e-10};
237 271 : coeffmap_p[1424]={ 1.000, -1.462e-3, 8.23e-7, -1.83e-10};
238 271 : coeffmap_p[1488]={ 1.000, -1.455e-3, 7.92e-7, -1.63e-10};
239 271 : coeffmap_p[1552]={ 1.000, -1.435e-3, 7.54e-7, -1.49e-10};
240 271 : coeffmap_p[1680]={ 1.000, -1.443e-3, 7.74e-7, -1.57e-10};
241 271 : coeffmap_p[1744]={ 1.000, -1.462e-3, 8.02e-7, -1.69e-10};
242 271 : coeffmap_p[1808]={ 1.000, -1.488e-3, 8.38e-7, -1.83e-10};
243 271 : coeffmap_p[1872]={ 1.000, -1.486e-3, 8.26e-7, -1.75e-10};
244 271 : coeffmap_p[1936]={ 1.000, -1.459e-3, 7.93e-7, -1.62e-10};
245 271 : coeffmap_p[2000]={ 1.000, -1.508e-3, 8.31e-7, -1.68e-10};
246 : ////////S
247 271 : coeffmap_p[2052]={1.000,-1.429e-3,7.52e-7,-1.47e-10};
248 271 : coeffmap_p[2180]={1.000,-1.389e-3,7.06e-7,-1.33e-10};
249 271 : coeffmap_p[2436]={1.000,-1.377e-3,6.90e-7,-1.27e-10};
250 271 : coeffmap_p[2564]={1.000,-1.381e-3,6.92e-7,-1.26e-10};
251 271 : coeffmap_p[2692]={1.000,-1.402e-3,7.23e-7,-1.40e-10};
252 271 : coeffmap_p[2820]={1.000,-1.433e-3,7.62e-7,-1.54e-10};
253 271 : coeffmap_p[2948]={1.000,-1.433e-3,7.46e-7,-1.42e-10};
254 271 : coeffmap_p[3052]={1.000,-1.467e-3,8.05e-7,-1.70e-10};
255 271 : coeffmap_p[3180]={1.000,-1.497e-3,8.38e-7,-1.80e-10};
256 271 : coeffmap_p[3308]={1.000,-1.504e-3,8.37e-7,-1.77e-10};
257 271 : coeffmap_p[3436]={1.000,-1.521e-3,8.63e-7,-1.88e-10};
258 271 : coeffmap_p[3564]={1.000,-1.505e-3,8.37e-7,-1.75e-10};
259 271 : coeffmap_p[3692]={1.000,-1.521e-3,8.51e-7,-1.79e-10};
260 271 : coeffmap_p[3820]={1.000,-1.534e-3,8.57e-7,-1.77e-10};
261 271 : coeffmap_p[3948]={1.000,-1.516e-3,8.30e-7,-1.66e-10};
262 : ///C
263 271 : coeffmap_p[4052]={1.000,-1.406e-3,7.41e-7,-1.48e-10};
264 271 : coeffmap_p[4180]={1.000,-1.385e-3,7.09e-7,-1.36e-10};
265 271 : coeffmap_p[4308]={1.000,-1.380e-3,7.08e-7,-1.37e-10};
266 271 : coeffmap_p[4436]={1.000,-1.362e-3,6.95e-7,-1.35e-10};
267 271 : coeffmap_p[4564]={1.000,-1.365e-3,6.92e-7,-1.31e-10};
268 271 : coeffmap_p[4692]={1.000,-1.339e-3,6.56e-7,-1.17e-10};
269 271 : coeffmap_p[4820]={1.000,-1.371e-3,7.06e-7,-1.40e-10};
270 271 : coeffmap_p[4948]={1.000,-1.358e-3,6.91e-7,-1.34e-10};
271 271 : coeffmap_p[5052]={1.000,-1.360e-3,6.91e-7,-1.33e-10};
272 271 : coeffmap_p[5180]={1.000,-1.353e-3,6.74e-7,-1.25e-10};
273 271 : coeffmap_p[5308]={1.000,-1.359e-3,6.82e-7,-1.27e-10};
274 271 : coeffmap_p[5436]={1.000,-1.380e-3,7.05e-7,-1.37e-10};
275 271 : coeffmap_p[5564]={1.000,-1.376e-3,6.99e-7,-1.31e-10};
276 271 : coeffmap_p[5692]={1.000,-1.405e-3,7.39e-7,-1.47e-10};
277 271 : coeffmap_p[5820]={1.000,-1.394e-3,7.29e-7,-1.45e-10};
278 271 : coeffmap_p[5948]={1.000,-1.428e-3,7.57e-7,-1.57e-10};
279 271 : coeffmap_p[6052]={1.000,-1.445e-3,7.68e-7,-1.50e-10};
280 271 : coeffmap_p[6148]={1.000,-1.422e-3,7.38e-7,-1.38e-10};
281 271 : coeffmap_p[6308]={1.000,-1.463e-3,7.94e-7,-1.62e-10};
282 271 : coeffmap_p[6436]={1.000,-1.478e-3,8.22e-7,-1.74e-10};
283 271 : coeffmap_p[6564]={1.000,-1.473e-3,8.00e-7,-1.62e-10};
284 271 : coeffmap_p[6692]={1.000,-1.455e-3,7.76e-7,-1.53e-10};
285 271 : coeffmap_p[6820]={1.000,-1.487e-3,8.22e-7,-1.72e-10};
286 271 : coeffmap_p[6948]={1.000,-1.472e-3,8.05e-7,-1.67e-10};
287 271 : coeffmap_p[7052]={1.000,-1.470e-3,8.01e-7,-1.64e-10};
288 271 : coeffmap_p[7180]={1.000,-1.503e-3,8.50e-7,-1.84e-10};
289 271 : coeffmap_p[7308]={1.000,-1.482e-3,8.19e-7,-1.72e-10};
290 271 : coeffmap_p[7436]={1.000,-1.498e-3,8.22e-7,-1.66e-10};
291 271 : coeffmap_p[7564]={1.000,-1.490e-3,8.18e-7,-1.66e-10};
292 271 : coeffmap_p[7692]={1.000,-1.481e-3,7.98e-7,-1.56e-10};
293 271 : coeffmap_p[7820]={1.000,-1.474e-3,7.94e-7,-1.57e-10};
294 271 : coeffmap_p[7948]={1.000,-1.448e-3,7.69e-7,-1.51e-10};
295 : //////X
296 271 : coeffmap_p[8052]={1.000,-1.403e-3,7.21e-7,-1.37e-10};
297 271 : coeffmap_p[8180]={1.000,-1.398e-3,7.10e-7,-1.32e-10};
298 271 : coeffmap_p[8308]={1.000,-1.402e-3,7.16e-7,-1.35e-10};
299 271 : coeffmap_p[8436]={1.000,-1.400e-3,7.12e-7,-1.32e-10};
300 271 : coeffmap_p[8564]={1.000,-1.391e-3,6.95e-7,-1.25e-10};
301 271 : coeffmap_p[8692]={1.000,-1.409e-3,7.34e-7,-1.49e-10};
302 271 : coeffmap_p[8820]={1.000,-1.410e-3,7.36e-7,-1.45e-10};
303 271 : coeffmap_p[8948]={1.000,-1.410e-3,7.34e-7,-1.43e-10};
304 271 : coeffmap_p[9052]={1.000,-1.403e-3,7.20e-7,-1.36e-10};
305 271 : coeffmap_p[9180]={1.000,-1.396e-3,7.09e-7,-1.31e-10};
306 271 : coeffmap_p[9308]={1.000,-1.432e-3,7.68e-7,-1.55e-10};
307 271 : coeffmap_p[9436]={1.000,-1.414e-3,7.43e-7,-1.47e-10};
308 271 : coeffmap_p[9564]={1.000,-1.416e-3,7.45e-7,-1.47e-10};
309 271 : coeffmap_p[9692]={1.000,-1.406e-3,7.26e-7,-1.39e-10};
310 271 : coeffmap_p[9820]={1.000,-1.412e-3,7.36e-7,-1.43e-10};
311 271 : coeffmap_p[9948]={1.000,-1.409e-3,7.29e-7,-1.39e-10};
312 271 : coeffmap_p[10052]={1.000,-1.421e-3,7.46e-7,-1.45e-10};
313 271 : coeffmap_p[10180]={1.000,-1.409e-3,7.25e-7,-1.36e-10};
314 271 : coeffmap_p[10308]={1.000,-1.402e-3,7.13e-7,-1.31e-10};
315 271 : coeffmap_p[10436]={1.000,-1.399e-3,7.09e-7,-1.29e-10};
316 271 : coeffmap_p[10564]={1.000,-1.413e-3,7.37e-7,-1.43e-10};
317 271 : coeffmap_p[10692]={1.000,-1.412e-3,7.34e-7,-1.41e-10};
318 271 : coeffmap_p[10820]={1.000,-1.401e-3,7.12e-7,-1.31e-10};
319 271 : coeffmap_p[10948]={1.000,-1.401e-3,7.12e-7,-1.31e-10};
320 271 : coeffmap_p[10052]={1.000,-1.401e-3,7.12e-7,-1.31e-10};
321 271 : coeffmap_p[11180]={1.000,-1.394e-3,6.99e-7,-1.24e-10};
322 271 : coeffmap_p[11308]={1.000,-1.394e-3,7.01e-7,-1.26e-10};
323 271 : coeffmap_p[11436]={1.000,-1.391e-3,6.94e-7,-1.22e-10};
324 271 : coeffmap_p[11564]={1.000,-1.389e-3,6.92e-7,-1.22e-10};
325 271 : coeffmap_p[11692]={1.000,-1.386e-3,6.80e-7,-1.15e-10};
326 271 : coeffmap_p[11820]={1.000,-1.391e-3,6.88e-7,-1.19e-10};
327 271 : coeffmap_p[11948]={1.000,-1.399e-3,6.97e-7,-1.22e-10};
328 : ///U
329 271 : coeffmap_p[12052]={1.000,-1.399e-3,7.17e-7,-1.34e-10};
330 271 : coeffmap_p[12180]={1.000,-1.392e-3,7.07e-7,-1.31e-10};
331 271 : coeffmap_p[12308]={1.000,-1.393e-3,7.19e-7,-1.38e-10};
332 271 : coeffmap_p[12436]={1.000,-1.393e-3,7.20e-7,-1.40e-10};
333 271 : coeffmap_p[12564]={1.000,-1.395e-3,7.19e-7,-1.38e-10};
334 271 : coeffmap_p[12692]={1.000,-1.397e-3,7.20e-7,-1.37e-10};
335 271 : coeffmap_p[12820]={1.000,-1.388e-3,7.06e-7,-1.32e-10};
336 271 : coeffmap_p[12948]={1.000,-1.397e-3,7.18e-7,-1.36e-10};
337 271 : coeffmap_p[13052]={1.000,-1.400e-3,7.27e-7,-1.40e-10};
338 271 : coeffmap_p[13180]={1.000,-1.406e-3,7.44e-7,-1.50e-10};
339 271 : coeffmap_p[13308]={1.000,-1.403e-3,7.37e-7,-1.47e-10};
340 271 : coeffmap_p[13436]={1.000,-1.392e-3,7.08e-7,-1.31e-10};
341 271 : coeffmap_p[13564]={1.000,-1.384e-3,6.94e-7,-1.24e-10};
342 271 : coeffmap_p[13692]={1.000,-1.382e-3,6.95e-7,-1.25e-10};
343 271 : coeffmap_p[13820]={1.000,-1.376e-3,6.88e-7,-1.24e-10};
344 271 : coeffmap_p[13948]={1.000,-1.384e-3,6.98e-7,-1.28e-10};
345 271 : coeffmap_p[14052]={1.000,-1.400e-3,7.36e-7,-1.48e-10};
346 271 : coeffmap_p[14180]={1.000,-1.397e-3,7.29e-7,-1.45e-10};
347 271 : coeffmap_p[14308]={1.000,-1.399e-3,7.32e-7,-1.45e-10};
348 271 : coeffmap_p[14436]={1.000,-1.396e-3,7.25e-7,-1.42e-10};
349 271 : coeffmap_p[14564]={1.000,-1.393e-3,7.20e-7,-1.39e-10};
350 271 : coeffmap_p[14692]={1.000,-1.384e-3,7.03e-7,-1.31e-10};
351 271 : coeffmap_p[14820]={1.000,-1.388e-3,7.06e-7,-1.32e-10};
352 271 : coeffmap_p[14948]={1.000,-1.393e-3,7.16e-7,-1.37e-10};
353 271 : coeffmap_p[15052]={1.000,-1.402e-3,7.38e-7,-1.48e-10};
354 271 : coeffmap_p[15180]={1.000,-1.407e-3,7.47e-7,-1.53e-10};
355 271 : coeffmap_p[15308]={1.000,-1.406e-3,7.41e-7,-1.48e-10};
356 271 : coeffmap_p[15436]={1.000,-1.399e-3,7.31e-7,-1.44e-10};
357 271 : coeffmap_p[15564]={1.000,-1.397e-3,7.28e-7,-1.43e-10};
358 271 : coeffmap_p[15692]={1.000,-1.401e-3,7.35e-7,-1.46e-10};
359 271 : coeffmap_p[15820]={1.000,-1.402e-3,7.34e-7,-1.45e-10};
360 271 : coeffmap_p[15948]={1.000,-1.399e-3,7.30e-7,-1.44e-10};
361 271 : coeffmap_p[16052]={1.000,-1.419e-3,7.59e-7,-1.54e-10};
362 271 : coeffmap_p[16180]={1.000,-1.419e-3,7.59e-7,-1.52e-10};
363 271 : coeffmap_p[16308]={1.000,-1.412e-3,7.40e-7,-1.44e-10};
364 271 : coeffmap_p[16436]={1.000,-1.407e-3,7.32e-7,-1.40e-10};
365 271 : coeffmap_p[16564]={1.000,-1.408e-3,7.32e-7,-1.41e-10};
366 271 : coeffmap_p[16692]={1.000,-1.410e-3,7.34e-7,-1.40e-10};
367 271 : coeffmap_p[16820]={1.000,-1.407e-3,7.27e-7,-1.38e-10};
368 271 : coeffmap_p[16948]={1.000,-1.423e-3,7.63e-7,-1.55e-10};
369 271 : coeffmap_p[17052]={1.000,-1.437e-3,7.87e-7,-1.66e-10};
370 271 : coeffmap_p[17180]={1.000,-1.438e-3,7.84e-7,-1.64e-10};
371 271 : coeffmap_p[17308]={1.000,-1.445e-3,7.98e-7,-1.71e-10};
372 271 : coeffmap_p[17436]={1.000,-1.452e-3,8.10e-7,-1.77e-10};
373 271 : coeffmap_p[17564]={1.000,-1.458e-3,8.13e-7,-1.70e-10};
374 271 : coeffmap_p[17692]={1.000,-1.456e-3,8.06e-7,-1.72e-10};
375 271 : coeffmap_p[17820]={1.000,-1.453e-3,8.00e-7,-1.68e-10};
376 271 : coeffmap_p[17948]={1.000,-1.452e-3,7.99e-7,-1.69e-10};
377 : /////K
378 271 : coeffmap_p[19052]={1.000,-1.419e-3,7.56e-7,-1.53e-10};
379 271 : coeffmap_p[19180]={1.000,-1.426e-3,7.70e-7,-1.59e-10};
380 271 : coeffmap_p[19308]={1.000,-1.433e-3,7.82e-7,-1.64e-10};
381 271 : coeffmap_p[19436]={1.000,-1.429e-3,7.73e-7,-1.60e-10};
382 271 : coeffmap_p[19564]={1.000,-1.427e-3,7.70e-7,-1.59e-10};
383 271 : coeffmap_p[19692]={1.000,-1.425e-3,7.65e-7,-1.56e-10};
384 271 : coeffmap_p[19820]={1.000,-1.430e-3,7.76e-7,-1.62e-10};
385 271 : coeffmap_p[19948]={1.000,-1.434e-3,7.81e-7,-1.63e-10};
386 271 : coeffmap_p[21052]={1.000,-1.448e-3,8.05e-7,-1.73e-10};
387 271 : coeffmap_p[21180]={1.000,-1.436e-3,7.84e-7,-1.63e-10};
388 271 : coeffmap_p[21308]={1.000,-1.441e-3,7.94e-7,-1.68e-10};
389 271 : coeffmap_p[21436]={1.000,-1.439e-3,7.89e-7,-1.66e-10};
390 271 : coeffmap_p[21564]={1.000,-1.442e-3,7.96e-7,-1.69e-10};
391 271 : coeffmap_p[21692]={1.000,-1.435e-3,7.81e-7,-1.61e-10};
392 271 : coeffmap_p[21820]={1.000,-1.442e-3,7.92e-7,-1.66e-10};
393 271 : coeffmap_p[21948]={1.000,-1.439e-3,7.82e-7,-1.61e-10};
394 271 : coeffmap_p[23052]={1.000,-1.401e-3,7.21e-7,-1.37e-10};
395 271 : coeffmap_p[23180]={1.000,-1.408e-3,7.31e-7,-1.41e-10};
396 271 : coeffmap_p[23308]={1.000,-1.407e-3,7.28e-7,-1.39e-10};
397 271 : coeffmap_p[23436]={1.000,-1.407e-3,7.31e-7,-1.41e-10};
398 271 : coeffmap_p[23564]={1.000,-1.419e-3,7.47e-7,-1.47e-10};
399 271 : coeffmap_p[23692]={1.000,-1.395e-3,7.10e-7,-1.33e-10};
400 271 : coeffmap_p[23820]={1.000,-1.413e-3,7.36e-7,-1.42e-10};
401 271 : coeffmap_p[23948]={1.000,-1.402e-3,7.21e-7,-1.36e-10};
402 271 : coeffmap_p[25052]={1.000,-1.402e-3,7.17e-7,-1.31e-10};
403 271 : coeffmap_p[25180]={1.000,-1.432e-3,7.73e-7,-1.58e-10};
404 271 : coeffmap_p[25308]={1.000,-1.407e-3,7.22e-7,-1.33e-10};
405 271 : coeffmap_p[25436]={1.000,-1.417e-3,7.43e-7,-1.45e-10};
406 271 : coeffmap_p[25564]={1.000,-1.422e-3,7.52e-7,-1.48e-10};
407 271 : coeffmap_p[25692]={1.000,-1.427e-3,7.59e-7,-1.52e-10};
408 271 : coeffmap_p[25820]={1.000,-1.416e-3,7.42e-7,-1.44e-10};
409 271 : coeffmap_p[25948]={1.000,-1.422e-3,7.46e-7,-1.45e-10};
410 : ///A
411 271 : coeffmap_p[28052]={1.000,-1.444e-3,7.61e-7,-1.44e-10};
412 271 : coeffmap_p[28180]={1.000,-1.439e-3,7.54e-7,-1.42e-10};
413 271 : coeffmap_p[28308]={1.000,-1.457e-3,7.87e-7,-1.58e-10};
414 271 : coeffmap_p[28436]={1.000,-1.457e-3,7.90e-7,-1.60e-10};
415 271 : coeffmap_p[28564]={1.000,-1.455e-3,7.87e-7,-1.59e-10};
416 271 : coeffmap_p[28692]={1.000,-1.458e-3,7.88e-7,-1.58e-10};
417 271 : coeffmap_p[28820]={1.000,-1.453e-3,7.81e-7,-1.56e-10};
418 271 : coeffmap_p[28948]={1.000,-1.460e-3,7.98e-7,-1.64e-10};
419 271 : coeffmap_p[31052]={1.000,-1.415e-3,7.44e-7,-1.44e-10};
420 271 : coeffmap_p[31180]={1.000,-1.408e-3,7.26e-7,-1.37e-10};
421 271 : coeffmap_p[31308]={1.000,-1.413e-3,7.28e-7,-1.36e-10};
422 271 : coeffmap_p[31436]={1.000,-1.394e-3,7.07e-7,-1.30e-10};
423 271 : coeffmap_p[31564]={1.000,-1.404e-3,7.23e-7,-1.37e-10};
424 271 : coeffmap_p[31692]={1.000,-1.427e-3,7.48e-7,-1.44e-10};
425 271 : coeffmap_p[31820]={1.000,-1.418e-3,7.48e-7,-1.48e-10};
426 271 : coeffmap_p[31948]={1.000,-1.413e-3,7.37e-7,-1.42e-10};
427 271 : coeffmap_p[34052]={1.000,-1.42e-3,7.28e-7,-1.34e-10};
428 271 : coeffmap_p[34180]={1.000,-1.46e-3,7.77e-7,-1.53e-10};
429 271 : coeffmap_p[34308]={1.000,-1.42e-3,7.41e-7,-1.42e-10};
430 271 : coeffmap_p[34436]={1.000,-1.42e-3,7.36e-7,-1.39e-10};
431 271 : coeffmap_p[34564]={1.000,-1.46e-3,7.76e-7,-1.52e-10};
432 271 : coeffmap_p[34692]={1.000,-1.42e-3,7.34e-7,-1.38e-10};
433 271 : coeffmap_p[34820]={1.000,-1.42e-3,7.34e-7,-1.39e-10};
434 271 : coeffmap_p[34948]={1.000,-1.45e-3,7.68e-7,-1.49e-10};
435 271 : coeffmap_p[37152]={1.000,-1.42e-3,7.47e-7,-1.44e-10};
436 271 : coeffmap_p[37280]={1.000,-1.41e-3,7.35e-7,-1.40e-10};
437 271 : coeffmap_p[37408]={1.000,-1.45e-3,7.65e-7,-1.46e-10};
438 271 : coeffmap_p[37536]={1.000,-1.41e-3,7.13e-7,-1.29e-10};
439 271 : coeffmap_p[37664]={1.000,-1.41e-3,7.30e-7,-1.38e-10};
440 271 : coeffmap_p[37792]={1.000,-1.45e-3,7.75e-7,-1.50e-10};
441 271 : coeffmap_p[37820]={1.000,-1.45e-3,7.68e-7,-1.49e-10};
442 271 : coeffmap_p[38048]={1.000,-1.41e-3,7.38e-7,-1.43e-10};
443 : //Q
444 271 : coeffmap_p[41052]={1.000,-1.453e-3,7.69e-7,-1.47e-10};
445 271 : coeffmap_p[41180]={1.000,-1.479e-3,8.03e-7,-1.61e-10};
446 271 : coeffmap_p[41308]={1.000,-1.475e-3,7.97e-7,-1.58e-10};
447 271 : coeffmap_p[41436]={1.000,-1.451e-3,7.73e-7,-1.51e-10};
448 271 : coeffmap_p[41564]={1.000,-1.450e-3,7.71e-7,-1.51e-10};
449 271 : coeffmap_p[41692]={1.000,-1.465e-3,7.79e-7,-1.49e-10};
450 271 : coeffmap_p[41820]={1.000,-1.460e-3,7.73e-7,-1.47e-10};
451 271 : coeffmap_p[41948]={1.000,-1.434e-3,7.47e-7,-1.40e-10};
452 271 : coeffmap_p[43052]={1.000,-1.428e-3,7.40e-7,-1.38e-10};
453 271 : coeffmap_p[43180]={1.000,-1.418e-3,7.29e-7,-1.34e-10};
454 271 : coeffmap_p[43308]={1.000,-1.433e-3,7.49e-7,-1.43e-10};
455 271 : coeffmap_p[43436]={1.000,-1.438e-3,7.55e-7,-1.45e-10};
456 271 : coeffmap_p[43564]={1.000,-1.419e-3,7.36e-7,-1.40e-10};
457 271 : coeffmap_p[43692]={1.000,-1.397e-3,7.13e-7,-1.33e-10};
458 271 : coeffmap_p[43820]={1.000,-1.423e-3,7.39e-7,-1.40e-10};
459 271 : coeffmap_p[43948]={1.000,-1.452e-3,7.68e-7,-1.47e-10};
460 :
461 :
462 :
463 :
464 :
465 :
466 :
467 271 : }
468 :
469 :
470 :
471 :
472 :
473 :
474 : } //# NAMESPACE CASA - END
|
