GridNavRetriever.GridNavRetriever.__init__() - Code Metrics - Unidata/python-awips - Measure and Improve Code Quality continuously with Scrutinizer

GridNavRetriever.GridNavRetriever.init() A
last analyzed 2020-10-19 17:29 UTC

↳ Parent: GridNavRetriever

Complexity

Conditions

Size

Total Lines	7
Code Lines	7

Duplication

Lines	0
Ratio	0 %

Importance

Changes

Metric	Value
eloc	7
dl	0
loc	7
rs	10
c	0
b	0
f	0
cc	1
nop	5

import os
import math
from awips import ThriftClient
from dynamicserialize.dstypes.gov.noaa.nws.ncep.common.dataplugin.gempak.request import GetGridNavRequest
from ctypes import *

EARTH_RADIUS = 6371200.0
DEG_TO_RAD = math.pi / 180.0
RAD_TO_DEG = 180.0 / math.pi
TWOPI = math.pi * 2.0
HALFPI = math.pi / 2.0
PI4TH = math.pi / 4.0
PI3RD = math.pi / 3.0


def createPolar(nsflag, clon, lat1, lon1, dx, dy, unit, nx, ny):
    clonr = clon * DEG_TO_RAD
    latr = lat1 * DEG_TO_RAD
    lonr = lon1 * DEG_TO_RAD
    if nsflag == 'N':
        x1 = EARTH_RADIUS * math.tan(PI4TH - latr/2.0) * math.sin(lonr-clonr)
        y1 = -1 * EARTH_RADIUS * math.tan(PI4TH - latr/2.0) * math.cos(lonr-clonr)
    else:
        x1 = EARTH_RADIUS * math.tan(PI4TH + latr/2.0) * math.sin(lonr-clonr)
        y1 = EARTH_RADIUS * math.tan(PI4TH + latr/2.0) * math.cos(lonr-clonr)

    if unit == 'm':
        tdx = dx / (1 + math.sin(PI3RD))
        tdy = dy / (1 + math.sin(PI3RD))
    else:
        tdx = (dx*1000.0) / (1 + math.sin(PI3RD))
        tdy = (dy*1000.0) / (1 + math.sin(PI3RD))

    x2 = x1 + tdx * (nx-1)
    y2 = y1 + tdy * (ny-1)
    xll = min(x1, x2)
    yll = min(y1, y2)
    xur = max(x1, x2)
    yur = max(y1, y2)

    if nsflag == 'N':
        latll = (HALFPI - 2*math.atan2(math.hypot(xll, yll), EARTH_RADIUS)) * RAD_TO_DEG
        rtemp = clonr + math.atan2(xll, -yll)
    else:
        latll = -1 * (HALFPI - 2*math.atan2(math.hypot(xll, yll), EARTH_RADIUS)) * RAD_TO_DEG
        rtemp = clonr + math.atan2(xll, yll)

    if rtemp > math.pi:
        lonll = (rtemp-TWOPI) * RAD_TO_DEG
    elif rtemp < -math.pi:
        lonll = (rtemp+TWOPI) * RAD_TO_DEG
    else:
        lonll = rtemp * RAD_TO_DEG

    if nsflag == 'N':
        latur = (HALFPI - 2*math.atan2(math.hypot(xur, yur), EARTH_RADIUS)) * RAD_TO_DEG
        rtemp = clonr + math.atan2(xur, -yur)
    else:
        latur = -1 * (HALFPI - 2*math.atan2(math.hypot(xur, yur), EARTH_RADIUS)) * RAD_TO_DEG
        rtemp = clonr + math.atan2(xur, yur)

    if rtemp > math.pi:
        lonur = (rtemp-TWOPI) * RAD_TO_DEG
    elif rtemp < -math.pi:
        lonur = (rtemp+TWOPI) * RAD_TO_DEG
    else:
        lonur = rtemp * RAD_TO_DEG

    return [latll, lonll, latur, lonur]


def createConic(nsflag, clon, lat1, lon1, dx, dy, unit, nx, ny, ang1, ang3):
    clonr = clon * DEG_TO_RAD
    latr = lat1 * DEG_TO_RAD
    lonr = lon1 * DEG_TO_RAD

    angle1 = HALFPI - (math.fabs(ang1) * DEG_TO_RAD)
    angle2 = HALFPI - (math.fabs(ang3) * DEG_TO_RAD)

    if ang1 == ang3:
        cc = math.cos(angle1)
    else:
        cc = (math.log(math.sin(angle2)) - math.log(math.sin(angle1))) \
             / (math.log(math.tan(angle2/2.0)) - math.log(math.tan(angle1/2.0)))

    er = EARTH_RADIUS / cc

    if nsflag == 'N':
        x1 = er * math.pow(math.tan((HALFPI-latr)/2.0), cc) * math.sin(cc*(lonr-clonr))
        y1 = -1.0 * er * math.pow(math.tan((HALFPI-latr)/2.0), cc) * math.cos(cc*(lonr-clonr))
    else:
        x1 = er * math.pow(math.tan((HALFPI+latr)/2.0), cc) * math.sin(cc*(lonr-clonr))
        y1 = er * math.pow(math.tan((HALFPI+latr)/2.0), cc) * math.cos(cc*(lonr-clonr))

    alpha = math.pow(math.tan(angle1/2.0), cc) / math.sin(angle1)

    if unit == 'm':
        x2 = x1 + (nx-1) * alpha * dx
        y2 = y1 + (ny-1) * alpha * dy
    else:
        x2 = x1 + (nx-1) * alpha * (dx*1000.0)
        y2 = y1 + (ny-1) * alpha * (dy*1000.0)

    xll = min(x1, x2)
    yll = min(y1, y2)
    xur = max(x1, x2)
    yur = max(y1, y2)

    if nsflag == 'N':
        latll = (HALFPI - 2.0 * math.atan(math.pow(math.hypot(xll, yll)/er, (1/cc)))) * RAD_TO_DEG
        rtemp = math.atan2(xll, -yll) * (1/cc) + clonr
    else:
        latll = (-1.0 * (HALFPI - 2.0 * math.atan(math.pow(math.hypot(xll, yll)/er, (1/cc))))) * RAD_TO_DEG
        rtemp = math.atan2(xll, yll) * (1/cc) + clonr

    if rtemp > math.pi:
        lonll = (rtemp-TWOPI) * RAD_TO_DEG
    elif rtemp < -math.pi:
        lonll = (rtemp+TWOPI) * RAD_TO_DEG
    else:
        lonll = rtemp * RAD_TO_DEG

    if nsflag == 'N':
        latur = (HALFPI - 2.0 * math.atan(math.pow(math.hypot(xur, yur)/er, (1/cc)))) * RAD_TO_DEG
        rtemp = math.atan2(xur, -yur) * (1/cc) + clonr
    else:
        latur = (-1.0 * (HALFPI - 2.0 * math.atan(math.pow(math.hypot(xur, yur)/er, (1/cc))))) * RAD_TO_DEG
        rtemp = math.atan2(xur, yur) * (1/cc) + clonr

    if rtemp > math.pi:
        lonur = (rtemp-TWOPI) * RAD_TO_DEG
    elif rtemp < -math.pi:
        lonur = (rtemp+TWOPI) * RAD_TO_DEG
    else:
        lonur = rtemp * RAD_TO_DEG

    return [latll, lonll, latur, lonur]


class StringConverter(Union):
    _fields_ = [("char", c_char*4), ("int", c_int), ("float", c_float)]



class GridNavRetriever:

    def __init__(self, server, pluginName, modelId, arrayLen):
        self.pluginName = pluginName
        self.modelId = modelId
        self.arrayLen = arrayLen
        self.host = os.getenv("DEFAULT_HOST", server)
        self.port = os.getenv("DEFAULT_PORT", "9581")
        self.client = ThriftClient.ThriftClient(self.host, self.port)

    def getNavBlk(self):
        """ Sends ThriftClient request and writes out received files."""
        req = GetGridNavRequest()
        req.setPluginName(self.pluginName)
        req.setModelId(self.modelId)
        resp = self.client.sendRequest(req)

        for i, rec in enumerate(resp):
            resp[i] = {
                key.decode() if isinstance(key, bytes) else key:
                    val.decode() if isinstance(val, bytes) else val
                for key, val in rec.items()
            }

        nav = []

        for record in resp:
            unit = record['spacingunit']
            sk = record['spatialkey']
            skarr = sk.split('/')

            nx = float(skarr[1])
            ny = float(skarr[2])
            dx = float(skarr[3])
            dy = float(skarr[4])

            sc = StringConverter()
            if record['projtype'] == 'LatLon':
                sc.char = 'CED '
                gemproj = 2.0
                ang1 = 0.0
                ang2 = 0.0
                ang3 = 0.0

                lllat = float(record['lowerleftlat'])
                lllon = float(record['lowerleftlon'])
                urlat = lllat + (dy * (ny-1))
                urlon = lllon + (dx * (nx-1))
                if lllon > 180:
                    lllon -= 360.0
                if urlon > 180:
                    urlon -= 360.0

            if record['projtype'] == 'Polar Stereographic':
                sc.char = 'STR '
                gemproj = 2.0
                if float(record['standard_parallel_1']) < 0.0:
                    ang1 = -90.0
                    nsflag = 'S'
                else:
                    ang1 = 90.0
                    nsflag = 'N'
                ang2 = float(record['central_meridian'])
                ang3 = 0.0

                lat1 = float(record['lowerleftlat'])
                lon1 = float(record['lowerleftlon'])
                coords = createPolar(nsflag, ang2, lat1, lon1, dx, dy, unit, nx, ny)
                lllat = coords[0]
                lllon = coords[1]
                urlat = coords[2]
                urlon = coords[3]

            if record['projtype'] == 'Lambert Conformal':
                sc.char = 'LCC '
                gemproj = 2.0

                ang1 = float(skarr[7])
                ang2 = float(record['central_meridian'])
                ang3 = float(skarr[8])
                if ang1 < 0.0:
                    nsflag = 'S'
                else:
                    nsflag = 'N'

                lat1 = float(record['lowerleftlat'])
                lon1 = float(record['lowerleftlon'])
                coords = createConic(nsflag, ang2, lat1, lon1, dx, dy, unit, nx, ny, ang1, ang3)
                lllat = coords[0]
                lllon = coords[1]
                urlat = coords[2]
                urlon = coords[3]

            # Fill up the output array of floats
            nav.append(gemproj)

            nav.append(sc.float)
            nav.append(1.0)
            nav.append(1.0)
            nav.append(nx)
            nav.append(ny)
            nav.append(lllat)

            nav.append(lllon)

            nav.append(urlat)

            nav.append(urlon)

            nav.append(ang1)

            nav.append(ang2)

            nav.append(ang3)


        for i in range(13, int(self.arrayLen)):
            nav.append(0.0)
        return nav

    def getAnlBlk(self):
        anl = []
        # Type
        anl.append(2.0)
        # Delta
        anl.append(1.0)
        # Extend area
        anl.append(0.0)
        anl.append(0.0)
        anl.append(0.0)
        anl.append(0.0)
        # Grid area
        anl.append(-90.0)
        anl.append(-180.0)
        anl.append(90.0)
        anl.append(180.0)
        # Data area
        anl.append(-90.0)
        anl.append(-180.0)
        anl.append(90.0)
        anl.append(180.0)
        for i in range(18, int(self.arrayLen)):
            anl.append(0.0)
        return anl


def getnavb(server, table, model, arrlen):
    gnr = GridNavRetriever(server, table, model, arrlen)
    return gnr.getNavBlk()


def getanlb(server, table, model, arrlen):
    gnr = GridNavRetriever(server, table, model, arrlen)
    return gnr.getAnlBlk()


# This is the standard boilerplate that runs this script as a main
if __name__ == '__main__':
    # Run Test
    srv = 'edex-cloud.unidata.ucar.edu'
    tbl = 'grid_info'
    mdl = 'NAM40'
    navlen = '256'
    print(getnavb(srv, tbl, mdl, navlen))
    anllen = '128'
    print(getanlb(srv, tbl, mdl, anllen))


1			import os
2			import math
3			from awips import ThriftClient
4			from dynamicserialize.dstypes.gov.noaa.nws.ncep.common.dataplugin.gempak.request import GetGridNavRequest
5			from ctypes import *
6
7			EARTH_RADIUS = 6371200.0
8			DEG_TO_RAD = math.pi / 180.0
9			RAD_TO_DEG = 180.0 / math.pi
10			TWOPI = math.pi * 2.0
11			HALFPI = math.pi / 2.0
12			PI4TH = math.pi / 4.0
13			PI3RD = math.pi / 3.0
14
15
16			def createPolar(nsflag, clon, lat1, lon1, dx, dy, unit, nx, ny):
17			clonr = clon * DEG_TO_RAD
18			latr = lat1 * DEG_TO_RAD
19			lonr = lon1 * DEG_TO_RAD
20			if nsflag == 'N':
21			x1 = EARTH_RADIUS * math.tan(PI4TH - latr/2.0) * math.sin(lonr-clonr)
22			y1 = -1 * EARTH_RADIUS * math.tan(PI4TH - latr/2.0) * math.cos(lonr-clonr)
23			else:
24			x1 = EARTH_RADIUS * math.tan(PI4TH + latr/2.0) * math.sin(lonr-clonr)
25			y1 = EARTH_RADIUS * math.tan(PI4TH + latr/2.0) * math.cos(lonr-clonr)
26
27			if unit == 'm':
28			tdx = dx / (1 + math.sin(PI3RD))
29			tdy = dy / (1 + math.sin(PI3RD))
30			else:
31			tdx = (dx*1000.0) / (1 + math.sin(PI3RD))
32			tdy = (dy*1000.0) / (1 + math.sin(PI3RD))
33
34			x2 = x1 + tdx * (nx-1)
35			y2 = y1 + tdy * (ny-1)
36			xll = min(x1, x2)
37			yll = min(y1, y2)
38			xur = max(x1, x2)
39			yur = max(y1, y2)
40
41			if nsflag == 'N':
42			latll = (HALFPI - 2math.atan2(math.hypot(xll, yll), EARTH_RADIUS)) RAD_TO_DEG
43			rtemp = clonr + math.atan2(xll, -yll)
44			else:
45			latll = -1 * (HALFPI - 2math.atan2(math.hypot(xll, yll), EARTH_RADIUS)) RAD_TO_DEG
46			rtemp = clonr + math.atan2(xll, yll)
47
48			if rtemp > math.pi:
49			lonll = (rtemp-TWOPI) * RAD_TO_DEG
50			elif rtemp < -math.pi:
51			lonll = (rtemp+TWOPI) * RAD_TO_DEG
52			else:
53			lonll = rtemp * RAD_TO_DEG
54
55			if nsflag == 'N':
56			latur = (HALFPI - 2math.atan2(math.hypot(xur, yur), EARTH_RADIUS)) RAD_TO_DEG
57			rtemp = clonr + math.atan2(xur, -yur)
58			else:
59			latur = -1 * (HALFPI - 2math.atan2(math.hypot(xur, yur), EARTH_RADIUS)) RAD_TO_DEG
60			rtemp = clonr + math.atan2(xur, yur)
61
62			if rtemp > math.pi:
63			lonur = (rtemp-TWOPI) * RAD_TO_DEG
64			elif rtemp < -math.pi:
65			lonur = (rtemp+TWOPI) * RAD_TO_DEG
66			else:
67			lonur = rtemp * RAD_TO_DEG
68
69			return [latll, lonll, latur, lonur]
70
71
72			def createConic(nsflag, clon, lat1, lon1, dx, dy, unit, nx, ny, ang1, ang3):
73			clonr = clon * DEG_TO_RAD
74			latr = lat1 * DEG_TO_RAD
75			lonr = lon1 * DEG_TO_RAD
76
77			angle1 = HALFPI - (math.fabs(ang1) * DEG_TO_RAD)
78			angle2 = HALFPI - (math.fabs(ang3) * DEG_TO_RAD)
79
80			if ang1 == ang3:
81			cc = math.cos(angle1)
82			else:
83			cc = (math.log(math.sin(angle2)) - math.log(math.sin(angle1))) \
84			/ (math.log(math.tan(angle2/2.0)) - math.log(math.tan(angle1/2.0)))
85
86			er = EARTH_RADIUS / cc
87
88			if nsflag == 'N':
89			x1 = er * math.pow(math.tan((HALFPI-latr)/2.0), cc) * math.sin(cc*(lonr-clonr))
90			y1 = -1.0 * er * math.pow(math.tan((HALFPI-latr)/2.0), cc) * math.cos(cc*(lonr-clonr))
91			else:
92			x1 = er * math.pow(math.tan((HALFPI+latr)/2.0), cc) * math.sin(cc*(lonr-clonr))
93			y1 = er * math.pow(math.tan((HALFPI+latr)/2.0), cc) * math.cos(cc*(lonr-clonr))
94
95			alpha = math.pow(math.tan(angle1/2.0), cc) / math.sin(angle1)
96
97			if unit == 'm':
98			x2 = x1 + (nx-1) * alpha * dx
99			y2 = y1 + (ny-1) * alpha * dy
100			else:
101			x2 = x1 + (nx-1) * alpha * (dx*1000.0)
102			y2 = y1 + (ny-1) * alpha * (dy*1000.0)
103
104			xll = min(x1, x2)
105			yll = min(y1, y2)
106			xur = max(x1, x2)
107			yur = max(y1, y2)
108
109			if nsflag == 'N':
110			latll = (HALFPI - 2.0 * math.atan(math.pow(math.hypot(xll, yll)/er, (1/cc)))) * RAD_TO_DEG
111			rtemp = math.atan2(xll, -yll) * (1/cc) + clonr
112			else:
113			latll = (-1.0 * (HALFPI - 2.0 * math.atan(math.pow(math.hypot(xll, yll)/er, (1/cc))))) * RAD_TO_DEG
114			rtemp = math.atan2(xll, yll) * (1/cc) + clonr
115
116			if rtemp > math.pi:
117			lonll = (rtemp-TWOPI) * RAD_TO_DEG
118			elif rtemp < -math.pi:
119			lonll = (rtemp+TWOPI) * RAD_TO_DEG
120			else:
121			lonll = rtemp * RAD_TO_DEG
122
123			if nsflag == 'N':
124			latur = (HALFPI - 2.0 * math.atan(math.pow(math.hypot(xur, yur)/er, (1/cc)))) * RAD_TO_DEG
125			rtemp = math.atan2(xur, -yur) * (1/cc) + clonr
126			else:
127			latur = (-1.0 * (HALFPI - 2.0 * math.atan(math.pow(math.hypot(xur, yur)/er, (1/cc))))) * RAD_TO_DEG
128			rtemp = math.atan2(xur, yur) * (1/cc) + clonr
129
130			if rtemp > math.pi:
131			lonur = (rtemp-TWOPI) * RAD_TO_DEG
132			elif rtemp < -math.pi:
133			lonur = (rtemp+TWOPI) * RAD_TO_DEG
134			else:
135			lonur = rtemp * RAD_TO_DEG
136
137			return [latll, lonll, latur, lonur]
138
139
140			class StringConverter(Union):
141			_fields_ = [("char", c_char*4), ("int", c_int), ("float", c_float)]
			0 ignored issues – show Comprehensibility Best Practice introduced 2018-10-30 15:45 UTC by Report Bug Copy Issue Report The variable `c_int` does not seem to be defined. Loading history... Comprehensibility Best Practice introduced 2018-10-30 15:45 UTC by Report Bug Copy Issue Report The variable `c_float` does not seem to be defined. Loading history... Comprehensibility Best Practice introduced 2018-10-30 15:45 UTC by Report Bug Copy Issue Report The variable `c_char` does not seem to be defined. Loading history...
142
143
144			class GridNavRetriever:
145
146			def __init__(self, server, pluginName, modelId, arrayLen):
147			self.pluginName = pluginName
148			self.modelId = modelId
149			self.arrayLen = arrayLen
150			self.host = os.getenv("DEFAULT_HOST", server)
151			self.port = os.getenv("DEFAULT_PORT", "9581")
152			self.client = ThriftClient.ThriftClient(self.host, self.port)
153
154			def getNavBlk(self):
155			""" Sends ThriftClient request and writes out received files."""
156			req = GetGridNavRequest()
157			req.setPluginName(self.pluginName)
158			req.setModelId(self.modelId)
159			resp = self.client.sendRequest(req)
160
161			for i, rec in enumerate(resp):
162			resp[i] = {
163			key.decode() if isinstance(key, bytes) else key:
164			val.decode() if isinstance(val, bytes) else val
165			for key, val in rec.items()
166			}
167
168			nav = []
169
170			for record in resp:
171			unit = record['spacingunit']
172			sk = record['spatialkey']
173			skarr = sk.split('/')
174
175			nx = float(skarr[1])
176			ny = float(skarr[2])
177			dx = float(skarr[3])
178			dy = float(skarr[4])
179
180			sc = StringConverter()
181			if record['projtype'] == 'LatLon':
182			sc.char = 'CED '
183			gemproj = 2.0
184			ang1 = 0.0
185			ang2 = 0.0
186			ang3 = 0.0
187
188			lllat = float(record['lowerleftlat'])
189			lllon = float(record['lowerleftlon'])
190			urlat = lllat + (dy * (ny-1))
191			urlon = lllon + (dx * (nx-1))
192			if lllon > 180:
193			lllon -= 360.0
194			if urlon > 180:
195			urlon -= 360.0
196
197			if record['projtype'] == 'Polar Stereographic':
198			sc.char = 'STR '
199			gemproj = 2.0
200			if float(record['standard_parallel_1']) < 0.0:
201			ang1 = -90.0
202			nsflag = 'S'
203			else:
204			ang1 = 90.0
205			nsflag = 'N'
206			ang2 = float(record['central_meridian'])
207			ang3 = 0.0
208
209			lat1 = float(record['lowerleftlat'])
210			lon1 = float(record['lowerleftlon'])
211			coords = createPolar(nsflag, ang2, lat1, lon1, dx, dy, unit, nx, ny)
212			lllat = coords[0]
213			lllon = coords[1]
214			urlat = coords[2]
215			urlon = coords[3]
216
217			if record['projtype'] == 'Lambert Conformal':
218			sc.char = 'LCC '
219			gemproj = 2.0
220
221			ang1 = float(skarr[7])
222			ang2 = float(record['central_meridian'])
223			ang3 = float(skarr[8])
224			if ang1 < 0.0:
225			nsflag = 'S'
226			else:
227			nsflag = 'N'
228
229			lat1 = float(record['lowerleftlat'])
230			lon1 = float(record['lowerleftlon'])
231			coords = createConic(nsflag, ang2, lat1, lon1, dx, dy, unit, nx, ny, ang1, ang3)
232			lllat = coords[0]
233			lllon = coords[1]
234			urlat = coords[2]
235			urlon = coords[3]
236
237			# Fill up the output array of floats
238			nav.append(gemproj)
			0 ignored issues – show introduced 2018-10-30 15:45 UTC by Report Bug Copy Issue Report The variable `gemproj` does not seem to be defined for all execution paths. Loading history...
239			nav.append(sc.float)
240			nav.append(1.0)
241			nav.append(1.0)
242			nav.append(nx)
243			nav.append(ny)
244			nav.append(lllat)
			0 ignored issues – show introduced 2018-10-30 15:45 UTC by Report Bug Copy Issue Report The variable `lllat` does not seem to be defined for all execution paths. Loading history...
245			nav.append(lllon)
			0 ignored issues – show introduced 2018-10-30 15:45 UTC by Report Bug Copy Issue Report The variable `lllon` does not seem to be defined for all execution paths. Loading history...
246			nav.append(urlat)
			0 ignored issues – show introduced 2018-10-30 15:45 UTC by Report Bug Copy Issue Report The variable `urlat` does not seem to be defined for all execution paths. Loading history...
247			nav.append(urlon)
			0 ignored issues – show introduced 2018-10-30 15:45 UTC by Report Bug Copy Issue Report The variable `urlon` does not seem to be defined for all execution paths. Loading history...
248			nav.append(ang1)
			0 ignored issues – show introduced 2018-10-30 15:45 UTC by Report Bug Copy Issue Report The variable `ang1` does not seem to be defined for all execution paths. Loading history...
249			nav.append(ang2)
			0 ignored issues – show introduced 2018-10-30 15:45 UTC by Report Bug Copy Issue Report The variable `ang2` does not seem to be defined for all execution paths. Loading history...
250			nav.append(ang3)
			0 ignored issues – show introduced 2018-10-30 15:45 UTC by Report Bug Copy Issue Report The variable `ang3` does not seem to be defined for all execution paths. Loading history...
251
252			for i in range(13, int(self.arrayLen)):
253			nav.append(0.0)
254			return nav
255
256			def getAnlBlk(self):
257			anl = []
258			# Type
259			anl.append(2.0)
260			# Delta
261			anl.append(1.0)
262			# Extend area
263			anl.append(0.0)
264			anl.append(0.0)
265			anl.append(0.0)
266			anl.append(0.0)
267			# Grid area
268			anl.append(-90.0)
269			anl.append(-180.0)
270			anl.append(90.0)
271			anl.append(180.0)
272			# Data area
273			anl.append(-90.0)
274			anl.append(-180.0)
275			anl.append(90.0)
276			anl.append(180.0)
277			for i in range(18, int(self.arrayLen)):
278			anl.append(0.0)
279			return anl
280
281
282			def getnavb(server, table, model, arrlen):
283			gnr = GridNavRetriever(server, table, model, arrlen)
284			return gnr.getNavBlk()
285
286
287			def getanlb(server, table, model, arrlen):
288			gnr = GridNavRetriever(server, table, model, arrlen)
289			return gnr.getAnlBlk()
290
291
292			# This is the standard boilerplate that runs this script as a main
293			if __name__ == '__main__':
294			# Run Test
295			srv = 'edex-cloud.unidata.ucar.edu'
296			tbl = 'grid_info'
297			mdl = 'NAM40'
298			navlen = '256'
299			print(getnavb(srv, tbl, mdl, navlen))
300			anllen = '128'
301			print(getanlb(srv, tbl, mdl, anllen))
302

Unidata / python-awips

GridNavRetriever.GridNavRetriever.__init__() A last analyzed 2020-10-19 17:29 UTC

Complexity

Size

Duplication

Importance

Duplication Side-by-Side

Filter issues like

GridNavRetriever.GridNavRetriever.init() A
last analyzed 2020-10-19 17:29 UTC