Changeset 4

Timestamp:

03/03/11 16:41:35 (13 years ago)

Author:

jbrlod

Message:

rst write functionnality

Location:

trunk/include

Files:

• meshmask.h (modified) (2 diffs)
• ncutil.cpp (modified) (1 diff)
• ncutil.h (modified) (1 diff)

trunk/include/meshmask.h

@@ -104 +104 @@
 void xinit_mesh_mask_nc();// lire le fichier mesh_mask.nc
 void xistate_rest_file();
+
+
+//Ecriture du netcdf
+  void write_rst_global_att(int ncid); //Ecrit les attributs globaux du restart standard
+
+
 void    xinit_mesh_mask_nc()
 {
@@ -342 +348 @@
     }
 }
+
+void write_rst_global_att(int ncid) {
+
+  //:DOMAIN_number_total = 1 ;
+  static int DOMAIN_number_total=1;
+  nc_put_att_int (ncid,NC_GLOBAL,"DOMAIN_number_total",NC_INT,1,&DOMAIN_number_total);
+
+  //:DOMAIN_number = 0 ;
+  static int DOMAIN_number=0;
+  nc_put_att_int (ncid,NC_GLOBAL,"DOMAIN_number",NC_INT,1,&DOMAIN_number);
+
+  //:DOMAIN_dimensions_ids = 1, 2 ;
+  static int DOMAIN_dimensions_ids[]={1,2};
+  nc_put_att_int (ncid,NC_GLOBAL,"DOMAIN_dimensions_ids",NC_INT,2,DOMAIN_dimensions_ids);
+
+  //:DOMAIN_size_global = 32, 22 ;
+  static int DOMAIN_size_global[]={NX,NY};
+  nc_put_att_int (ncid,NC_GLOBAL,"DOMAIN_size_global",NC_INT,2,DOMAIN_size_global);
+
+  //:DOMAIN_size_local = 32, 22 ;
+  static int DOMAIN_size_local[]={NX,NY};
+  nc_put_att_int (ncid,NC_GLOBAL,"DOMAIN_size_local",NC_INT,2,DOMAIN_size_local);
+
+  //:DOMAIN_position_first = 1, 1 ;
+  static int DOMAIN_position_first[]={1,1};
+  nc_put_att_int (ncid,NC_GLOBAL,"DOMAIN_position_first",NC_INT,2,DOMAIN_position_first);
+
+  //:DOMAIN_position_last = 32, 22 ;
+ static int DOMAIN_position_last[]={NX,NY};
+ nc_put_att_int (ncid,NC_GLOBAL,"DOMAIN_position_last",NC_INT,2,DOMAIN_position_last);
+
+ //:DOMAIN_halo_size_start = 0, 0 ;
+ static int DOMAIN_halo_size_start[]={0,0};
+ nc_put_att_int (ncid,NC_GLOBAL,"DOMAIN_halo_size_start",NC_INT,2,DOMAIN_halo_size_start);
+
+ //:DOMAIN_halo_size_end = 0, 0 ;
+ static int DOMAIN_halo_size_end[]={0,0};
+  nc_put_att_int (ncid,NC_GLOBAL,"DOMAIN_halo_size_end",NC_INT,2,DOMAIN_halo_size_end);
+
+  //DOMAIN_type = "BOX" ;
+  nc_put_att_text(ncid,NC_GLOBAL,"DOMAIN_type",3,"BOX");
+
+}

trunk/include/ncutil.cpp

@@ -6 +6 @@
 inline void     handle_error (int status, char *txt)
 {
-        if (status != NC_NOERR)2
+        if (status != NC_NOERR)
         { printf ("netcdf error: %s: netcdf statut: ", txt);
         cout << nc_strerror (status) << endl;
         exit (-1);
         }
+}
+
+int Ouvre_nc_write(char const *var_file)
+{
+  ///* ouvre un fichier nc pour écriture et renvoie son id var_file_id
+  int var_file_id;
+  int status = nc_create(var_file, NC_CLOBBER, &var_file_id);
+  if (status != NC_NOERR) {
+  char file_error[100];
+    sprintf(file_error,"open %s",var_file);
+    handle_error (status, file_error);
+  }
+  return(var_file_id);
 }
 

trunk/include/ncutil.h

@@ -11 +11 @@
 void    handle_error (int status, char *txt);///*fonction de message d'erreur
 int Ouvre_nc(char const *var_file );///* Ouvre un fichier nc et renvoie son id var_file_id
+int Ouvre_nc_write(char const *var_file);   ///* ouvre un fichier nc pour écriture et renvoie son id var_file_id
 int Var_id(int  var_file_id,char* varn); ///*initial traceurs actifs
+
 
 void vonclire(int t,int znbval, int ynbval, int xnbval, int ncid, int varid, double vect[]);