2 '''Calculate convergence angle from UTM coordinates in degrees
3 Requires `pip install docopt numpy utm km3astro`
6 get_convergence_angle.py <easting> <northing> <zone> <letter>
9 easting UTM Easting of the site
10 northing UTM Northing of the site
14 from docopt
import docopt
15 from km3astro.coord
import convergence_angle
19 args = docopt(__doc__, version=
'0.1')
20 easting = args[
'<easting>']
21 northing = args[
'<northing>']
23 letter = args[
'<letter>']
25 latitude, longitude = utm.to_latlon(float(easting), float(northing), int(zone), letter)
26 conv_angle = convergence_angle(latitude/180.*numpy.pi,longitude/180.*numpy.pi)*180./numpy.pi
27 print(
"lat, long", latitude, longitude)
28 print(
"convergence angle",conv_angle)
std::ostream & print(std::ostream &out, const JTestSummary &summary, const char delimiter=' ', const bool useColors=true)
Print test summary.