!
!  $Log: sofie_forward.pct,v $
!  Revision 1.6  2008/06/26 19:07:15  deaver
!  Added the NonLTE forward model.
!
!  Revision 1.5  2008/02/01 13:53:06  deaver
!  FOV offsets added to fov convolution routine in order to offset the lockdown position.
!
!  Revision 1.4  2008/01/28 13:07:28  deaver
!  Removed some of the old commented routines from HALOE days
!
!  Revision 1.3  2008/01/18 18:22:45  deaver
!  Added extend_tran_profile to  the EOL
!
!  Revision 1.2  2008/01/16 15:47:08  deaver
!  Added fov atmosphere solar limb darkening curve convolution
!
!  Revision 1.1  2007/12/19 18:17:56  deaver
!  Initial revision
!
!  Revision 1.1  2007/06/28 20:16:45  deaver
!  Initial Checkin
!
EOL =   H,     ! 0
        I,     ! 1
        J,     ! 2
        K,     ! 3
        a,     ! 4    !Determine which model to use
         b,    ! 5      !Call to MEGA
         c,    !  6
         d,    !  7    !Call to MLRAY
         A,    !  8
         e,    ! 9    !Convolve the FOV with the extended tranmission profile
         f,    ! 10   !create simulated signals  (apply non-linearity)
         z;
!
EOI =    0;
!
!
H: if_then_else: tm_inlte, 1,4;
I: incr: isc,1;

J:limbcalc    : ml, tm_ml, mgs, tm_mgs, tm_mgb, tm_mt, tm_mdet,
                tm_nch, tm_ntb, tm_npt,  tm_ntt,  tm_nst, nl, nl,
                ng, idg , ng, idg, tm_ngb, tm_idgb, re, wave, isc, 
		        tm_irfn, tm_irad, tm_icor, tm_itab,
                zt, pt, tt, gr, qmix, qg,              ! input atmosphere
		        za, tm_zt0, tm_pt0, tm_tt0, tm_gr0, tm_qmix, tm_qg, 
		        tm_qmixb, tm_qgb, tm_rlt, tm_rut, tearth,
                emr, rmr, esr, rsr, bsr,
                tm_ndx1, tm_ndx2,
                tm_tbl, pout, sout, tout, isl, iel, tran;


K: set_eoi: "",8 ;

a: if_then_else	:tran_opt, 5, 7;        	!0

b: sofie_mega_model:                            ! 1
	iblk,		nch,
	idg,		ng,
	isl,		iel,		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,
	tran, zt, tapr;   !,   zt, pt, tt,  xlat,za , iexo ;
						
c: set_eoi: "",8;				!2

d: sofie_lbl_model:                             !3
	lfilt,		lhvy,		loutput,	FALSE,
	iblk,		nch,		isc,
	idg,		mgs,		ng,
	ml,		nl,		isl,		iel,
	npt,		pout,		rlt,		rut,
	re,		tearth,		albedo,
	zt,		pt,		tt,
	qmix,		qg,		tran, xlat,za;

!Add fov convolution

A: extend_tran_profile:
		iblk,		nch,
		isl,		iel,		nl,
		nsr,		ner,		bottom_up,
		za,	     	tran,		fovflag,
		ifov_loop;


e:  fov_atm_convolution:
		nch,		iblk,	isl,	iel,
		nl,		lock,	appang,	trueang,
		tran,		conv_tran,	simexo,
		fovflag, idiff_ext, za ; 


f: create_sim_signal:
	ichando,	nchn,		iblk,		nch,
	isl,		iel,	   	conv_tran, 
    sim_signals,    idiff_ext;