//
//  $Id$
//-----------------------------------------------------------------------
//
//            (c) Copyright 2006 by GATS, Inc.,
//     11864 Canon Blvd, Suite 101, Newport News VA 23606
//
//  All Rights Reserved. No part of this software or publication may be
//  reproduced, stored in a retrieval system, or transmitted, in any form
//  or by any means, electronic, mechanical, photocopying, recording, or
//  otherwise without the prior written permission of GATS, Inc.
//
//-----------------------------------------------------------------------
//
// Module:	RefractionAngleTPRetrieval.cpp
//
// Author:	John Burton	<john@gats-inc.com>
//
// Date:  	Thu May  4 17:01:31 2006
//
//-----------------------------------------------------------------------
//
// Modification History:
//
//     $Log$
//
//-----------------------------------------------------------------------
//
//-----------------------------------------------------------------------
// Include Files:
//-----------------------------------------------------------------------
//

#include "ConfigFile.h"
#include "EventVar.h"
#include "EventPopul.hpp"
#include "RefractionAngleTPRetrieval.h"
#include <iostream>
#include <exception>
//
//-----------------------------------------------------------------------
// Defines and Macros:
//-----------------------------------------------------------------------
//



//
//-----------------------------------------------------------------------
// Global Variables:
//-----------------------------------------------------------------------
//



//
//-----------------------------------------------------------------------
// Utility Routines:
//-----------------------------------------------------------------------
//

int RefractionAngleTPRetrieval(Event& L0, Event& L1, Event& Tmp, Event& SD, ConfigFile& cf)
{

//  Currently this is a stub

int retCode = 0;
EventVar <double> zVarArray; 
EventVar <double> pVarArray; 
EventVar <double> tVarArray; 
EventVar <double> workArray;  
EventVar <double> angVarArray(zVarArray.size()); 
EventVar <double> errVarArray(zVarArray.size()); 

try {    
    L1.getEventVar("Z_Merged", zVarArray);
    L1.getEventVar("P_Merged", pVarArray);
    L1.getEventVar("T_Merged", tVarArray);

    angVarArray = zVarArray*0.001; 
    errVarArray = zVarArray*0.0; 

    EventPopul(zVarArray, "L1_RefractAlt", L1);
    EventPopul(pVarArray, "L1_RefractPress", L1);
    EventPopul(tVarArray, "L1_RefractTem", L1);
    EventPopul(angVarArray, "L1_RefracRetrievedAngles", L1);
    EventPopul(errVarArray, "L1_RefractRetrievedError", L1);

    std::cerr << "In RefractionAngleTPRetrieval \n";

    } catch (exception &ex) {
    
    cerr << "Exception occurred in RefractionAngleTPRetrieval: ";
    cerr << ex.what() << endl;
    
    retCode = -100;
    }    
    
    return retCode;
}