WRF & MPAS-A Support Forum

I do not how to use the pyvprm tool. The problem I see is that your evi and lswi have many values >10^36, which I guess results from some "blank" values. You may initial all values with zeros instead of leaving them blanks when running your scripts.

Or just try to modify your vprm_input file by changing unreasonable values to zeros manually (both evi and lswi cannot be >1), and then change the max and min values accordingly.

I just created my new input file now it has value between 0 to 1. Please check this file

You then need to check the variable names asked by WRF-Chem in Registry/registry.chem file. I think some of your variables are not correctly named. Then you can try to run wrf.exe and check other i/o issues, if any.

Hi
I checked the variables and set the vprm file accordingly but still i have an issue when i check in GEE and RES variables. I am attaching my vprm file and screenshot of wrfout file please have a look
thanks

The variable names should be in upper case. Note also that the first step of GEE and RES are normally zeros, and you should check the latter steps.

yes i already also checked with upper case but still same issue for all steps..I am sharing my namelist please check it whether i activated vprm functionality correctly or not

I do not see your new namelist and from your old namelist I can find no problem. The VEGFRA_VPRM need Time dimension (I am not sure about EVI_MAX or so). For further debugging, you may need to add some simple code in the VPRM module (chem/module_ghg_fluxes.f90) to print GEE/RES or vprm_in.

I got the values in both variables but in EBIO_GEE i have negative values and in EBIO_RES has +ve values is it right?
If its right then for now NEE should do this one.. NEE = -GEE + RES?

Yes, it is right since GEE is sink and RES is source. NEE=GEE+RES, not -GEE+RES.

My both variables (GEE and RES) has mol km-2/hr but CO2_BCK and ANT is in ppmv so can i change my VPRM variables units like this

EBIO_GEE * 0.04401 / (1e6 * 3600)
EBIO_RES * 0.04401 / (1e6 * 3600)

No. EBIO_GEE is flux and CO2_BCK is concentration. They should have different units and you cannot compare concentration with flux.

i want to calculate co2_total after that i will validate my results with observational data so what should i do?

The common approach is to add CO2_BIO, CO2_ANT, and CO2_OCE (and others if any) and subtracting CO2_BCK for proper times (CO2_BIO or so includes background to avoid negative concentrations during simulation).

Hi
I have an issue regarding IC BC.
My input MOZBC file contains CO2 values ranging from 403 to 500+ ppm, and in my study area the expected value is approximately 412 to 430 ppm. However, after running MOZBC, it populates CO2 background data into the wrfinput and wrfbdy files and the resulting CO2 is much narrower between 404.2 and 404.8 ppm.
Could you kindly guide me on what might be causing this issue?

I'm not familiar with MOZBC and did not met similar problems before. The following is just some conjecture. If the values vary in a narrow interval, it may result from a wrong grid mapping that map the values in a small areas to your domain. This can be checked by seeing the spatial distribution of the output values. You then may check your dimension variables in your input files.

Another possibility is that the MOZBC failed somewhat and some default values are placed. But since your output values are not strictly constant, this is not very likely the case. Did you try a normal WRF-Chem run for other species than CO2? Was the MOZBC output normal?

no i did not try ok so any other way to populate wrfinput and wrfbdy files with co2 parameter?

When I used WRF-VPRM, I used MOZBC and succeed in generating ICBC. So I did not have experience in debugging or so. I currently use (VPRM-)CMAQ which has its own tools, and I do not know other tools for wrfchem ICBC.

can you please share your input sample file of mozbc in which co2 parameters exists. File can be only 1 or 2 days because i just want to check it