Is it possible to extract the value of a variable from a template?

Dear all;

I am trying to simulate non-Newtonian fluid flow in porous media with the help of Palabos software. Palabos is a C+++ open-source library that solves flow problems with the lattice Boltzmann method. I am not familiar with object-oriented programming and I want to extract a value from a template that is used in every iteration. To be more precise; I want to print Omega (which is an (Nx-1)*(Ny-1) dimensional matrix) in every 250 iterations in the form of gif or data to validate my result. I added pcout to print variable in this function and it works but it is repeated every iteration which generates too much data ( imagine 200*200 lattice and 10000 iterations!).

I will be really glad if you could help me through this. the code is provided below. also, I uploaded

.h and .hh file to this message

Many thanks;

Elahe;

template<typename T, template<typename U> class Descriptor,int N>

T CarreauDynamics<T,Descriptor,N>::

getOmegaFromPiAndRhoBar(Array<T,SymmetricTensor<T,Descriptor>::n> const& PiNeq, T rhoBar) const

{

T nu0_nuInfoverCs2 = (global::CarreauParameters().getNu0()-global::CarreauParameters().getNuInf())*Descriptor<T>::invCs2;

T nuInfoverCs2 = global::CarreauParameters().getNuInf()*Descriptor<T>::invCs2;

T nMinusOneOverTwo = (global::CarreauParameters().getExponent() - (T)1)/(T)2;

T lambdaOverCs2sqr = global::CarreauParameters().getLambda()*Descriptor<T>::invCs2;

lambdaOverCs2sqr *= lambdaOverCs2sqr;

T piNeqNormSqr = SymmetricTensor<T,Descriptor>::tensorNormSqr(PiNeq);

T alpha = lambdaOverCs2sqr * piNeqNormSqr *(T)0.5

*Descriptor<T>::invRho(rhoBar)*Descriptor<T>::invRho(rhoBar);

T omega = carreauDynamicsTemplates<T,N>::fromPiAndRhoToOmega(alpha, nu0_nuInfoverCs2, nuInfoverCs2, nMinusOneOverTwo, this->getOmega());

T tau=(T)1/ omega;

T visco=((T)2 * tau - (T)1)/(T)6;

pcout << "" << visco << std::endl;

return omega;

}