java MakeSimpleContours > ! conts.gen
java MakeSimpleContours bound >! bound.gen
generate conts
input conts.gen
lines
q
generate bounds
input bound.gen
polygons
q
build conts lines
build bounds poly
AE;ec conts ;ef arc;sel all;densify 5;save;QUIT
tables
additem bounds.pat spacing 2 2 i
sel bounds.pat
calc spacing = 13 /* values imbedded in MakeSimpleContours.java
resel BOUNDS# < 12;calc spacing = 19
resel BOUNDS# < 10;calc spacing = 25
resel BOUNDS# < 8;calc spacing = 35
resel BOUNDS# < 6;calc spacing = 50
resel BOUNDS# < 4;calc spacing = 70
quit


&do ttol  &list 25 125 250
  &do it  &list 5 10 20 30 40 50 60
     &do enforce  &list on off
        &type elv%ttol%_%it%_%enforce%
        &s tol [calc %ttol% * 0.1]
        topogrid elv%ttol%_%it%_%enforce% 30
        contour conts conts-id
        datatype contour
        tolerance %tol%
        iterations %it%
        enforce %enforce%
        end
     &end
  &end
&end

GRID
setmask off
setcell 30
setwindow bounds
BOUNDmask = polygrid (bounds,spacing)
trueslope = atan(25.0 / boundmask) * deg
&s file_unit = [open autosummary openstatus -write]
&if %openstatus% ne 0 &then; &return openstatus was %openstatus%
&do ttol  &list 25 125 250
  &do it  &list 5 10 20 30 40 50 60
     &do enforce  &list on off
        setwindow  maxof
        slp%ttol%_%it%_%enforce% = slope(elv%ttol%_%it%_%enforce%,DEGREE) - trueslope
/*        asp_%ttol%_%it%_%enforce% = aspect(elv%ttol%_%it%_%enforce%)
        setwindow BOUNDmask
        slpdif = int(100.0 * (zonalmax(boundmask,slp%ttol%_%it%_%enforce%) - zonalmin(boundmask,slp%ttol%_%it%_%enforce%))+ 0.5)
        slpmean = int(100.0 * zonalmean(boundmask,slp%ttol%_%it%_%enforce%) + 0.5)
        slpstd = int(100.0 * zonalstd(boundmask,slp%ttol%_%it%_%enforce%) + 0.5)
/*        aspmin = ...
        &describe slpdif;&s meandif  %grd$mean%
        &describe slpstd;&s meanstd  %grd$mean%
        &type Write status was [write %file_unit% [quote %ttol%_%it%_%enforce% %meandif% %meanstd%]]
	stt%ttol%_%it%_%enforce% = combine(trueslope,slpdif,slpmean,slpstd)
        kill slpdif;kill slpstd;kill slpmean
     &end
  &end
&end
&ty close status = [close -all]
QUIT
tables
&do ttol  &list 25 125 250
  &do it  &list 5 10 20 30 40 50 60
     &do enforce  &list on off
       sel stt%ttol%_%it%_%enforce%.vat
       unload dmp%ttol%_%it%_%enforce%
     &end
  &end  
&end
setwindow  maxof
/* for web illustrations
setmask boundmask
tempg = ELV25_30_ON
slpx100p = int(100.0 * (slp25_30_ON + 10.0) + 0.5)


&GETPOINT;id conts arc %PNT$X% %PNT$Y%;cellv elv %PNT$X% %PNT$Y%