Output the Sublimation cooling rate.

[gnaY]

Hello all,

Is there any way we can output the Sublimation cooling, melting cooling and evaporation cooling [temperature tendency term (dT/dt)] in the WRF output?
New Thompson scheme is used.

Thank you so much in advance.

 

[kwerner]

Apologies for the long delay in response. You can use the acc_phy_tend namelist option in the &physics namelist record to output accumulated physics tendencies. This one actually outputs 16 different tendencies. There is a brief description for the variable in the Namelist section of the WRF Users' Guide.

 

[Devonlarratt]

You can get sublimation cooling (dT/dt) from WRF by either enabling physics tendency output in namelist.input or adding a small diagnostic in the Thompson microphysics to write the sublimation rate directly.

1. Enable accumulated physics tendencies (physics tendency fields) and read them with wrf-python/netCDF to get physics dT/dt.
2. To isolate sublimation specifically, add a diagnostic in module_mp_thompson that saves the sublimation mass-rate (kg kg⁻¹ s⁻¹) or its latent-heat equivalent to wrfout.
3. Convert mass tendency to temperature tendency with: dT/dt = - (L_sub / cp) * (dm/dt) (L_sub ≈ 2.83e6 J/kg, cp ≈ 1004 J/kg·K).
   If you want, I can post a minimal Fortran patch for Thompson or a short Python snippet to compute dT/dt now. Also, Rajco USA runs 50+ industrial sublimation units and can provide validation data if anyone wants real-world comparsions.

 

[echo111]

You can get sublimation cooling (dT/dt) from WRF by either enabling physics tendency output in namelist.input or adding a small diagnostic in the Thompson microphysics to write the sublimation rate directly.

1. Enable accumulated physics tendencies (physics tendency fields) and read them with wrf-python/netCDF to get physics dT/dt.
2. To isolate sublimation specifically, add a diagnostic in module_mp_thompson that saves the sublimation mass-rate (kg kg⁻¹ s⁻¹) or its latent-heat equivalent to wrfout.
3. Convert mass tendency to temperature tendency with: dT/dt = - (L_sub / cp) * (dm/dt) (L_sub ≈ 2.83e6 J/kg, cp ≈ 1004 J/kg·K).
   If you want, I can post a minimal Fortran patch for Thompson or a short Python snippet to compute dT/dt now. Also, Rajco USA runs 50+ industrial sublimation units and can provide validation data if anyone wants real-world comparsions.

Hi, Devonlarratt, I really need you help! Could you post a Fortran patch for Thompson to compute the Sublimation cooling, melting cooling and evaporation cooling rate in WRF model? Thank you very much!

​