;creates seltotal_h and seltotal_l, which contain the numbers of the ali/sar*****
;with high or low occupancy, respectively

x11=0
x12=0

do lb1 i=1,83
x50=x0

x20=(x50-1)*2+1
ud ic x20,x22,x55,x15,x16,x17,x18
picstat002

x20=x20+1
ud ic x20,x22,x55,x15,x16,x27,x18
picstat002

IF (x17.ge.x27) THEN
	x51=1
	x52=2
ELSE
	x51=2
	x52=1
ENDIF



ud n,x20
groups/grp{***x50}/sel{***x51}

do lb2 i=1,x20
x55=x0
x11=x11+1

ud ic x55,x21
groups/grp{***x50}/sel{***x51}

ud ic x21,x22
views/docview/view{***x50} 

sd x11,x22
seltotal_h

lb2
ud ice
groups/grp{***x50}/sel{***x51}
ud ice 
views/docview/view{***x50}

ud n,x20
groups/grp{***x50}/sel{***x52}

do lb3 i=1,x20
x55=x0
x12=x12+1

ud ic x55,x21
groups/grp{***x50}/sel{***x52}

ud ic x21,x22
views/docview/view{***x50}

sd x12,x22
seltotal_l

lb3
ud ice
groups/grp{***x50}/sel{***x52}
ud ice 
views/docview/view{***x50}

lb1

en d