Commit c91a6308 authored by Pierre Cazenave's avatar Pierre Cazenave

Fix the coordintes convert script.

parent 03ffd074
......@@ -3,7 +3,8 @@
# 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 | \
echo "lonDD,latDD" > grid.csv
awk '{if (NF == 4 && NR > 3) print $2,$3 }' tst_grd.dat | \
invproj -f %.6f \
+proj=tmerc \
+datum=NAD83 \
......@@ -11,5 +12,6 @@ awk '{if (NF == 4 && NR > 3) print , }' tst_grd.dat | \
+lat_0=42d50 \
+k=.9999666666666667 \
+x_0=900000 \
+y_0=0 \
> grid.csv
+y_0=0 |
tr " " "," \
>> 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