Changeset 67


Ignore:
Timestamp:
08/09/11 17:46:09 (13 years ago)
Author:
berrada
Message:

filtre lateral explicit fichier dta_lfexp.h

Location:
branches/branche-mb/src
Files:
1 added
5 edited

Legend:

Unmodified
Added
Removed

  • branches/branche-mb/src/cost.dat

    r63 r67  
    22 0.0000000000000000e+00 
    33 0.0000000000000000e+00 
    4  2.1474823164582097e+02 
    5  1.2958901480355472e+01 
    6  2.4208867493170015e+01 
    7  9.5310667871496211e+00 
    8  6.1790369949928481e+00 
    9  4.0542139095518230e+00 
    10  2.3555976045385250e+00 
    11  1.8747388242251850e+00 
    12  1.6582963346427140e+00 
    13  1.4291909107170053e+00 
    14  1.2994902018124086e+00 
    15  1.3334911336228903e+00 
    16  1.2859361879963762e+00 
    17  1.5496856239892498e+00 
    18  1.2418585591044098e+00 
    19  2.0442371997362545e+00 
    20  1.2971715810562954e+00 
    21  1.2975780013916505e+00 
    22  1.2986727516755405e+00 
    23  1.3057984268538101e+00 
    24  1.2852046369970458e+00 
    25  1.2266854183064928e+00 
    26  1.2247496446634496e+00 
    27  1.2151832426274585e+00 
    28  1.2154570853596762e+00 
    29  1.2858728960761006e+00 
    30  1.2152605999261445e+00 
    31  1.2152356639508144e+00 
    32  1.2858909584114377e+00 
    33  1.2859036070737340e+00 
    34  1.2859124488276841e+00 

  • branches/branche-mb/src/dta_lfil.h

    r62 r67  
    1818 
    1919  if(Yt==1){ 
    20     double zdep; 
    21     if (Yk==0)  
    22       zdep=fsdept( Yi , Yj , 1 ) - fsdept( Yi , Yj , 0 ); 
    23     else     
    24       zdep=fsdept( Yi , Yj , Yk ) - fsdept( Yi , Yj , Yk-1 ); 
    25     double wm12=1./(106000*sqrt(zdep)); //1/sqrt(ds1*ds2*ds3) 
     20    double zbtr = zbtr2(Yi,Yj); 
     21    double wm12=sqrt(zbtr);//1./(106000.); //1/sqrt(ds1*ds2) 
    2622    YS1=wm12*x1; 
    2723  } 
    … …  
    3632  } 
    3733  // 
    38     if(Yt==1 && Yk==0 && Yi<26 && Yi>22 && Yj<17 && Yj>13)  
     34  /*    if(Yt==1 && Yk==0 && Yi<26 && Yi>22 && Yj<17 && Yj>13)  
    3935      printf("nul1(%d,%d) = %24.16e \n",Yi+1,Yj+1,YS1); 
    4036    if(Yt==40 && Yk==0 && Yi<26 && Yi>22 && Yj<17 && Yj>13)  
    41       printf("nul40(%d,%d) = %24.16e \n",Yi+1,Yj+1,YS1); 
     37    printf("nul40(%d,%d) = %24.16e \n",Yi+1,Yj+1,YS1);*/ 
    4238} 
    4339 
    … …  
    5147 
    5248  if(Yt==1){ 
    53     double zdep; 
    54     if (Yk==0)  
    55       zdep=fsdept( Yi , Yj , 1 ) - fsdept( Yi , Yj , 0 ); 
    56     else     
    57       zdep=fsdept( Yi , Yj , Yk ) - fsdept( Yi , Yj , Yk-1 ); 
    58     double wm12=1./(106000*sqrt(zdep)); //1/sqrt(ds1*ds2*ds3) 
     49    double zbtr = zbtr2(Yi,Yj); 
     50    double wm12=sqrt(zbtr);//(106000.); //1/sqrt(ds1*ds2) 
     51    //    double wm12=1./(106000.); //1/sqrt(ds1*ds2) 
    5952    YJ1I1=wm12; 
    6053  } 
    … …  
    7063  } 
    7164  // 
    72   if(Ycurward==BACKWARD){ 
     65  /*  if(Ycurward==BACKWARD){ 
    7366    if(Yt==10 && Yk==0 && Yi<26 && Yi>22 && Yj<17 && Yj>13)  
    7467      printf("YG_1L(%d,%d) = %24.16e \n",Yi+1,Yj+1,YG1); 
    75   } 
     68  }*/ 
    7669} 

  • branches/branche-mb/src/sinobad.d

    r63 r67  
    11defval K_NEU 
     2#define FZIMP 
    23#define FILTER 
    3 #define FILZNOL 
    4 //#define FILNOZL 
     4//#define FILZNOL 
     5#define FILLNOZ 
    56//#define SOLSORYAO 
    67//#define OPTIMORDER 
    … …  
    2122defval TU       2 | start time step 
    2223#ifdef FILTER 
    23 defval OFTL      0 
    24 defval OFTZ      25 
    25 defval OFT       25 
     24defval OFTL      10 
     25defval OFTZ      0 
     26defval OFT       10 
    2627#endif 
    2728#ifndef FILTER 
    … …  
    2930#endif 
    3031//#exec disp_valdef 
     32#ifdef FZIMP 
     33defval K_FZIMP 
     34#endif 
    3135#ifdef FILTER 
    3236defval K_FILTER 
    … …  
    3539defval K_FILZNOL 
    3640#endif 
    37 #ifdef FILNOZL 
    38 defval K_FILNOZL 
     41#ifdef FILLNOZ 
     42defval K_FILLNOZ 
    3943#endif 
    4044#ifdef SOLSORYAO 
    … …  
    6266#ifdef FILTER 
    6367traj  Tcstf      M   0  1 
    64 #ifndef FILNOZL 
    6568#ifndef FILZNOL 
    6669traj  Tlfil      M   1  0    1  OFTL 
    6770#endif 
     71#ifndef FILLNOZ 
    6872traj  Tzfil      M   1  OFTL 1  OFTZ 
    6973#endif 
    … …  
    7579#ifdef FILTER 
    7680space  S3df      M    NX NY NZ    Tcstf 
    77 #ifndef FILNOZL 
    7881#ifndef FILZNOL 
    7982space  S3dtlf    M    NX NY NZ    Tlfil 
    8083#endif 
     84#ifndef FILLNOZ 
    8185space  S3dtzf    M    NX NY NZ    Tzfil 
    8286#endif 
    … …  
    112116#ifdef FILTER 
    113117modul    nu                     space   S3df                    noward            output 1      target 
    114 #ifndef FILNOZL 
    115118#ifndef FILZNOL 
    116 modul    dta_lfil               space   S3dtlf                  input 6           output 1      tempo 
    117 modul    ztuv_lfil              space   S3dtlf                  input 3           output 2      tempo 
    118 #endif 
    119 modul    dta_zfilt              space   S3dtzf                  input 4           output 1      tempo 
    120 //modul  dta_zfil               space   S3dtzf                  input 4           output 1      tempo 
    121 //modul  zwy_zfil               space   S3dtzf                  input 2           output 1      tempo 
    122 #endif 
     119modul    dta_lfexp              space   S3dtlf                  input 6           output 1      tempo 
     120#endif 
     121 
     122#ifndef FILLNOZ 
     123#ifdef FZIMP 
     124modul    dta_zfimp_t            space   S3dtzf                  input 2           output 1      tempo 
     125modul    dta_zfimp              space   S3dtzf                  input 3           output 1      tempo 
     126#endif 
     127#ifndef FZIMP 
     128modul    dta_zfexp              space   S3dtzf                  input 4           output 1      tempo 
     129#endif 
     130 
     131#endif 
     132 
    123133modul    ta_c                   space   S3d                     input 1           output 1       
    124134#endif 
    … …  
    306316 
    307317#ifdef FILTER 
    308 #ifndef FILNOZL 
     318 
    309319#ifndef FILZNOL 
    310 //#---------->dta_lfil 
    311 ctin     dta_lfil       1       from    nu                      1  i    j    k   
    312 ctin     dta_lfil       2       from    dta_lfil                1  i    j    k-1 t-1 
    313 ctin     dta_lfil       3       from    ztuv_lfil               1  i-1  j    k  t  
    314 ctin     dta_lfil       4       from    ztuv_lfil               1  i    j    k  t 
    315 ctin     dta_lfil       5       from    ztuv_lfil               2  i    j-1  k  t  
    316 ctin     dta_lfil       6       from    ztuv_lfil               2  i    j    k  t  
    317  
    318 //#---------->ztuv_lfil 
    319 ctin     ztuv_lfil      1       from    dta_lfil                1  i+1  j    k     t-1 
    320 ctin     ztuv_lfil      2       from    dta_lfil                1  i    j    k     t-1  
    321 ctin     ztuv_lfil      3       from    dta_lfil                1  i    j+1  k     t-1  
     320//#---------->dta_lfexp 
     321ctin     dta_lfexp      1       from    nu                      1  i    j    k 
     322ctin     dta_lfexp      2       from    dta_lfexp               1  i    j    k     t-1  
     323ctin     dta_lfexp      3       from    dta_lfexp               1  i-1  j    k     t-1  
     324ctin     dta_lfexp      4       from    dta_lfexp               1  i+1  j    k     t-1 
     325ctin     dta_lfexp      5       from    dta_lfexp               1  i    j-1  k     t-1  
     326ctin     dta_lfexp      6       from    dta_lfexp               1  i    j+1  k     t-1  
    322327#endif 
    323328 
    324329//#---------->dta_zfil 
     330#ifndef FILLNOZ 
     331#ifndef FILZNOL 
     332#ifndef FZIMP 
     333ctin     dta_zfexp      1       from    dta_lfexp               1  i    j    k     t 
     334#endif 
     335#ifdef FZIMP 
     336ctin     dta_zfimp      1       from    dta_lfexp               1  i    j    k     t  
     337#endif 
     338#endif 
     339 
    325340#ifdef FILZNOL 
    326 ctin     dta_zfilt      1       from    nu                      1  i    j    k   
    327 //ctin   dta_zfil       1       from    nu                      1  i    j    k   
    328 #endif 
    329 #ifndef FILZNOL 
    330 ctin     dta_zfil       1       from    dta_lfil                1  i    j    k     t 
    331 #endif 
    332 ctin     dta_zfilt      2       from    dta_zfilt               1  i    j    k-1     t-1  
    333 ctin     dta_zfilt      3       from    dta_zfilt               1  i    j    k   t-1  
    334 ctin     dta_zfilt      4       from    dta_zfilt               1  i    j    k+1     t-1 
    335 //ctin   dta_zfil       2       from    dta_zfil                1  i    j    k     t-1  
    336 //ctin   dta_zfil       3       from    zwy_zfil                1  i    j    k+1   t  
    337 //ctin   dta_zfil       4       from    zwy_zfil                1  i    j    k     t 
    338  
    339 //#---------->zwy_zfil 
    340 //ctin   zwy_zfil       1       from    dta_zfil                1  i    j    k-1   t-1 
    341 //ctin   zwy_zfil       2       from    dta_zfil                1  i    j    k     t-1  
     341#ifndef FZIMP 
     342ctin     dta_zfexp      1       from    nu                      1  i    j    k   
     343#endif 
     344#ifdef FZIMP 
     345ctin     dta_zfimp      1       from    nu                      1  i    j    k  
     346#endif 
     347#endif 
     348 
     349#ifdef FZIMP 
     350ctin     dta_zfimp      2       from    dta_zfimp               1  i    j    k+1   t  
     351ctin     dta_zfimp      3       from    dta_zfimp_t             1  i    j    k     t 
     352 
     353ctin     dta_zfimp_t    1       from    dta_zfimp_t             1  i    j    k-1 t 
     354ctin     dta_zfimp_t    2       from    dta_zfimp               1  i    j    k  t-1 
     355#endif 
     356#ifndef FZIMP 
     357ctin     dta_zfexp      2       from    dta_zfexp               1  i    j    k-1     t-1  
     358ctin     dta_zfexp      3       from    dta_zfexp               1  i    j    k   t-1  
     359ctin     dta_zfexp      4       from    dta_zfexp               1  i    j    k+1     t-1 
     360#endif 
    342361 
    343362//#---------->ta_c 
    344 ctin     ta_c           1       from    dta_zfilt               1  i    j    k   t 
    345 //ctin   ta_c           1       from    dta_zfil                1  i    j    k   t 
    346 #endif 
    347 #ifdef FILNOZL 
    348 ctin     ta_c           1       from    nu                      1  i    j    k    
    349 #endif 
    350 #endif 
     363#ifdef FZIMP 
     364ctin     ta_c           1       from    dta_zfimp               1  i    j    k   t 
     365#endif 
     366#ifndef FZIMP 
     367ctin     ta_c           1       from    dta_zfexp               1  i    j    k   t 
     368#endif 
     369#endif 
     370#ifdef FILLNOZ 
     371ctin     ta_c           1       from    dta_lfexp               1  i    j    k   t 
     372#endif 
     373 
     374 
     375#endif 
     376 
    351377//#ctin  ta_c                   1..NPCA from    pca_ta                  1..NPCA  1 
    352378 
    353379//#exec disp_ct_in 
    354 exec disp_modul 
     380//exec disp_modul 
    355381 
    356382 
    … …  
    916942#ifdef FILTER 
    917943// calcul du filtre 
    918 #ifndef FILNOZL 
    919944#ifndef FILZNOL 
    920945order modinspace S3dtlf 
    921       order YA3   
    922             order YA2 YA1 
    923                   ztuv_lfil 
    924             forder 
    925             order YA2 YA1 
    926                   dta_lfil 
    927             forder 
    928       forder 
    929 forder 
    930 #endif 
     946    order YA3 YA2 YA1 
     947          dta_lfexp 
     948    forder 
     949forder 
     950#endif 
     951#ifndef FILLNOZ 
    931952order modinspace S3dtzf 
    932 //    order YA3 YA2 YA1 
    933 //        zwy_zfil 
    934 //    forder 
     953#ifdef FZIMP 
     954    order YA2 YA1 
     955      order YA3 
     956          dta_zfimp_t 
     957      forder 
     958      order YB3 
     959          dta_zfimp 
     960      forder 
     961    forder 
     962#endif 
     963#ifndef FZIMP 
    935964    order YA3 YA2 YA1 
    936 //                dta_zfil 
    937           dta_zfilt 
     965          dta_zfexp 
    938966    forder 
    939 forder 
    940 #endif 
     967#endif 
     968forder 
     969#endif 
     970 
    941971order modinspace S3d 
    942972       order YA3 YA2 YA1 
    … …  
    11451175#ifdef FILTER 
    11461176 
    1147 #ifndef FILNOZL 
    11481177#ifndef FILZNOL 
    11491178order spaceintraj Tlfil 
    … …  
    11511180forder 
    11521181#endif 
     1182#ifndef FILLNOZ 
    11531183order spaceintraj Tzfil 
    11541184        S3dtzf 
    11551185forder 
    11561186#endif 
     1187 
    11571188order spaceintraj Tcst 
    11581189        S3d 

  • branches/branche-mb/src/sinobad.h

    r62 r67  
    1414 
    1515#ifdef K_FILTER 
     16#ifdef K_FZIMP 
     17void luzimp_init(); // initialise tab_luzimp 
     18#define tab_luzimp( i, j, k , t )  (tab_luzimp[(i)*(NY*NZ)+(j)*(NZ)+(k)][t]) // factorisation LU  
     19double tab_luzimp[NZ*NY*NX][2];   // tableau factorisation LU 
     20#endif 
    1621#define ta_eb( i, j, k )  (ta_eb[(i)*(NY*NZ)+(j)*(NZ)+(k)]) // ebauche 
    1722double ta_eb[NZ*NY*NX];   // ebauche 
    … …  
    8186#ifdef K_FILTER 
    8287double pdth_fil=100.;//7200.; 
    83 double pdtz_fil=100.;//7200.; 
    84 double fsaht_fil= 1.2e+2;//1000.; 
    85 double avt_fil=1.2e-5; 
    86 double lscale= 212000.; // length-scale 
     88double pdtz_fil=100.;//100.;//7200.; 
     89double fsaht_fil= 5.e+2;//1000.; 
     90double avt_fil=1.2e-05;//1.2e-5; 
     91//double lscale= 212000.; // length-scale 
     92double norm_fil; //coef de normalisation 
     93void normfil();// calcul du coef de normalisation norm_fil 
    8794#endif 
    8895 
    … …  
    146153  cal_ff(); 
    147154  xdom_init(); 
     155#ifdef K_FILTER 
     156#ifdef K_FZIMP 
     157  luzimp_init(); 
     158#endif 
     159  normfil(); 
     160#endif 
    148161  xsolmat_init(); 
    149162  xflt_rst(); 
    … …  
    151164  xdisplay(); 
    152165  xtraqsr_init(); 
    153    
    154166} 
    155167//____________________________________________________________________________ 
    … …  
    244256  normsup_alla_c(); 
    245257} 
     258 
     259//_____________________________________________________________________________ 
    246260#ifdef K_FILTER 
     261void normfil(){ 
     262  double norm_zfil=1.; 
     263  double norm_lfil=1.; 
     264#if defined (K_FILZNOL) /* && !defined (K_FILLNOZ)*/  
     265  norm_zfil= 2.*sqrt(PI*avt_fil*pdtz_fil*2.*OFTZ); 
     266#ifdef K_FZIMP 
     267  double fact=1.;  
     268  int M=2*(OFTZ-1); 
     269  for(int i=1;i<M;i++) 
     270      fact = 4.*fact * (double) i / (double) (M+i-1); 
     271  norm_zfil=2.*fact*sqrt(avt_fil*pdtz_fil); 
     272#endif 
     273#elif  defined (K_FILLNOZ) /*&& !defined (K_FILZNOL)*/   
     274  norm_lfil= 2.*sqrt(PI*fsaht_fil*pdth_fil*2.*OFTL); 
     275#endif 
     276  norm_fil= norm_zfil*norm_lfil; 
     277  printf("norm_fil = %23.16e\n",norm_fil); 
     278} 
     279//_____________________________________________________________________________ 
    247280void init_nu(){ 
    248281 
    … …  
    251284      for(int j=0;j<NY;j++) 
    252285        for(int i=0;i<NX;i++){ 
    253           YS_nu(0,i,j,k)=YS_nu(0,i,j,k)+1.*exp(-((i+1.-(3.*NX)/4.)*(i+1.-(3.*NX)/4.)+  
     286          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.)+  
    254287                                              (j+1.-(3.*NY)/4.)*(j+1.-(3.*NY)/4.))/4.) 
    255288            *(1.-exp(-double(k+1.-NZ)/NZ))*(1.-exp(-double(tt)/24.)); 
    256289            } 
    257290} 
     291//_____________________________________________________________________________ 
     292#ifdef K_FZIMP 
     293//_____________________________________________________________________________ 
     294void luzimp_init(){ 
     295  for(int j=0;j<NY;j++) 
     296    for(int i=0;i<NX;i++) 
     297      for(int k=0;k<NZ;k++){ // init a 0 
     298        tab_luzimp( i, j, k , 0 ) = 0.;// 0-->gam 
     299        tab_luzimp( i, j, k , 1 ) = 0.;// 1-->bet 
     300      } 
     301  double a0,a1; 
     302  for(int j=1;j<NY-1;j++) 
     303    for(int i=1;i<NX-1;i++){ 
     304      a1=avt_fil*pdtz_fil/fse3w(i,j,1); 
     305      tab_luzimp( i, j, 0 , 0 ) = 0.;// 0-->gam 
     306      tab_luzimp( i, j, 0 , 1 ) = fse3t(i,j,0)+a1;// 1-->bet 
     307#ifndef K_NEUM 
     308      a0=avt_fil*pdtz_fil/fse3w(i,j,0); 
     309      tab_luzimp( i, j, 0 , 1 ) += a0;// 1-->bet 
     310#endif 
     311      for(int k=1;k<NZ-2;k++){ 
     312        a0=avt_fil*pdtz_fil/fse3w(i,j,k); 
     313        a1=avt_fil*pdtz_fil/fse3w(i,j,k+1); 
     314        tab_luzimp( i, j, k , 0 ) = -a0/tab_luzimp( i, j, k-1 , 1); 
     315        tab_luzimp( i, j, k , 1 ) = fse3t(i,j,k)+a1+a0 + a0 * tab_luzimp( i, j, k , 0 ); 
     316      } 
     317      a0=avt_fil*pdtz_fil/fse3w(i,j,NZ-2); 
     318      tab_luzimp( i, j, NZ-2 , 0 ) = -a0/tab_luzimp( i, j, NZ-3 , 1); 
     319      tab_luzimp( i, j, NZ-2 , 1 ) = fse3t(i,j,NZ-2)+a0 + a0 * tab_luzimp( i, j, NZ-2 , 0 ); 
     320#ifndef K_NEUM 
     321      a1=avt_fil*pdtz_fil/fse3w(i,j,NZ-1); 
     322      tab_luzimp( i, j, NZ-2 , 1 ) += a1; 
     323#endif 
     324    } 
     325  /*  for(int k=0;k<NZ;k++){ // init a 0 
     326    printf("gam(%d)=%f \t",k+1,tab_luzimp( 5, 5, k , 0 )); 
     327    printf("bet(%d)=%f \n",k+1,tab_luzimp( 5, 5, k , 1 )); 
     328    }*/ 
     329} 
     330#endif 
    258331#endif 
    259332//_____________________________________________________________________________ 
    … …  
    678751 
    679752//---------------------------------------------- 
    680 #if defined (YE_ta_c) || defined (YE_pca_ta) 
     753#if defined (YE_ta_c) || defined (YE_pca_ta) || defined (YE_nu) 
    681754#define true_tac( i, j, k )   (true_tac[(k)*(NY*NX)+(j)*(NX)+(i)]) 
    682755double true_tac[NZ*NY*NX]; 
    … …  
    701774  int i,j,k; 
    702775  if(!strcmp(argv[1],"ta_c")){ 
    703 #if defined (YE_ta_c) || defined (YE_pca_ta) 
     776#if defined (YE_ta_c) || defined (YE_pca_ta) || defined (YE_nu) 
    704777    for(i=0;i<NX;i++) 
    705778      for(j=0;j<NY;j++) 
    … …  
    743816  int i,j,k; 
    744817   
    745 #if defined (YE_ta_c) || defined (YE_pca_ta) 
     818#if defined (YE_ta_c) || defined (YE_pca_ta) || defined (YE_nu) 
    746819  sum=0.;sumn=0.; 
    747820  for(i=0;i<NX;i++) 

  • branches/branche-mb/src/ta_c.h

    r62 r67  
    88forward (YREAL x1) 
    99{  
    10  
    11   //if(Yk==0 && Yi<14 && Yi>10 && Yj<14 && Yj>10) printf("nu(%d,%d) = %24.16e \n",Yi+1,Yj+1,x1); 
    12  
    13   double norm_fac; 
    14 #ifdef K_FILZNOL 
    15   norm_fac= 2.*sqrt(PI*avt_fil*pdtz_fil*2.*OFTZ);///1.e-01;//1./sqrt(2.*PI)/lscale; 
    16 #elif defined (K_FILNOZL)  
    17   norm_fac= 1.e+00; 
    18 #else 
    19   norm_fac= 1.e+00; 
    20 #endif 
    21   YS1 = ta_eb(Yi,Yj,Yk) +norm_fac*x1; 
     10  YS1 = ta_eb(Yi,Yj,Yk) +norm_fil*x1; 
    2211  if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1 || Yk==NZ-1){ 
    2312    YS1=0.; 
    … …  
    3019  //YS1=sum; 
    3120  //// 
    32   /* if(fabs(x1)>0.1)  
    33      printf("increm(%d,%d,%d) = %24.16e (%24.16e) \n",Yi+1,Yj+1,Yk+1,x1,YS1);*/ 
     21  //  printf("norm_fac= %24.16e\n",norm_fac); 
     22  if(fabs(norm_fil*x1)>5.)  
     23      printf("increm(%d,%d,%d) = %24.16e (%24.16e) \n",Yi+1,Yj+1,Yk+1,norm_fil*x1,YS1); 
    3424} 
    3525//=========================================================================== 
    … …  
    3828backward (YREAL x1) 
    3929{ 
    40   double norm_fac; 
    41 #ifdef K_FILZNOL 
    42   norm_fac= 2.*sqrt(PI*avt_fil*pdtz_fil*2.*OFTZ);///1.e-01;//1./sqrt(2.*PI)/lscale; 
    43 #elif defined (K_FILNOZL)  
    44   norm_fac= 1.e+00; 
    45 #else 
    46   norm_fac= 1.e+00; 
    47 #endif 
    48   YJ1I1=norm_fac;  
     30  YJ1I1=norm_fil;  
    4931  if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1 || Yk==NZ-1){ 
    5032    YJ1I1=0.; 
Note: See TracChangeset for help on using the changeset viewer.