/* Project the 10 meter dem from LAZ to UTM.
/* This is an example of how to use both the INPUT and the OUTPUT
/* subcommands to reproject a data set that does not have a defined projection
/* (assuming we are running the AML from Arc).

grid
dem10m_utm = projectgrid (dem10m, #, #, #, 1000)
input
projection lambert_azimuth
units meters
datum none
parameters
6370997
-100 0 0
45 0 0
~
~
output
projection utm
zone 12
datum nad83
units meters
parameters
end


/* And now we start an example of how to first define a projection for
/* a data set and then reproject it (assuming we are running the AML from Arc).
/* Use USGS parameters for projection info.
projectdefine grid dem30m
projection lambert_azimuth
units meters
datum none
parameters
6370997
-100 0 0
45 0 0
~
~

/* Project the 30 meter dem from LAZ to UTM.
grid
dem30m_utm = projectgrid (dem30m, #, cubic, #, 1000)
output
projection utm
zone 12
datum nad83
units meters
parameters
end