source: branches/branche-mb/src/ta_c.h @ 68

/***************************************************************************
                 module classe ta_c.h  -  description
***************************************************************************/
// mohamed.berrada@upmc.fr
// locean-ipsl.upmc, Paris, March 18, 2010
//===========================================================================
//                          methode forward
forward (YREAL x1)
{ 
  YS1 = ta_eb(Yi,Yj,Yk) +norm_fil*x1;
  if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1 || Yk==NZ-1){
    YS1=0.;
  }


  ///*  1..NPCA from    pca_ta                  1..NPCA  1*/
  //double sum=moy_ta(Yi,Yj,Yk);
  //for(int n=0;n<NPCA;n++) sum+=Yting[n]*stdev_ta[n]*shfs_ta(Yi,Yj,Yk,n);
  //YS1=sum;
  ////
  //  printf("norm_fac= %24.16e\n",norm_fac);
  if(fabs(norm_fil*x1)>5) 
      printf("increm(%d,%d,%d) = %24.16e (%24.16e) \n",Yi+1,Yj+1,Yk+1,norm_fil*x1,YS1);
}
//===========================================================================
//                         methode  backward

backward (YREAL x1)
{
  YJ1I1=norm_fil; 
  if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1 || Yk==NZ-1){
    YJ1I1=0.;
  }
   
  //for(int n=0;n<NPCA;n++)     
  //  Yjac[0][n]    = stdev_ta[n]*shfs_ta(Yi,Yj,Yk,n);
  ////
}

//===========================================================================
//***************************  FIN DU MODULE ta_c ***************************