!
!
!
EOL =    a,      !  0
         b,      !  1
         c,      !  2
	   M,N,               !call the tvib model   3,4
         d,      !  5
         e,      !  6
         f,      !  7
         g,      !  8
         h,      !  9
         i,      ! 10
         j,      ! 11
         k,      ! 12
         l,      ! 13
         m,      ! 14
         n,      ! 15
	 o,      ! 16
         p,      ! 17
         q,      ! 18
         r,      ! 19
         s,      ! 20
         t,      ! 21
         u,      ! 22
         z;
!
EOI =    0;
!
!

a: set_nvar: ifov_loop, 1;     		!  0  

b: set_var:				!1
	mgs,		ml,		ng,		nl,
	nreg,		nchn,		tfrac,		qfrac,
	papr,		tapr,		qmix,		vaprt,
	vapr;
!
! first_call is a flag that initilizes the convergence code
!
c: set_nvar: first_call, 1;   	! 2  Top of FOV Loop

!
!call the tvib model at the start of each pass
!
M: if_then_else: tm_inlte, 4, 5;           !3     
!N: nlteco2init      :tm_nlteco2init_eoi;     !4
N: nlteco2init;
!N: set_nvar: iray,5;


d: if_then_else: bottom_up,  6,  8;     !5
e: set_nvar:  iray, ner;     		!6
f: set_eoi: "",9    ;         		!7
g: set_nvar: iray, nsr;      		!8

h: set_nvar: num_iterations, 0;   	!9 

i: set_nvar: isl,iray;  			! 10
j: set_nvar: iel,iray;  			! 11

k: FORWARD;				! 12

l: derivative:				!13
	iray,		igdo,		ngdo,
	idg,		ng,		num_iterations,
	nl,		nchn,		itab,
	nreg,		ievent,		maxit,
	vapr,		vaprt,		tst,
	vary,		meas_signals,	
	evar,		zt,
	papr,		tapr,		zapr,
	qmix,		sim_signals,
	hydro_flag,	first_call,	converged,
	bottom_up,	nrays,		iscale_relax;	

m: if_then_else: hydro_flag, 15, 17   ;	! 14

n: hydro:  				!15
	nl, 		1, 		nreg, 
	nreg, 		re, 		xlat, 
	zapr, 		papr, 		tapr    ;

o: PATH_MASTER;				!16

p: if_then_else: converged, 18, 12;      !17

q: if_then_else: bottom_up,  19,  21;     ! 18

r: f_dse: iray, nsr, llll;		!19    !if (iray > nsr) return true.. decrement iray
s: if_then_else: llll, 9, 22;          !20   !if true   go to next layer else next fov loop

! Top town
t: ise: iray, ner, 9;			!21   increment iray 

! Next FOV LOOP
u: ise: ifov_loop, num_fov_loops, 2;    !22


!A: dump_initial: nl, 0, idg, zapr, papr, tapr, qmix,iblk;