Changeset 7 for trunk/include/ncutil.cpp

Timestamp:

03/07/11 11:49:01 (13 years ago)

Author:

jbrlod

Message:

rst write

File:

• trunk/include/ncutil.cpp (modified) (5 diffs)

trunk/include/ncutil.cpp

@@ -50 +50 @@
   return(var_id);
 }
- 
+
+void voncwrite(int t,int znbval,int ynbval, int xnbval, char const *var_name,int ncid,int dimids[], YREAL vect[])
+{ //Ecrrire le volume 3D au temps t
+  //le parametre t commence de 0.
+  // ecrit la variable var_name aux points (t, 1:znbval, 1:ynbval, 1:xnbval)
+  int status;
+
+  int var_id;
+  int local_dimids[4];
+  static size_t start[] = {0,0,0,0}; //3D+t
+  static size_t count[] = {1,1,1,1}; //3D+t
+  count[0]= (size_t) t+1;
+  count[1]= (size_t) znbval;
+  count[2]= (size_t) ynbval;
+  count[3]= (size_t) xnbval;
+  start[0]=t; //temps commence de 0;
+
+  local_dimids[0]=dimids[0]; //t
+  local_dimids[1]=dimids[1]; //z
+  local_dimids[2]=dimids[2]; //y
+  local_dimids[3]=dimids[3]; //x
+
+  //Define the netCDF variable
+  status=nc_redef(ncid);
+if (status != NC_NOERR)
+    {
+      handle_error (status, "Erreur change to def mode NetCDF ::voncwrite");
+    }
+ 
+ status=nc_def_var(ncid,var_name,NC_REAL,4,local_dimids,&var_id);
+ 
+  if (status != NC_NOERR)
+    {
+      printf("variable %s:",var_name);
+      handle_error (status, "Erreur define NetCDF variable ::voncwrite");
+    }
+ nc_enddef(ncid);
+  status=nc_put_vara_real(ncid,var_id,start,count,vect);
+ if (status != NC_NOERR)
+    {
+      handle_error (status, "Erreur writing NetCDF variable ::voncwrite");
+    }
+}
+
 void vonclire(int t,int znbval, int ynbval, int xnbval, int ncid, int varid, double vect[])
 {// lire volume 3D au temps t
@@ -76 +119 @@
 }
 
+void soncwrite(int t,int ynbval, int xnbval, char const *var_name,int ncid,int dimids[], YREAL vect[])
+{ //Ecrrire le volume 2D au temps t
+  //le parametre t commence de 0.
+  // ecrit la variable var_name aux points (t, 1:ynbval, 1:xnbval)
+  int status;
+
+  int var_id;
+  int local_dimids[3];
+  static size_t start[] = {0,0,0}; //2D+t
+  static size_t count[] = {1,1,1}; //2D+t
+  count[0]= (size_t) t+1;
+  count[1]= (size_t) ynbval;
+  count[2]= (size_t) xnbval;
+  start[0]=t; //temps commence de 0;
+
+  local_dimids[0]=dimids[0]; //t
+  local_dimids[1]=dimids[2]; //y
+  local_dimids[2]=dimids[3]; //x
+
+  //Define the netCDF variable
+  status=nc_redef(ncid);
+if (status != NC_NOERR)
+    {
+      handle_error (status, "Erreur change to def mode NetCDF ::voncwrite");
+    }
+ 
+ status=nc_def_var(ncid,var_name,NC_REAL,3,local_dimids,&var_id);
+ 
+  if (status != NC_NOERR)
+    {
+      printf("variable %s:",var_name);
+      handle_error (status, "Erreur define NetCDF variable ::voncwrite");
+    }
+ nc_enddef(ncid);
+  status=nc_put_vara_real(ncid,var_id,start,count,vect);
+ if (status != NC_NOERR)
+    {
+      handle_error (status, "Erreur writing NetCDF variable ::voncwrite");
+    }
+}
+
+
 void sonclire(int t,int ynbval, int xnbval, int ncid, int varid, double vect[])
 {// lire surface 2D au temps t
@@ -99 +184 @@
     }
 }
+
 void sonclire(int t,int ynbval, int xnbval, int ncid, int varid, float vect[])
 {// lire surface 2D au temps t
@@ -123 +209 @@
 }
 
+void snclire(int ynbval, int xnbval, int ncid, int varid, double vect[])
+{// lire surface 2D 
+  // (ynbval,xnbval,ncid, varid): 
+  // lit la variable varid aux points (1:ynbval,1:xnbval) et
+  // renvoie un vecteur de dim ynbval*xnbval
+ 
+  int           status;
+  static size_t   start[] = { 0, 0 };  // 2D
+  static size_t   count[] = { 1, 1 }; // 2D
+  count[0] =(size_t)  ynbval;
+  count[1] =(size_t)  xnbval;
+  
+  // Read an array of values
+  status = nc_get_vara_double (ncid, varid, start, count, vect);
+  
+ if (status != NC_NOERR)
+    {
+      handle_error (status, "Erreur lecture fichier NetCDF::snclire");
+    }
+}
+
+void snclire(int ynbval, int xnbval, int ncid, int varid, float vect[])
+{// lire surface 2D 
+  // (ynbval,xnbval,ncid, varid): 
+  // lit la variable varid aux points (1:ynbval,1:xnbval) et
+  // renvoie un vecteur de dim ynbval*xnbval
+ 
+  int           status;
+  static size_t   start[] = { 0, 0 };  // 2D
+  static size_t   count[] = { 1, 1 }; // 2D
+  count[0] =(size_t)  ynbval;
+  count[1] =(size_t)  xnbval;
+  
+  // Read an array of values
+  status = nc_get_vara_float (ncid, varid, start, count, vect);
+  
+ if (status != NC_NOERR)
+    {
+      handle_error (status, "Erreur lecture fichier NetCDF::snclire");
+    }
+}
+
+void zoncwrite(int t,int znbval,char const *var_name,int ncid,int dimids[], YREAL vect[])
+{ //Ecrrire la variable 1D au temps t
+  //le parametre t commence de 0.
+  // ecrit la variable var_name aux points (t, 1:znbval)
+  int status;
+
+  int var_id;
+  int local_dimids[2];
+  static size_t start[] = {0,0}; //1D+t
+  static size_t count[] = {1,1}; //1D+t
+  count[0]= (size_t) t+1;
+  count[1]= (size_t) znbval;
+  
+  start[0]=t; //temps commence de 0;
+
+  local_dimids[0]=dimids[0]; //t
+  local_dimids[1]=dimids[1]; //z
+  
+
+  //Define the netCDF variable
+  status=nc_redef(ncid);
+if (status != NC_NOERR)
+    {
+      handle_error (status, "Erreur change to def mode NetCDF ::voncwrite");
+    }
+ 
+ status=nc_def_var(ncid,var_name,NC_REAL,2,local_dimids,&var_id);
+ 
+  if (status != NC_NOERR)
+    {
+      printf("variable %s:",var_name);
+      handle_error (status, "Erreur define NetCDF variable ::voncwrite");
+    }
+ nc_enddef(ncid);
+  status=nc_put_vara_real(ncid,var_id,start,count,vect);
+ if (status != NC_NOERR)
+    {
+      handle_error (status, "Erreur writing NetCDF variable ::voncwrite");
+    }
+}
+
+
 void zonclire(int t,int znbval, int ncid, int varid, double vect[])
 {// lire surface 1D au temps t
@@ -218 +388 @@
     }
 }
+
+void znclire(int znbval, int ncid, int varid, double vect[])
+{// lire surface 1D 
+  // (znbval,ncid, varid): 
+  // lit la variable varid aux points (1:znbval) et
+  // renvoie un vecteur de dim znbval
+ 
+  int           status;
+  static size_t   start[] = { 0 };  //1D
+  static size_t   count[] = { 1 }; // 1D
+  count[0] =(size_t)  znbval;
+ 
+  
+  // Read an array of values
+  status = nc_get_vara_double (ncid, varid, start, count, vect);
+  
+ if (status != NC_NOERR)
+    {
+      handle_error (status, "Erreur lecture fichier NetCDF::zonclire");
+    }
+}
+
+
+void znclire(int znbval, int ncid, int varid, float vect[])
+{// lire surface 1D 
+  // (znbval,ncid, varid): 
+  // lit la variable varid aux points (1:znbval) et
+  // renvoie un vecteur de dim znbval
+ 
+  int           status;
+  static size_t   start[] = { 0 };  //1D
+  static size_t   count[] = { 1 }; // 1D
+  count[0] =(size_t)  znbval;
+ 
+  
+  // Read an array of values
+  status = nc_get_vara_float (ncid, varid, start, count, vect);
+  
+ if (status != NC_NOERR)
+    {
+      handle_error (status, "Erreur lecture fichier NetCDF::zonclire");
+    }
+}
+
+void sncwrite(int ynbval, int xnbval, char const *var_name,int ncid,int dimids[], YREAL vect[])
+{ //Ecrrire le volume 2D
+  // ecrit la variable var_name aux points (1:ynbval, 1:xnbval)
+  int status;
+
+  int var_id;
+  int local_dimids[2];
+  static size_t start[] = {0,0}; //2D
+  static size_t count[] = {1,1}; //2D
+  count[0]= (size_t) ynbval;
+  count[1]= (size_t) xnbval;
+ 
+  local_dimids[0]=dimids[2]; //y
+  local_dimids[1]=dimids[3]; //x
+
+  //Define the netCDF variable
+  status=nc_redef(ncid);
+if (status != NC_NOERR)
+    {
+      handle_error (status, "Erreur change to def mode NetCDF ::voncwrite");
+    }
+ 
+ status=nc_def_var(ncid,var_name,NC_REAL,2,local_dimids,&var_id);
+ 
+  if (status != NC_NOERR)
+    {
+      printf("variable %s:",var_name);
+      handle_error (status, "Erreur define NetCDF variable ::voncwrite");
+    }
+ nc_enddef(ncid);
+  status=nc_put_vara_real(ncid,var_id,start,count,vect);
+ if (status != NC_NOERR)
+    {
+      handle_error (status, "Erreur writing NetCDF variable ::voncwrite");
+    }
+}
+
+
+void zncwrite(int znbval,char const *var_name,int ncid,int dimids[], YREAL vect[])
+{ //Ecrrire la variable 1D 
+  // ecrit la variable var_name aux points (1:znbval)
+  int status;
+
+  int var_id;
+  int local_dimids[1];
+  static size_t start[] = {0}; //1D
+  static size_t count[] = {1}; //1D
+
+  count[0]= (size_t) znbval;
+  
+
+  local_dimids[0]=dimids[1]; //z
+  
+
+  //Define the netCDF variable
+  status=nc_redef(ncid);
+if (status != NC_NOERR)
+    {
+      handle_error (status, "Erreur change to def mode NetCDF ::voncwrite");
+    }
+ 
+ status=nc_def_var(ncid,var_name,NC_REAL,1,local_dimids,&var_id);
+ 
+  if (status != NC_NOERR)
+    {
+      printf("variable %s:",var_name);
+      handle_error (status, "Erreur define NetCDF variable ::voncwrite");
+    }
+ nc_enddef(ncid);
+  status=nc_put_vara_real(ncid,var_id,start,count,vect);
+ if (status != NC_NOERR)
+    {
+      handle_error (status, "Erreur writing NetCDF variable ::voncwrite");
+    }
+}