From bc0e6f97a796064cd52af9bfccd1607fecd7c45c Mon Sep 17 00:00:00 2001
From: Jessica Mejia <jzmejia990@gmail.com>
Date: Mon, 1 Apr 2024 14:38:04 -0400
Subject: [PATCH] bug fix in variable density implementation

---
 crevprop/__pycache__/crevasse.cpython-39.pyc | Bin 29735 -> 29752 bytes
 crevprop/__pycache__/fracture.cpython-39.pyc | Bin 21316 -> 21436 bytes
 crevprop/crevasse.py                         |   1 +
 crevprop/fracture.py                         |   5 +++++
 4 files changed, 6 insertions(+)

diff --git a/crevprop/__pycache__/crevasse.cpython-39.pyc b/crevprop/__pycache__/crevasse.cpython-39.pyc
index f870fb47db649a984c6691ad7c9c33e7c6d33296..4e8953e74558500620341fa7ada6b694183abe78 100644
GIT binary patch
delta 119
zcmZ4ff^o+SM&3kTUM>b8(D}!e*1eJUex{5DkXOR6fU$;gA!99L2~!DE4MP^Q1W1f!
z^S{ieOhT+hGC(y&EFgjpMA%OLkYmIF;xkBaF>cn#-O9wcfAhWkW+q0K$w@`Nj9)iz
MC@NuOj^Qu^0D0mdUjP6A

delta 102
zcmdn-f^qo^M&3kTUM>b8m><BNR=$z<ex{HdkXHhPH4Isd5)36w3z#=^WIbi#WhvqV
y$`sjx2${+MbBx%SfJ_F7&6c@anHcwNexKjW#K=6kxX73B%jSbcC5+6`9A*Gwyc`h#

diff --git a/crevprop/__pycache__/fracture.cpython-39.pyc b/crevprop/__pycache__/fracture.cpython-39.pyc
index 442bba03f717079344026b81fc09b679e2d10d07..fa429a4e30c8f5875241c873719125eb3be4d354 100644
GIT binary patch
delta 443
zcmX|-OH0E*6oq>yri9T(sSO3KI21)HA_WCss3_Fe!p7=GT*TU@G(KXT#QNNcJ9qOR
z)RhlLce)Y(fg7zGSAr}5fLb#Vyo-;^Ifu*PKAyte8K~J%2nW2$eBr{x;<XA%cxMgC
zTo4m)R@ibZDs_=OyQ+<hX6-zgB9=p}oJriGm0Pn-$MJ*N^>wmixka*VTdqSk3uMD|
zP3m+L-}b@n<OPBlZ|S)O=-^Mi=Lm+_VIKvZfA)RIu)tr2k0m(a##k3Xg`JGIH{PN^
zz)?g)q(THR&=P+f-$l`TVblpDtb!^9t8EB<(=9fiS_+mt)mv{;lhNn@7{bs{MF{_(
zDy2L;aU(IF(%9?tfc7Rr{u^Sbsv`W2WLoBlnGreUc&8GGu8NB#=B5WB&MwpY;oU#o
z@^yhdHl(-8hu+%B41E?{cA>Y&4uqw7sMHrMJ8$HuSvK-!$t@y@ed?WX#GB?AXn&1G
dm<p5yj@UzHv#~>oU`?sP{YnBWN;6h){1@4{bXWiY

delta 373
zcmdn9obkvqM!rN|UM>b8kn(0vI~TE$kI{fJZ?b?v^5)kDE14L@H%~R*%E+j++1qpl
z6QkMWFXlyz-kZ}cUa~U!Zw|C;WoB&OeBDupkuhpAi?gcm1||lEVh)gTOni(Yj9iRG
z`I~K=E11N0fFyW85<+~8e9UbBRhXFm78OiB?HbLN4-^G5b==&=FT<5FaWM+~XJX>|
z$MnCbaC4X27G}o$$@X5pqPIaxz%G_y`p?J2^_PjYsC4r}FFRJ=BA`K<Ohs8B^YSOZ
z@wZ}>ova>EWmXOm&jgEuhzgJZbCEBIEe;|=K!ga05Csv{Ktgr$^#DJ{;>l`(s*F{W
zeF7z!{xVHY57c68+B_+6GNTm4)gaY{AfgmRG=T__$$@UdlY@f8g_Sr&xI{TvxGXp{
Mxi~naxj49#07Mm782|tP

diff --git a/crevprop/crevasse.py b/crevprop/crevasse.py
index 573553e..d3fb1e2 100644
--- a/crevprop/crevasse.py
+++ b/crevprop/crevasse.py
@@ -1207,6 +1207,7 @@ def paterson_empirical_relation(self, z, C=0.02):
             density as a function of depth in kg/m^3 corresponding to the
             b array input to function.
         """
+        z=np.abs(z)
         return DENSITY_ICE-(DENSITY_ICE-self.surface_density)*np.exp(-C*z)
 
 
diff --git a/crevprop/fracture.py b/crevprop/fracture.py
index 550a685..168a6d4 100644
--- a/crevprop/fracture.py
+++ b/crevprop/fracture.py
@@ -513,6 +513,11 @@ def G(b, d, H):
     """Functional expression G in KI2 evaluation.
 
     G(gamma,lambda) where gamma=b/d, lambda=d/H
+    
+    b is depth below surface, d is crevasse depth, and H is ice thickness
+    all with units of meters. 
+    
+    
     from Tada et al., 1973 p 2.27, based on polynomial curve fitting to
     numerically derived solutions for the second component of the stress
     intensity factor KI2 for mode I crack opening. Implementiation also