!
!   $Log: sofie_path.pct,v $
!   Revision 1.7  2008/03/23 02:03:26  deaver
!   removed call to extrapolate_pmc_extinction
!
!   Revision 1.6  2008/03/10 11:10:47  deaver
!   added aero_refband to call to extrapolate_pmc_extinction
!
!   Revision 1.5  2008/02/01 14:00:28  deaver
!   Change the APR arrays to only have the gases that you are including for a
!   given channel .
!
!   Revision 1.4  2008/01/29 14:14:00  deaver
!   Added za to call to calc_limb_refraction
!
!   Revision 1.3  2008/01/14 14:37:48  deaver
!   had a bug affecting eoi after extrapolate_ref_aerosol was added
!
!   Revision 1.2  2008/01/11 16:28:23  deaver
!   Added extrapolate_pmc_extinction to produce a new Extinction profile
!   for the current band using a factor multiplied by a reference band.
!
!   Revision 1.1  2007/12/19 18:17:56  deaver
!   Initial revision
!
!
EOL   = 
           a,   !  0
           b,   !  1    IF(tran_opt ) THEN
           c,   !  2       Call mp_mass for mega
           d,   !  3
                      ! ELSE
           e,   !  4      ! Call lp_mass for Linepak
	       f,   !  5   IF( DV difference Extinction channel) THEN
	       g,   !  6       call make_dv_extinction to setup for difference extinction
           a,   ! 7
           c,   ! 8  call mass again for new setup
	       h,   !  9  calculate refraction angles

           z;
!
EOI    =   0;
!
! 


a: incr: isc, 1;

b: if_then_else	:	tran_opt,2,4; 

c: mp_mass:
	ml,		ml,		mgs,		mgs,
	npt,		nl,		nl,		irfn,
	1,		nl,		ng,		idg,
	ng,		idg,		isc,		re,
	wave,		iplat,		za,		zapr,
	papr,		tapr,		gr_apr,		qapr,
	qgapr,		zt,		pt,		tt,
	gt,		tearth,		qmix,		qg,
	tm_tbl,		pout;

d: set_eoi: "",5;

e: lp_mass:
	ml,		ml,		mgs,		mgs,
	npt,		nl,		nl,		irfn,
	1,		nl,		ng,		idg,
	ng,		idg,		isc,		re,
	wave,		iplat,		za,		zapr,
	papr,		tapr,		gr_apr,		qapr,
	qgapr,		zt,		pt,		tt,
	gt,		tearth,		qmix,		qg,
	pout;



f: if_then_else: idiff_ext, 6 , 9 ;   !5

g: make_dv_extinction2:            !6
	iblk,	nch,
	idg,  	 ng,	
	nl, 	 ml,
	idg_b,	irad_b,	icor_b,
	itab,	mgs,	
	tm_ntb,	ntt,	nst,
	pout,	rlt,	rut,
	re,	 tearth,	 albedo,
	qmix,	qg,	
	tm_tbl,	tout,	sout,
	emr,    rmr,  esr,
	rsr,	bsr, zt,
    zapr, qapr, qgapr,
	igdo, ngdo, ichando, nchn,
    za, meas_signals ;

!a and c again   !7,8

h: calc_limb_refraction:   !9
	1,		nl,		pout,		appang, 
	trueang,	za;

!A: dump_initial: nl,	ng,	idg,	zapr,	papr, tapr,	qapr;
!B: dump_initial: nl,    ng,     idg,    zt,   pt, tt,     qmix;