Skip to content

Commit

Permalink
conform code style
Browse files Browse the repository at this point in the history
  • Loading branch information
Osamu Miyawaki committed Nov 7, 2023
1 parent cefd5e5 commit d5fff71
Showing 1 changed file with 11 additions and 11 deletions.
22 changes: 11 additions & 11 deletions src/metpy/calc/thermo.py
Original file line number Diff line number Diff line change
Expand Up @@ -358,19 +358,19 @@ def dt_standard(p, t, params):
def dt_pseudoadiabatic(p, t, params):
rs = saturation_mixing_ratio._nounit(p, t)
frac = ((1 + rs) * (mpconsts.nounit.Rd * t + mpconsts.nounit.Lv * rs)
/ (mpconsts.nounit.Cp_d + rs * mpconsts.nounit.Cv_d
/ (mpconsts.nounit.Cp_d + rs * mpconsts.nounit.Cv_d
+ (mpconsts.nounit.Lv**2 * rs * (mpconsts.nounit.epsilon + rs)
/ (mpconsts.nounit.Rd * t**2))))
/ (mpconsts.nounit.Rd * t**2))))
return frac / p

def dt_reversible(p, t, params):
rs = saturation_mixing_ratio._nounit(p, t)
rl = params['rt'] - rs # assuming no ice content
frac = ((1 + params['rt']) * (mpconsts.nounit.Rd * t + mpconsts.nounit.Lv * rs)
/ (mpconsts.nounit.Cp_d + rs * mpconsts.nounit.Cv_d
+ rl * mpconsts.nounit.Cp_l + (mpconsts.nounit.Lv**2 * rs
* (mpconsts.nounit.epsilon + rs)
/ (mpconsts.nounit.Rd * t**2))))
/ (mpconsts.nounit.Cp_d + rs * mpconsts.nounit.Cv_d
+ rl * mpconsts.nounit.Cp_l + (mpconsts.nounit.Lv**2 * rs
* (mpconsts.nounit.epsilon + rs)
/ (mpconsts.nounit.Rd * t**2))))
return frac / p

def dt_so13(p, t, params):
Expand All @@ -382,11 +382,11 @@ def dt_so13(p, t, params):
rs = saturation_mixing_ratio._nounit(p, t)
qs = specific_humidity_from_mixing_ratio(rs)
frac = ((mpconsts.nounit.Rd * t + mpconsts.nounit.Lv * qs
+ ep * qs * mpconsts.nounit.Lv * (1 - params['rh0'])
* mpconsts.nounit.Rd * t / mpconsts.nounit.g)
+ ep * qs * mpconsts.nounit.Lv * (1 - params['rh0'])
* mpconsts.nounit.Rd * t / mpconsts.nounit.g)
/ (mpconsts.nounit.Cp_d
+ (mpconsts.nounit.Lv**2 * qs * mpconsts.nounit.epsilon
/ (mpconsts.nounit.Rd * t**2))))
+ (mpconsts.nounit.Lv**2 * qs * mpconsts.nounit.epsilon
/ (mpconsts.nounit.Rd * t**2))))
# cap lapse rate at dry adiabat (can be steeper with large entrainment rate)
frac = np.min([frac, mpconsts.nounit.Rd * t / mpconsts.nounit.Cp_d])
return frac / p
Expand All @@ -412,7 +412,7 @@ def dt_r14(p, t, params):
frac = (mpconsts.nounit.Rd * t / (mpconsts.nounit.g)
* mpconsts.nounit.Rv * t**2 / mpconsts.nounit.Lv
* ((-a2 + np.sqrt(a2**2 - 4 * a1 * a3)) / (2 * a1)
+ mpconsts.nounit.g / (mpconsts.nounit.Rd * t)))
+ mpconsts.nounit.g / (mpconsts.nounit.Rd * t)))
return frac / p

temperature = np.atleast_1d(temperature)
Expand Down

0 comments on commit d5fff71

Please sign in to comment.