diff --git a/src/oops/assimilation/BMatrix.h b/src/oops/assimilation/BMatrix.h
index c0c5e53a..74a58047 100644
--- a/src/oops/assimilation/BMatrix.h
+++ b/src/oops/assimilation/BMatrix.h
@@ -16,6 +16,7 @@
 #include "oops/assimilation/ControlIncrement.h"
 #include "oops/assimilation/CostFunction.h"
 #include "oops/base/ObsAuxCovariances.h"
+#include "oops/util/Logger.h"
 
 namespace oops {
 
@@ -33,7 +34,13 @@ template<typename MODEL, typename OBS> class BMatrix : private boost::noncopyabl
  public:
   explicit BMatrix(const CostFct_ & j): j_(j) {}
   void multiply(const CtrlInc_ & x, CtrlInc_ & bx) const {
+//tothink
+    Log::trace() << "BMatrix:: multiply  start" <<std::endl;
+    Log::trace() << "BMatrix:: multiply  start bx " <<bx<<std::endl;
+    Log::trace() << "BMatrix:: multiply  start const x " <<x<<std::endl;
     j_.jb().multiplyB(x, bx);
+    Log::trace() << "BMatrix:: multiply  start after bx " <<bx<<std::endl;
+    Log::trace() << "BMatrix:: multiply  start after const x " <<x<<std::endl;
   }
 
   /// Return ObsBias block of control variable error covariances
diff --git a/src/oops/assimilation/CostFunction.h b/src/oops/assimilation/CostFunction.h
index 2e84d1e3..d1e58e03 100644
--- a/src/oops/assimilation/CostFunction.h
+++ b/src/oops/assimilation/CostFunction.h
@@ -280,7 +280,8 @@ double CostFunction<MODEL, OBS>::evaluate(CtrlVar_ & fguess,
 
 template<typename MODEL, typename OBS>
 void CostFunction<MODEL, OBS>::computeGradientFG(CtrlInc_ & grad) const {
-  Log::trace() << "CostFunction::computeGradientFG start" << std::endl;
+  Log::trace() << "CostFunction::computeGradientFG start clt" << std::endl;
+  Log::trace() << "CostFunction::computeGradientFG start clt grad.var " << grad.state().variables()<<std::endl;
   PostProcessor<Increment_> pp;
   PostProcessorTLAD<MODEL> costad;
   this->zeroAD(grad);
@@ -290,7 +291,9 @@ void CostFunction<MODEL, OBS>::computeGradientFG(CtrlInc_ & grad) const {
     jterms_[jj]->computeCostAD(tmp, grad, costad);
   }
 
+  Log::trace() << "CostFunction::computeGradientFG start clt  beforerunADJ grad  " << grad.state().variables()<<std::endl;
   this->runADJ(grad, costad, pp);
+  Log::trace() << "CostFunction::computeGradientFG start clt afterunADJ grad  " << grad.state().variables()<<std::endl;
 
   for (size_t jj = 0; jj < jterms_.size(); ++jj) {
     jterms_[jj]->setPostProcAD();
diff --git a/src/oops/assimilation/DRPCGMinimizer.h b/src/oops/assimilation/DRPCGMinimizer.h
index 9cee7dc3..322e31cd 100644
--- a/src/oops/assimilation/DRPCGMinimizer.h
+++ b/src/oops/assimilation/DRPCGMinimizer.h
@@ -172,14 +172,17 @@ double DRPCGMinimizer<MODEL, OBS>::solve(CtrlInc_ & dx, CtrlInc_ & dxh, CtrlInc_
   scals.push_back(1.0/dotRr0);
 
   Log::info() << std::endl;
+  Log::trace() << " DRPCG Starting Iteration0 0.0dxh" <<dxh<<  std::endl;
   for (int jiter = 0; jiter < maxiter; ++jiter) {
     Log::info() << " DRPCG Starting Iteration " << jiter+1 << std::endl;
+    Log::trace() << " DRPCG Starting Iteration0 " << jiter << std::endl;
 
     if (jiter < 5 || (jiter + 1) % 5 == 0 || jiter + 1 == maxiter) {
       util::update_workflow_meter("iteration", jiter+1);
       util::printRunStats("DRPCG iteration " + std::to_string(jiter+1));
     }
 
+    Log::trace() << " DRPCG Starting Iteration0 1" << jiter << std::endl;
     if (jiter > 0) {
       // beta_{i} = r_{i+1}^T z_{i+1} / r_{i}^T z_{i}
       double beta = rdots/rdots_old;
@@ -193,6 +196,8 @@ double DRPCGMinimizer<MODEL, OBS>::solve(CtrlInc_ & dx, CtrlInc_ & dxh, CtrlInc_
       hh += pr;
     }
 
+   Log::trace() << " DRPCG Starting Iteration0 0.7pp" <<pp<< jiter << std::endl;
+   Log::trace() << " DRPCG Starting Iteration0 0.7qq" <<qq<< jiter << std::endl;
     // q_{i} = h_{i} + H^T R^{-1} H p_{i}
     HtRinvH.multiply(pp, qq);
     qq += hh;
@@ -202,10 +207,16 @@ double DRPCGMinimizer<MODEL, OBS>::solve(CtrlInc_ & dx, CtrlInc_ & dxh, CtrlInc_
     double alpha = rdots/rho;
 
     // dx_{i+1} = dx_{i} + alpha * p_{i}
+    Log::trace() << " DRPCG Starting Iteration0 1.9dx" <<dx<< jiter << std::endl;
+    Log::trace() << " DRPCG Starting Iteration0 1.9pp" <<pp<< jiter << std::endl;
     dx.axpy(alpha, pp);
     // dxh_{i+1} = dxh_{i} + alpha * h_{i} ! for diagnosing Jb
+    Log::trace() << " DRPCG Starting Iteration0 1.9hh" <<hh<< jiter << std::endl;
     dxh.axpy(alpha, hh);
+    Log::trace() << " DRPCG Starting Iteration0 1.9dxh" <<dxh<< jiter << std::endl;
     // r_{i+1} = r_{i} - alpha * q_{i}
+    Log::trace() << " DRPCG Starting Iteration0 1.9qq" <<qq<< jiter << std::endl;
+    Log::trace() << " DRPCG Starting Iteration0 1.9rr" <<rr<< jiter << std::endl;
     rr.axpy(-alpha, qq);
 
     // Compute the quadratic cost function
@@ -215,6 +226,7 @@ double DRPCGMinimizer<MODEL, OBS>::solve(CtrlInc_ & dx, CtrlInc_ & dxh, CtrlInc_
     double costJb = costJ0Jb + dot_product(dx, gradJb) + 0.5 * dot_product(dx, dxh);
     // Jo[dx_{i}] + Jc[dx_{i}] = J[dx_{i}] - Jb[dx_{i}]
     double costJoJc = costJ - costJb;
+    Log::trace() << " DRPCG Starting Iteration0 2.0rr" <<rr<< jiter << std::endl;
 
     // Re-orthogonalization
     for (int jj = 0; jj < jiter; ++jj) {
@@ -222,9 +234,19 @@ double DRPCGMinimizer<MODEL, OBS>::solve(CtrlInc_ & dx, CtrlInc_ & dxh, CtrlInc_
       rr.axpy(-proj, vvecs[jj]);
     }
 
+    Log::trace() << " DRPCG Starting Iteration0 2.5" << jiter << std::endl;
+    Log::trace() << " DRPCG Starting Iteration0 2.5pr" <<pr<< jiter << std::endl;
+    Log::trace() << " DRPCG Starting Iteration0 2.5rr" <<rr<< jiter << std::endl;
     // z_{i+1} = B LMP r_{i+1}
     lmp_.multiply(rr, pr);
+    Log::trace() << " DRPCG Starting Iteration0 3pr" << jiter << std::endl;
+    Log::trace() << pr << jiter << std::endl;
+    Log::trace() << " DRPCG Starting Iteration0 3zz" << jiter << std::endl;
+    Log::trace() << zz << jiter << std::endl;
+
+    Log::trace() << " DRPCG Starting Iteration0 3Before tothink B.multiply" << jiter << std::endl;
     B.multiply(pr, zz);
+    Log::trace() << " DRPCG Starting Iteration0 3after B.multiply" << jiter << std::endl;
 
     // r_{i}^T z_{i}
     rdots_old = rdots;
diff --git a/src/oops/base/EnsembleCovariance.h b/src/oops/base/EnsembleCovariance.h
index 8ca5e3ab..7230e8c3 100644
--- a/src/oops/base/EnsembleCovariance.h
+++ b/src/oops/base/EnsembleCovariance.h
@@ -89,6 +89,11 @@ EnsembleCovariance<MODEL>::EnsembleCovariance(const Geometry_ & resol, const Var
   ens_.reset(new IncrementSet<MODEL>(resol, vars, tmp, true));
   *ens_ -= ens_->ens_mean();
 
+ for (size_t i = 0; i < ens_->size(); ++i) {  //cltthink
+     auto member = (*ens_)[i];
+     Log::trace()<<"EnsembleCovariance::EnsembleCovariance start ith dx "<<i <<" "<<member<<std::endl;
+}
+
   // Setup ensemble transform
   if (conf.has("ensemble transform")) {
     const eckit::LocalConfiguration confEns(conf, "ensemble transform");
@@ -159,6 +164,7 @@ EnsembleCovariance<MODEL>::EnsembleCovariance(const Geometry_ & resol, const Var
     }
     loc_.reset(new Localization_(resol, locVars, locconf));
   }
+  Log::trace() << "EnsembleCovariance::EnsembleCovariance vars " <<vars<< std::endl;
 
   size_t current = eckit::system::ResourceUsage().maxResidentSetSize();
   this->setObjectSize(current - init);
diff --git a/src/oops/base/GetValuePosts.h b/src/oops/base/GetValuePosts.h
index 8e8dd2c8..138697e4 100644
--- a/src/oops/base/GetValuePosts.h
+++ b/src/oops/base/GetValuePosts.h
@@ -31,6 +31,8 @@ class GetValuesParameters : public oops::Parameters {
  public:
   Parameter<eckit::LocalConfiguration> variableChange{"variable change",
                                                       eckit::LocalConfiguration(), this};
+  Parameter<eckit::LocalConfiguration> linearvariableChange{"linear variable change",
+                                                      eckit::LocalConfiguration(), this};
 };
 
 /// \brief Fills GeoVaLs with requested variables at requested locations during model run
diff --git a/src/oops/base/GetValueTLADs.h b/src/oops/base/GetValueTLADs.h
index 4241191f..e36c6340 100644
--- a/src/oops/base/GetValueTLADs.h
+++ b/src/oops/base/GetValueTLADs.h
@@ -64,6 +64,7 @@ class GetValueTLADs : public PostBaseTLAD<MODEL> {
   std::vector<GetValPtr_> getvals_;
   std::map<util::DateTime, CVarPtr_> chvartlad_;
   const eckit::LocalConfiguration cvConf_;
+  const eckit::LocalConfiguration linearcvConf_;
 };
 
 // -----------------------------------------------------------------------------
@@ -71,8 +72,9 @@ template <typename MODEL, typename OBS>
 GetValueTLADs<MODEL, OBS>::GetValueTLADs(const GetValuesParameters<MODEL> & params,
               const util::DateTime & bgn, const util::DateTime & end)
   : PostBaseTLAD<MODEL>(bgn, end), geovars_(), linvars_(), getvals_(), chvartlad_(),
-    cvConf_(params.variableChange.value())
+    linearcvConf_(params.linearvariableChange.value()),cvConf_(params.variableChange.value())
 {
+  Log::trace() << "GetValueTLADs::GetValueTLADs clt linarcvConf_ " <<linearcvConf_ <<std::endl;
   Log::trace() << "GetValueTLADs::GetValueTLADs" << std::endl;
 }
 // -----------------------------------------------------------------------------
@@ -82,6 +84,8 @@ void GetValueTLADs<MODEL, OBS>::append(GetValPtr_ getval) {
   getvals_.push_back(getval);
   geovars_ += getval->requiredVariables();
   linvars_ += getval->linearVariables();
+  Log::trace() << "GetValuePosts::append start geovars_ " <<geovars_<< std::endl;
+  Log::trace() << "GetValuePosts::append start linvars_ " <<linvars_<< std::endl;
   Log::trace() << "GetValuePosts::append done" << std::endl;
 }
 
@@ -110,7 +114,7 @@ void GetValueTLADs<MODEL, OBS>::doProcessingTraj(const State_ & xx) {
 
   for (GetValPtr_ getval : getvals_) getval->process(zz);
 
-  CVarPtr_ cvtlad(new LinearVariableChange<MODEL>(xx.geometry(), cvConf_));
+  CVarPtr_ cvtlad(new LinearVariableChange<MODEL>(xx.geometry(), linearcvConf_));
   cvtlad->changeVarTraj(xx, linvars_);
   chvartlad_[xx.validTime()] = std::move(cvtlad);
 
@@ -170,6 +174,7 @@ void GetValueTLADs<MODEL, OBS>::doProcessingAD(Increment_ & dx) {
   Log::trace() << "GetValueTLADs::doProcessingAD start" << std::endl;
   const util::DateTime now = dx.validTime();
   ASSERT(chvartlad_.find(now) != chvartlad_.end());
+  Log::trace() << "GetValueTLADs::doProcessingAD linearvars" <<linvars_<<std::endl;
   Increment_ dz(dx.geometry(), linvars_, dx.validTime());
   dz.zero();
 
diff --git a/src/oops/base/GetValues.h b/src/oops/base/GetValues.h
index 2508b568..085ba7e0 100644
--- a/src/oops/base/GetValues.h
+++ b/src/oops/base/GetValues.h
@@ -690,6 +690,7 @@ template <typename MODEL, typename OBS>
 void GetValues<MODEL, OBS>::processAD(Increment_ & dx) {
   Log::trace() << "GetValues::processAD start" << std::endl;
   util::Timer timer("oops::GetValues", "processAD");
+  Log::trace() << "GetValues::processAD startclt variables " <<dx.variables()<< std::endl;
 
   for (size_t jtask = 0; jtask < ntasks_; ++jtask) {
 //  Mask obs outside time slot
diff --git a/src/oops/base/ObserversTLAD.h b/src/oops/base/ObserversTLAD.h
index c642096f..c02e917e 100644
--- a/src/oops/base/ObserversTLAD.h
+++ b/src/oops/base/ObserversTLAD.h
@@ -112,6 +112,7 @@ template <typename MODEL, typename OBS>
 void ObserversTLAD<MODEL, OBS>::initializeTraj(const Geometry_ & geom, const ObsAuxCtrls_ & ybias,
                                                PostProcTLAD_ & pp) {
   Log::trace() << "ObserversTLAD<MODEL, OBS>::initializeTraj start" << std::endl;
+  Log::trace() << "ObserversTLAD<MODEL, OBS>::initializeTraj clt params " << getValuesParams_<<std::endl;
   posts_.reset(new GetValueTLADs_(getValuesParams_, winbgn_, winend_));
   for (size_t jj = 0; jj < observers_.size(); ++jj) {
     if (observers_[jj]) {
diff --git a/src/oops/interface/Increment.h b/src/oops/interface/Increment.h
index 767b0418..426d582f 100644
--- a/src/oops/interface/Increment.h
+++ b/src/oops/interface/Increment.h
@@ -155,6 +155,7 @@ Increment<MODEL>::Increment(const Geometry_ & resol, const Variables & vars,
 {
   Log::trace() << "Increment<MODEL>::Increment starting" << std::endl;
   util::Timer timer(classname(), "Increment");
+  Log::trace() << "Increment<MODEL>::Increment vars " <<vars<< std::endl;
   increment_.reset(new Increment_(resol.geometry(), vars, time));
   this->setObjectSize(increment_->serialSize()*sizeof(double));
   Log::trace() << "Increment<MODEL>::Increment done" << std::endl;