isopleths for elements
Posted: 21 Nov 2015, 19:50
Hi,
I'm new in this forum. Recently, I started working with Theriak-Domino. The problem I'm facing is how to calculate isopleths for certain elements. It is simple task for end members of solid solution. It looks as if it should be also simple for elements like Si in phengite. I found in the UserGuide that there is a special statement for his task: El(site) Site occupancies, isopleths will be calculated. However, it doesn't work for me. I'm using scripts and the one which should calculate isopleths of Si I enclosed at the end of this post. I explored much wider space than between 3.0 to 3.49 of Si apfu. However, I always obtained empty PT diagram containing only the lower limit of phengite stability. What I'm doing wrong? Thank you in advance for for help.
Jacek
I'm new in this forum. Recently, I started working with Theriak-Domino. The problem I'm facing is how to calculate isopleths for certain elements. It is simple task for end members of solid solution. It looks as if it should be also simple for elements like Si in phengite. I found in the UserGuide that there is a special statement for his task: El(site) Site occupancies, isopleths will be calculated. However, it doesn't work for me. I'm using scripts and the one which should calculate isopleths of Si I enclosed at the end of this post. I explored much wider space than between 3.0 to 3.49 of Si apfu. However, I always obtained empty PT diagram containing only the lower limit of phengite stability. What I'm doing wrong? Thank you in advance for for help.
Jacek
Code: Select all
script.script21.txt
script21.plt
tcdb55c2_gig.txt
1.0000 300 0.010000 0.1000000E-08 0.1000000E-06 0.1000000E+01 0.1000000E-03 25 25 500
0 SI(67.9159)AL(15.8024)FE(3.3276)MG(2.7531)CA(0.9584)NA(5.7252)K(3.0434)TI(0.4414)MN(0.0326)H(50)O(?) * s25 core#1 H2O and q in excess
TC 500 650
P 2500 15000
PHNG SI(T1) 1 3.0 3.49 0.01
1
0.0000000E+00 0.0000000E+00
_script21_pix/