convert K to C for stomatal resistance (expected by Weseley 1989)

Author: SumilThakr

geoschem.go

@@ -694,7 +694,8 @@ func (gc *GEOSChem) W() NextData {
 			for j := 0; j < omega.Shape[1]; j++ {
 				for i := 0; i < omega.Shape[2]; i++ {
 					dz := -1 * math.Log(P.Get(k+1, j, i)/P.Get(k, j, i)) * rr * T.Get(k, j, i) / g // in meters
-					wVal := -omega.Get(k, j, i) * dz / (P.Get(k+1, j, i) - P.Get(k,))
+//					wVal := omega.Get(k, j, i) / (P.Get(k+1, j, i) - P.Get(k, j, i)) / dz
+					wVal := -omega.Get(k, j, i) * dz / (P.Get(k+1, j, i) - P.Get(k, j, i))
 					w.Set(wVal, k, j, i)
 				}
 			}
@@ -771,6 +772,7 @@ func (gc *GEOSChem) P() NextData {
 				for i := 0; i < PS.Shape[1]; i++ {
 					const hPa2Pa = 100.0                                      // Convert hPa to Pa.
 					p.Set(PS.Get(j, i)*bp.Get(k)+ap.Get(k)*hPa2Pa, k, j, i) // Pressure [Pa]
+//					p.Set((PS.Get(j, i)*bp.Get(k)+ap.Get(k))*hPa2Pa, k, j, i) // Pressure [Pa]
 				}
 			}
 		}

preproc.go

@@ -839,22 +839,22 @@ func stabilityMixingChemistry(LayerHeights *sparse.DenseArray, pblhFunc, ustarFu
 							To, p, ρparticle,
 							ρ, iSeasonP, seinfeldLU), 0, j, i)
 					SO2DryDep.AddVal(
-						seinfeld.DryDepGas(z, zo, u, L, To, ρ,
+						seinfeld.DryDepGas(z, zo, u, L, To-273.15, ρ,
 							G, Θsurface,
 							wesely1989.So2Data, iSeasonG,
 							weselyLU, rain, dew, true, false), 0, j, i)
 					NOxDryDep.AddVal(
-						seinfeld.DryDepGas(z, zo, u, L, To, ρ,
+						seinfeld.DryDepGas(z, zo, u, L, To-273.15, ρ,
 							G, Θsurface,
 							wesely1989.No2Data, iSeasonG,
 							weselyLU, rain, dew, false, false), 0, j, i)
 					NH3DryDep.AddVal(
-						seinfeld.DryDepGas(z, zo, u, L, To, ρ,
+						seinfeld.DryDepGas(z, zo, u, L, To-273.15, ρ,
 							G, Θsurface,
 							wesely1989.Nh3Data, iSeasonG,
 							weselyLU, rain, dew, false, false), 0, j, i)
 					VOCDryDep.AddVal(
-						seinfeld.DryDepGas(z, zo, u, L, To, ρ,
+						seinfeld.DryDepGas(z, zo, u, L, To-273.15, ρ,
 							G, Θsurface,
 							wesely1989.OraData, iSeasonG,
 							weselyLU, rain, dew, false, false), 0, j, i)