# -*- coding: utf-8 -*-
from __future__ import division
#Declarations
#The dictionary of parameters
#name,bname,type,family,unit,value,mode,description,group,min,max,list,enable,iscombocheckbox,isused
parameterDict = {}
try:
	if Parameter:
		pass
except NameError:
	class Parameter:
		def __init__(self, **d):
			pass
#Type:Variable
#BName:w_toc
#Family:Carbon Weight Fraction
#Unit:lbf/lbf
#Mode:In
#Description:Description
W_TOC = Variable("100093404108W500", "2WS_HRM", "w_toc", u"Carbon Weight Fraction", u"lbf/lbf")
parameterDict.update({'W_TOC' : Parameter(name='W_TOC',bname='w_toc',type='Variable',family='Carbon Weight Fraction',unit='lbf/lbf',value='100093404108W500.2WS_HRM.w_toc',mode='In',description='Description',group='',min='',max='',list='',enable='True',iscombocheckbox='False',isused='True')})
#Type:Number
#Unit:kg/m3
rhom_unit = u"kg/m3"
#Mode:In
#Description:Description
#Minimum:
#Maximum:
#List:
rhom = -9999.25
parameterDict.update({'rhom' : Parameter(name='rhom',bname='',type='Number',family='',unit='kg/m3',value='-9999.25',mode='In',description='Description',group='',min='',max='',list='',enable='True',iscombocheckbox='False',isused='True')})
#Type:Number
#Unit:kg/m3
DEN_KER_unit = u"kg/m3"
#Mode:In
#Description:Description
#Minimum:
#Maximum:
#List:
DEN_KER = -9999.25
parameterDict.update({'DEN_KER' : Parameter(name='DEN_KER',bname='',type='Number',family='',unit='kg/m3',value='-9999.25',mode='In',description='Description',group='',min='',max='',list='',enable='True',iscombocheckbox='False',isused='True')})
#Type:Variable
#BName:den
#Family:Bulk Density
#Unit:kg/m3
#Mode:In
#Description:Description
DEN = Variable("100093404108W500", "2WS_HRM", "DEN", u"Bulk Density", u"kg/m3")
parameterDict.update({'DEN' : Parameter(name='DEN',bname='den',type='Variable',family='Bulk Density',unit='kg/m3',value='100093404108W500.2WS_HRM.DEN',mode='In',description='Description',group='',min='',max='',list='',enable='True',iscombocheckbox='False',isused='True')})
#Type:Variable
#BName:phit_corr
#Family:Total Porosity
#Unit:v/v
#Mode:Out
#Description:Description
#Format:auto
PHIT_CORR = Variable("100093404108W500", "2WS_HRM", "PHIT_CORR", u"Total Porosity", u"v/v")
parameterDict.update({'PHIT_CORR' : Parameter(name='PHIT_CORR',bname='phit_corr',type='Variable',family='Total Porosity',unit='v/v',value='100093404108W500.2WS_HRM.PHIT_CORR',mode='Out',description='Description',group='',min='',max='',list='',enable='True',iscombocheckbox='False',isused='True')})
#Type:Number
#Unit:g/cm3
rhofl_unit = u"g/cm3"
#Mode:In
#Description:Description
#Minimum:
#Maximum:
#List:
rhofl = 1.05
parameterDict.update({'rhofl' : Parameter(name='rhofl',bname='',type='Number',family='',unit='g/cm3',value='1.1',mode='In',description='Description',group='',min='',max='',list='',enable='True',iscombocheckbox='False',isused='True')})
#DeclarationsEnd
### Begin Automatic Generation Loop ###
loopSize = W_TOC.referenceSize()
for loopIterator in xrange(loopSize):
	datasetIterator = loopIterator
	w_toc = W_TOC.value(loopIterator)
	den = DEN.value(loopIterator)
	phit_corr = MissingValue
	### Automatic Generation Loop End ###
	if w_toc and DEN!=MissingValue:
		phit_corr=(rhom-(den*(rhom*(w_toc/DEN_KER)-(w_toc)+1)))/(rhom-rhofl)
	phit_corr = limitValue(phit_corr, 0, 0.15)		
	print phit_corr
	### Begin Automatic Generation EndLoop ###
	PHIT_CORR.setValue(loopIterator, phit_corr)
PHIT_CORR.save(True)
### Automatic Generation EndLoop End ###

__author__ = """Kienan MARION (kmarion2)"""
__date__ = """2017-03-16"""
__version__ = """1.0"""
__group__ = """"""
__suffix__ = """"""
__prefix__ = """"""