%% Project using projection file from below

project grid gridname projected projectionfile

gridascii projected projected.xyz

%% Pit filling using ESRI fill method

grid

fill projected filled sink

quit

gridascii filled filled.xyz

%% Calculate flow direction

grid

flowdir = flowdirection( filled )

quit

gridascii flowdir flowdir.xyz

%% Calculate flow accumulation

grid

flowacc = flowaccumulation( flowdir )

quit

gridascii flowacc flowacc.xyz

---

%% To find original grid projection parameters use:

describe gridname

Sample projection parameters (projectionfile.prj):

INPUT

Projection GEOGRAPHIC

Datum NAD83

Zunits METERS

Units DD Spheroid GRS1980

PARAMETERS

OUTPUT

PROJECTION UTM

DATUM NAD83

ZONE 18

ZUNITS METERS

UNITS METERS

PARAMETERS

END