@@ -207,59 +207,63 @@ subroutine physconst_readnl(nlfile)
207
207
208
208
if (newg .or. newsday .or. newmwh2o .or. newcpwv .or. newmwdry .or. &
209
209
newrearth .or. newtmelt .or. newomega) then
210
+
211
+ ! Populate the new constants into module after mpi_bcast
212
+ if (newg) gravit = user_defined_gravit
213
+ if (newsday) sday = user_defined_sday
214
+ if (newmwh2o) mwh2o = user_defined_mwh2o
215
+ if (newcpwv) cpwv = user_defined_cpwv
216
+ if (newmwdry) mwdry = user_defined_mwdry
217
+ if (newcpair) cpair = user_defined_cpair
218
+ if (newrearth) rearth = user_defined_rearth
219
+ if (newtmelt) tmelt = user_defined_tmelt
220
+ if (newomega) omega = user_defined_omega
221
+
210
222
if (masterproc) then
211
223
write (iulog, * ) banner
212
224
write (iulog, * ) ' *** New Physical Constant Values set ' , &
213
225
' via namelist ***'
214
226
write (iulog, * ) bline
215
227
write (iulog, * ) ' *** Physical Constant Old Value New Value ***'
216
228
if (newg) then
217
- gravit = user_defined_gravit
218
229
field = ' GRAVIT'
219
230
write (iulog, 2000 ) field, shr_const_g, gravit
220
231
end if
221
232
if (newsday) then
222
- sday = user_defined_sday
223
233
field = ' SDAY'
224
234
write (iulog, 2000 ) field, shr_const_sday, sday
225
235
end if
226
236
if (newmwh2o) then
227
- mwh2o = user_defined_mwh2o
228
237
field = ' MWH20'
229
238
write (iulog, 2000 ) field, shr_const_mwwv, mwh2o
230
239
end if
231
240
if (newcpwv) then
232
- cpwv = user_defined_cpwv
233
241
field = ' CPWV'
234
242
write (iulog, 2000 ) field, shr_const_cpwv, cpwv
235
243
end if
236
244
if (newmwdry) then
237
- mwdry = user_defined_mwdry
238
245
field = ' MWDRY'
239
246
write (iulog, 2000 ) field, shr_const_mwdair, mwdry
240
247
end if
241
248
if (newcpair) then
242
- cpair = user_defined_cpair
243
249
field = ' CPAIR'
244
250
write (iulog, 2000 ) field, shr_const_cpdair, cpair
245
251
end if
246
252
if (newrearth) then
247
- rearth = user_defined_rearth
248
253
field = ' REARTH'
249
254
write (iulog, 2000 ) field, shr_const_rearth, rearth
250
255
end if
251
256
if (newtmelt) then
252
- tmelt = user_defined_tmelt
253
257
field = ' TMELT'
254
258
write (iulog, 2000 ) field, shr_const_tkfrz, tmelt
255
259
end if
256
260
if (newomega) then
257
- omega = user_defined_omega
258
261
field = ' OMEGA'
259
262
write (iulog, 2000 ) field, shr_const_omega, omega
260
263
end if
261
264
write (iulog,* ) banner
262
265
end if
266
+
263
267
rga = 1._kind_phys / gravit
264
268
rearth_recip = 1._kind_phys / rearth
265
269
if (.not. newomega) then
0 commit comments