Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
[mesh] update to MDAL 0.0.9 (NetCDF Support)
- Loading branch information
1 parent
e25d634
commit 9acc5c1
Showing
6 changed files
with
155 additions
and
18 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,79 @@ | ||
| /* | ||
| MDAL - Mesh Data Abstraction Library (MIT License) | ||
| Copyright (C) 2018 Peter Petrik (zilolv at gmail dot com) | ||
| */ | ||
|
|
||
| #include "mdal_gdal_netcdf.hpp" | ||
| #include "mdal_utils.hpp" | ||
|
|
||
| MDAL::LoaderGdalNetCDF::LoaderGdalNetCDF( const std::string &netCDFFile ) | ||
| : MDAL::LoaderGdal( netCDFFile, "GRIB" ) | ||
| , mTimeDiv( 1.0 ) | ||
| { | ||
| } | ||
|
|
||
| std::string MDAL::LoaderGdalNetCDF::GDALFileName( const std::string &fileName ) | ||
| { | ||
| #ifdef WIN32 | ||
| // Force usage of the predefined GDAL driver | ||
| // http://gis.stackexchange.com/a/179167 | ||
| // on Windows, HDF5 driver is checked before NETCDF driver in GDAL | ||
| return "NETCDF:\"" + fileName + "\""; | ||
| #else | ||
| return fileName; | ||
| #endif | ||
| } | ||
|
|
||
| bool MDAL::LoaderGdalNetCDF::parseBandInfo( const MDAL::GdalDataset *cfGDALDataset, const MDAL::LoaderGdal::metadata_hash &metadata, std::string &band_name, double *time, bool *is_vector, bool *is_x ) | ||
| { | ||
| MDAL_UNUSED( cfGDALDataset ); | ||
|
|
||
| metadata_hash::const_iterator iter; | ||
|
|
||
| iter = metadata.find( "netcdf_dim_time" ); | ||
| if ( iter == metadata.end() ) return true; //FAILURE, skip no-time bands | ||
| *time = parseMetadataTime( iter->second ) / mTimeDiv; | ||
|
|
||
| // NAME | ||
| iter = metadata.find( "long_name" ); | ||
| if ( iter == metadata.end() ) | ||
| { | ||
| iter = metadata.find( "netcdf_varname" ); | ||
| if ( iter == metadata.end() ) return true; //FAILURE, should be always present | ||
| band_name = iter->second; | ||
| } | ||
| else | ||
| { | ||
| band_name = iter->second; | ||
| } | ||
|
|
||
| // Loop throught all additional dimensions but time | ||
| for ( iter = metadata.begin(); iter != metadata.end(); ++iter ) | ||
| { | ||
| std::string key = iter->first; | ||
| if ( MDAL::contains( key, "netcdf_dim_" ) ) | ||
| { | ||
| key = MDAL::replace( key, "netcdf_dim_", "" ); | ||
| if ( key != "time" ) | ||
| { | ||
| band_name += "_" + key + ":" + iter->second; | ||
| } | ||
| } | ||
| } | ||
|
|
||
| // Parse X, Y components if present | ||
| parseBandIsVector( band_name, is_vector, is_x ); | ||
|
|
||
| return false; // SUCCESS | ||
| } | ||
|
|
||
| void MDAL::LoaderGdalNetCDF::parseGlobals( const MDAL::LoaderGdal::metadata_hash &metadata ) | ||
| { | ||
| metadata_hash::const_iterator iter = metadata.find( "time#units" ); | ||
| if ( iter != metadata.end() ) | ||
| { | ||
| std::string units = iter->second; | ||
| mTimeDiv = MDAL::parseTimeUnits( units ); | ||
| // TODO store reference time from iter->second too, see crayfish_netcdf.cpp | ||
| } | ||
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,34 @@ | ||
| /* | ||
| MDAL - Mesh Data Abstraction Library (MIT License) | ||
| Copyright (C) 2018 Peter Petrik (zilolv at gmail dot com) | ||
| */ | ||
|
|
||
| #ifndef MDAL_GDAL_NETCDF_HPP | ||
| #define MDAL_GDAL_NETCDF_HPP | ||
|
|
||
| #include "mdal_gdal.hpp" | ||
| #include "mdal_data_model.hpp" | ||
| #include "mdal.h" | ||
| #include <string> | ||
|
|
||
| namespace MDAL | ||
| { | ||
|
|
||
| class LoaderGdalNetCDF: public LoaderGdal | ||
| { | ||
| public: | ||
| LoaderGdalNetCDF( const std::string &netCDFFile ); | ||
| private: | ||
| std::string GDALFileName( const std::string &fileName ) override; | ||
| bool parseBandInfo( const MDAL::GdalDataset *cfGDALDataset, | ||
| const metadata_hash &metadata, std::string &band_name, | ||
| double *time, bool *is_vector, bool *is_x | ||
| ) override; | ||
| void parseGlobals( const metadata_hash &metadata ) override; | ||
|
|
||
| //! delimiter to get time in hours | ||
| double mTimeDiv; | ||
| }; | ||
|
|
||
| } // namespace MDAL | ||
| #endif // MDAL_GDAL_NETCDF_HPP |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters