[ncl-talk] Slope and aspect

[Michael Notaro]

Guess I may have answered my own question about aspect....seems to work.


; a b c

; d e f

; g h i


; delta = 800.


a=1000.

b=800.

c=600.

d=1200.

e=1000.

f=800.

g=1400.

h=1200.

i=1000.


dzdx=((a+2.*d+g)-(c+2.*f+i))/(8.*800.)

dzdy=((g+2.*h+i)-(a+2.*b+c))/(8.*800.)


aspect=atan2(dzdx,dzdy)*45./atan(1.0)


print(aspect)



Michael Notaro
Associate Director
Nelson Institute Center for Climatic Research
University of Wisconsin-Madison
Phone: (608) 261-1503
Email: mnotaro at wisc.edu


________________________________
From: Michael Notaro
Sent: Friday, January 19, 2018 3:16 PM
To: ncl-talk at ucar.edu
Subject: Slope and aspect


Is there a way to compute slope and aspect

in NCL given an equal area netcdf dataset of elevation (m)?
Michael


Michael Notaro
Associate Director
Nelson Institute Center for Climatic Research
University of Wisconsin-Madison
Phone: (608) 261-1503
Email: mnotaro at wisc.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180119/4315c2fe/attachment.html>