Changeset 14 for trunk/include/meshmask.h

Timestamp:

03/18/11 17:35:53 (13 years ago)

Author:

jbrlod

Message:

experience 3D

File:

• trunk/include/meshmask.h (modified) (2 diffs)

trunk/include/meshmask.h

@@ -410 +410 @@
 }
 
-void write_rst_var(int t,   int ncid, int dimids[]) {
+void write_out_var_init (int ncid, int dimids[])
+{
+  //Ecrit les variables de navigation pour les sorties
 YREAL vect[NZ*NY*NX];
  
+
 //nav_lon
  for(int j=0;j<NY;j++)
@@ -432 +435 @@
 }
 zncwrite(NZ,"nav_lev",ncid,dimids,vect);
+}
+
+void write_out_var(int t, int ncid, char const *suff )
+{
+  YREAL vect [NZ*NY*NX];
+  char varname[200];
+ int dimids[4];
+//retrouve les dimensions
+ nc_inq_dimid(ncid,"x",&dimids[3]);
+ nc_inq_dimid(ncid,"y",&dimids[2]);
+ nc_inq_dimid(ncid,"z",&dimids[1]);
+ nc_inq_dimid(ncid,"t",&dimids[0]);
+
+  //u
+  if (t==0) sprintf(varname,"ua_c_%s",suff);
+  else sprintf (varname,"ub_%s",suff);
+  for(int k=0;k<NZ;k++)
+    for(int j=0;j<NY;j++)
+      for(int i=0;i<NX;i++){
+        if (t==0) vect[(k)*(NY*NX)+(j)*(NX)+(i)]=YS_ua_c(0,i,j,k);
+        else vect[(k)*(NY*NX)+(j)*(NX)+(i)]=YS_ub(0,i,j,k,t-1);
+     }
+  voncwrite(0,NZ,NY,NX,varname,ncid,dimids,vect);
+  
+  //u
+  if (t==0) sprintf(varname,"va_c_%s",suff);
+  else sprintf (varname,"vb_%s",suff);
+  for(int k=0;k<NZ;k++)
+    for(int j=0;j<NY;j++)
+      for(int i=0;i<NX;i++){
+        if (t==0) vect[(k)*(NY*NX)+(j)*(NX)+(i)]=YS_va_c(0,i,j,k);
+        else vect[(k)*(NY*NX)+(j)*(NX)+(i)]=YS_vb(0,i,j,k,t-1);
+      }
+  voncwrite(0,NZ,NY,NX,varname,ncid,dimids,vect);
+  
+  //t
+  if (t==0) sprintf(varname,"ta_c_%s",suff);
+  else sprintf (varname,"tb_%s",suff);
+  for(int k=0;k<NZ;k++)
+    for(int j=0;j<NY;j++)
+      for(int i=0;i<NX;i++){
+        if (t==0) vect[(k)*(NY*NX)+(j)*(NX)+(i)]=YS_ta_c(0,i,j,k);
+        else vect[(k)*(NY*NX)+(j)*(NX)+(i)]=YS_tb(0,i,j,k,t-1);
+      }
+  voncwrite(0,NZ,NY,NX,varname,ncid,dimids,vect);
+  
+  //s
+  if (t==0) sprintf(varname,"sa_c_%s",suff);
+  else sprintf (varname,"sb_%s",suff);
+  for(int k=0;k<NZ;k++)
+    for(int j=0;j<NY;j++)
+      for(int i=0;i<NX;i++){
+        if (t==0) vect[(k)*(NY*NX)+(j)*(NX)+(i)]=YS_sa_c(0,i,j,k);
+        else vect[(k)*(NY*NX)+(j)*(NX)+(i)]=YS_sb(0,i,j,k,t-1);
+      }
+  voncwrite(0,NZ,NY,NX,varname,ncid,dimids,vect);
+  
+  //ssh
+  if (t==0) sprintf(varname,"sshn_c_%s",suff);
+  else sprintf (varname,"sshb_%s",suff);
+  for(int j=0;j<NY;j++)
+    for(int i=0;i<NX;i++){
+      if (t==0)  vect[(j)*(NX)+(i)]=YS_sshn_c(0,i,j);
+      else vect[(j)*(NX)+(i)]=YS_sshb(0,i,j,t-1);
+    }
+  soncwrite(0,NY,NX,varname,ncid,dimids,vect);
+  
+}
+
+
+void write_rst_var(int t,   int ncid, int dimids[]) {
+YREAL vect[NZ*NY*NX];
+ 
+//nav_lon
+ for(int j=0;j<NY;j++)
+      for(int i=0;i<NX;i++){
+        vect[(j)*(NX)+(i)]=nav_lon( i, j );
+      }
+ sncwrite(NY,NX,"nav_lon",ncid,dimids,vect);
+ 
+ //nav_lat
+ for(int j=0;j<NY;j++)
+   for(int i=0;i<NX;i++){
+     vect[(j)*(NX)+(i)]=nav_lat( i, j );
+   }
+ sncwrite(NY,NX,"nav_lat",ncid,dimids,vect);
+ 
+ //nav_lev
+ for(int k=0;k<NZ;k++){
+  vect[k]=nav_lev[k];  
+}
+zncwrite(NZ,"nav_lev",ncid,dimids,vect);
 
 //time_counter