trunk/models/small_components.a4l

REQUIRE "atoms.a4l";
(* => atoms.a4l, measures.a4l, system.a4l, basemodel.a4l *)
PROVIDE "components.a4l";
(*
 *
 * WARNINGS:
 * - The contents of this file represent a physical
 * properties database of minimal functionality. Such a
 * database is just messy conceptually in a declarative language,
 * so we do not apologize for the modeling style found in this
 * file.
 *
 * - We, of course, would like to replace it with a small wrapper to a
 * thermodynamic information database with a much wider range
 * of chemical species.  We are aware of no such database of
 * significant size being placed in the public domain
 * in source form which we can distribute to all our users.
 *
 * 1:
 * Add species or new correlation coefficients as you need them.
 * Models you create in the ASCEND IV language are Not subject to
 * the GNU Public License (GPL) UNLESS you base those models on GNU Public
 * Licensed ASCEND Libraries. All ASCEND libraries distributed from
 * Carnegie Mellon are distributed under the GPL unless explicitly noted
 * as being in the public domain in the distributed source code.
 *
 * In the interest of promoting research, we make an EXCEPTION to the
 * above condition for proprietary physical property data integrated
 * with ASCEND library models for non-commercial research purposes only.
 *
 * If you base a commercial application on GNU Public Licensed ASCEND IV
 * libraries or modifications or extensions of those libraries, then the
 * models you create must be released in source code form per the GPL.
 *
 * 2: If you have such a properties database to donate, please let us know.
 * 3:
 * If you create a wrapper to a proprietary database for use with
 * ASCEND in a way that entangles your code with our sources, you
 * are required under the terms of the ASCEND GPL
 * to _give_ the wrapper code back to us and to make it otherwise
 * available for public use. This does not require making the
 * proprietary database public, just the interface. If releasing
 * such an interface violates the proprietary licensing, then do not
 * create it in an entangled fashion.
 * 4:
 * If you have questions about any of the above, please contact us
 * ascend+developers@cs.cmu.edu and aw0a@cs.cmu.edu. We will consider
 * alternative licensing arrangements on a case-by-case basis subject
 * to keeping the lawyers and accountants on all sides happy.
 *)

(*
 *  components.a4l
 *  by Joseph J. Zaher and Ben Allan
 *  Part of the ASCEND Library
 *  $Date: 1998/08/10 16:03:36 $
 *  $Revision: 1.1 $
 *  $Author: ballan $
 *  $Source: /afs/cs.cmu.edu/project/ascend/Repository/models/small_components.a4l,v $
 *
 *  This file is part of the ASCEND Modeling Library.
 *
 *  Copyright (C) 1994 Joseph J Zaher
 *  Copyright (C) 1997 Benjamin Andrew Allan
 *
 *  The ASCEND Modeling Library is free software; you can redistribute
 *  it and/or modify it under the terms of the GNU General Public
 *  License as published by the Free Software Foundation; either
 *  version 2 of the License, or (at your option) any later version.
 *
 *  The ASCEND Modeling Library is distributed in hope that it
 *  will be useful, but WITHOUT ANY WARRANTY; without even the implied
 *  warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
 *  See the GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with this program.  If not, see <http://www.gnu.org/licenses/>.
 *)

(**
 **
 **  C  O  M  P  O  N  E  N  T  S  .  A 4 L
 **  ----------------------------------------
 **
 **  AUTHOR:	  Joseph J. Zaher
 **
 **  DATES:	   07/91 - Original code.
 **		   02/92 - Made compatible with new set version of ASCEND.
 **			   Expanded data base, revised vapor pressure data,
 **			   and added UNIFAC group and subgroup sets with help
 **			   of Bob Huss.
 **		   07/92 - Expanded data base with help of Kay C. Dee.
 **		   08/92 - Replaced name attribute of each component with a
 **		   	   formula attribute.  Component names are free to be
 **			   specified by the user.
 **		   03/94 - Made compatible with gnu-ascend.
 **		   08/94 - Made compatible with H,G or H,S thermo library,
 **			   and with Wilson models written by Boyd Safrit.
 **		  02/96 - Made compatible with ASCEND IV. Ben Allan.
 **		  01/97 - Addition OF td_component_constants by _____
 **		  01/97 - All constants checked by Jennifer Perry
 **		   03/97 - Added container models for data groups. BAA.
 **		  03/98 - 99% of the code in this file should be reduced
 **		 		to a call on a database.
 **
 **
 **  CONTENTS:	ASCEND structure for component physical property constants.
 **		   All anticipated constants which may be used by the models
 **		   of "thermodynamics.a4l" are created in a general model where
 **		   a consistent reference state (298.15{K} and 1.0{atm}) is
 **		   chosen.  Specific refinements to actual chemical species are
 **		   made UNIVERSAL to ensure only one instantiation of the
 **		   constants for each component is maintained throughout a
 **		   user's simulation.  A reference enthalpy and entropy is
 **		   chosen by default to be those of formation of the component
 **		   from its elements at the reference state.
 **
 **
 **)

(* Wilson pairs checked by Boyd
i_propanol-water
i_propanol-ethylene_glycol
ethylene_glycol-water
acetone-chloroform
acetone-benzene
acetone-methanol
methanol-water
benzene-chloroform
propadiene-propylene	provided by Vince Verneulli OF Sim Sci
propadiene-propane	provided by Vince Verneulli OF Sim Sci
propylene-propane	provided by Vince Verneulli OF Sim Sci
*)

MODEL compmodel() REFINES cmumodel();
METHODS
METHOD check_self;
END check_self;
METHOD check_all;
END check_all;
METHOD default_self;
END default_self;
METHOD default_all;
END default_all;
METHOD specify;
END specify;
METHOD scale_self;
END scale_self;
METHOD scale_all;
END scale_all;
METHOD bound_self;
END bound_self;
METHOD bound_all;
END bound_all;
END compmodel;

UNIVERSAL MODEL UNIFAC_constants() REFINES compmodel;

(*
REFERENCE:
The reference for the bulk of this MODEL is
Reid, Prausnitz & Poling, Properties of Gases and Liquids,
1986, Chapter 8.
DISCLAIMER: DISCLAIMER: DISCLAIMER: DISCLAIMER: DISCLAIMER: DISCLAIMER:
The authors of ASCEND and Carnegie Mellon University make
absolutely NO WARRANTY about the accuracy of this transcription
of the RPP data or of the original data itself, nor do they
provide any guarantee that the data here represented is
suitable for any purpose academic or commercial.
DISCLAIMER: DISCLAIMER: DISCLAIMER: DISCLAIMER: DISCLAIMER: DISCLAIMER:
NOTES:
The frontiers of group contribution methods have moved on
well past the data presented here and the correlations
it is used with. Users with a penchant for thermodynamic
detail are advised to MODEL their own favorite correlations.
COMMON USAGE:
The local name of every instance of this MODEL
(since it is UNIVERSAL) should always be uc.
In the comments that follow, we will use uc where applicable.
ANNOTATED:
10/96 Ben Allan
*)


	groups		IS_A set OF integer_constant;
(*
* uc.groups is the main group list (the ones we support) from RPP
* Table 8.21 column 1.
* Currently supported groups are (in RPP Main index)
* 1,2,3,4,5,6,7,8,9,23, and 11
*)
	groups	   :== [1..47];
(* fix a 10, 12-22, 24-47
*)

(*
* uc.sub[i IN uc.groups] is the set of secondary groups for the ith
* main group. This is the column 2 data of RPP Table 8-21.
* The UNION of uc.sub[i] is uc.subgroups, basically the set of
* rows in Table 8-21 that we have entered.
*)
	sub[groups]	 IS_A set OF symbol_constant;
	sub[1]   :== ['CH3', 'CH2', 'CH', 'C'];
	sub[2]   :== ['CH2=CH', 'CH=CH', 'CH2=C', 'CH=C', 'C=C'];
	sub[3]   :== ['ACH', 'AC'];
	sub[4]   :== ['ACCH3', 'ACCH2', 'ACCH'];
	sub[5]   :== ['OH'];
	sub[6]   :== ['CH3OH'];
	sub[7]   :== ['H2O'];
	sub[8]   :== ['ACOH'];
	sub[9]   :== ['CH3CO', 'CH2CO'];
	sub[10]  :== ['CHO'];
	sub[11]  :== ['CH3COO', 'CH2COO'];
	sub[12]  :== ['HCOO'];
	sub[13]  :== ['CH3O', 'CH2O', 'CH-O', 'FCH2O'];
	sub[14]  :== ['CH3NH2', 'CH2NH2', 'CHNH2'];
	sub[15]  :== ['CH3NH', 'CH2NH', 'CHNH'];
	sub[16]  :== ['CH3N', 'CH2N'];
	sub[17]  :== ['ACNH2'];
	sub[18]  :== ['C5H5N', 'C5H4N', 'C5H3N'];
	sub[19]  :== ['CH3CN', 'CH2CN'];
	sub[20]  :== ['COOH', 'HCOOH'];
	sub[21]  :== ['CH2Cl', 'CHCl', 'CCl'];
	sub[22]  :== ['CH2Cl2', 'CHCl2', 'CCl2'];
	sub[23]  :== ['CHCl3', 'CCl3'];
	sub[24]  :== ['CCl4'];
	sub[25]  :== ['ACCl'];
	sub[26]  :== ['CH3NO2', 'CH2NO2', 'CHNO2'];
	sub[27]  :== ['ACNO2'];
	sub[28]  :== ['CS2'];
	sub[29]  :== ['CH3SH', 'CH2SH'];
	sub[30]  :== ['Furfural'];
	sub[31]  :== ['(CH2OH)2'];
	sub[32]  :== ['I'];
	sub[33]  :== ['Br'];
	sub[34]  :== ['CH-trip-C', 'C-trip-C'];
	sub[35]  :== ['Me2SO'];
	sub[36]  :== ['Acry'];
	sub[37]  :== ['Cl(C=C)'];
	sub[38]  :== ['ACF'];
	sub[39]  :== ['DMF-1', 'DMF-2'];
	sub[40]  :== ['CF3', 'CF2', 'CF'];
	sub[41]  :== ['COO'];
	sub[42]  :== ['SiH3', 'SiH2', 'SiH', 'Si'];
	sub[43]  :== ['SiH2O', 'SiHO', 'SiO'];
	sub[44]  :== ['NMP'];
	sub[45]  :== ['tert-N'];
	sub[46]  :== ['Amide'];
	sub[47]  :== ['CON(Me)2', 'CONMeCH2', 'CON(CH2)2'];
	subgroups	   IS_A set OF symbol_constant;

(*
 * uc.	subgroups is the total list of secondary groups we accomodate.
 *)
	subgroups			:== UNION[sub[i] | i IN groups];
	group[subgroups]		IS_A integer_constant;
(*
 * Next we have a horribly named array, group, that lets us look up the
 * main group given a particular subgroup. So, for example,
 * if we want to know the main group of 'CHCl3', then
 * group['CHCl3'] will give it to us: 'CCl3'.
 * Note, then, that the assignments which follow in most
 * cases expand to have several left-hand-sides.
 *)
	FOR i IN groups CREATE
	group[sub[i]]	:== i;
	END FOR;

	a[groups][groups] "binary group interaction parameters" IS_A UNIFAC_a;
(*
 * Now we have uc.a[m][n], a portion of RPP Table 8-22, the main group
 * interaction matrix which is asymmetric and in principle
 * dense but in our practice rather sparse.
 *   Adding a main group to this MODEL means adding a row and
 * column to uc.a, which is a non-trivial exercise until we
 * get a proper tabular assignment syntax implemented.
 * Note that such an syntax has been completely scoped but
 * we've just not yet got round to implementation.
 * When we do implement, this will be the first test MODEL.
 * a[i][j] should have i and j being integer_constant!
 *)

(* we would like the ascend TABLE syntax or a database
 * external lookup to handle this data entry.
 *)
(* row 1 *)
	a[1][1]		:==	0.0	{K};
	a[1][2]		:==	86.020	{K};
	a[1][3]		:==	61.13	{K};
	a[1][4]		:==	76.50	{K};
	a[1][5]		:==	986.5	{K};
	a[1][6]		:==	697.2	{K};
	a[1][7]		:==	1318.0	{K};
	a[1][8]		:==	1333.0	{K};
	a[1][9]		:==	476.4	{K};
	a[1][10]	:==	677	{K};
	a[1][11]	:==	232.1	{K};
	a[1][12]	:==	741.4	{K};
	a[1][13]	:==	251.5	{K};
	a[1][14]	:==	391.5	{K};
	a[1][15]	:==	225.7	{K};
	a[1][16]	:==	206.6	{K};
	a[1][17]	:==	920.7	{K};
	a[1][18]	:==	287.7	{K};
	a[1][19]	:==	597	{K};
	a[1][20]	:==	663.5	{K};
	a[1][21]	:==	35.93	{K};
	a[1][22]	:==	53.76	{K};
	a[1][23]	:==	24.90	{K};
	a[1][24]	:==	104.3	{K};
	a[1][25]	:==	321.5	{K};
	a[1][26]	:==	661.5	{K};
	a[1][27]	:==	543	{K};
	a[1][28]	:==	153.6	{K};
	a[1][29]	:==	184.4	{K};
	a[1][30]	:==	354.5	{K};
	a[1][31]	:==	3025	{K};
	a[1][32]	:==	335.8	{K};
	a[1][33]	:==	479.5	{K};
	a[1][34]	:==	298.9	{K};
	a[1][35]	:==	526.5	{K};
	a[1][36]	:==	689	{K};
	a[1][37]	:==	-4.189	{K};
	a[1][38]	:==	125.8	{K};
	a[1][39]	:==	485.3	{K};
	a[1][40]	:==	-2.859	{K};
	a[1][41]	:==	387.1	{K};
	a[1][42]	:==	-450.4	{K};
	a[1][43]	:==	252.7	{K};
	a[1][44]	:==	13.89	{K};
	a[1][45]	:==	383	{K};
	a[1][46]	:==	-1380	{K};
	a[1][47]	:==	729	{K};

(* row 2 *) (* the ?!? numbers in group 2 are bogus. probably entered by ka.
 * they have been replaced. part looks like a paste of row 4.
 *)
(* old:
	a[2][1]		:==	2520	{K}; (*?!?*)
	a[2][2]		:==	0.0	{K}; (*?!?*)
	a[2][3]		:==	340.7	{K}; (*?!?*)
	a[2][4]		:==	4102	{K}; (*?!?*)
	a[2][5]		:==	693.9	{K}; (*?!?*)
	a[2][6]		:==	1509.0	{K}; (*?!?*)
	a[2][7]		:==	634.2	{K}; (*?!?*)
	a[2][8]		:==	547.4	{K}; (*?!?*)
	a[2][9]		:==	524.5	{K}; (*?!?*)
	a[2][11]	:==	71.23	{K}; (*?!?*)
 END old *)
	a[2][1]		:==	-35.36	{K};
	a[2][2]		:==	0.0	{K};
	a[2][3]		:==	38.81	{K};
	a[2][4]		:==	74.15	{K};
	a[2][5]		:==	524.1	{K};
	a[2][6]		:==	787.6	{K};
	a[2][7]		:==	270.6	{K};
	a[2][8]		:==	526.1	{K};
	a[2][9]		:==	182.6	{K};
	a[2][10]	:==	448.8	{K};
	a[2][11]	:==	37.85	{K};
	a[2][12]	:==	449.1	{K};
	a[2][13]	:==	214.5	{K};
	a[2][14]	:==	240.9	{K};
	a[2][15]	:==	163.9	{K};
	a[2][16]	:==	61.11	{K};
	a[2][17]	:==	749.3	{K};
	a[2][18]	:==	0	{K};
	a[2][19]	:==	336.9	{K};
	a[2][20]	:==	318.9	{K};
	a[2][21]	:==	204.6	{K};
	a[2][22]	:==	5.892	{K};
(* old
	a[2][23]	:==	4584.0	{K}; (*?!?*)
   END old *)

	a[2][23]	:==	-13.99	{K};
	a[2][24]	:==	-109.7	{K};
	a[2][25]	:==	393.1	{K};
	a[2][26]	:==	357.5	{K};
	a[2][27]	:==	0	{K};
	a[2][28]	:==	76.3	{K};
	a[2][29]	:==	0	{K};
	a[2][30]	:==	0	{K};
	a[2][31]	:==	0	{K};
	a[2][32]	:==	0	{K};
	a[2][33]	:==	0	{K};
	a[2][34]	:==	31.14	{K};
	a[2][35]	:==	-137.4	{K};
	a[2][36]	:==	0	{K};
	a[2][37]	:==	-66.46	{K};
	a[2][38]	:==	0	{K};
	a[2][39]	:==	-70.45	{K};
	a[2][40]	:==	0	{K};
	a[2][41]	:==	48.33	{K};
	a[2][42]	:==	0	{K};
	a[2][43]	:==	0	{K};
	a[2][44]	:==	0	{K};
	a[2][45]	:==	0	{K};
	a[2][46]	:==	2340	{K};
	a[2][47]	:==	0	{K};
(* row 3 *)
	a[3][1]		:==	-11.12	{K};
	a[3][2]		:==	3.446	{K};
	a[3][3]		:==	0.0	{K};
	a[3][4]		:==	167.0	{K};
	a[3][5]		:==	636.10	{K};
	a[3][6]		:==	637.3	{K};
	a[3][7]		:==	903.8	{K};
	a[3][8]		:==	1329.0	{K};
	a[3][9]		:==	25.77	{K};
	a[3][11]	:==	5.994	{K};
	a[3][10]	:==	347.3	{K};
	a[3][12]	:==	-92.55	{K};
	a[3][13]	:==	32.14	{K};
	a[3][14]	:==	161.7	{K};
	a[3][15]	:==	122.8	{K};
	a[3][16]	:==	90.49	{K};
	a[3][17]	:==	648.2	{K};
	a[3][18]	:==	-4.449	{K};
	a[3][19]	:==	212.5	{K};
	a[3][20]	:==	537.4	{K};
	a[3][21]	:==	-18.81	{K};
	a[3][22]	:==	-144.4	{K};
	a[3][23]	:==	-231.9	{K};
	a[3][24]	:==	3	{K};
	a[3][25]	:==	538.2	{K};
	a[3][26]	:==	168	{K};
	a[3][27]	:==	194.9	{K};
	a[3][28]	:==	52.07	{K};
	a[3][29]	:==	-10.43	{K};
	a[3][30]	:==	-64.69	{K};
	a[3][31]	:==	210.4	{K};
	a[3][32]	:==	113.3	{K};
	a[3][33]	:==	-13.59	{K};
	a[3][34]	:==	0	{K};
	a[3][35]	:==	169.9	{K};
	a[3][36]	:==	0	{K};
	a[3][37]	:==	-259.1	{K};
	a[3][38]	:==	389.3	{K};
	a[3][39]	:==	245.6	{K};
	a[3][40]	:==	0	{K};
	a[3][41]	:==	103.5	{K};
	a[3][42]	:==	-432.3	{K};
	a[3][43]	:==	238.9	{K};
	a[3][44]	:==	-23.88	{K};
	a[3][45]	:==	109	{K};
	a[3][46]	:==	75.9	{K};
	a[3][47]	:==	784	{K};
(* row 4 *)
	a[4][1]		:==	-69.70	{K};
	a[4][2]		:==	-113.6	{K};
	a[4][3]		:==	-146.80	{K};
	a[4][4]		:==	0.0	{K};
	a[4][5]		:==	803.20	{K};
	a[4][6]		:==	603.2	{K};
	a[4][7]		:==	5695.00	{K};
	a[4][8]		:==	547.4	{K};
	a[4][9]		:==	-52.10	{K};
	a[4][10]	:==	586.6	{K};
	a[4][11]	:==	5688.0	{K};

	a[4][12]	:==	115.2	{K};
	a[4][13]	:==	213.1	{K};
	a[4][14]	:==	0	{K};
	a[4][15]	:==	49.29	{K};
	a[4][16]	:==	23.5	{K};
	a[4][17]	:==	664.2	{K};
	a[4][18]	:==	52.8	{K};
	a[4][19]	:==	6096	{K};
	a[4][20]	:==	603.8	{K};
	a[4][21]	:==	-114.1	{K};
	a[4][22]	:==	-111	{K};
	a[4][23]	:==	-12.14	{K};
	a[4][24]	:==	-141.3	{K};
	a[4][25]	:==	-126.9	{K};
	a[4][26]	:==	3629	{K};
	a[4][27]	:==	4448	{K};
	a[4][28]	:==	-9.451	{K};
	a[4][29]	:==	0	{K};
	a[4][30]	:==	-20.36	{K};
	a[4][31]	:==	4975	{K};
	a[4][32]	:==	0	{K};
	a[4][33]	:==	-171.3	{K};
	a[4][34]	:==	0	{K};
	a[4][35]	:==	4284	{K};
	a[4][36]	:==	0	{K};
	a[4][37]	:==	0	{K};
	a[4][38]	:==	101.4	{K};
	a[4][39]	:==	5629	{K};
	a[4][40]	:==	0	{K};
	a[4][41]	:==	69.26	{K};
	a[4][42]	:==	683.3	{K};
	a[4][43]	:==	355.5	{K};
	a[4][44]	:==	6.214	{K};
	a[4][45]	:==	1320	{K};
	a[4][46]	:==	482	{K};
	a[4][47]	:==	386	{K};
(* row 5 *)
	a[5][1]		:==	156.40	{K};
	a[5][2]		:==	475.000	{K};
	a[5][3]		:==	89.60	{K};
	a[5][4]		:==	25.82	{K};
	a[5][5]		:==	0.0	{K};
	a[5][6]		:==	-137.1	{K};
	a[5][7]		:==	353.50	{K};
	a[5][8]		:==	-259.7	{K};
	a[5][9]		:==	84.0	{K};
	a[5][10]	:==	441.8	{K};
	a[5][11]	:==	101.1	{K};
	a[5][12]	:==	193.1	{K};
	a[5][13]	:==	28.06	{K};
	a[5][14]	:==	83.02	{K};
	a[5][15]	:==	42.7	{K};
	a[5][16]	:==	-323	{K};
	a[5][17]	:==	-52.39	{K};
	a[5][18]	:==	170	{K};
	a[5][19]	:==	6.712	{K};
	a[5][20]	:==	199	{K};
	a[5][21]	:==	75.62	{K};
	a[5][22]	:==	-112.1	{K};
	a[5][23]	:==	-98.12	{K};
	a[5][24]	:==	143.1	{K};
	a[5][25]	:==	287.8	{K};
	a[5][26]	:==	256.5	{K};
	a[5][27]	:==	157.1	{K};
	a[5][28]	:==	477	{K};
	a[5][29]	:==	147.5	{K};
	a[5][30]	:==	-120.5	{K};
	a[5][31]	:==	-318.9	{K};
	a[5][32]	:==	313.5	{K};
	a[5][33]	:==	133.4	{K};
	a[5][34]	:==	0	{K};
	a[5][35]	:==	-202.1	{K};
	a[5][36]	:==	0	{K};
	a[5][37]	:==	225.8	{K};
	a[5][38]	:==	44.78	{K};
	a[5][39]	:==	-143.9	{K};
	a[5][40]	:==	0	{K};
	a[5][41]	:==	190.3	{K};
	a[5][42]	:==	-817.7	{K};
	a[5][43]	:==	202.7	{K};
	a[5][44]	:==	796.9	{K};
	a[5][45]	:==	0	{K};
	a[5][46]	:==	0	{K};
	a[5][47]	:==	0	{K};
(* row 6 *)
	a[6][1]		:==	16.51	{K};
	a[6][2]		:==	-12.520	{K};
	a[6][3]		:==	-50.00	{K};
	a[6][4]		:==	-44.50	{K};
	a[6][5]		:==	249.1	{K};
	a[6][6]		:==	0.0	{K};
	a[6][7]		:==	-181.0	{K};
	a[6][8]		:==	-101.7	{K};
	a[6][9]		:==	23.39	{K};
	a[6][10]	:==	306.4	{K};
	a[6][11]	:==	-10.72	{K};
	a[6][12]	:==	193.4	{K};
	a[6][13]	:==	-128.6	{K};
	a[6][14]	:==	359.3	{K};
	a[6][15]	:==	266	{K};
	a[6][16]	:==	53.9	{K};
	a[6][17]	:==	489.7	{K};
	a[6][18]	:==	580.5	{K};
	a[6][19]	:==	36.23	{K};
	a[6][20]	:==	-289.5	{K};
	a[6][21]	:==	-38.32	{K};
	a[6][22]	:==	-102.5	{K};
	a[6][23]	:==	-139.4	{K};
	a[6][24]	:==	-67.8	{K};
	a[6][25]	:==	17.12	{K};
	a[6][26]	:==	75.14	{K};
	a[6][27]	:==	0	{K};
	a[6][28]	:==	-31.09	{K};
	a[6][29]	:==	37.84	{K};
	a[6][30]	:==	0	{K};
	a[6][31]	:==	0	{K};
	a[6][32]	:==	0	{K};
	a[6][33]	:==	106.3	{K};
	a[6][34]	:==	0	{K};
	a[6][35]	:==	-399.3	{K};
	a[6][36]	:==	0	{K};
	a[6][37]	:==	33.47	{K};
	a[6][38]	:==	-48.25	{K};
	a[6][39]	:==	-172.4	{K};
	a[6][40]	:==	0	{K};
	a[6][41]	:==	165.7	{K};
	a[6][42]	:==	0	{K};
	a[6][43]	:==	0	{K};
	a[6][44]	:==	0	{K};
	a[6][45]	:==	214	{K};
	a[6][46]	:==	0	{K};
	a[6][47]	:==	0	{K};
(* row 7 *)
	a[7][1]		:==	300.00	{K};
	a[7][2]		:==	496.100	{K};
	a[7][3]		:==	362.30	{K};
	a[7][4]		:==	377.60	{K};
	a[7][5]		:==	-229.10	{K};
	a[7][6]		:==	289.6	{K};
	a[7][7]		:==	0.0	{K};
	a[7][8]		:==	324.5	{K};
	a[7][9]		:==	-195.40	{K};
	a[7][10]	:==	-257.3	{K};
	(* a[7][11]	:==	14.42	{K};  ?!?*)
	a[7][11]	:==	72.87	{K};
	a[7][12]	:==	0	{K};
	a[7][13]	:==	540.5	{K};
	a[7][14]	:==	48.89	{K};
	a[7][15]	:==	168	{K};
	a[7][16]	:==	304	{K};
	a[7][17]	:==	-59.29	{K};
	a[7][18]	:==	459	{K};
	a[7][19]	:==	112.6	{K};
	a[7][20]	:==	-14.09	{K};
	a[7][21]	:==	325.4	{K};
	a[7][22]	:==	370.4	{K};
	a[7][23]	:==	353.7	{K};
	a[7][24]	:==	497.5	{K};
	a[7][25]	:==	678.2	{K};
	a[7][26]	:==	220.6	{K};
	a[7][27]	:==	399.5	{K};
	a[7][28]	:==	887.1	{K};
	a[7][29]	:==	0	{K};
	a[7][30]	:==	188	{K};
	a[7][31]	:==	13.53	{K};
	a[7][32]	:==	0	{K};
	a[7][33]	:==	0	{K};
	a[7][34]	:==	0	{K};
	a[7][35]	:==	-139	{K};
	a[7][36]	:==	160.8	{K};
	a[7][37]	:==	0	{K};
	a[7][38]	:==	0	{K};
	a[7][39]	:==	319	{K};
	a[7][40]	:==	0	{K};
	a[7][41]	:==	-197.5	{K};
	a[7][42]	:==	0	{K};
	a[7][43]	:==	0	{K};
	a[7][44]	:==	832.2	{K};
	a[7][45]	:==	365	{K};
	a[7][46]	:==	0	{K};
	a[7][47]	:==	0	{K};
(* row 8 *)
	a[8][1]		:==	275.8	{K};
	a[8][2]		:==	217.5	{K};
	a[8][3]		:==	25.34	{K};
	a[8][4]		:==	244.2	{K};
	a[8][5]		:==	-451.6	{K};
	a[8][6]		:==	-265.2	{K};
	a[8][7]		:==	-601.8	{K};
	a[8][8]		:==	0.0	{K};
	a[8][9]		:==	-356.1	{K};
	a[8][10]	:==	0	{K};
	a[8][11]	:==	-449.4	{K};
	a[8][12]	:==	0	{K};
	a[8][13]	:==	0	{K};
	a[8][14]	:==	0	{K};
	a[8][15]	:==	0	{K};
	a[8][16]	:==	0	{K};
	a[8][17]	:==	119.9	{K};
	a[8][18]	:==	-305.5	{K};
	a[8][19]	:==	0	{K};
	a[8][20]	:==	0	{K};
	a[8][21]	:==	0	{K};
	a[8][22]	:==	0	{K};
	a[8][23]	:==	0.0	{K};
	a[8][24]	:==	1827	{K};
	a[8][25]	:==	0	{K};
	a[8][26]	:==	0	{K};
	a[8][27]	:==	0	{K};
	a[8][28]	:==	0	{K};
	a[8][29]	:==	0	{K};
	a[8][30]	:==	0	{K};
	a[8][31]	:==	-687.1	{K};
	a[8][32]	:==	0	{K};
	a[8][33]	:==	0	{K};
	a[8][34]	:==	0	{K};
	a[8][35]	:==	0	{K};
	a[8][36]	:==	0	{K};
	a[8][37]	:==	0	{K};
	a[8][38]	:==	0	{K};
	a[8][39]	:==	0	{K};
	a[8][40]	:==	0	{K};
	a[8][41]	:==	-494.2	{K};
	(* a[8][42]	:==	undefined in book. *)
	a[8][43]	:==	0	{K};
	a[8][44]	:==	0	{K};
	a[8][45]	:==	0	{K};
	a[8][46]	:==	0	{K};
	a[8][47]	:==	0	{K};
(* row 9 *)
	a[9][1]		:==	26.76	{K};
	a[9][2]		:==	42.920	{K};
	a[9][3]		:==	140.10	{K};
	a[9][4]		:==	365.80	{K};
	a[9][5]		:==	164.5	{K};
	a[9][6]		:==	108.7	{K};
	a[9][7]		:==	472.5	{K};
	a[9][8]		:==	-133.1	{K};
	a[9][9]		:==	0.0	{K};
	a[9][10]	:==	-37.36	{K};
	a[9][11]	:==	-213.7	{K};
	a[9][12]	:==	-38.47	{K};
	a[9][13]	:==	-103.6	{K};
	a[9][14]	:==	0	{K};
	a[9][15]	:==	0	{K};
	a[9][16]	:==	-169	{K};
	a[9][17]	:==	6201	{K};
	a[9][18]	:==	165.1	{K};
	a[9][19]	:==	481.7	{K};
	a[9][20]	:==	669.4	{K};
	a[9][21]	:==	-191.7	{K};
	a[9][22]	:==	-284	{K};
	a[9][23]	:==	-354.6	{K};
	a[9][24]	:==	-39.2	{K};
	a[9][25]	:==	174.5	{K};
	a[9][26]	:==	137.5	{K};
	a[9][27]	:==	0	{K};
	a[9][28]	:==	216.1	{K};
	a[9][29]	:==	-46.28	{K};
	a[9][30]	:==	-163.7	{K};
	a[9][31]	:==	0	{K};
	a[9][32]	:==	53.59	{K};
	a[9][33]	:==	245.2	{K};
	a[9][34]	:==	-246.2	{K};
	a[9][35]	:==	-44.58	{K};
	a[9][36]	:==	0	{K};
	a[9][37]	:==	-34.57	{K};
	a[9][38]	:==	0	{K};
	a[9][39]	:==	-61.7	{K};
	a[9][40]	:==	0	{K};
	a[9][41]	:==	-18.8	{K};
	a[9][42]	:==	0	{K};
	a[9][43]	:==	0	{K};
	a[9][44]	:==	0	{K};
	a[9][45]	:==	135	{K};
	a[9][46]	:==	-1680	{K};
	a[9][47]	:==	-58	{K};
(* row 10 *)
	a[10][1]	:==	505.7	{K};
	a[10][2]	:==	56.3	{K};
	a[10][3]	:==	23.39	{K};
	a[10][4]	:==	106.6	{K};
	a[10][5]	:==	-404.8	{K};
	a[10][6]	:==	-340.2	{K};
	a[10][7]	:==	232.7	{K};
	a[10][8]	:==	0	{K};
	a[10][9]	:==	128	{K};
	a[10][10]	:==	0	{K};
	a[10][11]	:==	-110.3	{K};
	a[10][12]	:==	11.31	{K};
	a[10][13]	:==	304.1	{K};
	a[10][14]	:==	0	{K};
	a[10][15]	:==	0	{K};
	a[10][16]	:==	0	{K};
	a[10][17]	:==	0	{K};
	a[10][18]	:==	0	{K};
	a[10][19]	:==	0	{K};
	a[10][20]	:==	0	{K};
	a[10][21]	:==	751.9	{K};
	a[10][22]	:==	0	{K};
	a[10][23]	:==	-483.7	{K};
	a[10][24]	:==	0	{K};
	a[10][25]	:==	0	{K};
	a[10][26]	:==	0	{K};
	a[10][27]	:==	0	{K};
	a[10][28]	:==	0	{K};
	a[10][29]	:==	0	{K};
	a[10][30]	:==	0	{K};
	a[10][31]	:==	0	{K};
	a[10][32]	:==	0	{K};
	a[10][33]	:==	0	{K};
	a[10][34]	:==	0	{K};
	a[10][35]	:==	0	{K};
	a[10][36]	:==	0	{K};
	a[10][37]	:==	0	{K};
	a[10][38]	:==	0	{K};
	a[10][39]	:==	0	{K};
	a[10][40]	:==	0	{K};
	a[10][41]	:==	0	{K};
	a[10][42]	:==	0	{K};
	a[10][43]	:==	0	{K};
	a[10][44]	:==	0	{K};
	a[10][45]	:==	-7.18	{K};
	a[10][46]	:==	333	{K};
	a[10][47]	:==	6810	{K};
(* row 11 *)
	a[11][1]	:==	114.8	{K};
	a[11][2]	:==	132.1	{K};
	a[11][3]	:==	85.84	{K};
	a[11][4]	:==	-170.0	{K};
	a[11][5]	:==	245.4	{K};
	a[11][6]	:==	249.6	{K};
	a[11][7]	:==	10000.0	{K};
	a[11][8]	:==	-36.72	{K};
	a[11][9]	:==	372.2	{K};
	a[11][10]	:==	185.1	{K};
	a[11][11]	:==	0.0	{K};
	a[11][12]	:==	372.9	{K};
	a[11][13]	:==	-235.7	{K};
	a[11][14]	:==	0	{K};
	a[11][15]	:==	-73.5	{K};
	a[11][16]	:==	0	{K};
	a[11][17]	:==	475.5	{K};
	a[11][18]	:==	0	{K};
	a[11][19]	:==	494.6	{K};
	a[11][20]	:==	660.2	{K};
	a[11][21]	:==	-34.74	{K};
	a[11][22]	:==	108.9	{K};
	a[11][23]	:==	-209.7	{K};
	a[11][24]	:==	54.47	{K};
	a[11][25]	:==	629	{K};
	a[11][26]	:==	-81.13	{K};
	a[11][27]	:==	0	{K};
	a[11][28]	:==	183	{K};
	a[11][29]	:==	0	{K};
	a[11][30]	:==	202.3	{K};
	a[11][31]	:==	-101.7	{K};
	a[11][32]	:==	148.3	{K};
	a[11][33]	:==	18.88	{K};
	a[11][34]	:==	0	{K};
	a[11][35]	:==	52.08	{K};
	a[11][36]	:==	-28.61	{K};
	a[11][37]	:==	-83.3	{K};
	a[11][38]	:==	0	{K};
	a[11][39]	:==	0	{K};
	a[11][40]	:==	0	{K};
	a[11][41]	:==	560.2	{K};
	a[11][42]	:==	0	{K};
	a[11][43]	:==	0	{K};
	a[11][44]	:==	0	{K};
	a[11][45]	:==	-54.6	{K};
	a[11][46]	:==	0	{K};
	a[11][47]	:==	6960	{K};
(* row 12 *)
	a[12][1]	:==	90.49	{K};
	a[12][2]	:==	-62.55	{K};
	a[12][3]	:==	1967	{K};
	a[12][4]	:==	2347	{K};
	a[12][5]	:==	191.2	{K};
	a[12][6]	:==	155.7	{K};
	a[12][7]	:==	0	{K};
	a[12][8]	:==	0	{K};
	a[12][9]	:==	70.42	{K};
	a[12][10]	:==	35.35	{K};
	a[12][11]	:==	-261.1	{K};
	a[12][12]	:==	0	{K};
	a[12][13]	:==	0	{K};
	a[12][14]	:==	0	{K};
	a[12][15]	:==	0	{K};
	a[12][16]	:==	0	{K};
	a[12][17]	:==	0	{K};
	a[12][18]	:==	0	{K};
	a[12][19]	:==	0	{K};
	a[12][20]	:==	-356.3	{K};
	a[12][21]	:==	0	{K};
	a[12][22]	:==	0	{K};
	a[12][23]	:==	-287.2	{K};
	a[12][24]	:==	36.84	{K};
	a[12][25]	:==	0	{K};
	a[12][26]	:==	0	{K};
	a[12][27]	:==	0	{K};
	a[12][28]	:==	0	{K};
	a[12][29]	:==	4.339	{K};
	a[12][30]	:==	0	{K};
	a[12][31]	:==	0	{K};
	a[12][32]	:==	0	{K};
	a[12][33]	:==	0	{K};
	a[12][34]	:==	0	{K};
	a[12][35]	:==	0	{K};
	a[12][36]	:==	0	{K};
	a[12][37]	:==	0	{K};
	a[12][38]	:==	0	{K};
	a[12][39]	:==	0	{K};
	a[12][40]	:==	0	{K};
	a[12][41]	:==	-70.24	{K};
	a[12][42]	:==	0	{K};
	a[12][43]	:==	0	{K};
	a[12][44]	:==	0	{K};
	a[12][45]	:==	0	{K};
	a[12][46]	:==	0	{K};
	a[12][47]	:==	0	{K};
(* row 13 *)
	a[13][1]	:==	83.36	{K};
	a[13][2]	:==	26.51	{K};
	a[13][3]	:==	52.13	{K};
	a[13][4]	:==	65.69	{K};
	a[13][5]	:==	237.7	{K};
	a[13][6]	:==	238.4	{K};
	a[13][7]	:==	-314.7	{K};
	a[13][8]	:==	0	{K};
	a[13][9]	:==	191.1	{K};
	a[13][10]	:==	-7.838	{K};
	a[13][11]	:==	461.3	{K};
	a[13][12]	:==	0	{K};
	a[13][13]	:==	0	{K};
	a[13][14]	:==	0	{K};
	a[13][15]	:==	141.7	{K};
	a[13][16]	:==	0	{K};
	a[13][17]	:==	0	{K};
	a[13][18]	:==	0	{K};
	a[13][19]	:==	-18.51	{K};
	a[13][20]	:==	664.6	{K};
	a[13][21]	:==	301.1	{K};
	a[13][22]	:==	137.8	{K};
	a[13][23]	:==	-154.3	{K};
	a[13][24]	:==	47.67	{K};
	a[13][25]	:==	66.15	{K};
	a[13][26]	:==	95.18	{K};
	a[13][27]	:==	0	{K};
	a[13][28]	:==	140.9	{K};
	a[13][29]	:==	-8.538	{K};
	a[13][30]	:==	0	{K};
	a[13][31]	:==	-20.11	{K};
	a[13][32]	:==	-149.5	{K};
	a[13][33]	:==	-202.3	{K};
	a[13][34]	:==	0	{K};
	a[13][35]	:==	172.1	{K};
	a[13][36]	:==	0	{K};
	a[13][37]	:==	240.2	{K};
	a[13][38]	:==	-273.9	{K};
	a[13][39]	:==	254.8	{K};
	a[13][40]	:==	0	{K};
	a[13][41]	:==	417	{K};
	a[13][42]	:==	0	{K};
	a[13][43]	:==	0	{K};
	a[13][44]	:==	0	{K};
	a[13][45]	:==	5780	{K};
	a[13][46]	:==	131	{K};
	a[13][47]	:==	0	{K};
(* row 14 *)
	a[14][1]	:==	-30.48	{K};
	a[14][2]	:==	1.163	{K};
	a[14][3]	:==	-44.850	{K};
	a[14][4]	:==	0	{K};
	a[14][5]	:==	-164	{K};
	a[14][6]	:==	-481.7	{K};
	a[14][7]	:==	-330.4	{K};
	a[14][8]	:==	0	{K};
	a[14][9]	:==	0	{K};
	a[14][10]	:==	0	{K};
	a[14][11]	:==	0	{K};
	a[14][12]	:==	0	{K};
	a[14][13]	:==	0	{K};
	a[14][14]	:==	0	{K};
	a[14][15]	:==	63.72	{K};
	a[14][16]	:==	-41.11	{K};
	a[14][17]	:==	-200.7	{K};
	a[14][18]	:==	0	{K};
	a[14][19]	:==	0	{K};
	a[14][20]	:==	0	{K};
	a[14][21]	:==	0	{K};
	a[14][22]	:==	0	{K};
	a[14][23]	:==	0	{K};
	a[14][24]	:==	-99.81	{K};
	a[14][25]	:==	68.81	{K};
	a[14][26]	:==	0	{K};
	a[14][27]	:==	0	{K};
	a[14][28]	:==	0	{K};
	a[14][29]	:==	-70.14	{K};
	a[14][30]	:==	0	{K};
	a[14][31]	:==	0	{K};
	a[14][32]	:==	0	{K};
	a[14][33]	:==	0	{K};
	a[14][34]	:==	0	{K};
	a[14][35]	:==	0	{K};
	a[14][36]	:==	0	{K};
	a[14][37]	:==	0	{K};
	a[14][38]	:==	0	{K};
	a[14][39]	:==	0	{K};
	a[14][40]	:==	0	{K};
	a[14][41]	:==	0	{K};
	a[14][42]	:==	0	{K};
	a[14][43]	:==	0	{K};
	a[14][44]	:==	0	{K};
	a[14][45]	:==	0	{K};
	a[14][46]	:==	0	{K};
	a[14][47]	:==	0	{K};
(* row 15 *)
	a[15][1]	:==	65.33	{K};
	a[15][2]	:==	-28.7	{K};
	a[15][3]	:==	-22.31	{K};
	a[15][4]	:==	223	{K};
	a[15][5]	:==	-150	{K};
	a[15][6]	:==	-500	{K};
	a[15][7]	:==	-448.2	{K};
	a[15][8]	:==	0	{K};
	a[15][9]	:==	0	{K};
	a[15][10]	:==	0	{K};
	a[15][11]	:==	136	{K};
	a[15][12]	:==	0	{K};
	a[15][13]	:==	-49.3	{K};
	a[15][14]	:==	108.8	{K};
	a[15][15]	:==	0	{K};
	a[15][16]	:==	-189.2	{K};
	a[15][17]	:==	0	{K};
	a[15][18]	:==	0	{K};
	a[15][19]	:==	0	{K};
	a[15][20]	:==	0	{K};
	a[15][21]	:==	0	{K};
	a[15][22]	:==	0	{K};
	a[15][23]	:==	0	{K};
	a[15][24]	:==	71.23	{K};
	a[15][25]	:==	4350	{K};
	a[15][26]	:==	0	{K};
	a[15][27]	:==	0	{K};
	a[15][28]	:==	0	{K};
	a[15][29]	:==	0	{K};
	a[15][30]	:==	0	{K};
	a[15][31]	:==	0	{K};
	a[15][32]	:==	0	{K};
	a[15][33]	:==	0	{K};
	a[15][34]	:==	0	{K};
	a[15][35]	:==	0	{K};
	a[15][36]	:==	0	{K};
	a[15][37]	:==	0	{K};
	a[15][38]	:==	0	{K};
	a[15][39]	:==	0	{K};
	a[15][40]	:==	0	{K};
	a[15][41]	:==	-38.77	{K};
	a[15][42]	:==	0	{K};
	a[15][43]	:==	0	{K};
	a[15][44]	:==	0	{K};
	a[15][45]	:==	0	{K};
	a[15][46]	:==	0	{K};
	a[15][47]	:==	0	{K};
(* row 16 *)
	a[16][1]	:==	-83.98	{K};
	a[16][2]	:==	-25.38	{K};
	a[16][3]	:==	-223.9	{K};
	a[16][4]	:==	109.9	{K};
	a[16][5]	:==	28.6	{K};
	a[16][6]	:==	-406.8	{K};
	a[16][7]	:==	-598.8	{K};
	a[16][8]	:==	0	{K};
	a[16][9]	:==	225.3	{K};
	a[16][10]	:==	0	{K};
	a[16][11]	:==	0	{K};
	a[16][12]	:==	0	{K};
	a[16][13]	:==	0	{K};
	a[16][14]	:==	38.89	{K};
	a[16][15]	:==	865.9	{K};
	a[16][16]	:==	0	{K};
	a[16][17]	:==	0	{K};
	a[16][18]	:==	0	{K};
	a[16][19]	:==	0	{K};
	a[16][20]	:==	0	{K};
	a[16][21]	:==	0	{K};
	a[16][22]	:==	-73.85	{K};
	a[16][23]	:==	-352.9	{K};
	a[16][24]	:==	-8.238	{K};
	a[16][25]	:==	-86.36	{K};
	a[16][26]	:==	0	{K};
	a[16][27]	:==	0	{K};
	a[16][28]	:==	0	{K};
	a[16][29]	:==	0	{K};
	a[16][30]	:==	0	{K};
	a[16][31]	:==	0	{K};
	a[16][32]	:==	0	{K};
	a[16][33]	:==	0	{K};
	a[16][34]	:==	0	{K};
	a[16][35]	:==	243.1	{K};
	a[16][36]	:==	0	{K};
	a[16][37]	:==	0	{K};
	a[16][38]	:==	-196.3	{K};
	a[16][39]	:==	0	{K};
	a[16][40]	:==	0	{K};
	a[16][41]	:==	0	{K};
	a[16][42]	:==	0	{K};
	a[16][43]	:==	0	{K};
	a[16][44]	:==	0	{K};
	a[16][45]	:==	0	{K};
	a[16][46]	:==	0	{K};
	a[16][47]	:==	0	{K};
(* row 17 *)
	a[17][1]	:==	1139	{K};
	a[17][2]	:==	2000	{K};
	a[17][3]	:==	247.5	{K};
	a[17][4]	:==	762.8	{K};
	a[17][5]	:==	-17.4	{K};
	a[17][6]	:==	-118.1	{K};
	a[17][7]	:==	-367.8	{K};
	a[17][8]	:==	-253.1	{K};
	a[17][9]	:==	-450.3	{K};
	a[17][10]	:==	0	{K};
	a[17][11]	:==	-294.8	{K};
	a[17][12]	:==	0	{K};
	a[17][13]	:==	0	{K};
	a[17][14]	:==	-15.07	{K};
	a[17][15]	:==	0	{K};
	a[17][16]	:==	0	{K};
	a[17][17]	:==	0	{K};
	a[17][18]	:==	0	{K};
	a[17][19]	:==	-281.6	{K};
	a[17][20]	:==	0	{K};
	a[17][21]	:==	287	{K};
	a[17][22]	:==	0	{K};
	a[17][23]	:==	0	{K};
	a[17][24]	:==	882	{K};
	a[17][25]	:==	287.9	{K};
	a[17][26]	:==	0	{K};
	a[17][27]	:==	-139.3	{K};
	a[17][28]	:==	0	{K};
	a[17][29]	:==	0	{K};
	a[17][30]	:==	0	{K};
	a[17][31]	:==	-136.9	{K};
	a[17][32]	:==	0	{K};
	a[17][33]	:==	0	{K};
	a[17][34]	:==	0	{K};
	a[17][35]	:==	0	{K};
	a[17][36]	:==	0	{K};
	a[17][37]	:==	0	{K};
	a[17][38]	:==	0	{K};
	a[17][39]	:==	-334.4	{K};
	a[17][40]	:==	0	{K};
	a[17][41]	:==	-89.42	{K};
	a[17][42]	:==	0	{K};
	a[17][43]	:==	0	{K};
	a[17][44]	:==	0	{K};
	a[17][45]	:==	0	{K};
	a[17][46]	:==	0	{K};
	a[17][47]	:==	0	{K};
(* row 18 *)
	a[18][1]	:==	-101.6	{K};
	a[18][2]	:==	0	{K};
	a[18][3]	:==	31.87	{K};
	a[18][4]	:==	49.8	{K};
	a[18][5]	:==	-132.3	{K};
	a[18][6]	:==	378.2	{K};
	a[18][7]	:==	-332.9	{K};
	a[18][8]	:==	-341.6	{K};
	a[18][9]	:==	-51.54	{K};
	a[18][10]	:==	0	{K};
	a[18][11]	:==	0	{K};
	a[18][12]	:==	0	{K};
	a[18][13]	:==	0	{K};
	a[18][14]	:==	0	{K};
	a[18][15]	:==	0	{K};
	a[18][16]	:==	0	{K};
	a[18][17]	:==	0	{K};
	a[18][18]	:==	0	{K};
	a[18][19]	:==	-169.7	{K};
	a[18][20]	:==	-153.7	{K};
	a[18][21]	:==	0	{K};
	a[18][22]	:==	-351.6	{K};
	a[18][23]	:==	-114.7	{K};
	a[18][24]	:==	-165.1	{K};
	a[18][25]	:==	0	{K};
	a[18][26]	:==	0	{K};
	a[18][27]	:==	0	{K};
	a[18][28]	:==	0	{K};
	a[18][29]	:==	0	{K};
	a[18][30]	:==	0	{K};
	a[18][31]	:==	0	{K};
	a[18][32]	:==	0	{K};
	a[18][33]	:==	0	{K};
	a[18][34]	:==	0	{K};
	a[18][35]	:==	0	{K};
	a[18][36]	:==	0	{K};
	a[18][37]	:==	0	{K};
	a[18][38]	:==	0	{K};
	a[18][39]	:==	0	{K};
	a[18][40]	:==	0	{K};
	a[18][41]	:==	0	{K};
	a[18][42]	:==	0	{K};
	a[18][43]	:==	0	{K};
	a[18][44]	:==	0	{K};
	a[18][45]	:==	0	{K};
	a[18][46]	:==	0	{K};
	a[18][47]	:==	0	{K};
(* row 19 *)
	a[19][1]	:==	24.82	{K};
	a[19][2]	:==	-40.62	{K};
	a[19][3]	:==	-22.97	{K};
	a[19][4]	:==	-138.4	{K};
	a[19][5]	:==	-185.4	{K};
	a[19][6]	:==	157.8	{K};
	a[19][7]	:==	242.8	{K};
	a[19][8]	:==	0	{K};
	a[19][9]	:==	-287.5	{K};
	a[19][10]	:==	0	{K};
	a[19][11]	:==	-266.6	{K};
	a[19][12]	:==	0	{K};
	a[19][13]	:==	38.81	{K};
	a[19][14]	:==	0	{K};
	a[19][15]	:==	0	{K};
	a[19][16]	:==	0	{K};
	a[19][17]	:==	777.4	{K};
	a[19][18]	:==	134.3	{K};
	a[19][19]	:==	0	{K};
	a[19][20]	:==	0	{K};
	a[19][21]	:==	88.75	{K};
	a[19][22]	:==	-152.7	{K};
	a[19][23]	:==	-15.62	{K};
	a[19][24]	:==	-54.86	{K};
	a[19][25]	:==	52.31	{K};
	a[19][26]	:==	-0.515	{K};
	a[19][27]	:==	0	{K};
	a[19][28]	:==	230.9	{K};
	a[19][29]	:==	21.37	{K};
	a[19][30]	:==	0	{K};
	a[19][31]	:==	0	{K};
	a[19][32]	:==	0	{K};
	a[19][33]	:==	0	{K};
	a[19][34]	:==	-203	{K};
	a[19][35]	:==	0	{K};
	a[19][36]	:==	81.57	{K};
	a[19][37]	:==	3.509	{K};
	a[19][38]	:==	0	{K};
	a[19][39]	:==	0	{K};
	a[19][40]	:==	0	{K};
	a[19][41]	:==	120.3	{K};
	a[19][42]	:==	0	{K};
	a[19][43]	:==	0	{K};
	a[19][44]	:==	0	{K};
	a[19][45]	:==	0	{K};
	a[19][46]	:==	0	{K};
	a[19][47]	:==	0	{K};
(* row 20 *)
	a[20][1]	:==	315.3	{K};
	a[20][2]	:==	1264	{K};
	a[20][3]	:==	62.32	{K};
	a[20][4]	:==	268.2	{K};
	a[20][5]	:==	-151	{K};
	a[20][6]	:==	1020	{K};
	a[20][7]	:==	-66.17	{K};
	a[20][8]	:==	0	{K};
	a[20][9]	:==	-297.8	{K};
	a[20][10]	:==	0	{K};
	a[20][11]	:==	-256.3	{K};
	a[20][12]	:==	312.5	{K};
	a[20][13]	:==	-338.5	{K};
	a[20][14]	:==	0	{K};
	a[20][15]	:==	0	{K};
	a[20][16]	:==	0	{K};
	a[20][17]	:==	0	{K};
	a[20][18]	:==	-313.5	{K};
	a[20][19]	:==	0	{K};
	a[20][20]	:==	0	{K};
	a[20][21]	:==	44.42	{K};
	a[20][22]	:==	120.2	{K};
	a[20][23]	:==	76.75	{K};
	a[20][24]	:==	212.7	{K};
	a[20][25]	:==	0	{K};
	a[20][26]	:==	0	{K};
	a[20][27]	:==	0	{K};
	a[20][28]	:==	0	{K};
	a[20][29]	:==	0	{K};
	a[20][30]	:==	0	{K};
	a[20][31]	:==	0	{K};
	a[20][32]	:==	0	{K};
	a[20][33]	:==	-95	{K};
	a[20][34]	:==	0	{K};
	a[20][35]	:==	-561.2	{K};
	a[20][36]	:==	0	{K};
	a[20][37]	:==	-11.16	{K};
	a[20][38]	:==	0	{K};
	a[20][39]	:==	-246.5	{K};
	a[20][40]	:==	0	{K};
	a[20][41]	:==	-337	{K};
	a[20][42]	:==	169.3	{K};
	a[20][43]	:==	127.2	{K};
	a[20][44]	:==	0	{K};
	a[20][45]	:==	0	{K};
	a[20][46]	:==	0	{K};
	a[20][47]	:==	0	{K};
(* row 21 *)
	a[21][1]	:==	91.46	{K};
	a[21][2]	:==	97.51	{K};
	a[21][3]	:==	4.68	{K};
	a[21][4]	:==	122.9	{K};
	a[21][5]	:==	562.2	{K};
	a[21][6]	:==	529	{K};
	a[21][7]	:==	698.2	{K};
	a[21][8]	:==	0	{K};
	a[21][9]	:==	286.3	{K};
	a[21][10]	:==	-47.51	{K};
	a[21][11]	:==	35.38	{K};
	a[21][12]	:==	0	{K};
	a[21][13]	:==	225.4	{K};
	a[21][14]	:==	0	{K};
	a[21][15]	:==	0	{K};
	a[21][16]	:==	0	{K};
	a[21][17]	:==	429.7	{K};
	a[21][18]	:==	0	{K};
	a[21][19]	:==	-62.41	{K};
	a[21][20]	:==	326.4	{K};
	a[21][21]	:==	0	{K};
	a[21][22]	:==	108.3	{K};
	a[21][23]	:==	249.2	{K};
	a[21][24]	:==	62.42	{K};
	a[21][25]	:==	464.4	{K};
	a[21][26]	:==	32.75	{K};
	a[21][27]	:==	0	{K};
	a[21][28]	:==	450.1	{K};
	a[21][29]	:==	59.02	{K};
	a[21][30]	:==	0	{K};
	a[21][31]	:==	0	{K};
	a[21][32]	:==	0	{K};
	a[21][33]	:==	-125.9	{K};
	a[21][34]	:==	0	{K};
	a[21][35]	:==	0	{K};
	a[21][36]	:==	0	{K};
	a[21][37]	:==	-245.4	{K};
	a[21][38]	:==	0	{K};
	a[21][39]	:==	0	{K};
	a[21][40]	:==	0	{K};
	a[21][41]	:==	63.67	{K};
	a[21][42]	:==	0	{K};
	a[21][43]	:==	0	{K};
	a[21][44]	:==	0	{K};
	a[21][45]	:==	0	{K};
	a[21][46]	:==	0	{K};
	a[21][47]	:==	0	{K};
(* row 22 *)
	a[22][1]	:==	34.01	{K};
	a[22][2]	:==	18.25	{K};
	a[22][3]	:==	121.3	{K};
	a[22][4]	:==	140.8	{K};
	a[22][5]	:==	747.7	{K};
	a[22][6]	:==	669.9	{K};
	a[22][7]	:==	708.7	{K};
	a[22][8]	:==	0	{K};
	a[22][9]	:==	423.2	{K};
	a[22][10]	:==	0	{K};
	a[22][11]	:==	-132.9	{K};
	a[22][12]	:==	0	{K};
	a[22][13]	:==	-197.7	{K};
	a[22][14]	:==	0	{K};
	a[22][15]	:==	0	{K};
	a[22][16]	:==	-141.4	{K};
	a[22][17]	:==	0	{K};
	a[22][18]	:==	587.3	{K};
	a[22][19]	:==	258.6	{K};
	a[22][20]	:==	339.6	{K};
	a[22][21]	:==	-84.53	{K};
	a[22][22]	:==	0	{K};
	a[22][23]	:==	0	{K};
	a[22][24]	:==	56.33	{K};
	a[22][25]	:==	0	{K};
	a[22][26]	:==	0	{K};
	a[22][27]	:==	0	{K};
	a[22][28]	:==	0	{K};
	a[22][29]	:==	0	{K};
	a[22][30]	:==	0	{K};
	a[22][31]	:==	0	{K};
	a[22][32]	:==	177.6	{K};
	a[22][33]	:==	0	{K};
	a[22][34]	:==	0	{K};
	a[22][35]	:==	215	{K};
	a[22][36]	:==	0	{K};
	a[22][37]	:==	0	{K};
	a[22][38]	:==	0	{K};
	a[22][39]	:==	0	{K};
	a[22][40]	:==	0	{K};
	a[22][41]	:==	-96.87	{K};
	a[22][42]	:==	0	{K};
	a[22][43]	:==	0	{K};
	a[22][44]	:==	0	{K};
	a[22][45]	:==	0	{K};
	a[22][46]	:==	0	{K};
	a[22][47]	:==	0	{K};
(* row 23 *)
	a[23][1]	:==	36.70	{K};
	a[23][2]	:==	51.060	{K};
	a[23][3]	:==	288.5	{K};
	a[23][4]	:==	33.61	{K};
	a[23][5]	:==	742.1	{K};
	a[23][6]	:==	649.1	{K};
	a[23][7]	:==	826.7	{K};
	a[23][8]	:==	0.0	{K};
	a[23][9]	:==	552.1	{K};
	a[23][10]	:==	242.8	{K};
	a[23][11]	:==	176.5	{K};
	a[23][12]	:==	488.9	{K};
	a[23][13]	:==	-20.93	{K};
	a[23][14]	:==	0	{K};
	a[23][15]	:==	0	{K};
	a[23][16]	:==	-293.7	{K};
	a[23][17]	:==	0	{K};
	a[23][18]	:==	18.98	{K};
	a[23][19]	:==	74.04	{K};
	a[23][20]	:==	1346	{K};
	a[23][21]	:==	-157.1	{K};
	a[23][22]	:==	0	{K};
	a[23][23]	:==	0.0	{K};
	a[23][24]	:==	-30.1	{K};
	a[23][25]	:==	0	{K};
	a[23][26]	:==	0	{K};
	a[23][27]	:==	0	{K};
	a[23][28]	:==	116.6	{K};
	a[23][29]	:==	0	{K};
	a[23][30]	:==	-64.38	{K};
	a[23][31]	:==	0	{K};
	a[23][32]	:==	86.4	{K};
	a[23][33]	:==	0	{K};
	a[23][34]	:==	0	{K};
	a[23][35]	:==	363.7	{K};
	a[23][36]	:==	0	{K};
	a[23][37]	:==	111.2	{K};
	a[23][38]	:==	0	{K};
	a[23][39]	:==	0	{K};
	a[23][40]	:==	0	{K};
	a[23][41]	:==	255.8	{K};
	a[23][42]	:==	0	{K};
	a[23][43]	:==	0	{K};
	a[23][44]	:==	0	{K};
	a[23][45]	:==	0	{K};
	a[23][46]	:==	0	{K};
	a[23][47]	:==	0	{K};

(* fixme
(* row 24 *)
	a[24][1]	:==
	a[24][2]	:==
	a[24][3]	:==
	a[24][4]	:==
	a[24][5]	:==
	a[24][6]	:==
	a[24][7]	:==
	a[24][8]	:==
	a[24][9]	:==
	a[24][10]	:==
	a[24][11]	:==
	a[24][12]	:==
	a[24][13]	:==
	a[24][14]	:==
	a[24][15]	:==
	a[24][16]	:==
	a[24][17]	:==
	a[24][18]	:==
	a[24][19]	:==
	a[24][20]	:==
	a[24][21]	:==
	a[24][22]	:==
	a[24][23]	:==
	a[24][24]	:==
	a[24][25]	:==
	a[24][26]	:==
	a[24][27]	:==
	a[24][28]	:==
	a[24][29]	:==
	a[24][30]	:==
	a[24][31]	:==
	a[24][32]	:==
	a[24][33]	:==
	a[24][34]	:==
	a[24][35]	:==
	a[24][36]	:==
	a[24][37]	:==
	a[24][38]	:==
	a[24][39]	:==
	a[24][40]	:==
	a[24][41]	:==
	a[24][42]	:==
	a[24][43]	:==
	a[24][44]	:==
	a[24][45]	:==
	a[24][46]	:==
	a[24][47]	:==

(* row 25 *)
	a[25][1]	:==
	a[25][2]	:==
	a[25][3]	:==
	a[25][4]	:==
	a[25][5]	:==
	a[25][6]	:==
	a[25][7]	:==
	a[25][8]	:==
	a[25][9]	:==
	a[25][10]	:==
	a[25][11]	:==
	a[25][12]	:==
	a[25][13]	:==
	a[25][14]	:==
	a[25][15]	:==
	a[25][16]	:==
	a[25][17]	:==
	a[25][18]	:==
	a[25][19]	:==
	a[25][20]	:==
	a[25][21]	:==
	a[25][22]	:==
	a[25][23]	:==
	a[25][24]	:==
	a[25][25]	:==
	a[25][26]	:==
	a[25][27]	:==
	a[25][28]	:==
	a[25][29]	:==
	a[25][30]	:==
	a[25][31]	:==
	a[25][32]	:==
	a[25][33]	:==
	a[25][34]	:==
	a[25][35]	:==
	a[25][36]	:==
	a[25][37]	:==
	a[25][38]	:==
	a[25][39]	:==
	a[25][40]	:==
	a[25][41]	:==
	a[25][42]	:==
	a[25][43]	:==
	a[25][44]	:==
	a[25][45]	:==
	a[25][46]	:==
	a[25][47]	:==

(* row 26 *)
	a[26][1]	:==
	a[26][2]	:==
	a[26][3]	:==
	a[26][4]	:==
	a[26][5]	:==
	a[26][6]	:==
	a[26][7]	:==
	a[26][8]	:==
	a[26][9]	:==
	a[26][10]	:==
	a[26][11]	:==
	a[26][12]	:==
	a[26][13]	:==
	a[26][14]	:==
	a[26][15]	:==
	a[26][16]	:==
	a[26][17]	:==
	a[26][18]	:==
	a[26][19]	:==
	a[26][20]	:==
	a[26][21]	:==
	a[26][22]	:==
	a[26][23]	:==
	a[26][24]	:==
	a[26][25]	:==
	a[26][26]	:==
	a[26][27]	:==
	a[26][28]	:==
	a[26][29]	:==
	a[26][30]	:==
	a[26][31]	:==
	a[26][32]	:==
	a[26][33]	:==
	a[26][34]	:==
	a[26][35]	:==
	a[26][36]	:==
	a[26][37]	:==
	a[26][38]	:==
	a[26][39]	:==
	a[26][40]	:==
	a[26][41]	:==
	a[26][42]	:==
	a[26][43]	:==
	a[26][44]	:==
	a[26][45]	:==
	a[26][46]	:==
	a[26][47]	:==

(* row 27 *)
	a[27][1]	:==
	a[27][2]	:==
	a[27][3]	:==
	a[27][4]	:==
	a[27][5]	:==
	a[27][6]	:==
	a[27][7]	:==
	a[27][8]	:==
	a[27][9]	:==
	a[27][10]	:==
	a[27][11]	:==
	a[27][12]	:==
	a[27][13]	:==
	a[27][14]	:==
	a[27][15]	:==
	a[27][16]	:==
	a[27][17]	:==
	a[27][18]	:==
	a[27][19]	:==
	a[27][20]	:==
	a[27][21]	:==
	a[27][22]	:==
	a[27][23]	:==
	a[27][24]	:==
	a[27][25]	:==
	a[27][26]	:==
	a[27][27]	:==
	a[27][28]	:==
	a[27][29]	:==
	a[27][30]	:==
	a[27][31]	:==
	a[27][32]	:==
	a[27][33]	:==
	a[27][34]	:==
	a[27][35]	:==
	a[27][36]	:==
	a[27][37]	:==
	a[27][38]	:==
	a[27][39]	:==
	a[27][40]	:==
	a[27][41]	:==
	a[27][42]	:==
	a[27][43]	:==
	a[27][44]	:==
	a[27][45]	:==
	a[27][46]	:==
	a[27][47]	:==

(* row 28 *)
	a[28][1]	:==
	a[28][2]	:==
	a[28][3]	:==
	a[28][4]	:==
	a[28][5]	:==
	a[28][6]	:==
	a[28][7]	:==
	a[28][8]	:==
	a[28][9]	:==
	a[28][10]	:==
	a[28][11]	:==
	a[28][12]	:==
	a[28][13]	:==
	a[28][14]	:==
	a[28][15]	:==
	a[28][16]	:==
	a[28][17]	:==
	a[28][18]	:==
	a[28][19]	:==
	a[28][20]	:==
	a[28][21]	:==
	a[28][22]	:==
	a[28][23]	:==
	a[28][24]	:==
	a[28][25]	:==
	a[28][26]	:==
	a[28][27]	:==
	a[28][28]	:==
	a[28][29]	:==
	a[28][30]	:==
	a[28][31]	:==
	a[28][32]	:==
	a[28][33]	:==
	a[28][34]	:==
	a[28][35]	:==
	a[28][36]	:==
	a[28][37]	:==
	a[28][38]	:==
	a[28][39]	:==
	a[28][40]	:==
	a[28][41]	:==
	a[28][42]	:==
	a[28][43]	:==
	a[28][44]	:==
	a[28][45]	:==
	a[28][46]	:==
	a[28][47]	:==

(* row 29 *)
	a[29][1]	:==
	a[29][2]	:==
	a[29][3]	:==
	a[29][4]	:==
	a[29][5]	:==
	a[29][6]	:==
	a[29][7]	:==
	a[29][8]	:==
	a[29][9]	:==
	a[29][10]	:==
	a[29][11]	:==
	a[29][12]	:==
	a[29][13]	:==
	a[29][14]	:==
	a[29][15]	:==
	a[29][16]	:==
	a[29][17]	:==
	a[29][18]	:==
	a[29][19]	:==
	a[29][20]	:==
	a[29][21]	:==
	a[29][22]	:==
	a[29][23]	:==
	a[29][24]	:==
	a[29][25]	:==
	a[29][26]	:==
	a[29][27]	:==
	a[29][28]	:==
	a[29][29]	:==
	a[29][30]	:==
	a[29][31]	:==
	a[29][32]	:==
	a[29][33]	:==
	a[29][34]	:==
	a[29][35]	:==
	a[29][36]	:==
	a[29][37]	:==
	a[29][38]	:==
	a[29][39]	:==
	a[29][40]	:==
	a[29][41]	:==
	a[29][42]	:==
	a[29][43]	:==
	a[29][44]	:==
	a[29][45]	:==
	a[29][46]	:==
	a[29][47]	:==

(* row 30 *)
	a[30][1]	:==
	a[30][2]	:==
	a[30][3]	:==
	a[30][4]	:==
	a[30][5]	:==
	a[30][6]	:==
	a[30][7]	:==
	a[30][8]	:==
	a[30][9]	:==
	a[30][10]	:==
	a[30][11]	:==
	a[30][12]	:==
	a[30][13]	:==
	a[30][14]	:==
	a[30][15]	:==
	a[30][16]	:==
	a[30][17]	:==
	a[30][18]	:==
	a[30][19]	:==
	a[30][20]	:==
	a[30][21]	:==
	a[30][22]	:==
	a[30][23]	:==
	a[30][24]	:==
	a[30][25]	:==
	a[30][26]	:==
	a[30][27]	:==
	a[30][28]	:==
	a[30][29]	:==
	a[30][30]	:==
	a[30][31]	:==
	a[30][32]	:==
	a[30][33]	:==
	a[30][34]	:==
	a[30][35]	:==
	a[30][36]	:==
	a[30][37]	:==
	a[30][38]	:==
	a[30][39]	:==
	a[30][40]	:==
	a[30][41]	:==
	a[30][42]	:==
	a[30][43]	:==
	a[30][44]	:==
	a[30][45]	:==
	a[30][46]	:==
	a[30][47]	:==

(* row 31 *)
	a[31][1]	:==
	a[31][2]	:==
	a[31][3]	:==
	a[31][4]	:==
	a[31][5]	:==
	a[31][6]	:==
	a[31][7]	:==
	a[31][8]	:==
	a[31][9]	:==
	a[31][10]	:==
	a[31][11]	:==
	a[31][12]	:==
	a[31][13]	:==
	a[31][14]	:==
	a[31][15]	:==
	a[31][16]	:==
	a[31][17]	:==
	a[31][18]	:==
	a[31][19]	:==
	a[31][20]	:==
	a[31][21]	:==
	a[31][22]	:==
	a[31][23]	:==
	a[31][24]	:==
	a[31][25]	:==
	a[31][26]	:==
	a[31][27]	:==
	a[31][28]	:==
	a[31][29]	:==
	a[31][30]	:==
	a[31][31]	:==
	a[31][32]	:==
	a[31][33]	:==
	a[31][34]	:==
	a[31][35]	:==
	a[31][36]	:==
	a[31][37]	:==
	a[31][38]	:==
	a[31][39]	:==
	a[31][40]	:==
	a[31][41]	:==
	a[31][42]	:==
	a[31][43]	:==
	a[31][44]	:==
	a[31][45]	:==
	a[31][46]	:==
	a[31][47]	:==

(* row 32 *)
	a[32][1]	:==
	a[32][2]	:==
	a[32][3]	:==
	a[32][4]	:==
	a[32][5]	:==
	a[32][6]	:==
	a[32][7]	:==
	a[32][8]	:==
	a[32][9]	:==
	a[32][10]	:==
	a[32][11]	:==
	a[32][12]	:==
	a[32][13]	:==
	a[32][14]	:==
	a[32][15]	:==
	a[32][16]	:==
	a[32][17]	:==
	a[32][18]	:==
	a[32][19]	:==
	a[32][20]	:==
	a[32][21]	:==
	a[32][22]	:==
	a[32][23]	:==
	a[32][24]	:==
	a[32][25]	:==
	a[32][26]	:==
	a[32][27]	:==
	a[32][28]	:==
	a[32][29]	:==
	a[32][30]	:==
	a[32][31]	:==
	a[32][32]	:==
	a[32][33]	:==
	a[32][34]	:==
	a[32][35]	:==
	a[32][36]	:==
	a[32][37]	:==
	a