!
!
EOL =    a,  !  0 
         b,  !  1
         c,  !  2
         d,  !  3
         e,  !  4
         f,  !  5
         g,  !  6 
         h,  !  7
  	 i,  !  8
	 j,  !  9
         k,  ! 10
        B, ! 11
         l,  ! 12
        A, !13
         m,  ! 14
	 n,  ! 15
	 o,  ! 16
         z;
!
EOI =    0;
!
!

a: set_nvar: ievent_no, istart_event;   ! Set ievent_no equal to starting event number

b: read_l1_event:
	ievent_no, 	re,		    xlat, 		xlong,
	dateloc,	nl_hi,		iorbit,		no_data,
	exolock,	l1_ver,		zoffsun,	rse,
	fovflag,        poorman;


c: if_then_else:	no_data, 16, 3;   ! if no_data goto the next event

d: init_aerosols;  !: InputOffsets, NonLinCoeff;

e: get_l1_zpt:
	ztm,		ptm,		ttm,		ntm,
	mr,		grm,		 pout, npt, zoff;

f: readfasatm2: 
	latmos,		atmos,		re,		xlat,
	mr,		ntm,		ngg,		nam,
	igs,		ztm,		ptm,		ttm,
	grm,		qm,		qgm,		FALSE,    !this false is for IHYDRO
	inatmos, 	TRUE,  		FALSE, 30.0;     !True for PPMV conversion FALSE for reading in ZPT

g: get_waccmc2_climatology:   use_waccm, 	waccm_dir,	ptm, ztm,	qm,
	igs,	ntm,	ngg;

h:  get_msis_data: use_msis, use_saber, O_x2, naltsab, saberO_alt, saberO_NH, saberO_SH, ztm, qm, ntm, igs, ngg;

!
! create a high resolution library that is used to initialize retrievals
!  This is zr, pr, tr, gr, q and qgg,  which has nr levels containing ngg gases 
!  with GAS IDs   igs
!

i: create_z_grid:  za,  nr, mr, 150.0 , 2.9, 0.2 ;

j: set_nvar: isc, 1;

k: lp_mass:
        mr,             ml,             mgg,            mgg,
        npt,            nr,             ntm,           FALSE,    ! Not using refraction here
        1,              nr,             ngg,            igs,
        ngg,            igs,            isc,            re,
        wave,           iplat,          za,             ztm,
        ptm,            ttm,            grm,             qm,
        qgm,            zr,             pr,             tr,
        gr,             tearth,         q,              qgg,
        pout;
!
! finished with climatology
!
B: set_nvar: pmcflag, FALSE;

l: AUX_DATA_OUTPUT;

A: CORRECT_PMC;
!A: set_eoi:"",14;

m: CHANNEL_LOOP;

n: insert_l2_summary: l1_ver, l2_ver, exolock, pmcflag, cycle, z_abelout;

o: ise: ievent_no, iend_event, 1;