c**d54    process first 1H-1H planes (real and imag in 15n)
;
cl
;
; nt=4, np=640, ni=60, ni2=16, ni3=14
;
;
;
def input fidc
ty -----------------------------------------
   ty Building matrix...
   bld ec283_4d_ccnoesy 4 512 128 64 64 0 y
   mat ec283_4d_ccnoesy w
   ty Matrix  opened and write-enabled...
   ty
;end
;
; 				D1 processing
: This data collected with array = phase, phase2, phase3
;
;
  esc out
  if &out ne 0 quit
	def phase0 110
	def phase1 0
 for cube 1 14
     eva icube (2*&cube)
     eva rcube (&icube-1)
     ty  Processing cubes &rcube and &icube 
 for tier 1 16
     eva itier (2*&tier)
     eva rtier (&itier-1)
     ty  Processing tiers &rtier and &itier $
     esc out
;     if &out ne 0 scram
     for row 1 60
       eva irow (2*&row)
       eva rrow (&irow-1)
       ;
       ; Process phase=1, phase2=1, phase3=1
       re &input
       cnv 1 50
       qss 320 90 1
       bc .2
       zf 512
       ft
       ph
       def datsiz 512
       ;pol 4
       ;abl 4 40
       red
       sto 0 &rrow &rtier &rcube
       ;
       ;process phase=1, phase2=1, phase3=2
       re &input
       cnv 1 50
       qss 320 90 1
       bc .2
       zf 512
       ft
       ph
       def datsiz 512
       ;pol 4
       ;abl 4 40
       red
       sto 0 &rrow &rtier &icube
       ;
       ; Process phase=1, phase2=2, phase3=1
       re &input
       cnv 1 50
       qss 320 90 1
       bc .2
       zf 512
       ft
       ph
       def datsiz 512
       ;pol 4
       ;abl 4 40
       red
       sto 0 &rrow &itier &rcube
       ;
       ; Process phase=1, phase2=2, phase3=2
       re &input
       cnv 1 50
       qss 320 90 1
       bc .2
       zf 512
       ft
       ph
       def datsiz 512
       ;pol 4
       ;abl 4 40
       red
       sto 0 &rrow &itier &icube
       ; Process phase=2, phase2=1, phase3=1
       re &input
       cnv 1 50
       qss 320 90 1
       bc .2
       zf 512
       ft
       ph
       def datsiz 512
       ;pol 4
       ;abl 4 40
       red
       sto 0 &irow &rtier &rcube
       ;
       ;process phase=2, phase2=1, phase3=2
       re &input
       cnv 1 50
       qss 320 90 1
       bc .2
       zf 512
       ft
       ph
       def datsiz 512
       ;pol 4
       ;abl 4 40
       red
       sto 0 &irow &rtier &icube
       ;
       ; Process phase=2, phase2=2, phase3=1
       re &input
       cnv 1 50
       qss 320 90 1
       bc .2
       zf 512
       ft
       ph
       def datsiz 512
       ;pol 4
       ;abl 4 40
       red
       sto 0 &irow &itier &rcube
       ;
       ; Process phase=2, phase2=2, phase3=2
       re &input
       cnv 1 50
       qss 320 90 1
       bc .2
       zf 512
       ft
       ph
       def datsiz 512
       ;pol 4
       ;abl 4 40
       red
       sto 0 &irow &itier &icube
     nex
   nex
   nex
   tim r t3

;end
	def phase0 0
	def phase1 0
;	bun 2
for tier 1 28 
for plane 1 32 
	for col 1 512
	  loa &col 0 &plane &tier
	esc out
	  if &out ne 0 quit
	  def datype 1
	  def datsiz 60
	  qsb 60 90 1
          bc .2
	  zf 128
gv 1 reel imagin
eva reel (&reel*0.5)
eva imagin (&imagin*0.5)
pv 1 &reel &imagin
	  ft
	  ph
	  red
	  rev
	  sto &col 0 &plane &tier
	ty col=&col plane=&plane tier=&tier$
	nex
nex
nex
;end
def phase0 90
def phase1 -180
    for plane 1 32
    for row 1 128
    for col 1 512
    loa &col &row 0 &plane
     def datype 1
     def datsiz 16
    bc .2
    qsb 16 90 1
    zf 64
     ft
     ph
     rev
     red
     sto &col &row 0 &plane
     ty  D3-plane=&plane $
    nex
    nex
    nex
def phase0 90
def phase1 -180
    for plane 1 64
    for row 1 128
    for col 1 512
    loa &col &row &plane 0
     def datype 1
     def datsiz 14
    bc .2
;    lpl 14 5 5 15 32
;    qsb 32 90 1
    qsb 14 90 1
    zf 64
     ft
     ph
     rev
     red
     sto &col &row &plane 0
     ty  D3-plane=&plane $
    nex
    nex
    nex
quit:
end