def test_abp2latlon_2(): ''' abp2latlon(): check (ne=120, ei=84, ej=79, panel=2) ''' from cs_ll import ij2ab, abp2latlon, latlon2abp ne, ngq = 120, 4 panel = 2 gi1, gj1, ei1, ej1 = 4, 4, 84, 79 gi2, gj2, ei2, ej2 = 1, 4, 85, 79 gi3, gj3, ei3, ej3 = 1, 1, 85, 80 gi4, gj4, ei4, ej4 = 4, 1, 84, 80 a1, b1 = ij2ab(ne, ngq, ei1, ej1, gi1, gj1) a2, b2 = ij2ab(ne, ngq, ei2, ej2, gi2, gj2) a3, b3 = ij2ab(ne, ngq, ei3, ej3, gi3, gj3) a4, b4 = ij2ab(ne, ngq, ei4, ej4, gi4, gj4) lat1, lon1 = abp2latlon(a1,b1,panel) lat2, lon2 = abp2latlon(a2,b2,panel) lat3, lon3 = abp2latlon(a3,b3,panel) lat4, lon4 = abp2latlon(a4,b4,panel) aa_equal([lat1,lon1], [lat2,lon2], 15) aa_equal([lat3,lon3], [lat2,lon2], 15) aa_equal([lat3,lon3], [lat4,lon4], 15) aa_equal([lat1,lon1], [lat4,lon4], 15)
def test_abp2latlon(): ''' abp2latlon(): MVP low accuracy ''' from cs_ll import ij2ab, abp2latlon ne, ngq = 120, 4 panel = 2 gi1, gj1, ei1, ej1 = 4, 4, 84, 79 gi2, gj2, ei2, ej2 = 1, 4, 85, 79 gi3, gj3, ei3, ej3 = 1, 1, 85, 80 gi4, gj4, ei4, ej4 = 4, 1, 84, 80 a1, b1 = ij2ab(ne, ngq, ei1, ej1, gi1, gj1) a2, b2 = ij2ab(ne, ngq, ei2, ej2, gi2, gj2) a3, b3 = ij2ab(ne, ngq, ei3, ej3, gi3, gj3) a4, b4 = ij2ab(ne, ngq, ei4, ej4, gi4, gj4) lat1, lon1 = abp2latlon(a1,b1,panel) lat2, lon2 = abp2latlon(a2,b2,panel) lat3, lon3 = abp2latlon(a3,b3,panel) lat4, lon4 = abp2latlon(a4,b4,panel) ''' print('') print('{:.18f}, {:.18f}'.format(a1, b1)) print('{:.18f}, {:.18f}'.format(a2, b2)) #print('{:.18f}, {:.18f}'.format(a3, b3)) #print('{:.18f}, {:.18f}'.format(a4, b4)) print('') print('{:.18f}, {:.18f}'.format(lat1, lon1)) print('{:.18f}, {:.18f}'.format(lat2, lon2)) #print('{:.18f}, {:.18f}'.format(lat3, lon3)) #print('{:.18f}, {:.18f}'.format(lat4, lon4)) ''' aa_equal([lat1,lon1], [lat2,lon2], 15) aa_equal([lat2,lon2], [lat3,lon3], 15) aa_equal([lat3,lon3], [lat4,lon4], 15) aa_equal([lat1,lon1], [lat4,lon4], 15)