Commit d4252121 authored by Pierre Cazenave's avatar Pierre Cazenave

Add script to convert the grid file to lat/long.

parent a4d12cc4
#!/usr/bin/env bash
# Convert the input to lat/long (assuming the proj line in the nml is
# right...).
awk '{if (NF == 4 && NR > 3) print , }' tst_grd.dat | \
invproj -f %.6f \
+proj=tmerc \
+datum=NAD83 \
+lon_0=-70d10 \
+lat_0=42d50 \
+k=.9999666666666667 \
+x_0=900000 \
+y_0=0 \
> grid.csv
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment