/******************************************************************
 * FUNCTION: prntad - ideal Adiabatic printout of the var matrix  *
 *                    derivatives                                 *
 * FILENAME prntad.c                                              *
 * PROGRAMMER I.Urieli    (FORTRAN)                               *
 *            E. Malroy   (C Translation)                         *
 * DATE:      11/30/94  (FORTRAN), 03/13/96 (C Translation)       *
 * LAST MODIFIED: 08/18/96 by I. Urieli                           *
 * INCLUDE:  <stdio.h>,<math.h>,"adiabatic.h","define.h"          *
 * GLOBAL VARIABLES: infile, outfile, printfile                   *
 * PROTOTYPES: void prntad(double var[][COL], double dvar[][COL]);*
 * PRE:   the var and dvar matrices have all their values         *
 * POST:  various values have been printed out                    *
 ******************************************************************/
#include <stdio.h>
#include <math.h>
#include "define.h"
#include "adiabatic.h"

void prntad(double var[ROWV][COL], /* matrix of variables        */
            double dvar[ROWD][COL]) /* matrix of derivatives     */
{
    /* variables                                                 */
   double eff;     /* thermal efficency of refrigerator          */
   double cop;     /* coefficient of performance (refrigeration) */
   double vtotcc;  /* vtotcc: total volume(cc)                   */
   double pbar;    /* pbar: pressure(bar)                        */
   double m;       /* m: total mass of gas in machine(kg)        */
   double gak;     /* mass flow in cooler (kg/s)                 */
   double gar;     /* mass flow in regen  (kg/s)                 */
   double gah;     /* mass flow in heater (kg/s)                 */

   int i;          /* index                                      */
   int j;          /* index                                      */
   int deg;        /* deg:  cycle angle (degrees)                */

   eff = var[W][COL-1]/var[QH][COL-1];
   cop = var[QH][COL-1]/var[W][COL-1];
   printf("\nIdeal Adiabatic simulation results:\n");
   printf(" Qk (W) = %10.2f, Qr (W) = %10.3f, Qh (W) = %10.2f\n",
            var[QK][COL-1]*freq,var[QR][COL-1]*freq,var[QH][COL-1]*freq);
   printf(" Wc (W) = %10.2f, We (W) = %10.2f, W (W) = %10.2f\n",
            var[WC][COL-1]*freq,var[WE][COL-1]*freq,var[W][COL-1]*freq);
   printf(" eff (W/Qh) = %6.3lf, COP(Qh/W) = %6.3lf\n", eff, cop);
   fprintf(printfile,"\nIdeal Adiabatic simulation results:\n");
   fprintf(printfile," Qk (W) = %10.2f, Qr (W) = %10.3f, Qh (W) = %10.2f\n",
            var[QK][COL-1]*freq,var[QR][COL-1]*freq,var[QH][COL-1]*freq);
   fprintf(printfile," Wc (W) = %10.2f, We (W) = %10.2f, W (W) = %10.2f\n",
            var[WC][COL-1]*freq,var[WE][COL-1]*freq,var[W][COL-1]*freq);
   fprintf(printfile," eff (W/Qh) = %6.3f, COP(Qh/W) = %6.3f\n", eff, cop);
   printf("\n theta  V(cc)  P(bar)"
     "    Qk(J)    Qr(J)    Qh(J)    Wc(J)    We(J)     W(J)\n");
   fprintf(printfile,"\n theta  V(cc)  P(bar)"
     "    Qk(J)    Qr(J)    Qh(J)    Wc(J)    We(J)     W(J)\n");

   for(i = 0;i<37;i++) {
      deg = (i)*10;
      vtotcc = (var[VC][i] + vk + vr + vh + var[VE][i])*1e6;
      pbar = var[P][i]*1e-5;
      printf(" %4d%8.1f%8.2f",deg,vtotcc,pbar);
      fprintf(printfile," %4d%8.1f%8.2f", deg, vtotcc,pbar);
      for(j=QK;j<=W;j++) {
         printf("%9.1f", var[j][i]);
         fprintf(printfile,"%9.1f", var[j][i]);
      }
      printf("\n");
      fprintf(printfile, "\n");
   }
   printf("\n Vk(cc) =%7.1f Vr(cc) =%7.1f Vh(cc) =%7.1f\n",
          vk*1e6, vr*1e6, vh*1e6);
   fprintf(printfile,"\n Vk(cc) =%7.1f Vr(cc) =%7.1f Vh(cc) =%7.1f\n",
          vk*1e6, vr*1e6, vh*1e6);
   printf(" Tk(K) =%7.1f   Tr(K) =%7.1f   Th(K) =%7.1f\n",
          tk, tr, th);
   fprintf(printfile," Tk(K) =%7.1f   Tr(K) =%7.1f   Th(K) =%7.1f\n",
          tk, tr, th);
   printf("\n theta  Vc(cc)  Ve(cc)   V(cc)   Tc(K)   Te(K)\n");
   fprintf(printfile,"\n theta  Vc(cc)  Ve(cc)   V(cc)   Tc(K)   Te(K)\n");

   for(i = 0;i<COL;i++) {
      deg = (i)*10;
      vtotcc = (var[VC][i] + vk + vr + vh + var[VE][i])*1e6;
      printf(" %4d%8.1f%8.1f%8.1f%8.1f%8.1f\n",
      deg, var[VC][i]*1e6,var[VE][i]*1e6,vtotcc,var[TC][i],var[TE][i]);
      fprintf(printfile," %4d%8.1f%8.1f%8.1f%8.1f%8.1f\n",
      deg, var[VC][i]*1e6,var[VE][i]*1e6,vtotcc,var[TC][i],var[TE][i]);
   }
   printf("\n theta mc(gm) mk(gm) mr(gm) mh(gm) me(gm) "
       "M(gm) gAk(gm/s) gAr(gm/s) gAh(gm/s)\n");
   fprintf(printfile,"\n theta mc(gm) mk(gm) mr(gm) mh(gm) me(gm) "
       "M(gm) gAk(gm/s) gAr(gm/s) gAh(gm/s)\n");

   for(i = 0;i<COL;i++) {
      deg = (i)*10;
      m = 0.0;
      for(j = MC;j<=ME;j++)
         m = m + var[j][i];
      gak = (var[GACK][i] + var[GAKR][i])*omega*1e3/2.0;
      gar = (var[GAKR][i] + var[GARH][i])*omega*1e3/2.0;
      gah = (var[GARH][i] + var[GAHE][i])*omega*1e3/2.0;
      printf(" %5d", deg);
      fprintf(printfile," %5d", deg);
      for(j = MC;j<=ME;j++) {
         printf("%7.3f", var[j][i]*1e3);
         fprintf(printfile,"%7.3f", var[j][i]*1e3);
      }
      printf("%7.3f%10.3f%10.3f%10.3f\n", m*1e3, gak, gar, gah);
      fprintf(printfile,"%7.3f%10.3f%10.3f%10.3f\n", m*1e3, gak, gar, gah);
   }
}