c**d54    process first 1H-1H planes (real and imag in 15n)
;
cl
cmx
;
;ran on glacier
;processed with grad_sort_nd_mak origfid fid 32 50 0 1024
;
;
def input fidc
ty -----------------------------------------
   ty Building matrix...
   bld ir24_hncacb_baker_nolp 3 512 128 128 0 y
   mat ir24_hncacb_baker_nolp w
   ty Matrix opened and write-enabled...
   ty
;
; 				D1 processing
;
;
  esc out
  if &out ne 0 quit
	def phase0 65
	def phase1 0
 for tier 1 32
     eva itier (2*&tier)
     eva rtier (&itier-1)
     ty  Processing tiers &rtier and &itier $
     esc out
;     if &out ne 0 scram
     for row 1 50
       eva irow (2*&row)
       eva rrow (&irow-1)
       ;
       ; Process phase2=1, phase=1
       re &input
       cnv 1 25
       qss 512 90 1
       zf 1024
       ft
       ph
       def datsiz 512
       ;pol 4
       ;abl 4 40
       red
       sto 0 &rrow &rtier
       ;
       ;process phase2=1, phase=2
       re &input
       cnv 1 25
       qss 512 90 1
       zf 1024
       ft
       ph
       def datsiz 512
       ;pol 4
       ;abl 4 40
       red
       sto 0 &irow &rtier
       ;
       ; Process phase2=2, phase=1
       re &input
       cnv 1 25
       qss 512 90 1
       zf 1024
       ft
       ph
       def datsiz 512
       ;pol 4
       ;abl 4 40
       red
       sto 0 &rrow &itier
       ;
       ; Process phase2=2, phase=2
       re &input
       cnv 1 25
       qss 512 90 1
       zf 1024
       ft
       ph
       def datsiz 512
       ;pol 4
       ;abl 4 40
       red
       sto 0 &irow &itier
     nex
   nex
   tim r t3

;end
	def phase0 0
	def phase1 0
;	bun 2
for plane 1 64 
	for col 1 512
	  loa &col 0 &plane
	esc out
	  if &out ne 0 quit
	  def datype 1
	  def datsiz 50
          lpl 50 8 8 51 96
	  qsb 50 90 1
	  zf 128
	  ft
	  ph
	  red
          ; rev ; needed on glacier, but not on baker
          csl 0
	  sto &col 0 &plane
	ty col=&col plane=&plane$
	nex
nex
;end
bun 0
   bun 3
   ty There are &vector vectors to process.
   tim z
def phase0 0
def phase1 0
   for ivec 1 &vector
     lwb
     def datype 1
     def datsiz 32
;
   def datsiz 32 ;=27
   eva datso (&datsiz) ;=27
   eva shift (&datsiz-1) ;=26
   eva fpnt (&datsiz+1) ;=28
   eva lpnt (&shift+&datsiz) ;=53
   eva lpf (&lpnt+1) ;=54
   eva lpl (&lpnt+&datso) ;=80
   eva hafpnt (nint(&lpnt/2)) ;=27
;
; conduct mirror image LP
   ;  def datsiz &lpnt
   ;  srv &fpnt &lpnt 0 0
   ;  shr &shift
  ;   stb 1
  ;   shl &shift
  ;   def datsiz &datso
  ;   cnj
  ;   rev
  ;   def datsiz &lpnt
  ;   srv &datso &lpnt 0 0
  ;   adb 1
  ;   ldb 1
     ;
     ;The following LP setup looks good but defies std setup...
  ;   lpl &lpnt 4 4 &lpf &lpl
  ;   shl &shift
  ;;   eva datsiz (&datso+&datso)
  ;   ss &datsiz 90
  ;   zf &d3size
  qsb 32 90 1
  zf 128
     ft
     ph
     rev
     red
     swb
     ty  D3-vector=&ivec $
   next
   tim r t2
   bun 0
quit:
end