Evgenii B. Rudnyi
Blog http://blog.rudnyi.ru
Name Last modified Size
Parent Directory 25-Mar-2009 15:34 -
dhscor1.for 31-Dec-2004 22:47 15k
dkcl.for 31-Dec-2004 22:47 8k
dkcl_ml1.for 31-Dec-2004 22:48 17k
dkcl_ml1.q 16-Jun-1992 00:40 1k
dkcl_ml2.for 31-Dec-2004 22:48 16k
dkcl_ml3.for 31-Dec-2004 22:48 17k
kcl.cp 31-Dec-2004 19:25 1k
kcl.lst 31-Dec-2004 20:40 6k
kcl.old.lst 31-Dec-2004 22:18 6k
kcl_ml1.lst 31-Dec-2004 22:25 40k
kcl_ml1.old.lst 16-Jun-1992 09:33 37k
kcl_ml2.lst 31-Dec-2004 22:36 10k
kcl_ml2.old.lst 16-Jun-1992 01:45 10k
kcl_ml3.lst 31-Dec-2004 22:41 9k
kcl_ml3.old.lst 16-Jun-1992 23:23 11k
kcltotal.dat 31-Dec-2004 19:45 9k
makefile 31-Dec-2004 22:42 1k
testPrint.for 31-Dec-2004 21:09 1k
http://Evgenii.Rudnyi.Ru/
Evgenii B. Rudnyi, 2004
Use at your own risk.
Vaporization thermodynamics of KCl
==================================
In directory, there is the code that has been written to make
computations for papers [1] and [2]. I have made it in 1992 for Lahey
compiler at my PC.
I was able to fix the code to work with g77 3.4 on my Mac. It happened
not to be too difficult. The most puzzling is that when nowadays one
prints an integer with F format, the result is zero. Look at
testPrint.for, in the old time one would receive 1.000 as output.
Strange, isn't it? Another strange thing is with puzzling warnings
from IMSL but it looks like that results are okay.
Requirements: IMSL subroutines ZXSSQ and LINV3P. The latter can be
change to a LAPACK equivalent (see my eqconst) but still one needs
ZXSSQ. Well, some replacement is in TDLIB'00.
*.old.lst are old listings. The new listings are quite close (use diff
-b). There is some difference but it probably can be ignored. The most
funny is zeros on the diagonal of the correlation matrix. The reason
is in testPrint.for.
Files:
Experimental Data
-----------------
kcl.cp contains the heat capacity;
kcltotal.dat contains all the experimental data to process;
dhscor1.for contains code that is used by all other routines.
It reads the files above.
The statistical hypothesis for ML
---------------------------------
dkcl_ml1.q contains the expert conclusion. This file is read by *_ml? code.
Main code
---------
dkcl.for computes the OLS solution in [1] and [2]
dkcl_ml1.for computes the ML(II) solution in [1] and [2]
dkcl_ml2.for computes the ML(III-1) solution in [1] and [2]
dkcl_ml3.for computes the ML(III-2) solution in [1] and [2]
[1] E.B. Rudnyi, D.W. Bonnell, D.W. Hastie. Vaporization
thermodynamics of KCl. Combined processing of total vapor pressure and
evaporation rates. Vestn. Mosk. Univ., ser. Khimiya, 35, 291-308
(1994) (in Russian). Moscow Univ. Chem. Bul., 1994, N 2 (English
transl).
[2] E.B. Rudnyi, D.W. Bonnell, D.W. Hastie. Vaporization
thermodynamics of KCl. Combining vapor pressure and gravimetrics data.
Unpublished manuscript, 1998.