def set_retrieval_discretization(self, dx=[2., 1., 0.5], dy=[2., 1., 0.5]):
l = []
p = []
for x, y in zip(dx, dy):
d = discretize( self.border, x, y, self.link['Trx'], self.link['Try'],
self.link['Recx'], self.link['Recy'])
l.append(d[0])
p.append(d[1])
self.nesting_lij = l
self.nesting_pij = p
def set_rsl_discretization( self, resolx=0.25, resoly=0.25 ):
'''
Sets discretization for RSL at a given discretization step.
Inputs:
resolx - scalar. Discretization in x axis, [km].
By default, resolx = 0.25
resoly - scalar. Discretization in y axis, [km].
By default, resoly = 0.25
'''
# discretize the area with (resolx,resoly) sq.km resolution
self.lij, self.pij = discretize( self.border, resolx, resoly,
self.link['Trx'], self.link['Try'],
self.link['Recx'], self.link['Recy'])