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

WRF_VPRM have no negative value in CO2_BIO and CO2_OCE

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.


New member
Hi there,
I ran WRF_VPRM recently but got all positive values in CO2_BIO and CO2_OCE. It should have some negative values but seemed to just output the values larger than zero.
Any suggestion is appreciated.

namelist.input is attached.


  • namelist.input
    12.1 KB · Views: 127

I believe these quantities are concentrations, so cannot be zero. Are you interested in the fluxes? e.g., see this line in module_ghg_fluxes.F

chem(i,1,j,p_co2_bio)= chem(i,1,j,p_co2_bio) + conv_rho* (eghg_bio(i,1,j,p_ebio_gee) + eghg_bio(i,1,j,p_ebio_res)) ! both uptake and release

Dear Pwang
Mike is right.
I believe what really happened is that the advection scheme is positive definite. So the model indeed only keeps positive concentrations.
So to solve this issue, either adding a constant offset to both initial and boundary conditions, e.g. 400 ppm, as Mike said.
or using global CO2 simulation outputs (e.g., CarbonTracker) as initial and boundary conditions, as we did in our regional simulations:

Xiao-Ming Hu