/***************************************************************************
                 module classe ta_c1z.h  -  description
***************************************************************************/
// Mohamed Berrada
// locean-ipsl.upmc, Paris, December 15, 2009
//===========================================================================
//                          methode forward
forward ()
{ /*	1..NPCA	from	pca			1..NPCA  1*/
#ifdef YE_pca
  if(Yk!=NZ-1){
    double sum=ta_c_moy[Yk];
    if(Yi==28 && Yj==12) { 
      for(int j=0;j<NPCA;j++) sum+=Yting[j]*stdev_pca[j]*shfs[Yk][j];
      YS1=sum;
      //     printf("taz1=%23.16e\n",sum);
      //     printf("tac=%23.16e\n\n",Yting[NPCA]);
    }
    else
      YS1=Yting[NPCA];//YS_ta_c(0,Yi,Yj,Yk);
  }
#else
  YS1=Yting[NPCA];//YS_ta_c(0,Yi,Yj,Yk);
#endif
  //
}
//===========================================================================
//                         methode  backward

backward ()
{
#ifdef YE_pca
  if(Yk!=NZ-1){
    if(Yi==28 && Yj==12){  
      for(int j=0;j<NPCA;j++) 	
	Yjac[0][j]    = stdev_pca[j]*shfs[Yk][j];
    }
    else{
      Yjac[0][NPCA]=1.;//YG_ta_c(0,Yi,Yj,Yk)+=YG1;
    }
  }
#else
  Yjac[0][NPCA]=1.;//YG_ta_c(0,Yi,Yj,Yk)+=YG1;
#endif

  //
}

//===========================================================================
//**************************  FIN DU MODULE ta_c1z **************************