glife/locations/lact_lib.qsrc

# lact_lib

!!	---- new lactation string index array for all lactation functions ----------
!!	lactation['active'] = old "lactate" variable
!!	lactation['induced'] = old lactateinduced variable
!!	lactation['pc_aware'] = old pcs_knowsaboutlactation variable
!!	lactation['lactaterate'] = old lactaterate variable now in ml*1000
!!	lactation['breastmv'] = old lactatemv variable now in ml*1000
!!	lactation['breastmm'] = old lactatemm variable now in ml*1000
!!	lactation['prolactinlvl'] = old prolactinlvl variable
!!	lactation['milkedvolume'] = old milkedvolume variable
!!	lactation['lactatemess'] = old lactatemess variable
!!	lactation['engorgement'] = old lact_engorgement variable
!!	lactation['dailyoverdemand'] = old milkoverdemand variable
!!	lactation['inducegrowth'] = old milkinduceday variable
!!	lactation['nipgrowth'] = variable to count nipplegrowth with lactation start nipples will grow by 2 points.
!!	lactation['milkgrowday'] = old milkgrowday variable.
!!	lactation['breastdensitytype'] = variable to define breast type. 0 = not initialized, 1 = 17% gland tissue, 2 = 36% gland tissue, 3 = 55% gland tissue, 4 = 74% gland tissue, 5 = 93% gland tissue
!!	lactation['breastglandtissue'] = the amount of glandular tissue per breast
!!	lactation['alveoliexpandlvl'] = variable for alveoli expansion 0 = not expanded 5 fully expanded.
!!	lactation['breastcount'] = how many breasts does sveta have? Usually should be set to 2.
!!	lactation['maturebreast'] = Did Svetas breast fully change during pregnancy yet? 0 = no, 1 = yes
!!	lactation['breastpumped'] = old breastpumped variable.
!!	---- old variables ---------------------------------------------------------
!!	- These Variables still need to be replaced/removed or substituted. -
!!	breastcounter
!!	di_i
!!	emb_i
!!	energytomilkpump
!!	lactaterate_use
!!	max_nip_flow_rate
!!	max_sucktion_demand
!!	mbarrcool
!!	mbarrfill
!!	mbarrloca
!!	mbarrmage
!!	mbarrtype
!!	milk_demand
!!	nip_flow
!!	nip_flow_mod
!!	nipple_cream
!!	pcs_ask_mward_help
!!	pcs_know_mward
!!	pcs_knows_electric_pump
!!	pcs_massagedherbreasts
!!	pcs_usedbreastpumponherself
!!	piits_i
!!	pump_bottletype
!!	pump_stop_brake
!!	pump_timestamp
!!	pumpinggrowth
!!	pumptime
!!	pumptime_left
!!	pumptype
!!	rbc_index
!!	reg_bottlem_count
!!	reg_bottles_count
!!	soreness_mod
!!	suck_flow_rate
!!	suck_until_pain
!!	temp_breastsize
!!	temp_lactaterate
!!	temp_timestamp
!!	temp_var
!!	temp_var2
!!	ts_i
!!	umbc_i
!!	usbc_i
!!	vmb_i

!!	This function is called only once in the begining of the game. It is calculating svetas mammary gland density and alveolicount and sets the breastcount which will be 2 naturally
if $ARGS[0] = 'init_breasttissue' and lactation['breasttissueinitiated'] <= 0:
	!!	use func('lact_lib','init_breasttissue')
	lactation['breasttissueinitiated'] = 1
	!! breast density type. Density of glands.
	!! 0 = not initialized, 1 = 17% gland tissue, 2 = 36% gland tissue, 3 = 55% gland tissue, 4 = 74% gland tissue, 5 = 93% gland tissue. This is based on studies I found on the internet.
	!! The logic behind the below is, when losing fat, breasts usually volunteer first. So the more bodyfat sveta has, the more she has in her breasts, and therefore the percentage of gland tissue goes down.
	!! Breasts cannot lose glandular tissue (if not cut out, yikes), therefore this will define the maximum they can shrink if fat is removed. At the moment this is done by Genbsize. The below will replace this in the future.
	!! Genbsize is the breastsize that is given to Sveta by her parents. During puberty the breasts swell up to this size with mostly being fat tissue. Due to genetic variability and hormones, breasts will develop more or
	!! less mammary gland tissue during puberty. As the game starts after Svetas puberty, this is impossible to simulate. That is why I came up with the below part. It will calculate the breastdensity using Svetas body fat.
	!! During pregnancy most of the breast fat is getting replaced by conjuctive tissue and mammary glands. But the genetical breastsize also affects the total mammary gland growth. So for example, if Sveta is a chubby,
	!! and has D cups with a breastdensitytype of 1, and she looses weight, her breast will shrink down to B cups. If her Genbsize is D, then she will grow D cups or bigger during pregnancy. But that is part of another function.
	lactation['breastcount'] = 2
	if lactation['breastdensitytype'] <= 0:
		if pregtimes > 0:
			lactation['breastdensitytype'] = 5
		elseif salo <= 20:
			lactation['breastdensitytype'] = 5
		elseif salo <= 40:
			lactation['breastdensitytype'] = 4
		elseif salo <= 60:
			lactation['breastdensitytype'] = 3
		elseif salo <= 80:
			lactation['breastdensitytype'] = 2
		elseif salo <= 100:
			lactation['breastdensitytype'] = 1
		end
	end
	!!	breast tissue that can have mammary glands is nbsize + magicf2b
	!!	Silicone is obviously not a good material to grow milkglands in
	!!	full breastccm would be nbsize + magicf2b + silicone.
	!!	Getting percentage tissue from breastccm.
	lactation['breastccm'] = func('lact_lib','bsizetoccm',(nbsize + magicf2b + silicone))
	lactation['nbsizepercent'] = (nbsize/(nbsize + magicf2b + silicone))*100
	lactation['magicf2bpercent'] = (magicf2b/(nbsize + magicf2b + silicone))*100
	lactation['siliconepercent'] = (silicone/(nbsize + magicf2b + silicone))*100
	lactation['nbsizepercent_wos'] = (nbsize/(nbsize + magicf2b))*100
	lactation['magicf2bpercent_wos'] = (magicf2b/(nbsize + magicf2b))*100
	lactation['nbsizechange'] = nbsize
	lactation['magicf2bchange'] = magicf2b
	lactation['siliconechange'] = silicone
	!!	Getting the usuable amount of breast tissue which can actually have mammary glands (in ccm)
	lactation['breastuseabletissue'] = ((lactation['nbsizepercent'] + lactation['magicf2bpercent'])*lactation['breastccm'])/100
	lactation['breastsiliconeccm'] = (lactation['siliconepercent']*lactation['breastccm'])/100
	!!	Getting mammary gland tissue per breast, based on the breastdensitytype. Explaination on the percentages is in the head of the function.
	if lactation['breastdensitytype'] = 1:
		lactation['breastglandtissue'] = (lactation['breastuseabletissue']*17)/100
		lactation['breastdensitypercent'] = 17
	elseif lactation['breastdensitytype'] = 2:
		lactation['breastglandtissue'] = (lactation['breastuseabletissue']*36)/100
		lactation['breastdensitypercent'] = 36
	elseif lactation['breastdensitytype'] = 3:
		lactation['breastglandtissue'] = (lactation['breastuseabletissue']*55)/100
		lactation['breastdensitypercent'] = 55
	elseif lactation['breastdensitytype'] = 4:
		lactation['breastglandtissue'] = (lactation['breastuseabletissue']*74)/100
		lactation['breastdensitypercent'] = 74
	elseif lactation['breastdensitytype'] = 5:
		lactation['breastglandtissue'] = (lactation['breastuseabletissue']*93)/100
		lactation['breastdensitypercent'] = 93
	end
	!!	calculating breast fat per breast
	lactation['breastfat'] = lactation['breastuseabletissue'] - lactation['breastglandtissue']
	!!	calculating alveolicount per breast Alveoli volume 2.18mm² (flat), milk volume 0.18mm³ (flat) and 4mm³ (inflated), milk volume 2mm³ (inflated) 1000 mm³ = 1 cm³
	if lactation['active'] <= 0:
		lactation['alveolicount'] = ((lactation['breastglandtissue']*100000)/218)
		lactation['alveoliexpandlvl'] = 0
		lactation['nipgrowth'] = 0
	else
		lactation['alveolicount'] = ((lactation['breastglandtissue']*1000)/4)
		lactation['alveoliexpandlvl'] = 10
		lactation['nipgrowth'] = 2
	end
	!!	storing start volumes for later use
	lactation['breastccm_start'] = lactation['breastccm']
	lactation['breastfat_start'] = lactation['breastfat']
	lactation['breastglandtissue_start'] = lactation['breastglandtissue']
	lactation['alveolicount_start'] = lactation['alveolicount']
	lactation['nipgrowth_start'] = lactation['nipgrowth']
	lactation['breastdensitytype_start'] = lactation['breastdensitytype']
	lactation['usablebtissue_start'] = lactation['usablebtissue']
	lactation['usablebtissue_change'] = lactation['usablebtissue']
end

if $ARGS[0] = 'bsizetoccm':
	!!	function to calculate volume per one breast from npc bust size.
	!!	Getting breast ccm per breast function based on the bra size chart in here https://en.wikipedia.org/wiki/Bra_size
	!!	I tried to create a math polynom for this, but ended up with x^16 monsters and inaccurate numbers, so I am going this approach.
	!!	This is the most accurate function I was able to make, and is doing the job really good. So please do not touch it.
	!!	usage func('lact_lib','bsizetoccm', (nbsize + magicf2b + silicone))
	if ARGS[1] < 0:
		result = 100
	elseif ARGS[1] <= 15:
		result = ((ARGS[1]*14)+100)
	elseif ARGS[1] <= 20:
		result = (((ARGS[1]-15)*16)+310)
	elseif ARGS[1] <= 25:
		result = (((ARGS[1]-20)*18)+390)
	elseif ARGS[1] <= 30:
		result = (((ARGS[1]-25)*22)+480)
	elseif ARGS[1] <= 35:
		result = (((ARGS[1]-30)*24)+590)
	elseif ARGS[1] <= 40:
		result = (((ARGS[1]-35)*28)+710)
	elseif ARGS[1] <= 45:
		result = (((ARGS[1]-40)*30)+850)
	elseif ARGS[1] <= 50:
		result = (((ARGS[1]-45)*36)+1000)
	elseif ARGS[1] <= 55:
		result = (((ARGS[1]-50)*38)+1180)
	elseif ARGS[1] <= 60:
		result = (((ARGS[1]-55)*42)+1370)
	elseif ARGS[1] <= 65:
		result = (((ARGS[1]-60)*46)+1580)
	elseif ARGS[1] <= 70:
		result = (((ARGS[1]-65)*50)+1810)
	elseif ARGS[1] <= 75:
		result = (((ARGS[1]-70)*56)+2060)
	elseif ARGS[1] <= 80:
		result = (((ARGS[1]-75)*60)+2340)
	elseif ARGS[1] <= 85:
		result = (((ARGS[1]-80)*72)+2640)
	elseif ARGS[1] > 85:
		result = (((ARGS[1]-85)*84)+3000)
	else
		result = 100
	end
end

if $ARGS[0] = 'bccmtosize':
	!!	backwardsfunction of bsizetoccm. It just reverts breast ccm to "nbsize + magicf2b + silicone" 
	!!	func('lact_lib','bccmtosize',lactation['breastccm'])
	if ARGS[1] < 0:
		result = 0
	elseif ARGS[1] <= 310:
		result = ((ARGS[1] - 100)/14)
	elseif ARGS[1] <= 390:
		result = (((ARGS[1] - 310)/16)+15)
	elseif ARGS[1] <= 480:
		result = (((ARGS[1] - 390)/18)+20)
	elseif ARGS[1] <= 590:
		result = (((ARGS[1] - 480)/22)+25)
	elseif ARGS[1] <= 710:
		result = (((ARGS[1] - 590)/24)+30)
	elseif ARGS[1] <= 850:
		result = (((ARGS[1] - 710)/28)+35)
	elseif ARGS[1] <= 1000:
		result = (((ARGS[1] - 850)/30)+40)
	elseif ARGS[1] <= 1180:
		result = (((ARGS[1] - 1000)/36)+45)
	elseif ARGS[1] <= 1370:
		result = (((ARGS[1] - 1180)/38)+50)
	elseif ARGS[1] <= 1580:
		result = (((ARGS[1] - 1370)/42)+55)
	elseif ARGS[1] <= 1810:
		result = (((ARGS[1] - 1580)/46)+60)
	elseif ARGS[1] <= 2060:
		result = (((ARGS[1] - 1810)/50)+65)
	elseif ARGS[1] <= 2340:
		result = (((ARGS[1] - 2060)/56)+70)
	elseif ARGS[1] <= 2640:
		result = (((ARGS[1] - 2340)/60)+75)
	elseif ARGS[1] <= 3000:
		result = (((ARGS[1] - 2640)/72)+80)
	elseif ARGS[1] > 3000:
		result = (((ARGS[1] - 3000)/84)+85)
	else
		result = 0
	end
end

!!	function to recalculate breastccm size. Basically an update function to the init_breasttissue, but without the initialization. Used in breastgrowth function and can be used for debugging or ... whatever comes up in the future.
if $ARGS[0] = 'breastrecalc':
	!!	use func('lact_lib','breastrecalc')
	lactation['breastccm'] = func('lact_lib','bsizetoccm',(nbsize + magicf2b + silicone))
	lactation['nbsizepercent'] = (nbsize/(nbsize + magicf2b + silicone))*100
	lactation['magicf2bpercent'] = (magicf2b/(nbsize + magicf2b + silicone))*100
	lactation['siliconepercent'] = (silicone/(nbsize + magicf2b + silicone))*100
	lactation['magicf2bpercent_wos'] = (magicf2b/(nbsize + magicf2b))*100
	lactation['breastuseabletissue'] = ((lactation['nbsizepercent'] + lactation['magicf2bpercent'])*lactation['breastccm'])/100
	lactation['breastsiliconeccm'] = (lactation['siliconepercent']*lactation['breastccm'])/100
	if lactation['alveolicount'] > 0:
		lactation['breastglandtissue']= (lactation['alveolicount'] * (218 + (lactation['alveoliexpandlvl'] * 20)))/100000
	else
		if lactation['breastdensitytype'] = 1:
			lactation['breastglandtissue'] = (lactation['breastuseabletissue']*17)/100
			lactation['breastdensitypercent'] = 17
		elseif lactation['breastdensitytype'] = 2:
			lactation['breastglandtissue'] = (lactation['breastuseabletissue']*36)/100
			lactation['breastdensitypercent'] = 36
		elseif lactation['breastdensitytype'] = 3:
			lactation['breastglandtissue'] = (lactation['breastuseabletissue']*55)/100
			lactation['breastdensitypercent'] = 55
		elseif lactation['breastdensitytype'] = 4:
			lactation['breastglandtissue'] = (lactation['breastuseabletissue']*74)/100
			lactation['breastdensitypercent'] = 74
		elseif lactation['breastdensitytype'] = 5:
			lactation['breastglandtissue'] = (lactation['breastuseabletissue']*93)/100
			lactation['breastdensitypercent'] = 93
		end
		if lactation['active'] <= 0:
			lactation['alveolicount'] = ((lactation['breastglandtissue']*100000)/218)
			lactation['alveoliexpandlvl'] = 0
			lactation['breastmm'] = 0
			lactation['breastmv'] = 0
			lactation['lactaterate'] = 0			
			if pcs_nips > lactation['nipgrowth']: pcs_nips -= lactation['nipgrowth']
			lactation['nipgrowth'] = 0
		else
			lactation['alveolicount'] = ((lactation['breastglandtissue']*100000)/418)
			lactation['alveoliexpandlvl'] = 10
			lactation['nipgrowth'] = 2
			func('lact_lib','set_breastmm')
			if lactation['prolactinlvl'] < 200: lactation['prolactinlvl'] = 200
			lactation['lactaterate'] = 1000*(lactation['prolactinlvl']/200)
			lactation['breastmv'] = 0
			lactation['nipgrowth'] = 2
			pcs_nips += lactation['nipgrowth']
		end
	end
	lactation['breastfat'] = lactation['breastuseabletissue'] - lactation['breastglandtissue']
end

!!	This function is used for breastgrowth. When nbsize, magicf2b, silicone, alveolicount or the alveoliexpandlvl are changed the breastccm changes. Everything will be written back to the specific variable
if $ARGS[0] = 'BreastGrowth':
	!!	usage func('lact_lib','BreastGrowth')
	if lactation['breastccm'] <> func('lact_lib','bsizetoccm',(nbsize + magicf2b + silicone)) or lactation['alveolicount_change'] <> lactation['alveolicount'] or lactation['alveoliexpandlvl_change'] <> lactation['alveoliexpandlvl']:
		!!	There are general breast growth events that will increase nbsize, magicf2b or silicone. So the breast "contents" have to be recalculated. This is happening here with the breastrecalc function.
		!!	Now the fun part. Warning, crazy math ahead. I got headaches while doing this, so better not mess with it or you will break everything.
		!!	first - check for alveoligrowth: yes, do the crazy stuff
		if lactation['alveolicount_change'] <> lactation['alveolicount'] or lactation['alveoliexpandlvl_change'] <> lactation['alveoliexpandlvl']:
			!!	calculating the ccm all mammarygland tissue will use up in the breast, based on alveolicount and alveoliexpandlvl. With every level, cmm will rise by 20 per alveoli.
			lactation['breastglandtissue']= (lactation['alveolicount'] * (218 + (lactation['alveoliexpandlvl'] * 20)))/100000
			!!	Crazy Math going on here. When I was coding this, only god and I knew what I was doing. Now only god knows.
			!!	Well, this one recalculates all the CCM stuff to normal nbsize, magicf2b and silicone. It works, no idea why, but it works.
			lactation['breastccm'] = (lactation['breastglandtissue'] + lactation['breastfat'] + lactation['breastsiliconeccm'])
			lactation['siliconepercent'] = (lactation['breastsiliconeccm']/lactation['breastccm'])*100
			lactation['magicf2bccm'] = ((lactation['breastccm'] - ((lactation['breastccm']*lactation['siliconepercent'])/100))*lactation['magicf2bpercent_wos'])/100
			lactation['nbsizebccm'] = lactation['breastccm'] - lactation['breastsiliconeccm'] - lactation['magicf2bccm']
			lactation['nbsizepercent'] = (lactation['nbsizebccm']/lactation['breastccm'])*100
			lactation['magicf2bpercent'] = (lactation['magicf2bccm']/lactation['nbsizebccm'])*100
			nbsize = (func('lact_lib','bccmtosize',lactation['breastccm'])*lactation['nbsizepercent'])/100
			magicf2b = (func('lact_lib','bccmtosize',lactation['breastccm'])*lactation['magicf2bpercent'])/100
			lactation['nbsizechange'] = nbsize
			lactation['magicf2bchange'] = magicf2b
			lactation['siliconechange'] = silicone
		elseif (lactation['nbsizechange'] + lactation['siliconechange'] + lactation['magicf2bchange']) <> (nbsize + magicf2b + silicone):
		!!	 and (lactation['alveolicount_change'] = lactation['alveolicount'] or lactation['alveoliexpandlvl_change'] = lactation['alveoliexpandlvl'])
			func('lact_lib','breastrecalc')
		end
	end
end

!!	new lactation['breastmm'] function. Warning! lactation['breastmm'] is now in microliter (ml*1000) for more acuracy. Alveoli count per breast times 2mm³ (2mm³ = 0,002ml) times breastcount.
!!	Breastcount can be changed later on if someone comes up with multibreasts or breast amputation.
if $ARGS[0] = 'set_breastmm':
	lactation['breastmm'] = (lactation['alveolicount']*2)*lactation['breastcount']
end

!!	Just a function used for debugging. Can be deleted if not needed anymore.
if $ARGS[0] = 'show_breast_stat':
	!!	func('lact_lib','show_breast_stat')
	'lactation[''breastcount''] = <<lactation[''breastcount'']>>'
	'lactation[''breasttissueinitiated''] = <<lactation[''breasttissueinitiated'']>>'
	'lactation[''breastccm''] = <<lactation[''breastccm'']>>'
	'lactation[''maturebreast''] = <<lactation[''maturebreast'']>>'
	'lactation[''alveolicount''] = <<lactation[''alveolicount'']>>'
	'lactation[''prolactinlvl''] = <<lactation[''prolactinlvl'']>>'
	'lactation[''dailyoverdemand''] = <<lactation[''dailyoverdemand'']>>'
	'lactation[''alveoliexpandlvl''] = <<lactation[''alveoliexpandlvl'']>>'
	'lactation[''active''] = <<lactation[''active'']>>'
	'lactation[''lactaterate''] = <<lactation[''lactaterate'']>>'
	'lactation[''breastmm''] = <<lactation[''breastmm'']>>'
	'lactation[''breastmv''] = <<lactation[''breastmv'']>>'
end

!!	milk production function
if $ARGS[0] = 'prod_milk':
	!!	function for Svetas milk production. Milk production will cost stamina, hydration, and hunger. this function is used in breastcycle
	!!	func('lact_lib','prod_milk',milkamount)
	!!	result = produced breast milk amount in ml*1000
	!!	ARGS[1] = Requested milk production amount produced in ml*1000. Usually should be the lactaterate.
	!!	If Sveta has salo <= 15, she will not produce breast milk as she is basically starving and her body will shut down everything that is causing too much energy consumption.
	if lactation['active'] > 0 and salo > 15:
		!!	checking how much hydration, stamina and hunger is needed to produce x amount of breast milk
		!!	checking if enough stamina and hydra is available
		if (ARGS[1]/10000) < pcs_stam and (ARGS[1]/10000) < pcs_hydra:
			!!	full requested milk amount can be produced. substracting from stamina, hydration energy.
			if (ARGS[1]/10000) < 1:
				pcs_stam -= 1
				pcs_hydra -= 1
			else
				pcs_stam -= (ARGS[1]/10000)
				pcs_hydra -= (ARGS[1]/10000)
			end
			if pcs_stam < 1: pcs_stam = 1
			if pcs_hydra < 1: pcs_hydra = 1
			if (ARGS[1]/9000) < pcs_energy:
				if (ARGS[1]/9000) < 1:
					pcs_energy -= 1
				else
					pcs_energy -= (ARGS[1]/9000)
					if pcs_energy < 1: pcs_energy = 1
				end
				result = ARGS[1]
			else
				!!	If energy is too low, milk will be produced but there is a chance that salo will shrink.
				pcs_energy = 1
				if rand(0,100) > 95: salo -= 1
				if pain['breasts'] < 20: pain['breasts'] += 1
				result = ARGS[1]
			end
		elseif pcs_stam <= 0 or pcs_hydra <= 0:
			!!	if sveta has not enough hydration or stamina, no milk is produced.
			result = 0
		else
			if (pcs_stam*10000) < (pcs_hydra*10000):
				temp_var = (pcs_stam*10000)
				pcs_stam = 1
				if (ARGS[1]/10000) < pcs_hydra:
					pcs_hydra -= (ARGS[1]/10000)
				else
					pcs_hydra = 1
				end
			else
				temp_var = (pcs_hydra*10000)
				pcs_hydra = 1
				if (ARGS[1]/10000) < pcs_stam:
					pcs_stam -= (ARGS[1]/10000)
				else
					pcs_stam = 1
				end
			end
			if (temp_var/9000) < pcs_energy:
				if (temp_var/9000) < 1:
					pcs_energy -= 1
				else
					pcs_energy -= (temp_var/9000)
					if pcs_energy < 1: pcs_energy = 1
				end
			else
				!!	If energy is too low, milk will be produced but there is a chance that salo will shrink.
				pcs_energy = 1
				if rand(0,100) > 95: salo -= 1
				if pain['breasts'] < 20: pain['breasts'] += 1
			end
			if pcs_energy < 1: pcs_energy = 1
			if pcs_hydra < 1: pcs_hydra = 1
			if pcs_stam < 1: pcs_stam = 1
			result = temp_var			
		end
	else
		result = 0
	end	
end

!!	Lactation On/Off switch function
if $ARGS[0] = 'lact_switch':
	!!	Usage
	!!	func('lact_lib','lact_switch')
	if lactation['active'] <= 0:	
		!!	Sveta is not aware of the lactation yet. She will notice when playing with her breasts or applying pressure to them. Or during pumping or when she starts to leak.
		lactation['pc_aware'] = 0
		!!	setting lactation to active
		lactation['active'] = 1
		!!	setting breast maximum milk volume via function
		func('lact_lib','set_breastmm')
		!!	resetting prolactinlvl so that the initial production rate is not too high
		if lactation['prolactinlvl'] < 200: lactation['prolactinlvl'] = 200
		!!	setting the lactation rate. Warning! lactation rate is in ml*1000 now for higher accuracy.
		lactation['lactaterate'] = 1000*(lactation['prolactinlvl']/200)
		!!	setting milk volume in breasts to 0. Will fill up with lactaterate in another part of the breastcycle
		lactation['breastmv'] = 0
		!!	nipples will grow a bit with extended milk channels. If lactation stops, this will be substracted from the nipples again and nipples will shrink.
		lactation['nipgrowth'] = 2
		pcs_nips += lactation['nipgrowth']
		if lactation['alveoliexpandlvl'] < 10: lactation['alveoliexpandlvl'] = 10
	else
		!!	If lactation is on, the stuff below will switch it off again.
		lactation['active'] = 0
		lactation['lactaterate'] = 0
		lactation['breastmv'] = 0
		lactation['breastmm']
		if pcs_nips > lactation['nipgrowth']: pcs_nips -= lactation['nipgrowth']
		lactation['nipgrowth'] = 0	
	end
end

!!	new breast cycle function, will replace the monster lactation function in femcyc
if $ARGS[0] = 'breastcycle':
	if lactation['breasttissueinitiated'] <= 0: func('lact_lib','init_breasttissue')
	!!	usage func('lact_lib','breastcycle')
	!!	if lactating or not, breasts will grow during pregnancy. This will only happen during first full pregnancy
	if pregchem => 2191 and lactation['maturebreast'] <= 0:
		!!	Alveoligrowth/breastgrowth + 78876 new aveoli per breast average during pregnancy
		!!	Checking pregchem last. It should not be small than 0.
		if lactation['pregchemlast'] <= 0:
			lactation['pregchemlast'] = 2191
		end
		if lactation['preggrowth'] <= 78876 and (pregchem - lactation['pregchemlast']) > 0:
			lactation['preggrowth'] = (pregchem - 2191) * 36
			lactation['alveolicount'] = lactation['alveolicount_start'] + lactation['preggrowth']
			lactation['growthsoreness_on'] = 1
			if pain['breasts'] <= 20: pain['breasts'] = 20
		else
			lactation['growthsoreness_on'] = 0
			lactation['maturebreast'] = 1
		end
		lactation['pregchemlast'] = pregchem
	elseif PregChem > 4383 and lactation['maturebreast'] > 0:
	!!	pregnancy prolactinlvl rise
		if lactation['prolactinlvl'] < 200 and lactation['active'] <= 0:
			if rand(0,100) > 25:
				lactation['prolactinlvl'] +=1
			end
		end
	end
	!!	checking for prolactin and if svetas breasts were worked (e.g. breast stimulation via massage, pumps, by mouth, or whatever)
	if lactation['active'] <= 0:
		!!	prolactinlvl calculation. If sveta never was pregnant, prolactinlvl will grow slower
		if lactation['milkgrowday'] < daystart:
			if (lactation['dailyoverdemand']/200) <= 0 and lactation['dailyoverdemand'] > 0:
				lactation['prolactinlvl'] += 1
			elseif (lactation['dailyoverdemand']/200) >= 10 and pregtimes <= 0:
				lactation['prolactinlvl'] += 10
			elseif (lactation['dailyoverdemand']/200) >= 50 and pregtimes > 1:
				lactation['prolactinlvl'] += 50
			else
				if lactation['prolactinlvl'] > 0 and pregchem < 2191:
					if ((lactation['prolactinlvl']/2) + (lactation['prolactinlvl']/5) - (lactation['prolactinlvl']/8)) < 1:
						lactation['prolactinlvl'] -= 1
					else
						lactation['prolactinlvl'] -= ((lactation['prolactinlvl']/2) + (lactation['prolactinlvl']/5) - (lactation['prolactinlvl']/8))
					end
				end
			end
		end
		if lactation['prolactinlvl'] > 200: lactation['prolactinlvl'] = 200
		!!	expand level cannot be over 5
		if lactation['alveoliexpandlvl'] > 10: lactation['alveoliexpandlvl'] = 10
		!!	Breast growth due to mammarygland expansion. Happens if prolactin level is high enough and will cause breastsoreness.
		!!	If prolactinlvl lowers, milkglands will shrink again.
		if lactation['milkgrowday'] < daystart and lactation['prolactinlvl'] >= 100 and lactation['alveoliexpandlvl'] < 10:
			lactation['alveoliexpandlvl_change'] = lactation['alveoliexpandlvl']
			lactation['alveoliexpandlvl'] += 1
			lactation['growthsoreness_on'] = 1
			if pain['breasts'] <= 20: pain['breasts'] = 20
		elseif lactation['milkgrowday'] < daystart and lactation['prolactinlvl'] < 100 and lactation['alveoliexpandlvl'] > 0:
			lactation['alveoliexpandlvl_change'] = lactation['alveoliexpandlvl']
			lactation['alveoliexpandlvl'] -= 1
			lactation['growthsoreness_on'] = 1
			if pain['breasts'] <= 20: pain['breasts'] = 20
		end
		!!	If sveta meets alveoliexpandlvl 10 and her prolactinlvl is between random 200 and 250, her lactation will set in. This is happening during induction or during pregnancy
		if rand(1,50) + lactation['prolactinlvl']  >= rand(200,250) and lactation['alveoliexpandlvl'] >= 10:
			func('lact_lib','lact_switch')
		end
	else
	!!	things that happen when Sveta is lactating lactation['active'] > 1
		if lactation['milkgrowday'] < daystart:
			!!	Lactation rate growth. The higher the lactationrate is, the slower it will grow. I have not found a study on this, so this is pure asumption and common sense.
			if lactation['dailyoverdemand'] > 0 and lactation['breastpumped'] >= 0:
				if lactation['lactaterate'] > 350000:
					lactation['lactaterate'] += (5*lactation['prolactinlvl'])/2
				elseif lactation['lactaterate'] > 600000:
					lactation['lactaterate'] += lactation['prolactinlvl']/2
				else
					lactation['lactaterate'] += (10*lactation['prolactinlvl'])/2
				end
				prolactinlvl = 100
			elseif lactation['dailyoverdemand'] <= 0 and lactation['breastpumped'] > 0:
			!!	If there was no milk overdemand, sveta will keep her lactationrate level.
				lactation['prolactinlvl'] = 100
				lactation['dailyoverdemand'] = 0
				lactation['breastpumped'] = 0
			else
			!!	If there was no demand at all, svetas lactation rate will drop. Weaning off will take 7-10 days usually. If there is no engorgment happening, it will take longer (up to 30 days)
				lactation['dailyoverdemand'] = 0
				lactation['breastpumped'] = 0
				if lactation['lactaterate'] > 0:
					!!	lactation rate is dropping with below function.
					lactation['prolactinlvl'] = 100
					if ((lactation['lactaterate']/2) + (lactation['lactaterate']/5) - (lactation['lactaterate']/8)) < 1:
						lactation['lactaterate'] -= 1
					else
						lactation['lactaterate'] = (lactation['lactaterate'] - (lactation['lactaterate']/2) + (lactation['lactaterate']/5) - (lactation['lactaterate']/8))
					end
				elseif lactation['lactaterate'] <= 0 and lactation['maturebreast'] > 0 and pregchem >= 4383:
					!!	during pregnancy lactation is not switched off. Svetas breasts will always produce a bit off milk
					lactation['lactaterate'] = 10
				elseif lactation['lactaterate'] <= 0:
					!!	switching off lactation production and making sure prolactin is dropping. Nipples shrink too. There are some changes that are permanent, especially during pregnancy.
					func('lact_lib','lact_switch')
				end
			end
		end
		!!	Filling lactaterate into svetas breast every hour
		if lactation['breastmv'] < lactation['breastmm'] and lactation['lactaterate'] > 50:
			lactation['breastmv'] += func('lact_lib','prod_milk',lactation['lactaterate'])
			if lactation['breastmv'] > lactation['breastmm']: lactation['breastmv'] = lactation['breastmm']
		end
	end
	if lactation['milkgrowday'] < daystart:
		func('lact_lib','BreastGrowth')
		lactation['dailyoverdemand'] = 0
		lactation['alveolicount_change'] = lactation['alveolicount']
		lactation['milkgrowday'] = daystart
	end	
end

!!	this needs to be reworked
if $ARGS[0] = 'lactate_optout':
	!! func('lact_lib','lactate_optout')
	!! use gs 'lact_lib','lactate_optout'
	!! Used for switching off all lactation related content
	lactation['active'] = 0
	prolactinlvl = 0
	lactatemess = 0
	lact_engorgement = 0
	lactation['breastmv'] = 0
	lactation['breastmm'] = 0
	lactaterate = 0
	milk_demand = 0
end

!!	get suckflow rate function
if $ARGS[0] = 'get_suckflowrate':
	!!	lactation['suck_flow_rate'] = func('lact_lib','get_suckflowrate', sucktion_type)
	!!	The suck flow rate is the demand that is applied to Svetas breasts per minute. Be it breast massaging, direct nipple sucking or the simple use of a breast pump.
	!!	ARGS[1] = sucktion_type: the type of sucktion. adult_mouth [0], child_mouth [1], baby_mouth [2], manual_breastpump [3] , electric_breastpump [4], hand_expressing [5], cow milker [6], weak baby_mouth [7], stimulation during sex/sex play [8]
	!!	what type of sucktion is happening? different sucktion can cause a different demand of breastmilk.
	!!	Result is base rate of sucktion 10ml/minute x 1000. So if the sucktion rate is 10ml/min then suck_flow_rate = 100.
	if ARGS[1] = 0:
		!!	'adult sucking'
		result = 12000
	elseif ARGS[1] = 1:
		!!	'child sucking'
		result = 11000
	elseif ARGS[1] = 2:
		!!	'baby sucking'
		result = 8000
	elseif ARGS[1] = 3:
		!!	'manual_breastpump sucking'
		result = 18000
	elseif ARGS[1] = 4:
		!!	'electric_breastpump sucking'
		result = 40000
	elseif ARGS[1] = 5:
		!!	'hand expression'
		result = 10500
	elseif ARGS[1] = 6:
		!!	'cow milker sucking'
		result = 100000
	elseif ARGS[1] = 7:
		!!	'weak baby sucking'
		result = 5000
	elseif ARGS[1] = 8:
		!!	'manual stimulation during sex/sex play'
		result = 6000
	else
		!!	'suck flow else branch'
		result = 10000
	end	
end

!!	get max nipple flow rate function
if $ARGS[0] = 'get_maxnipflowrate':
	!!	lactation['max_nip_flow_rate'] = func('lact_lib','get_maxnipflowrate')
	!!	How much milk can be sucked out through svetas nipples per minute.
	!!	Result is base rate is 10ml/minute x 1000.
	!!	Calculating the nip flow mod. The more svetas breasts were worked/sucked, the more milk can go through her nipples.
	lactation['nip_flow_mod'] = breastcounter * 100
	!!	IT IS OVER 9000! 
	if lactation['nip_flow_mod'] > 9000: lactation['nip_flow_mod'] = 9000
	if pcs_nips > 0:
		result = (6000 + lactation['nip_flow_mod'])
	elseif pcs_nips > 20:
		result = (8000 + lactation['nip_flow_mod'])
	elseif pcs_nips > 40:
		result = (10000 + lactation['nip_flow_mod'])
	elseif pcs_nips > 60:
		result = (12000 + lactation['nip_flow_mod'])
	elseif pcs_nips > 80:
		result = (14000 + lactation['nip_flow_mod'])
	elseif pcs_nips > 100:
		result = (16000 + lactation['nip_flow_mod'])
	else
		result = (10000 + lactation['nip_flow_mod'])
	end
end

!!	get mood modification
if $ARGS[0] = 'get_moodmod':
	!!	during breastfeeding or pumping or whatever, svetas mood can be pushed up.
	!!	pcs_mood += func('lact_lib','get_maxnipflowrate', sucktion type, sucktion time)
	!!	ARGS[1] = sucktion_type: the type of sucktion. adult_mouth [0], child_mouth [1], baby_mouth [2], manual_breastpump [3] , electric_breastpump [4], hand_expressing [5], cow milker [6], weak baby_mouth [7], stimulation during sex/sex play [8]
	!!	Result = additional mood
	if ARGS[1] = 0:
		!!	'adult sucking'
		lactation['moodmod'] = ARGS[2]/8
	elseif ARGS[1] = 1:
		!!	'child sucking'
		lactation['moodmod'] = ARGS[2]/2
	elseif ARGS[1] = 2:
		!!	'baby sucking'
		lactation['moodmod'] = ARGS[2]
	elseif ARGS[1] = 7:
		!!	'weak baby sucking'
		lactation['moodmod'] = ARGS[2]
	else
		lactation['moodmod'] = 0
	end
	if lactation['active'] <= 0: lactation['moodmod'] = (lactation['moodmod']/2)
	if (lactation['moodmod'] + pcs_mood) > 100:
		result = (100 - pcs_mood)
	else
		result = lactation['moodmod']
	end
end

!!	get soreness_mod
if $ARGS[0] = 'get_sorenessmod':
	!!	lactation['soreness_mod'] = func('lact_lib','get_sorenessmod')
	if lactation['nipple_cream_applied'] > 0:
		lactation['nipple_cream_applied'] = 0
		if (breastcounter/10) > 20:
			result = 40
		elseif (breastcounter/10) <= 0:
			result = 20
		else
			result = ((breastcounter/10) + 20)
		end
	else
		if (breastcounter/10) > 20:
			result = 20
		elseif (breastcounter/10) <= 0:
			result = 1
		else
			result = (breastcounter/10)
		end
	end
end

!!	get breastmilktime function
if $ARGS[0] = '$get_breastmilk_time':
	!!	func('lact_lib','$get_breastmilk_time', sucktion_type, milk_demand)
	!!	ARGS[1] = sucktion_type: the type of sucktion. adult_mouth [0], child_mouth [1], baby_mouth [2], manual_breastpump [3] , electric_breastpump [4], hand_expressing [5], cow milker [6], weak baby_mouth [7], stimulation during sex/sex play [8]
	!!	ARGS[2] = milk_demand: how much milk is demanded. Warning: Should be ml and not in ml*10 for ease of use. :D
	!!	RESULT = Time in Minutes
	!!	Usage: breastmilktimetoget = func('lact_lib','$get_breastmilk_time', 2, 150)
	!!	This function is just calculating the amount of time needed to get a specific amount of milk from svetas breasts. It is not pumping milk from her or changes anything else.
	!!	If you want to milk sveta, use get_breastmilk.
	if lactation['active'] > 0:
		!! suck_flow_rate initialized Base Rate of sucking is 10ml/minute.
		lactation['suck_flow_rate'] = func('lact_lib','get_suckflowrate', ARGS[1])
		!!	max nip flow rate initialized
		lactation['max_nip_flow_rate'] = func('lact_lib','get_maxnipflowrate')
		!!	getting temp lactaterate
		if (lactaterate/60) < 1 and lactaterate > 0:
			temp_lactaterate = 1
		else
			temp_lactaterate = lactaterate/60
		end

		if lactation['breastmv'] >= (ARGS[2]*10):
			if lactation['max_nip_flow_rate'] <= lactation['suck_flow_rate']:
				temp_var = (ARGS[2]*10)/lactation['max_nip_flow_rate']
			else
				temp_var = (ARGS[2]*10)/lactation['suck_flow_rate']
			end
			RESULT = temp_var
		elseif lactation['breastmv'] > 0 and lactation['breastmv'] < (ARGS[2]*10):
			if lactation['max_nip_flow_rate'] <= lactation['suck_flow_rate']:
				temp_var = lactation['breastmv']/lactation['max_nip_flow_rate']
				temp_var2 = (ARGS[2]*10) - lactation['breastmv']
				if lactation['max_nip_flow_rate'] <= temp_lactaterate:
					temp_var += temp_var2/lactation['max_nip_flow_rate']
				else
					temp_var += temp_var2/(temp_lactaterate*10)
				end
			else
				temp_var = lactation['breastmv']/lactation['suck_flow_rate']
				temp_var2 = (ARGS[2]*10) - lactation['breastmv']
				if lactation['suck_flow_rate'] <= temp_lactaterate:
					temp_var += temp_var2/lactation['suck_flow_rate']
				else
					temp_var += temp_var2/(temp_lactaterate*10)
				end
			end
			RESULT = temp_var
		elseif lactation['breastmv'] <=0:
			if lactation['max_nip_flow_rate'] <= lactation['suck_flow_rate']:
				if lactation['max_nip_flow_rate'] <= temp_lactaterate:
					temp_var = (ARGS[2]*10)/lactation['max_nip_flow_rate']
				else
					temp_var = (ARGS[2]*10)/(temp_lactaterate*10)
				end
			else
				if lactation['suck_flow_rate'] <= temp_lactaterate:
					temp_var = (ARGS[2]*10)/lactation['suck_flow_rate']
				else
					temp_var = (ARGS[2]*10)/(temp_lactaterate*10)
				end
			end
			RESULT = temp_var
		else
			RESULT = 0
		end
	else
		RESULT = 0
	end
end

if $ARGS[0] = '$get_breastmilk':
	!!	func('lact_lib','$get_breastmilk', sucktion_type, sucktion_time)
	!!	ARGS[1] = sucktion_type: the type of sucktion. adult_mouth [0], child_mouth [1], baby_mouth [2], manual_breastpump [3] , electric_breastpump [4], hand_expressing [5], cow milker [6], weak baby_mouth [7], stimulation during sex/sex play [8]
	!!	ARGS[2] = sucktion_time: how many minutes the sucking/pumping takes place
	!!	Usage: breastmilkvolume = func('lact_lib','$get_breastmilk', 4, 20)
	!! 	RESULT = Milk in ml*10
	!!	Sveta, someone or something tries to milk Svetas breasts, so the daily variable breastpumped is switched on. This is used for the breastcycle/femcycle. If this is not switched on, possible changes like prolactinlvl etc go down.
	lactation['breastpumped'] = 1
	!! suck_flow_rate initialized Base Rate of sucking is 10ml/minute.
	lactation['suck_flow_rate'] = func('lact_lib','get_suckflowrate', ARGS[1])
	!!	max nip flow rate initialized
	lactation['max_nip_flow_rate'] = func('lact_lib','get_maxnipflowrate')
	!!	will there be a mood change?
	pcs_mood += func('lact_lib','get_maxnipflowrate',ARGS[1],ARGS[2])
	!!	get soreness that will result from pumping/sucking/massaging
	lactation['soreness_mod'] = func('lact_lib','get_sorenessmod')
	!!	If Sveta is not lactating we still need to do some actions. This is handled in here. Nipple pain + breast pain are applied as well as stamina is used.
	if lactation['active']<= 0:
		if lactation['dailyoverdemand'] <= 0: lactation['dailyoverdemand'] = 0
		pump_stop_brake = 0
		!!	This is how much milk demand is applied to the breasts. In this function it is needed to calculate the nipple pain and breast pain.
		max_sucktion_demand = lactation['suck_flow_rate']*ARGS[2]
		if pain['nipples'] <= 60:
			suck_until_pain = (((60 - pain['nipples'])*lactation['soreness_mod'])*1000)/(lactation['suck_flow_rate']/18)
		end
		temp_var = pain['nipples'] + max_sucktion_demand/(lactation['soreness_mod']*10000)
		if temp_var > 60:
			pain['nipples'] = 65
			pump_stop_brake = 1
			pumptime = suck_until_pain
			max_sucktion_demand = lactation['suck_flow_rate']*pumptime
		else
			pain['nipples'] += max_sucktion_demand/(lactation['soreness_mod']*100000)
		end
		temp_var = pain['breasts'] + max_sucktion_demand/10000
		if temp_var > 20:
			pain['breasts'] = 25
		else
			pain['breasts'] += max_sucktion_demand/10000
		end
		!!	Pumping breasts is exhausting even if no milk comes. Pumping will impact the stamina of Sveta if she is not lactating.
		temp_var = pcs_stam - (max_sucktion_demand/30000)
		if temp_var < 0:
			pcs_stam = 0
			!!	energy usage flag on.
			energytomilkpump = 1
			pcs_sweat += pumptime/30
		else
			pcs_stam -= max_sucktion_demand/30000
			pcs_sweat += pumptime/30
			!!	energy usage flag on.
			energytomilkpump = 1
		end
		if pumptime < 10:
			breastcounter += 1
		else
			breastcounter += pumptime/10
		end
		lactation['dailyoverdemand'] += max_sucktion_demand
		result = 0
	elseif lactation['active'] > 0:
		!!	How much milk demand is applied to the breasts.
		max_sucktion_demand = lactation['suck_flow_rate']*ARGS[2]
		!!	We want to know if we can just suck out what is already stored or if additional milk production is needed to go through the whole pumping/sucking.
		!!	If we only need 20ml from a 100ml storage, we would only substract 20ml, leaving 80ml inside the breasts.
		!!	Also lets make sure we start at zero from the stuff we get from Sveta.
		milkedvolume = 0
		!!	Case: Enough milk in Svetas breasts, covering the demand
		if lactation['breastmv'] >= max_sucktion_demand:
			!!	Checking if svetas nipple flow is enough to serve the demand. If yes - get the demand; If no - get what she is able to give us.
			if lactation['max_nip_flow_rate'] >= lactation['suck_flow_rate']:
				milkedvolume = max_sucktion_demand
				lactation['breastmv'] -= max_sucktion_demand
			else:
				milkedvolume = lactation['max_nip_flow_rate']*ARGS[2]
				lactation['breastmv'] -= milkedvolume
			end
		!!	Case: Not enough milk in Svetas breasts.
		elseif lactation['breastmv'] < max_sucktion_demand:
			!!	We have two sub cases in this one. Sveta has no lactation['breastmv'] and Sveta has not enough lactation['breastmv']. This means she has to produce more milk, but first we need to know how much is needed.
			!!	temp_var will be the amount of milk Sveta has to produce to serve the demand. pump time left is the time left after lactation['breastmv'] is drained.
			if lactation['breastmv'] <= 0:
				temp_var = max_sucktion_demand
				pumptime_left = ARGS[2]
			else
				temp_var = max_sucktion_demand - lactation['breastmv']
				if lactation['max_nip_flow_rate'] >= lactation['suck_flow_rate']:
					pumptime_left = ARGS[2] - (lactation['breastmv']/lactation['suck_flow_rate'])
				else
					pumptime_left = ARGS[2] - (lactation['breastmv']/lactation['max_nip_flow_rate'])
				end
				if pumptime_left < 0: pumptime_left = 0
			end
			!!	Now we know how much time we have and how much demand we have. Time to check if Sveta can do it.
			!!	For the case, that svetas milk flow is not big enough pumptime left will be 0, therefore this is handled here.
			if pumptime_left >0:
				if lactation['max_nip_flow_rate'] >= lactation['suck_flow_rate']:
					!!	the lactaterate is per hour, so it is possible that lactaterate will end up being zero when devided by 60. therefore everything here is multiplied by 100.
					if (lactation['suck_flow_rate']*100) <= ((lactaterate*100)/60):
						milkedvolume += lactation['suck_flow_rate']*pumptime_left
					else
						milkedvolume += ((((lactaterate*100)/60)*pumptime_left)/100)
						breastcounter += 2
					end
				else
					!!	Same here as above.
					if (lactation['max_nip_flow_rate']*100) <= ((lactaterate*100)/60):
						milkedvolume += lactation['max_nip_flow_rate']*pumptime_left
						breastcounter += 2
					else
						milkedvolume += ((((lactaterate*100)/60)*pumptime_left)/100)
						breastcounter += 2
					end
				end
				!!	Produced milk needs to be subtracted from the hourly lactate rate, so femcycle is only adding the milk that was not produced yet in the next hour
				if pumptime_left > 60:
					if ((pumptime_left*100)/60) > 100:
						lactaterate_use -= ((((lactaterate*100)/60)*(pumptime_left mod 60))/100)
					else
						lactaterate_use = lactaterate
					end
				else
					if lactaterate >= lactaterate_use:
						lactaterate_use -= milkedvolume
					else
						lactaterate_use = lactaterate
					end
				end
				if lactaterate_use <= 0: lactaterate_use = 0
				!!	Sveta had to make milk, so she needs to use some of her body resources to make the milk.
				if milkedvolume > 0:
					!!	energy usage flag on. Can be used in events to show that sveta is getting hungry or thirsty
					energytomilkpump = 1
					pcs_stam -= max_sucktion_demand/30000
					if pcs_stam <= 0: pcs_stam = 0
					pcs_sweat += pumptime/30
					pcs_hydra -= milkedvolume/120
					if pcs_hydra < 1: pcs_hydra = 1
					pcs_energy -= milkedvolume/1200
				end
			end
			!!	Now we check if sveta was able to make it or not. If not we have milk overdemand which will be used in femcycle to make her produce more milk next day.
			if temp_var > milkedvolume:
				lactation['dailyoverdemand'] += (temp_var - milkedvolume)
			end
			!!	time to add the pumped lactation['breastmv'] to the milkedvolume and empty lactation['breastmv'].
			milkedvolume += lactation['breastmv']
			lactation['breastmv'] = 0
			!!	Here we have nipple and breastpain. The idea is - if there is more sucktion than the nipples can take more pain is applied. If the breasts make less milk then demanded pain to the breasts is applied.
			!!	Else - there will only be nipple pain.
			if lactation['max_nip_flow_rate'] < lactation['suck_flow_rate']:
				temp_var = pain['nipples'] + ((lactation['suck_flow_rate'] - lactation['max_nip_flow_rate'])*ARGS[2])/2000 - (lactation['soreness_mod']*ARGS[2])*15
				if temp_var >= 60:
					pain['nipples'] = 65
				elseif temp_var <= 0:
					pain['nipples'] += 1
				else
					pain['nipples'] += (((lactation['suck_flow_rate'] - lactation['max_nip_flow_rate'])*ARGS[2])/2000 - (lactation['soreness_mod']*ARGS[2])*15)
				end
			else
				if ARGS[2] < 10:
					pain['nipples'] += 1
				elseif ARGS[2] >= 10 and lactation['soreness_mod'] > ARGS[2]:
					pain['nipples'] += 1
				else
					pain['nipples'] += (ARGS[2]/10 - lactation['soreness_mod']/10)
				end
			end
			if lactation['max_nip_flow_rate'] >= lactation['suck_flow_rate']:
				if (lactation['suck_flow_rate']*100) > ((lactaterate*100)/60):
					temp_var = pain['breasts'] + (((lactation['suck_flow_rate']*100)-((lactaterate*100)/60))*pumptime_left)/1000
					if temp_var > 20:
						pain['breasts'] = 20
					else
						pain['breasts'] += (((lactation['suck_flow_rate']*100)-((lactaterate*100)/60))*pumptime_left)/1000
					end
				end
			else
				if (lactation['max_nip_flow_rate']*100) > ((lactaterate*100)/60):
					temp_var = pain['breasts'] + (((lactation['max_nip_flow_rate']*100)-((lactaterate*100)/60))*pumptime_left)/1000
					if temp_var > 20:
						pain['breasts'] = 20
					else
						pain['breasts'] += (((lactation['max_nip_flow_rate']*100)-((lactaterate*100)/60))*pumptime_left)/1000
					end
				end
			end
		end
		if pumptime < 10:
			breastcounter += 1
		else
			breastcounter += pumptime/10
		end
		!!	Done - Result is the milked volume.
		RESULT = milkedvolume
	end
end

--- lact_lib ---------------------------------