Changeset 24256 for issm/trunk-jpl/src/m/mech/steadystateiceshelftemp.py

Timestamp:

10/18/19 06:52:19 (5 years ago)

Author:

bdef

Message:

BUG: fixing some negative number syntax

File:

• issm/trunk-jpl/src/m/mech/steadystateiceshelftemp.py (modified) (1 diff)

issm/trunk-jpl/src/m/mech/steadystateiceshelftemp.py

@@ -50 +50 @@
     except FloatingPointError:
         print('steadystateiceshelf warning: overflow encountered in multipy / divide / exp, trying another formulation.')
-        temperature[pos] = -(((Tb[pos] - Ts[pos]) * ki / wi[pos] + Hi[pos] * Tb[pos]) / np.exp(Hi[pos] * wi[pos] / ki) - Hi[pos] * Ts[pos] + (Tb[pos] - Ts[pos]) * ki / wi[pos]) / (Hi[pos] * (1 - np.exp(- Hi[pos] * wi[pos] / ki)))
+        temperature[pos] = -(((Tb[pos] - Ts[pos]) * ki / wi[pos] + Hi[pos] * Tb[pos]) / np.exp(Hi[pos] * wi[pos] / ki) - Hi[pos] * Ts[pos] + (Tb[pos] - Ts[pos]) * ki / wi[pos]) / (Hi[pos] * (1 - np.exp(-Hi[pos] * wi[pos] / ki)))
 
     #temperature should not be less than surface temp