integer NGridT
      double precision deltaT,GridT0,sumW(0:MaxGridT,maxprot),
     & sumE(0:MaxGridT,maxprot),sumEsq(0:MaxGridT,maxprot),
     & sumEbis(0:MaxGridT,maxprot),
     & sumQ(0:MaxGridT,maxprot),sumQsq(0:MaxGridT,maxprot),
     & sumRMS(0:MaxGridT,maxprot),sumRGY(0:MaxGridT,maxprot),
     & sumEQ(0:MaxGridT,maxprot),sumEEQ(0:MaxGridT,maxprot),
     & sumEEE(0:MaxGridT,maxprot),Rgas
      common /thermal_var/ deltaT,GridT0,sumW,sumE,sumEbis,sumEsq,
     & sumQ,sumRMS,sumRGY,sumQsq,sumEQ,sumEEE,sumEEQ,Rgas,nGridT