Changeset 68 for branches/branche-mb/src/sinobad.h

Timestamp:

08/30/11 16:13:04 (13 years ago)

Author:

berrada

Message:

commit final avec filtre

File:

• branches/branche-mb/src/sinobad.h (modified) (4 diffs)

branches/branche-mb/src/sinobad.h

@@ -87 +87 @@
 double pdth_fil=100.;//7200.;
 double pdtz_fil=100.;//100.;//7200.;
-double fsaht_fil= 5.e+2;//1000.;
+double fsaht_fil= 5.e+02;//1000.;
 double avt_fil=1.2e-05;//1.2e-5;
 //double lscale= 212000.; // length-scale
 double norm_fil; //coef de normalisation
 void normfil();// calcul du coef de normalisation norm_fil
+void xset_fdt_fco();// set fdt fcoef
 #endif
 
@@ -259 +260 @@
 //_____________________________________________________________________________
 #ifdef K_FILTER
+//_____________________________________________________________________________
 void normfil(){
   double norm_zfil=1.;
   double norm_lfil=1.;
-#if defined (K_FILZNOL) /* && !defined (K_FILLNOZ)*/ 
-  norm_zfil= 2.*sqrt(PI*avt_fil*pdtz_fil*2.*OFTZ);
-#ifdef K_FZIMP
+  norm_lfil= 2.*sqrt(PI*fsaht_fil*pdth_fil*2.*(OFTL-1.));// schema explicite
+  norm_zfil= 2.*sqrt(PI*avt_fil*pdtz_fil*2.*(OFTZ-1.));// schema explicite 
+#ifdef K_FZIMP  // schema implicite 
   double fact=1.; 
   int M=2*(OFTZ-1);
@@ -271 +273 @@
   norm_zfil=2.*fact*sqrt(avt_fil*pdtz_fil);
 #endif
-#elif  defined (K_FILLNOZ) /*&& !defined (K_FILZNOL)*/  
-  norm_lfil= 2.*sqrt(PI*fsaht_fil*pdth_fil*2.*OFTL);
-#endif
+
+#if defined (K_FILZNOL)
+  norm_lfil=1.;
+#endif
+#if  defined (K_FILLNOZ)
+  norm_zfil=1.;
+#endif
+
   norm_fil= norm_zfil*norm_lfil;
   printf("norm_fil = %23.16e\n",norm_fil);
 }
 //_____________________________________________________________________________
-void init_nu(){
-
-  for(int tt=0;tt<24;tt++)
+void xset_fdt_fco(int argc, char *argv[]){// set fdt fcoef
+  if (argc!=3){
+    printf("\n inappropriate arguments choice in command xset_fdt_fco\n");
+    exit(99);
+  }
+  if(!strcmp(argv[1],"th")){
+    pdth_fil=atof(argv[2]);
+  }
+  else if(!strcmp(argv[1],"tz")){
+    pdtz_fil=atof(argv[2]);
+  }
+  else if(!strcmp(argv[1],"kh")){
+    fsaht_fil=atof(argv[2]);
+  }
+  else if(!strcmp(argv[1],"kz")){
+    avt_fil=atof(argv[2]);
+  }
+  else{
+    printf("\n inappropriate first argument choice in command xset_fdt_fco\n");
+    exit(99);
+  }
+}
+//_____________________________________________________________________________
+void init_nu(int argc, char *argv[]){
+  if (argc==1){
+    for(int tt=0;tt<24;tt++)
+      for(int k=0;k<NZ;k++)
+        for(int j=0;j<NY;j++)
+          for(int i=0;i<NX;i++){
+            YS_nu(0,i,j,k)=YS_nu(0,i,j,k)+0.1*exp(-((i+1.-(3.*NX)/4.)*(i+1.-(3.*NX)/4.)+ 
+                                                    (j+1.-(3.*NY)/4.)*(j+1.-(3.*NY)/4.))/4.)
+              *(1.-exp(-double(k+1.-NZ)/NZ))*(1.-exp(-double(tt)/24.));
+          }
+  }
+  else if  (argc==5){
     for(int k=0;k<NZ;k++)
       for(int j=0;j<NY;j++)
         for(int i=0;i<NX;i++){
-          YS_nu(0,i,j,k)=YS_nu(0,i,j,k)+0.1*exp(-((i+1.-(3.*NX)/4.)*(i+1.-(3.*NX)/4.)+ 
-                                              (j+1.-(3.*NY)/4.)*(j+1.-(3.*NY)/4.))/4.)
-            *(1.-exp(-double(k+1.-NZ)/NZ))*(1.-exp(-double(tt)/24.));
-            }
+          YS_nu(0,i,j,k)=0.;
+          }
+
+    float v=atof(argv[1]);
+    int ji=atoi(argv[2])-1;
+    int jj=atoi(argv[3])-1;
+    int jk=atoi(argv[4])-1;
+    YS_nu(0,ji,jj,jk)=v;
+  }
+  else{
+    printf("\n inappropriate arguments choice in command init_nu\n");
+    exit(99);
+  }
 }
 //_____________________________________________________________________________
@@ -300 +348 @@
       }
   double a0,a1;
+  double a2,b2,c2;
   for(int j=1;j<NY-1;j++)
     for(int i=1;i<NX-1;i++){
       a1=avt_fil*pdtz_fil/fse3w(i,j,1);
-      tab_luzimp( i, j, 0 , 0 ) = 0.;// 0-->gam
-      tab_luzimp( i, j, 0 , 1 ) = fse3t(i,j,0)+a1;// 1-->bet
+      b2= 1.-a1/fse3t(i,j,0);
 #ifndef K_NEUM
       a0=avt_fil*pdtz_fil/fse3w(i,j,0);
-      tab_luzimp( i, j, 0 , 1 ) += a0;// 1-->bet
-#endif
+      b2 -= (a0/fse3t(i,j,0));// 1-->bet
+#endif
+      tab_luzimp( i, j, 0 , 0 ) = 0.;// 0-->gam
+      tab_luzimp( i, j, 0 , 1 ) = b2;// 1-->bet
       for(int k=1;k<NZ-2;k++){
         a0=avt_fil*pdtz_fil/fse3w(i,j,k);
         a1=avt_fil*pdtz_fil/fse3w(i,j,k+1);
-        tab_luzimp( i, j, k , 0 ) = -a0/tab_luzimp( i, j, k-1 , 1);
-        tab_luzimp( i, j, k , 1 ) = fse3t(i,j,k)+a1+a0 + a0 * tab_luzimp( i, j, k , 0 );
+        a2 = a0/fse3t(i,j,k);
+        b2 =  1.-a1/fse3t(i,j,k)-a0/fse3t(i,j,k);
+        c2 = a0/fse3t(i,j,k-1);
+        tab_luzimp( i, j, k , 0 ) = c2/tab_luzimp( i, j, k-1 , 1);
+        tab_luzimp( i, j, k , 1 ) = b2 - a2 * tab_luzimp( i, j, k , 0 );
       }
       a0=avt_fil*pdtz_fil/fse3w(i,j,NZ-2);
-      tab_luzimp( i, j, NZ-2 , 0 ) = -a0/tab_luzimp( i, j, NZ-3 , 1);
-      tab_luzimp( i, j, NZ-2 , 1 ) = fse3t(i,j,NZ-2)+a0 + a0 * tab_luzimp( i, j, NZ-2 , 0 );
+      a2 = a0/fse3t(i,j,NZ-2);
+      b2 =  1.-a0/fse3t(i,j,NZ-2);
 #ifndef K_NEUM
       a1=avt_fil*pdtz_fil/fse3w(i,j,NZ-1);
-      tab_luzimp( i, j, NZ-2 , 1 ) += a1;
-#endif
-    }
-  /*  for(int k=0;k<NZ;k++){ // init a 0
-    printf("gam(%d)=%f \t",k+1,tab_luzimp( 5, 5, k , 0 ));
-    printf("bet(%d)=%f \n",k+1,tab_luzimp( 5, 5, k , 1 ));
-    }*/
+      b2 -= a1/fse3t(i,j,NZ-2);
+#endif
+      c2 = a0/fse3t(i,j,NZ-3);
+     
+      tab_luzimp( i, j, NZ-2 , 0 ) = c2/tab_luzimp( i, j, NZ-3 , 1);
+      tab_luzimp( i, j, NZ-2 , 1 ) = b2 - a2 * tab_luzimp( i, j, NZ-2 , 0 );
+    }
+  for(int k=0;k<NZ;k++){ // init a 0
+    printf("gam(%d)=%f \t",k+1,tab_luzimp( 23, 15, k , 0 ));
+    printf("bet(%d)=%f \n",k+1,tab_luzimp( 23, 15, k , 1 ));
+    }
 }
 #endif