@@ -115,8 +115,8 @@ inline void init_scavimptbl(View2DHost scavimptblvol,
115
115
// clang-format on
116
116
template <typename FUNC, typename VectorType>
117
117
KOKKOS_INLINE_FUNCTION void
118
- calculate_cloudy_volume (const int nlev, const VectorType& cld, FUNC lprec,
119
- const bool is_tot_cld, const View1D& cldv) {
118
+ calculate_cloudy_volume (const int nlev, const VectorType & cld, FUNC lprec,
119
+ const bool is_tot_cld, const View1D & cldv) {
120
120
// BAD CONSTANT
121
121
const Real small_value_30 = 1 .e -30 ;
122
122
const Real small_value_36 = 1 .e -36 ;
@@ -1136,13 +1136,13 @@ void wetdepa_v2(const Real deltat, const Real pdel, const Real cmfdqr,
1136
1136
*/
1137
1137
1138
1138
KOKKOS_INLINE_FUNCTION
1139
- void clddiag (const int nlev, haero::ConstColumnView temperature, haero::ConstColumnView pmid,
1140
- haero::ConstColumnView pdel, const View1D &cmfdqr, const View1D &evapc,
1139
+ void clddiag (const int nlev, haero::ConstColumnView temperature,
1140
+ haero::ConstColumnView pmid, haero::ConstColumnView pdel,
1141
+ const View1D &cmfdqr, const View1D &evapc,
1141
1142
const haero::ConstColumnView &cldt, const View1D &cldcu,
1142
- const View1D &cldst,
1143
- const haero::ConstColumnView &evapr, const haero::ConstColumnView & prain,
1144
- const View1D &cldv, const View1D &cldvcu,
1145
- const View1D &cldvst, const View1D &rain) {
1143
+ const View1D &cldst, const haero::ConstColumnView &evapr,
1144
+ const haero::ConstColumnView &prain, const View1D &cldv,
1145
+ const View1D &cldvcu, const View1D &cldvst, const View1D &rain) {
1146
1146
// Calculate local precipitation production rate
1147
1147
// In src/chemistry/aerosol/wetdep.F90, (prain + cmfdqr) is used for
1148
1148
// source_term
@@ -1231,27 +1231,23 @@ void cloud_diagnostics(const ThreadTeam &team,
1231
1231
// NOTE: The k loop inside clddiag cannot be converted to parallel_for
1232
1232
// because precabs requires values from the previous elevation (k-1).
1233
1233
Kokkos::single (Kokkos::PerTeam (team), [=]() {
1234
- wetdep::clddiag (nlev, temperature, pmid, pdel,
1235
- cmfdqr, evapc, cldt, cldcu,
1234
+ wetdep::clddiag (nlev, temperature, pmid, pdel, cmfdqr, evapc, cldt, cldcu,
1236
1235
cldst, evapr, prain,
1237
1236
// outputs
1238
1237
cldv, cldvcu, cldvst, rain);
1239
1238
});
1240
1239
}
1241
1240
template <typename VectorIntType>
1242
- KOKKOS_INLINE_FUNCTION
1243
- void set_f_act (const ThreadTeam &team, VectorIntType& isprx,
1244
- const View1D &f_act_conv_coarse,
1245
- const View1D &f_act_conv_coarse_dust,
1246
- const View1D &f_act_conv_coarse_nacl,
1247
- haero::ConstColumnView pdel, haero::ConstColumnView prain,
1248
- const View1D &cmfdqr, const ConstView1D &evapr,
1249
- const View2D &state_q, const View2D &ptend_q, const Real dt,
1250
- const int nlev) {
1241
+ KOKKOS_INLINE_FUNCTION void
1242
+ set_f_act (const ThreadTeam &team, VectorIntType &isprx,
1243
+ const View1D &f_act_conv_coarse, const View1D &f_act_conv_coarse_dust,
1244
+ const View1D &f_act_conv_coarse_nacl, haero::ConstColumnView pdel,
1245
+ haero::ConstColumnView prain, const View1D &cmfdqr,
1246
+ const ConstView1D &evapr, const View2D &state_q,
1247
+ const View2D &ptend_q, const Real dt, const int nlev) {
1251
1248
1252
1249
Kokkos::parallel_for (Kokkos::TeamVectorRange (team, nlev), [&](int k) {
1253
- isprx[k] = aero_model::examine_prec_exist (k, pdel, prain,
1254
- cmfdqr, evapr);
1250
+ isprx[k] = aero_model::examine_prec_exist (k, pdel, prain, cmfdqr, evapr);
1255
1251
1256
1252
aero_model::set_f_act_coarse (k, state_q, ptend_q, dt, f_act_conv_coarse[k],
1257
1253
f_act_conv_coarse_dust[k],
@@ -1261,13 +1257,11 @@ void set_f_act(const ThreadTeam &team, VectorIntType& isprx,
1261
1257
1262
1258
// Computes lookup table for aerosol impaction/interception scavenging rates
1263
1259
template <typename VectorIntType>
1264
- KOKKOS_INLINE_FUNCTION
1265
- void modal_aero_bcscavcoef_get (const ThreadTeam &team, const Diagnostics &diags,
1266
- const VectorIntType& isprx, const View2D &scavimptblvol,
1267
- const View2D &scavimptblnum,
1268
- const View1D &scavcoefnum,
1269
- const View1D &scavcoefvol, const int imode,
1270
- const int nlev) {
1260
+ KOKKOS_INLINE_FUNCTION void modal_aero_bcscavcoef_get (
1261
+ const ThreadTeam &team, const Diagnostics &diags,
1262
+ const VectorIntType &isprx, const View2D &scavimptblvol,
1263
+ const View2D &scavimptblnum, const View1D &scavcoefnum,
1264
+ const View1D &scavcoefvol, const int imode, const int nlev) {
1271
1265
Kokkos::parallel_for (Kokkos::TeamVectorRange (team, nlev), [&](int k) {
1272
1266
scavcoefnum[k] = scavcoefvol[k] = 0 ;
1273
1267
const bool let_it_rain = (isprx[k] == 1 );
@@ -1284,14 +1278,11 @@ void modal_aero_bcscavcoef_get(const ThreadTeam &team, const Diagnostics &diags,
1284
1278
1285
1279
// Computes lookup table for aerosol impaction/interception scavenging rates
1286
1280
template <typename VectorIntType>
1287
- KOKKOS_INLINE_FUNCTION
1288
- void modal_aero_bcscavcoef_get (const ThreadTeam &team,
1289
- const View2D &wet_geometric_mean_diameter_i,
1290
- const VectorIntType & isprx, const View2D &scavimptblvol,
1291
- const View2D &scavimptblnum,
1292
- const View1D &scavcoefnum,
1293
- const View1D &scavcoefvol, const int imode,
1294
- const int nlev) {
1281
+ KOKKOS_INLINE_FUNCTION void modal_aero_bcscavcoef_get (
1282
+ const ThreadTeam &team, const View2D &wet_geometric_mean_diameter_i,
1283
+ const VectorIntType &isprx, const View2D &scavimptblvol,
1284
+ const View2D &scavimptblnum, const View1D &scavcoefnum,
1285
+ const View1D &scavcoefvol, const int imode, const int nlev) {
1295
1286
Kokkos::parallel_for (Kokkos::TeamVectorRange (team, nlev), [&](int k) {
1296
1287
scavcoefnum[k] = scavcoefvol[k] = 0 ;
1297
1288
const bool let_it_rain = (isprx[k] == 1 );
@@ -1317,9 +1308,8 @@ void define_act_frac(const ThreadTeam &team, const View1D &sol_facti,
1317
1308
});
1318
1309
}
1319
1310
template <typename VectorType>
1320
- KOKKOS_INLINE_FUNCTION
1321
- void compute_q_tendencies_phase_1 (
1322
- Real &scavt, Real &bcscavt, Real &rcscavt, const VectorType& rtscavt_sv,
1311
+ KOKKOS_INLINE_FUNCTION void compute_q_tendencies_phase_1 (
1312
+ Real &scavt, Real &bcscavt, Real &rcscavt, const VectorType &rtscavt_sv,
1323
1313
const Real f_act_conv, const Real scavcoefnum, const Real scavcoefvol,
1324
1314
const Real totcond, const Real cmfdqr, const Real conicw, const Real evapc,
1325
1315
const Real evapr, const Real prain, const Real dlf, const Real cldt,
@@ -1360,9 +1350,8 @@ void compute_q_tendencies_phase_1(
1360
1350
rtscavt_sv);
1361
1351
}
1362
1352
template <typename VectorType>
1363
- KOKKOS_INLINE_FUNCTION
1364
- void compute_q_tendencies_phase_2 (
1365
- Real &scavt, Real &bcscavt, Real &rcscavt, const VectorType& rtscavt_sv,
1353
+ KOKKOS_INLINE_FUNCTION void compute_q_tendencies_phase_2 (
1354
+ Real &scavt, Real &bcscavt, Real &rcscavt, const VectorType &rtscavt_sv,
1366
1355
const Real qqcw_tmp, const Real tracer,
1367
1356
1368
1357
// const Prognostics &progs,
@@ -1843,8 +1832,8 @@ void aero_model_wetdep(
1843
1832
mam4::water_uptake::modal_aero_water_uptake_dr (
1844
1833
// inputs
1845
1834
calcsizedata.nspec_amode , calcsizedata.specdens_amode ,
1846
- calcsizedata.spechygro , calcsizedata.lspectype_amode ,
1847
- state_q_kk, temperature (kk), pmid (kk), cldt (kk), dgnumdry_m_kk,
1835
+ calcsizedata.spechygro , calcsizedata.lspectype_amode , state_q_kk,
1836
+ temperature (kk), pmid (kk), cldt (kk), dgnumdry_m_kk,
1848
1837
// outputs
1849
1838
dgnumwet_m_kk, qaerwat_m_kk, wetdens_kk);
1850
1839
}
0 commit comments