Python mf_opx Examples

Example #1

File: katz_mf.py Project: ap4909/geochemistry_data

def katz_mf(props_name,props,cond,h2o):	#calculates melt fraction as in Katz 2003. props_name=mineralogy of mantle 
	from scipy.optimize import fsolve
	import katz_melt_frac
	def fcpx_out(P,comp):
			
			#print comp
			Mcpx=comp["cpx"]
				
			#r0=0.446
			
			r0=0.45	
			r1=0.08
			return Mcpx/(r0+r1*P)
	if cond[0]!=0.0 or cond[1]!=0.0:# if not an empty P,T,H2O estimate in file
		if props_name=="Fo=92":
				A_val=1159.6#value from kelley 2010
		elif props_name=="Fo=91":
				A_val=1122.65#average of Fo=90,Fo=92
		elif props_name=="Fo=90":
				A_val=1085.7
		
		x0=fsolve(katz_melt_frac.mf_root_finding,0,(cond[1],cond[0],h2o,A_val))[0]


		if x0>fcpx_out(cond[1],props[props_name]):
			x0=katz_melt_frac.mf_opx(cond[1],cond[0],fcpx_out(cond[1],props[props_name]),h2o,A_val)
	else:
		return 0.0
	return x0

Example #2

def katz_mf_new(
    mg_num, cpx_mode, cond, h2o
):  #calculates melt fraction as in Katz 2003. props_name=mineralogy of mantle
    cpx_mode = 0.15
    from scipy.optimize import fsolve
    import katz_melt_frac
    mg_nums = [0.893419654, 0.9010621194, 0.9184344261]
    A_vals = [1085.7, 1122.65, 1159.6]

    if cond[0] != 0.0 or cond[
            1] != 0.0:  # if not an empty P,T,H2O estimate in file
        A_val = 1085.7

        x0 = fsolve(katz_melt_frac.mf_root_finding, 0,
                    (cond[1], cond[0], h2o, A_val))[0]

        if x0 == 0:
            print "Katz_melt_frac_error", cond[1], cond[0], h2o, mg_num

        if x0 > fcpx_out(cond[1], {"cpx": cpx_mode}):

            x0 = katz_melt_frac.mf_opx(cond[1], cond[0],
                                       fcpx_out(cond[1], {"cpx": cpx_mode}),
                                       h2o, A_val)
    else:
        return 0.0
    return x0

Example #3

        def katz_mf_kelley_solidus(
            self, mg_num, cpx_mode, cond, h2o
        ):  #calculates melt fraction as in Katz 2003. props_name=mineralogy of mantle
            from scipy.optimize import fsolve
            import katz_melt_frac
            if cond[0] != 0.0 or cond[
                    1] != 0.0:  # if not an empty P,T,H2O estimate in file

                A_val = 1159.66061

                self.x0 = fsolve(katz_melt_frac.mf_root_finding, 0,
                                 (cond[1], cond[0], h2o, A_val))[0]
                if self.x0 == 0:
                    print "Katz_melt_frac_error", self.name, cond[1], cond[
                        0], h2o, props_name

                if self.x0 > self.fcpx_out(cond[1], {"cpx": cpx_mode}):

                    self.x0 = katz_melt_frac.mf_opx(
                        cond[1], cond[0],
                        self.fcpx_out(cond[1], {"cpx": cpx_mode}), h2o, A_val)

            else:
                return 0.0
            return self.x0

Example #4

File: pet_read_in_with_labels.py Project: ap4909/geochemistry_data

		def katz_mf_kelley_solidus(self,mg_num,cpx_mode,cond,h2o):	#calculates melt fraction as in Katz 2003. props_name=mineralogy of mantle 
			from scipy.optimize import fsolve
			import katz_melt_frac
			if cond[0]!=0.0 or cond[1]!=0.0:# if not an empty P,T,H2O estimate in file
				
				A_val=1159.66061	
				
				self.x0=fsolve(katz_melt_frac.mf_root_finding,0,(cond[1],cond[0],h2o,A_val))[0]
				if self.x0==0:
					print "Katz_melt_frac_error",self.name,cond[1],cond[0],h2o,props_name
					
		
				if self.x0>self.fcpx_out(cond[1],{"cpx":cpx_mode}):
					
					self.x0=katz_melt_frac.mf_opx(cond[1],cond[0],self.fcpx_out(cond[1],{"cpx":cpx_mode}),h2o,A_val)
				
			else:
				return 0.0
			return self.x0

Example #5

File: compare_kelley_katz_mf_pressure.py Project: ap4909/geochemistry_data

def kelley_mf_arc_with_cpx_out(props_name,cond,h2o):#finds melt fraction using eqn 4 from kelley 2010. cond in format [T,P], using the more "refractory" constants
			import math
			
			import katz_melt_frac
			a=-5.1404654
			b=132.899012
			
			
				
			c=1159.66061#value from kelley 2010
			x=-136.88
			y=332.01
		
			print  cond[1],math.log(cond[1])
			x0= (-60.0*h2o**(1./1.85)-cond[0]+(a*cond[1]**2+b*cond[1]+c))/(-1*(x*math.log(cond[1])+y))
			if x0>0.1:
					
				x0=katz_melt_frac.mf_opx(cond[1],cond[0],0.1,h2o,1159)
			
			return x0

Example #6

File: pet_read_in_with_labels.py Project: xinzhou1006/geochemistry_data

		def katz_mf(self,props):	#calculates melt fraction as in Katz 2003. props_name=mineralogy of mantle 
			from scipy.optimize import fsolve
			import katz_melt_frac
			self.fracs={}
			for props_name in props:
				
				for self.con,self.est in self.pts.iteritems():#looping through conditions(fo,h2o,fo2) - self.con and corresponding P,T,Estimate - self.est
					if self.est!=0.0:# if not an empty P,T,H2O estimate in file
						if props_name=="Fo=92":
							A_val=1159.6#value from kelley 2010
						else:
							A_val=1085.7
			
						self.x0=fsolve(katz_melt_frac.mf_root_finding,0,(self.est[0],self.est[1],self.est[2],A_val))[0]
		
		
						if self.x0>self.fcpx_out(self.est[0],props[props_name]):
			
							self.x0=katz_melt_frac.mf_opx(self.est[0],self.est[1],self.fcpx_out(self.est[0],props[props_name]),self.est[2],A_val)
						self.fracs[self.con]=self.x0

Example #7

File: compare_kelley_katz_mf.py Project: xinzhou1006/geochemistry_data

def kelley_mf_arc_with_cpx_out(props_name,cond,h2o):#finds melt fraction using eqn 4 from kelley 2010. cond in format [T,P], using the more "refractory" constants
			import math
			
			import katz_melt_frac
			a=-5.1404654
			b=132.899012
			
			
				
			c=1159.66061#value from kelley 2010
			x=-136.88
			y=332.01
		
			print  cond[1],math.log(cond[1])
			x0= (-60.0*h2o**(1./1.85)-cond[0]+(a*cond[1]**2+b*cond[1]+c))/(-1*(x*math.log(cond[1])+y))
			if x0>0.1:
					
				x0=katz_melt_frac.mf_opx(cond[1],cond[0],0.1,h2o,1159)
			
			return x0

Example #8

File: pet_read_in_with_labels.py Project: ap4909/geochemistry_data

		def katz_mf_new(self,mg_num,cpx_mode,cond,h2o):	#calculates melt fraction as in Katz 2003. props_name=mineralogy of mantle 
			from scipy.optimize import fsolve
			import katz_melt_frac
			mg_nums=[0.893419654,0.9010621194,0.9184344261]
			A_vals=[1085.7,1122.65,1159.6]
			
			if cond[0]!=0.0 or cond[1]!=0.0:# if not an empty P,T,H2O estimate in file
				A_val=np.interp(mg_num,mg_nums,A_vals)
				
				self.x0=fsolve(katz_melt_frac.mf_root_finding,0,(cond[1],cond[0],h2o,A_val))[0]
				
				if self.x0==0:
					print "Katz_melt_frac_error",self.name,cond[1],cond[0],h2o,mg_num
					
		
				if self.x0>self.fcpx_out(cond[1],{"cpx":cpx_mode}):
					
					self.x0=katz_melt_frac.mf_opx(cond[1],cond[0],self.fcpx_out(cond[1],{"cpx":cpx_mode}),h2o,A_val)
			else:
				return 0.0
			return self.x0

Example #9

File: pet_read_in_with_labels.py Project: xinzhou1006/geochemistry_data

		def katz_mf(self,props_name,props,cond,h2o):	#calculates melt fraction as in Katz 2003. props_name=mineralogy of mantle 
			from scipy.optimize import fsolve
			import katz_melt_frac
			if cond[0]!=0.0 or cond[1]!=0.0:# if not an empty P,T,H2O estimate in file
				if props_name=="Fo=92":
						A_val=1159.6#value from kelley 2010
				elif props_name=="Fo=91":
						A_val=1122.65#average of Fo=90,Fo=92
				elif props_name=="Fo=90":
						A_val=1085.7	
				
				self.x0=fsolve(katz_melt_frac.mf_root_finding,0,(cond[1],cond[0],h2o,A_val))[0]
				if self.x0==0:
					print "Katz_melt_frac_error",self.name,cond[1],cond[0],h2o,props_name
					
		
				if self.x0>self.fcpx_out(cond[1],props[props_name]):
					self.x0=katz_melt_frac.mf_opx(cond[1],cond[0],self.fcpx_out(cond[1],props[props_name]),h2o,A_val)
			else:
				return 0.0
			return self.x0

Example #10

def katz_mf(
    props_name, props, cond, h2o
):  #calculates melt fraction as in Katz 2003. props_name=mineralogy of mantle
    from scipy.optimize import fsolve
    import katz_melt_frac

    def fcpx_out(P, comp):

        #print comp
        Mcpx = comp["cpx"]

        #r0=0.446

        r0 = 0.45
        r1 = 0.08
        return Mcpx / (r0 + r1 * P)

    if cond[0] != 0.0 or cond[
            1] != 0.0:  # if not an empty P,T,H2O estimate in file
        if props_name == "Fo=92":
            A_val = 1159.6  #value from kelley 2010
        elif props_name == "Fo=91":
            A_val = 1122.65  #average of Fo=90,Fo=92
        elif props_name == "Fo=90":
            A_val = 1085.7

        x0 = fsolve(katz_melt_frac.mf_root_finding, 0,
                    (cond[1], cond[0], h2o, A_val))[0]

        if x0 > fcpx_out(cond[1], props[props_name]):
            x0 = katz_melt_frac.mf_opx(cond[1], cond[0],
                                       fcpx_out(cond[1], props[props_name]),
                                       h2o, A_val)
    else:
        return 0.0
    return x0