/**

    @class TruncRefracIn.cpp
    @author Brian Magill
    @datecreated 11/01/2007

    $Date:$
    $Revision:$

    @copyright
     (©) 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.

    @brief: Removes the missing values from the solar extents and corresponding time         

*/
//----------------------------------------------------------------------
//

#include <iostream>
#include <vector>
#include <cassert>
#include <cmath>
#include "TruncRefracIn.h"
#include "Vec2Val.hpp"

using namespace std;

void TruncRefracIn::operator()(std::valarray<double> const &InTime,
                    std::valarray<double> const &InValue,
                    std::valarray<double> &TruncTime,
                    std::valarray<double> &TruncValue) const   
{                                    
vector<double> workTime;
vector<double> workArray;
unsigned long indx;

assert(InTime.size() == InValue.size());

for(indx = 0; indx < InValue.size(); indx++)
    {
    if (startTime <= InTime[indx] && InTime[indx] <= endTime)
        {
        workTime.push_back(InTime[indx]);
        workArray.push_back(InValue[indx]);
        }
    }
    
Vec2Val(workTime, TruncTime);
Vec2Val(workArray, TruncValue);
}