Wed Jan 28 14:39:30 2004
e100n40.shade
mbasa
mbasa
cell


generated by r.mapcalc
eval(x = (e100n40[-1,-1] + 2 * e100n40[0,-1] + e100n40[1,-1] - 
e100n40[-1,1] - 2 * e100n40[0,1] - e100n40[1,1]) / (8.000000 * 
ewres()), y = (e100n40[-1,-1] + 2 * e100n40[-1,0] + e100n40[-1,1] - 
e100n40[1,-1] - 2 * e100n40[1,0] - e100n40[1,1]) / (8.000000 * 
nsres()), slope = 90.000000 - atan(sqrt(x * x + y * y)), a = 
round(atan(x, y)), a = if(isnull(a), 1, a), aspect = if(x != 0 || y 
!= 0, if(a, a, 360.000000)), cang = sin(45) * sin(slope) + cos(45) * 
cos(slope) * cos(180 - aspect), if(cang < 0.000000, 0.000000, 
100.000000 * cang), if(isnull(cang), null(), 100.000000 * cang))
