Sensitivity test on QVAPOR

[Chetan]

Hi everyone,

I want to do the sensitivity test on QVAPOR by modifying it for the first few vertical layers in wrfrst file and restarting the simulation from that point. To do this, do I also need to modify the QVAPOR in the wrfbdy file at the same time? If so, how to modify the QVAPOR_B[][] in the boundary file, which seems much larger compared to QVAPOR in the wrfrst file?.



Thank you in advance.

Best,
Chetan

 

[kwerner]

Hi Chetan,
Apologies for the delay. Yes, you will need to modify the boundary conditions, as well, since the restart program uses both the wrfbdy file, and the wrfrst file(s). The wrfbdy_d01 file is in netCDF format, so you would modify it, as you would any other netCDF file. Unfortunately I'm not aware of a simple way to do this, but it's possible others have done something similar, and may have a script they could share. It's possible you can find some helpful suggestions by searching the content in this forum and/or a web search.

 

[Chetan]

Thank you @kwerner for your response.

 

[ztx]

Hi everyone,

I want to do the sensitivity test on QVAPOR by modifying it for the first few vertical layers in wrfrst file and restarting the simulation from that point. To do this, do I also need to modify the QVAPOR in the wrfbdy file at the same time? If so, how to modify the QVAPOR_B[][] in the boundary file, which seems much larger compared to QVAPOR in the wrfrst file?.



Thank you in advance.

Best,
Chetan

Hi Chetan,
This post might answer your question about why there is a significant difference of the value of QVAPOR between wrfbdy and wrfrst: Definitions of _BTXE and _BXE in wrfbdy output
In a word, the reason is as the member davegill said in that post: "Each of these 3d fields is "coupled" or "mass weighted" with the 3d value of d(Pd)/d(eta). These d(Pd)/d(eta) values range from 0 to about 10^5, and are effectively used as scale factors. They DRAMATICALLY change the interpretation of the value of the original field."

About how to modify the wrfbdy netCDF file, I think the most hard part maybe how to get the 3d value d(Pd)/d(eta), so you can transform the original value to the value in wrfbdy. But I have not yet find a script that can do this.