Example #1
0
 def test_Fstdc_ezgetlalo_KnownValues2(self):
     """Fstdc_ezgetlalo corners should give known result with known input"""
     (ni,nj) = self.la.shape
     grtyp='L'
     grref='L'
     (ig1,ig2,ig3,ig4) =  Fstdc.cxgaig(grtyp,-69.5,180.0,0.5,0.5)
     hasAxes = 0
     doCorners = 1
     (i0,j0) = (0,0)
     (la2,lo2,cla2,clo2) = Fstdc.ezgetlalo((ni,nj),grtyp,(grref,ig1,ig2,ig3,ig4),(None,None),hasAxes,(i0,j0),doCorners)
     if numpy.any(self.la!=la2):
         print("\nLAT Expected:\n",self.la)
         print("LAT Got:\n",la2)
     if numpy.any(self.lo!=lo2):
         print("\nLON Expected:\n",self.lo)
         print("LON Got:\n",lo2)
     self.assertFalse(numpy.any(self.la!=la2))
     self.assertFalse(numpy.any(self.lo!=lo2))
     for ic in range(0,4):
         if numpy.any(self.cla[ic,...]!=cla2[ic,...]):
             print("\n",ic,'cla')
             print("LAT Expected:\n",self.cla[ic,...])
             print("LAT Got:\n",cla2[ic,...])
         self.assertFalse(numpy.any(self.cla[ic,...]!=cla2[ic,...]))
         if numpy.any(self.clo[ic,...]!=clo2[ic,...]):
             print("\n",ic,'clo')
             print("\nLON Expected:\n",self.clo[ic,...])
             print("LON Got:\n",clo2[ic,...])
         self.assertFalse(numpy.any(self.clo[ic,...]!=clo2[ic,...]))
Example #2
0
 def test_Fstdc_ezgetlalo_KnownValues2(self):
     """Fstdc_ezgetlalo corners should give known result with known input"""
     (ni,nj) = self.la.shape
     grtyp='L'
     grref='L'
     (ig1,ig2,ig3,ig4) =  Fstdc.cxgaig(grtyp,-69.5,180.0,0.5,0.5)
     hasAxes = 0
     doCorners = 1
     (i0,j0) = (0,0)
     (la2,lo2,cla2,clo2) = Fstdc.ezgetlalo((ni,nj),grtyp,(grref,ig1,ig2,ig3,ig4),(None,None),hasAxes,(i0,j0),doCorners)
     if numpy.any(self.la!=la2):
         print "\nLAT Expected:\n",self.la
         print "LAT Got:\n",la2
     if numpy.any(self.lo!=lo2):
         print "\nLON Expected:\n",self.lo
         print "LON Got:\n",lo2
     self.assertFalse(numpy.any(self.la!=la2))
     self.assertFalse(numpy.any(self.lo!=lo2))
     for ic in range(0,4):
         if numpy.any(self.cla[ic,...]!=cla2[ic,...]):
             print "\n",ic,'cla'
             print "LAT Expected:\n",self.cla[ic,...]
             print "LAT Got:\n",cla2[ic,...]
         self.assertFalse(numpy.any(self.cla[ic,...]!=cla2[ic,...]))
         if numpy.any(self.clo[ic,...]!=clo2[ic,...]):
             print "\n",ic,'clo'
             print "\nLON Expected:\n",self.clo[ic,...]
             print "LON Got:\n",clo2[ic,...]
         self.assertFalse(numpy.any(self.clo[ic,...]!=clo2[ic,...]))
Example #3
0
 def test_Fstdc_ezgetlalo_KnownValues2(self):
     """Fstdc_ezgetlalo corners should give known result with known input"""
     (ni, nj) = self.la.shape
     grtyp = "L"
     grref = "L"
     (ig1, ig2, ig3, ig4) = Fstdc.cxgaig(grtyp, -89.5, 180.0, 0.5, 0.5)
     hasAxes = 0
     doCorners = 1
     (i0, j0) = (0, 0)
     (la2, lo2, cla2, clo2) = Fstdc.ezgetlalo(
         (ni, nj), grtyp, (grref, ig1, ig2, ig3, ig4), (None, None), hasAxes, (i0, j0), doCorners
     )
     if numpy.any(self.la != la2):
         print self.la
         print la2
     if numpy.any(self.lo != lo2):
         print self.lo
         print lo2
     self.assertFalse(numpy.any(self.la != la2))
     self.assertFalse(numpy.any(self.lo != lo2))
     for ic in range(0, 4):
         if numpy.any(self.cla[ic, ...] != cla2[ic, ...]):
             print ic, "cla"
             print self.cla[ic, ...]
             print cla2[ic, ...]
         self.assertFalse(numpy.any(self.cla[ic, ...] != cla2[ic, ...]))
         if numpy.any(self.clo[ic, ...] != clo2[ic, ...]):
             print ic, "clo"
             print self.clo[ic, ...]
             print clo2[ic, ...]
         self.assertFalse(numpy.any(self.clo[ic, ...] != clo2[ic, ...]))
Example #4
0
 def gridL(self,dlalo=0.5,nij=10):
     """provide grid and rec values for other tests"""
     grtyp='L'
     grref=grtyp
     la0 = 0.-dlalo*(nij/2.)
     lo0 = 180.-dlalo*(nij/2.)
     ig14 = (ig1,ig2,ig3,ig4) =  Fstdc.cxgaig(grtyp,la0,lo0,dlalo,dlalo)
     axes = (None,None)
     hasAxes = 0
     ij0 = (1,1)
     doCorners = 0
     (la,lo) = Fstdc.ezgetlalo((nij,nij),grtyp,(grref,ig1,ig2,ig3,ig4),axes,hasAxes,ij0,doCorners)
     return (grtyp,ig14,(nij,nij),la,lo)
Example #5
0
 def gridL(self,dlalo=0.5,nij=10):
     """provide grid and rec values for other tests"""
     grtyp='L'
     grref=grtyp
     la0 = 0.-dlalo*(nij/2.)
     lo0 = 180.-dlalo*(nij/2.)
     ig14 = (ig1,ig2,ig3,ig4) =  Fstdc.cxgaig(grtyp,la0,lo0,dlalo,dlalo)
     axes = (None,None)
     hasAxes = 0
     ij0 = (1,1)
     doCorners = 0
     (la,lo) = Fstdc.ezgetlalo((nij,nij),grtyp,(grref,ig1,ig2,ig3,ig4),axes,hasAxes,ij0,doCorners)
     return (grtyp,ig14,(nij,nij),la,lo)
Example #6
0
 def test_Fstdc_ezgetlalo_KnownValues(self):
     """Fstdc_ezgetlalo should give known result with known input"""
     (ni,nj) = self.la.shape
     grtyp='L'
     grref='L'
     (ig1,ig2,ig3,ig4) =  Fstdc.cxgaig(grtyp,-69.5,180.0,0.5,0.5)
     hasAxes = 0
     doCorners = 0
     (i0,j0) = (0,0)
     (la2,lo2) = Fstdc.ezgetlalo((ni,nj),grtyp,(grref,ig1,ig2,ig3,ig4),(None,None),hasAxes,(i0,j0),doCorners)
     if numpy.any(self.la!=la2):
         print "\nLAT Expected:\n",self.la
         print "LAT Got:\n",la2
     if numpy.any(self.lo!=lo2):
         print "\nLON Expected:\n",self.lo
         print "LON Got:\n",lo2
     self.assertFalse(numpy.any(self.la!=la2))
     self.assertFalse(numpy.any(self.lo!=lo2))
Example #7
0
 def test_Fstdc_ezgetlalo_Dieze_KnownValues(self):
     """Fstdc_ezgetlalo with #-grid should give known result with known input"""
     (ni,nj) = self.la.shape
     grtyp='#'
     grref='L'
     (ig1,ig2,ig3,ig4) =  Fstdc.cxgaig(grref,0.,0.,1.,1.)
     xaxis = self.lo[:,0].reshape((self.lo.shape[0],1)).copy('FORTRAN')
     yaxis = self.la[0,:].reshape((1,self.la.shape[1])).copy('FORTRAN')
     hasAxes = 1
     doCorners = 0
     (i0,j0) = (2,2)
     (la2,lo2) = Fstdc.ezgetlalo((ni-1,nj-1),grtyp,(grref,ig1,ig2,ig3,ig4),(xaxis,yaxis),hasAxes,(i0,j0),doCorners)
     if numpy.any(self.la[1:,1:]!=la2):
         print "\nLAT Expected:\n",self.la[1:,1:]
         print "LAT Got:\n",la2
     if numpy.any(self.lo[1:,1:]!=lo2):
         print "\nLON Expected:\n",self.lo[1:,1:]
         print "LON Got:\n",lo2
     self.assertFalse(numpy.any(self.la[1:,1:]!=la2))
     self.assertFalse(numpy.any(self.lo[1:,1:]!=lo2))
Example #8
0
 def test_Fstdc_ezgetlalo_Dieze_KnownValues(self):
     """Fstdc_ezgetlalo with #-grid should give known result with known input"""
     (ni,nj) = self.la.shape
     grtyp='#'
     grref='L'
     (ig1,ig2,ig3,ig4) =  Fstdc.cxgaig(grref,0.,0.,1.,1.)
     xaxis = self.lo[:,0].reshape((self.lo.shape[0],1)).copy('FORTRAN')
     yaxis = self.la[0,:].reshape((1,self.la.shape[1])).copy('FORTRAN')
     hasAxes = 1
     doCorners = 0
     (i0,j0) = (2,2)
     (la2,lo2) = Fstdc.ezgetlalo((ni-1,nj-1),grtyp,(grref,ig1,ig2,ig3,ig4),(xaxis,yaxis),hasAxes,(i0,j0),doCorners)
     if numpy.any(self.la[1:,1:]!=la2):
         print("\nLAT Expected:\n",self.la[1:,1:])
         print("LAT Got:\n",la2)
     if numpy.any(self.lo[1:,1:]!=lo2):
         print("\nLON Expected:\n",self.lo[1:,1:])
         print("LON Got:\n",lo2)
     self.assertFalse(numpy.any(self.la[1:,1:]!=la2))
     self.assertFalse(numpy.any(self.lo[1:,1:]!=lo2))
Example #9
0
 def test_Fstdc_ezgetlalo_KnownValues(self):
     """Fstdc_ezgetlalo should give known result with known input"""
     (ni, nj) = self.la.shape
     grtyp = 'L'
     grref = 'L'
     (ig1, ig2, ig3, ig4) = Fstdc.cxgaig(grtyp, -89.5, 180.0, 0.5, 0.5)
     hasAxes = 0
     doCorners = 0
     (i0, j0) = (0, 0)
     (la2, lo2) = Fstdc.ezgetlalo((ni, nj), grtyp,
                                  (grref, ig1, ig2, ig3, ig4), (None, None),
                                  hasAxes, (i0, j0), doCorners)
     if numpy.any(self.la != la2):
         print self.la
         print la2
     if numpy.any(self.lo != lo2):
         print self.lo
         print lo2
     self.assertFalse(numpy.any(self.la != la2))
     self.assertFalse(numpy.any(self.lo != lo2))
Example #10
0
 def test_Fstdc_ezgetlalo_Dieze_KnownValues(self):
     """Fstdc_ezgetlalo with #-grid should give known result with known input"""
     (ni, nj) = self.la.shape
     grtyp = "#"
     grref = "L"
     (ig1, ig2, ig3, ig4) = Fstdc.cxgaig(grref, 0.0, 0.0, 1.0, 1.0)
     xaxis = self.lo[:, 0].reshape((self.lo.shape[0], 1)).copy("FORTRAN")
     yaxis = self.la[0, :].reshape((1, self.la.shape[1])).copy("FORTRAN")
     hasAxes = 1
     doCorners = 0
     (i0, j0) = (2, 2)
     (la2, lo2) = Fstdc.ezgetlalo(
         (ni - 1, nj - 1), grtyp, (grref, ig1, ig2, ig3, ig4), (xaxis, yaxis), hasAxes, (i0, j0), doCorners
     )
     if numpy.any(self.la[1:, 1:] != la2):
         print self.la[1:, 1:]
         print la2
     if numpy.any(self.lo[1:, 1:] != lo2):
         print self.lo[1:, 1:]
         print lo2
     self.assertFalse(numpy.any(self.la[1:, 1:] != la2))
     self.assertFalse(numpy.any(self.lo[1:, 1:] != lo2))
Example #11
0
 def test_Fstdc_ezgetlalo_Z_KnownValues(self):
     """Fstdc_ezgetlalo with Z grid should give known result with known input"""
     (ni, nj) = self.la.shape
     grtyp = 'Z'
     grref = 'L'
     (ig1, ig2, ig3, ig4) = Fstdc.cxgaig(grref, 0., 0., 1., 1.)
     xaxis = self.lo[:, 0].reshape((self.lo.shape[0], 1)).copy('FORTRAN')
     yaxis = self.la[0, :].reshape((1, self.la.shape[1])).copy('FORTRAN')
     hasAxes = 1
     doCorners = 0
     (i0, j0) = (0, 0)
     (la2, lo2) = Fstdc.ezgetlalo(
         (ni, nj), grtyp, (grref, ig1, ig2, ig3, ig4), (xaxis, yaxis),
         hasAxes, (i0, j0), doCorners)
     if numpy.any(self.la != la2):
         print self.la
         print la2
     if numpy.any(self.lo != lo2):
         print self.lo
         print lo2
     self.assertFalse(numpy.any(self.la != la2))
     self.assertFalse(numpy.any(self.lo != lo2))