@@ -72,7 +72,7 @@ function initializeFireParcel(
7272) : Parcel {
7373 // Start with parcel props from ambient air
7474 const z = background . z [ 0 ] ;
75- let thetal = background . theta [ 0 ] ; // CLASS assumes RH_near-surf < 100%
75+ let thetal = background . theta [ 0 ] ; // CLASS assumes RH_near-surf < 100%
7676 const thetavAmbient = background . thetav [ 0 ] ;
7777 let qt = background . qt [ 0 ] ;
7878 const rho = background . rho [ 0 ] ;
@@ -92,7 +92,7 @@ function initializeFireParcel(
9292 ( 2 * rho * cp * thetavAmbient * ( 1 + plumeConfig . bW ) ) ;
9393 const w = Math . cbrt ( fac_w * fire . h0 ) ;
9494
95- // Add excess temperature/humidity
95+ // Add excess temperature/humidity
9696 const dtheta = FFire / ( rho * cp * w ) ;
9797 const dqv = FqFire / ( rho * w ) ;
9898 thetal += dtheta ;
@@ -102,7 +102,7 @@ function initializeFireParcel(
102102 const T = saturationAdjustment ( thetal , qt , p , exner ) ;
103103 const qsat = qsatLiq ( p , T ) ;
104104 const ql = Math . max ( qt - qsat , 0 ) ;
105- const theta = thetal + Lv / cp / exner * ql
105+ const theta = thetal + ( Lv / cp / exner ) * ql ;
106106 const thetav = virtualTemperature ( theta , qt , ql ) ;
107107 const rh = ( ( qt - ql ) / qsat ) * 100 ;
108108 const Td = dewpoint ( qt , p / 100 ) ;
@@ -167,7 +167,7 @@ export function calculatePlume(
167167 const T = saturationAdjustment ( thetal , qt , bg . p [ i ] , bg . exner [ i ] ) ;
168168 const qsat = qsatLiq ( bg . p [ i ] , T ) ;
169169 const ql = Math . max ( qt - qsat , 0 ) ;
170- const theta = thetal + Lv / cp / bg . exner [ i ] * ql
170+ const theta = thetal + ( Lv / cp / bg . exner [ i ] ) * ql ;
171171 const thetav = virtualTemperature ( theta , qt , ql ) ;
172172 const rh = ( ( qt - ql ) / qsat ) * 100 ;
173173 const Td = dewpoint ( qt , bg . p [ i ] / 100 ) ;
0 commit comments