//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
// Copyright 2014 California Institute of Technology. ALL RIGHTS RESERVED.
// 
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
// 
// http://www.apache.org/licenses/LICENSE-2.0
// 
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.
// 
// United States Government Sponsorship acknowledged. This software is subject to
// U.S. export control laws and regulations and has been classified as 'EAR99 NLR'
// (No [Export] License Required except when exporting to an embargoed country,
// end user, or in support of a prohibited end use). By downloading this software,
// the user agrees to comply with all applicable U.S. export laws and regulations.
// The user has the responsibility to obtain export licenses, or other export
// authority as may be required before exporting this software to any 'EAR99'
// embargoed foreign country or citizen of those countries.
//
// Author: Piyush Agram
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~





#ifndef WATERMASK_H
#define WATERMASK_H

#ifndef ERR_MESSAGE
#define ERR_MESSAGE cout << "Error in file " << __FILE__ << " at line " << __LINE__  << " Exiting" <<  endl; exit(1);
#endif


class Polygon
{
    private:
        int npoints;
        double **xy;
    public:
        Polygon();
        ~Polygon();
        void allocate(int n);
        void setPoint(int i, double x, double y);
        int isInside(double testx, double testy);
        void print();
};


class WaterBody
{
    private:
        Polygon *shapes;
        int nshapes;
        double x0, y0;
        double dx, dy;
        int width, height;

    public:
        WaterBody(int n);
        ~WaterBody();
        int isWater(double x, double y);
        void allocate(int ind, int n);
        void setShapeData(int ind, int i, double x, double y);
        void setDimensions(int ww, int ll);
        void setTopLeft(double xx, double yy);
        void setSpacing(double ddx, double ddy);
        void fillGrid(char* filename);
        void makemask(char* lonfile, char* latfile, char* outfile);
        void printShape(int i);
};

#endif //WATERMASK_H