#!/usr/bin/env python
import numpy as np
import meta, pyproj, sord

# parameters
extent = meta.extent
projection = meta.projection
transform = meta.transform
proj = sord.coord.Transform( pyproj.Proj( **projection ), **transform )

# topo
topo, extent = sord.data.topo( extent, path='mapdata' )
ddeg = 0.5 / 60.0
lon, lat = extent

# mesh
x, y = np.loadtxt( 'src/grid-1.0x1.0-gp1_h0.125.ll', usecols=(0,1), unpack=True )
z = sord.coord.interp2( lon[0], lat[0], ddeg, ddeg, topo, x, y )
np.array( x, 'f' ).tofile( 'lon' )
np.array( y, 'f' ).tofile( 'lat' )
np.array( z, 'f' ).tofile( 'topo' )