Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

How to run WRF_VPRM in WRF4.3?

zhengzich

New member
Hi guys, I'm currently learning to use WRF_VPRM.
I carefully read the WRF_GHG manual and the WRF-CHEM user guide 4.4.
But I have a question now, I use WRF-CHEM version 4.3, I hope to be able to run WRF-VPRM with chem_opt=16, use VPRM to provide fluxes and combine with other lists to calculate changes in CO 2 concentration.

Do I need to just use WRF4.3, or do I need to take the WRF-GHG code and combine it with WRF4.3?

I know this question may be a little silly, but it's important to me and I look forward to your answers.
 
Last edited:
I assume you refer to Beck et al (2011) when talking about the WRF-GHG manual? This is unfortunally an outdated document as since then (2012) the VPRM module was implemented in the WRF-Chem code. So generally, you don't need to adjust the model code.

The main task that you need to do in order to 'activate' the VPRM module is provide the correct input fields.
WRF-Chem expects a file (for every domain and time step) containing VEGFRA_VPRM, EVI, EVI_MIN, EVI_MAX, LSWI, LSWI_MIN and LSWI_MAX, which then should be refered to by using auxinput15 in the namelist.

Additionally, make sure that you set bio_emiss_opt=16 or 17 in the namelist.
 
I assume you refer to Beck et al (2011) when talking about the WRF-GHG manual? This is unfortunally an outdated document as since then (2012) the VPRM module was implemented in the WRF-Chem code. So generally, you don't need to adjust the model code.

The main task that you need to do in order to 'activate' the VPRM module is provide the correct input fields.
WRF-Chem expects a file (for every domain and time step) containing VEGFRA_VPRM, EVI, EVI_MIN, EVI_MAX, LSWI, LSWI_MIN and LSWI_MAX, which then should be refered to by using auxinput15 in the namelist.

Additionally, make sure that you set bio_emiss_opt=16 or 17 in the namelist.
Hello, I just started learning WRF-VPRM, and I would like to ask you how to generate files containing VEGFRA_VPRM, EVI, EVI_MIN, EVI_MAX, LSWI, LSWI_MIN and LSWI_MAX data.

I am going to try VPRM_Preprocessor, but LPDAAC no longer provides the Linux version of MRT; then I tried pyVPRM, but I also encountered some problems; so I would like to ask if you have any recommended methods for generating VPRM module data.

Your help is very important to me, and I look forward to your reply.
 
Dear zengchang fan

I personally only have experience with the (old) VPRM Preprocessor as it was installed together with the other required tools a few years ago on the computing system that I use. I would assume the best approach nowadays is to try using pyVPRM, which I have not tried yet myself.
Maybe you can try contacting the creators of the tool? Or you can be more specific about the problems you encounter there and hope that some other users on this forum can help you with that. Good luck!
 
Dear zengchang fan

I personally only have experience with the (old) VPRM Preprocessor as it was installed together with the other required tools a few years ago on the computing system that I use. I would assume the best approach nowadays is to try using pyVPRM, which I have not tried yet myself.
Maybe you can try contacting the creators of the tool? Or you can be more specific about the problems you encounter there and hope that some other users on this forum can help you with that. Good luck!
Dear callewaerts

Thanks for your reply, it will help me to clarify the focus of the next step.
Specifically, in the process of using pyVPRM, the packages in the environment have errors one after another, but I am trying to solve these problems one by one.
Thanks again for your advice, I wish you all the best!
 
Dear callewaerts

Thanks for your reply, it will help me to clarify the focus of the next step.
Specifically, in the process of using pyVPRM, the packages in the environment have errors one after another, but I am trying to solve these problems one by one.
Thanks again for your advice, I wish you all the best!
I haven't used py processing, but VPRMpreproc does require some specific debugging. In this calculation process, it is mainly necessary to calculate a percentage based on the local class in order to enter the WRF, which is a more critical step.
 
I haven't used py processing, but VPRMpreproc does require some specific debugging. In this calculation process, it is mainly necessary to calculate a percentage based on the local class in order to enter the WRF, which is a more critical step.
pyVPRM also needs debugging, but more problems occur in the packages in the environment. Thank you for sharing your experience, I will pay attention to the debugging of VPRM parameters.
 
Hello, I just started learning WRF-VPRM, and I would like to ask you how to generate files containing VEGFRA_VPRM, EVI, EVI_MIN, EVI_MAX, LSWI, LSWI_MIN and LSWI_MAX data.

I am going to try VPRM_Preprocessor, but LPDAAC no longer provides the Linux version of MRT; then I tried pyVPRM, but I also encountered some problems; so I would like to ask if you have any recommended methods for generating VPRM module data.

Your help is very important to me, and I look forward to your reply.
你好,我也在研究VPRM,请问可以加个联系方式交流一下吗
 
i've been trying to use PyVPRM too, but have got stuck. The examples read in NetCDF files from MODIS, but the files downloaded from MODIS are in HDF4, which causes it not to succeed. Has anyone else encountered this, and how have they solved it?
 
i've been trying to use PyVPRM too, but have got stuck. The examples read in NetCDF files from MODIS, but the files downloaded from MODIS are in HDF4, which causes it not to succeed. Has anyone else encountered this, and how have they solved it?
There is no need to convert the files into Netcdf. You can simply use the HDF files that you download from the MODIS server. And if you are doing wrf-vprm, we can communicate about that!!!
 
I have now run PyVPRM successfully, with some help from the author (Theo Glauch) also; I noticed he has edited the readme to include this information too. For additional reference, I had to edit the vprm_preprocessor_new.py file on l.68 to target *.hdf files rather than *.nc files, and also install gdal and rasterio in my conda installation, which I hadn't done.
 
Top