# Metview Macro

# **************************** LICENSE START ***********************************
#
# Copyright 2012 ECMWF. This software is distributed under the terms
# of the Apache License version 2.0. In applying this license, ECMWF does not
# waive the privileges and immunities granted to it by virtue of its status as
# an Intergovernmental Organization or submit itself to any jurisdiction.
#
# ***************************** LICENSE END ************************************


tuv_data = read("../TUV_Data")
tuv_data = tuv_data[1,3]

lat = 50
lon = 10

p = nearest_gridpoint(tuv_data, lat, lon)
print('p: ', p)

n = interpolate(tuv_data, lat, lon)
print('n: ', n)

i = nearest_gridpoint_info(tuv_data, lat, lon)
print('i: ', i)

print('first latitude: ', i[1].latitude)

