Changeset 67 for branches/branche-mb/src/dta_lfil.h

Timestamp:

08/09/11 17:46:09 (13 years ago)

Author:

berrada

Message:

filtre lateral explicit fichier dta_lfexp.h

File:

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

branches/branche-mb/src/dta_lfil.h

@@ -18 +18 @@
 
   if(Yt==1){
-    double zdep;
-    if (Yk==0) 
-      zdep=fsdept( Yi , Yj , 1 ) - fsdept( Yi , Yj , 0 );
-    else    
-      zdep=fsdept( Yi , Yj , Yk ) - fsdept( Yi , Yj , Yk-1 );
-    double wm12=1./(106000*sqrt(zdep)); //1/sqrt(ds1*ds2*ds3)
+    double zbtr = zbtr2(Yi,Yj);
+    double wm12=sqrt(zbtr);//1./(106000.); //1/sqrt(ds1*ds2)
     YS1=wm12*x1;
   }
@@ -36 +32 @@
   }
   //
-    if(Yt==1 && Yk==0 && Yi<26 && Yi>22 && Yj<17 && Yj>13) 
+  /*    if(Yt==1 && Yk==0 && Yi<26 && Yi>22 && Yj<17 && Yj>13) 
       printf("nul1(%d,%d) = %24.16e \n",Yi+1,Yj+1,YS1);
     if(Yt==40 && Yk==0 && Yi<26 && Yi>22 && Yj<17 && Yj>13) 
-      printf("nul40(%d,%d) = %24.16e \n",Yi+1,Yj+1,YS1);
+    printf("nul40(%d,%d) = %24.16e \n",Yi+1,Yj+1,YS1);*/
 }
 
@@ -51 +47 @@
 
   if(Yt==1){
-    double zdep;
-    if (Yk==0) 
-      zdep=fsdept( Yi , Yj , 1 ) - fsdept( Yi , Yj , 0 );
-    else    
-      zdep=fsdept( Yi , Yj , Yk ) - fsdept( Yi , Yj , Yk-1 );
-    double wm12=1./(106000*sqrt(zdep)); //1/sqrt(ds1*ds2*ds3)
+    double zbtr = zbtr2(Yi,Yj);
+    double wm12=sqrt(zbtr);//(106000.); //1/sqrt(ds1*ds2)
+    //    double wm12=1./(106000.); //1/sqrt(ds1*ds2)
     YJ1I1=wm12;
   }
@@ -70 +63 @@
   }
   //
-  if(Ycurward==BACKWARD){
+  /*  if(Ycurward==BACKWARD){
     if(Yt==10 && Yk==0 && Yi<26 && Yi>22 && Yj<17 && Yj>13) 
       printf("YG_1L(%d,%d) = %24.16e \n",Yi+1,Yj+1,YG1);
-  }
+  }*/
 }