/home/pascal/depot/filedata/src/cloudsatfiledata.h

/***************************************************************************
 *   Copyright (C) 2005 by Nicolas PASCAL   *
 *   nicolas.pascal@icare.univ-lille1.fr   *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 *   This program is distributed in the hope that it will be useful,       *
 *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
 *   GNU General Public License for more details.                          *
 *                                                                         *
 *   You should have received a copy of the GNU General Public License     *
 *   along with this program; if not, write to the                         *
 *   Free Software Foundation, Inc.,                                       *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
 ***************************************************************************/
#ifndef CLOUDSATFILEDATA_H
#define CLOUDSATFILEDATA_H

#include "hdffiledata.h"
#include "pixel.h"
#include "satellitefiledata.h"

class CLOUDSATFileData : public SatelliteFileData, public HDFFileData {
    static const float colocation_tolerance;

    /***** CLOUDSAT constants *****/
    string latitude_name;
    string longitude_name;
    static const string time_name;
    static const string time_origin_name;

    /***** Product characteristics *****/
    string product;
    string version;
    int lat_lon_index_max;
    short orbit_nb;

    void init();
    void set_time_coverage();
    void parse_filename(const string & short_filename);
    bool check_filename(const string & short_filename);
    void set_lat_lon_index_max();
public:
    static const int nb_vbins=125;
    static const int nb_profile_max=40000;
    static const int nb_layers=5;
    CLOUDSATFileData(const string &name, const string &mode="r");
    virtual ~CLOUDSATFileData();
    void load_geolocation_data();
    void free_geolocation_data();

    const bool contain_location(const float &lat,const float &lon, const double &tolerance=0.01);
    const bool contain_data(const float &lat, const float &lon, const double &time, const double &colocation_tolerance=0.01) ;
    const float get_nearest_point_distance(const float &lat,const float &lon,const float coloc_tolerance=CLOUDSATFileData::colocation_tolerance);
    string get_product() const { return product; }
    string get_version() const { return version; }
     const bool get_index(const float &lat, const float& lon, int &nearest_pix_idx, const float colocation_tolerance=0.01);
     virtual void get_vindex(vector < vector < int > > &v_index, const float &lat, const float& lon,
                                            const float colocation_tolerance=0.01 ) {
        UnimplementedMethod e(__FILE__,__LINE__,__PRETTY_FUNCTION__);
        throw e;
     };

    const int get_nb_geo_points() { return lat_lon_index_max + 1; }
    virtual void load_v_pixel() ;
    virtual void close_data_file() {
        free_hdf_file();
    };
    virtual void open_data_file() {
        load_hdf_file();
    };
    virtual void get_pixel_coord ( const vector < int > & ipix, float &lat, float &lon, double &time ) {
        if (!is_geolocation_data_loaded())
            load_geolocation_data();
        lat  = lat_data  [ ipix[0] ];
        lon  = lon_data  [ ipix[0] ];
        time = time_data [ ipix[0] ];
    };
};

#endif