Changeset 75 for trunk/src/sinobad.h

Timestamp:

03/02/12 09:38:08 (12 years ago)

Author:

jbrlod

Message:

new option for xchangmode

File:

• trunk/src/sinobad.h (modified) (8 diffs)

trunk/src/sinobad.h

@@ -46 +46 @@
 #define hdivn_at_TU( i, j, k) (hdivn_at_TU[(k)*(NY*NX)+(j)*(NX)+(i)]) 
 #define gcx_at_TU( i, j) (gcx_at_TU[(j)*(NX)+(i)]) 
+#define masque_obs( i, j) (masque_obs[(j)*(NX)+(i)]) 
+
+
 double  tb_neuler1[NZ*NY*NX]; //
 double  sb_neuler1[NZ*NY*NX]; //
@@ -54 +57 @@
 double  hdivn_at_TU[NZ*NY*NX]; //
 double  gcx_at_TU[NY*NX]; //
+
+short masque_obs[NY*NX]; // 
 
 void true_target_in_tab(int argc, char *argv[]); // charge YS_*a_c(0,i,j,k) dans true_*a_c(i,j,k)
@@ -107 +112 @@
 void xistate_init(int argc, char *argv[]);
 void xrst_save(int argc, char *argv[]);
-
+void xinit_masque_obs();
 //////////////////////////////////////////////////////////////
 // Les fonctions OBLIGATOIRES
@@ -163 +168 @@
   normfil();
 #endif
+  xinit_masque_obs();
   xsolmat_init();
   xflt_rst();
@@ -195 +201 @@
 void after_it (int nit)
 {
-  xdisplay();
+  // xdisplay();
 }
 //_____________________________________________________________________________
@@ -253 +259 @@
        return(0);
    }
+
+   if(indic==YIO_SAVESTATE && (!strcmp(nmmod,"tb") || !strcmp(nmmod,"sb")) && savemode==3) return(masque_obs(iaxe,jaxe));
 
     
@@ -647 +655 @@
 }
 
+void xinit_masque_obs(){
+  int i,j;
+  for (i=0;i<NX;i++)
+    for(j=0;j<NY;j++)
+    masque_obs(i,j)=0;
+}
+
 void xchangesavemode(int argc, char *argv[]){
   short newmode=atoi(argv[1]);
@@ -656 +671 @@
     printf("\n savemode : 1 point (Yi=%d,Yj=%d,Yk=%d)\n",isave,jsave,ksave);
   }
-
-}
+  if (newmode==3) { //sauver N profil
+    int Np; // Nombre d'obs
+    int ix,iy,comp=0; //indices de boucles
+    float dx,dy; // pas dans le tableau masque 
+    int nx,ny ; //nombre d'obs selon x et y
+    Np=atoi(argv[2]);
+    nx=(int)sqrt(Np);
+    ny=Np/nx;
+    dx=(float)NX/(nx+1);
+    dy=(float)NY/(ny+1);
+    for(ix=1;ix<=nx;ix++)
+      for(iy=1;iy<=ny;iy++) {
+        masque_obs((int)(ix*dx),(int)(iy*dy))=1;
+        comp++;
+      }
+    printf("\n savemode : %d profils observés(dx=%f,dy=%f)\n",comp,dx,dy);
+  } 
+}
+
 
 //_____________________________________________________________________________