Changeset 5605

Timestamp:

Aug 5, 2008, 11:12:01 AM (16 years ago)

Author:

jack

Message:

Extract the height data from the NetCDF file.

Location:

anuga_work/development/anugavis/src

Files:

• height_quantity.cc (modified) (2 diffs)
• height_quantity.hh (modified) (1 diff)

anuga_work/development/anugavis/src/height_quantity.cc

@@ -14 +14 @@
 #  include <OpenGL/glu.h>
 #endif
+#include <netcdf.h>
 #include "height_quantity.hh"
 
@@ -46 +47 @@
 
 void HeightQuantity::set_SWWFile(const shared_ptr<SWWFile> &file){
+  int ncstatus;
   this->sww_file = file;
   this->points = \
     shared_array<float>(new float[this->sww_file->number_of_points]);
   this->dynamic = this->sww_file->nc_inq_varndims_by_name(this->name) == 2;
+  if((ncstatus = nc_inq_varid(this->sww_file->netcdf_id, this->name.c_str(),
+                              &this->varid)) != NC_NOERR)
+    throw nc_strerror(ncstatus);
 }
 
 void HeightQuantity::compile(int frame){
-  
+  size_t start_1[1] = {0};
+  size_t start_2[2] = {0, 0};
+  size_t span_1[1];
+  size_t span_2[2] = {1, 0};
+  int ncstatus;
+  if(this->dynamic){
+    span_2[1] = this->sww_file->number_of_points;
+    ncstatus = nc_get_vara_float(this->sww_file->netcdf_id, this->varid,
+                                 start_2, span_2, this->points.get());
+  }else{
+    span_1[0] = this->sww_file->number_of_points;
+    ncstatus = nc_get_vara_float(this->sww_file->netcdf_id, this->varid,
+                                 start_1, span_1, this->points.get());
+  }
+  if(ncstatus != NC_NOERR) throw nc_strerror(ncstatus);
+  /* TODO: Compile the triangles into an OGL displaylist. */
 }

anuga_work/development/anugavis/src/height_quantity.hh

@@ -46 +46 @@
   int frame_number; /* The most recently drawn frame (i.e., the one
                        cached in the display list. */
+  int varid;
   bool dynamic;
   GLuint display_list;