Comments (8)
turning back on shedding seems to fix this... will close this issue once I confirm by completing the PalEON millennial runs
from ed2.
In the above diagnostics where it gives the D(CAN_CO2) error it appears that CAN_CO2 is 30 ppm, which is impossible. You might want to check your drivers to make sure there are not zeros in there somewhere (e.g. accidentally padding the end of a month). That could easily explain the D(CAN_CO2), if the model's trying to do a step change from 270 ppm to 0 ppm
from ed2.
After digging into the role of shedding in triggering the CO2 error a bit, it looks like the error may be stemming from how we constraining the amount of wood/leaf evaporation (wflxwc, wflxlc) around lines 1409 & 1749 in rk4_derivs.f90.
Examples shown below use wflxwc as example, but all of it applies to wflxlc too.
Currently constrain on wflx to canopy:
wflxwc = min(wflxwc,max_dwdt+wood_intercepted-wshed)
The CO2 error appears to be tied to the second option: max_dwdt+wood_intercepted-wshed
where
max_dwdt = initp%wood_water(ico)/dt
wshed = max(0.d0,( (initp%wood_water(ico) + wood_intercepted*dt) - max_wood_water) / dt)
If we get rid of wood_intercepted from the wflcwc the problem seems to be resolved (at least I've gotten a couple decades past where I was so far). My rationale for trying this was that the shed calculation takes into account wood_intercepted already and so we might be double-accounting here and incorrectly constraining the leaf/wood water fluxes to the canopy.
Thoughts?
from ed2.
New Update:
It looks like the draw-down of Canopy CO2 is from GPP being extremely high as a result of the the vm0 I specified for my dominant PFT based on estimates from Pecan.
I think the high vm0 wouldn't be a problem if I knew how to tweak the respiration parameters accordingly. This looks like it links in to some of the other conversations happening to add some sort of sanity or balance check with the PFT parameters to make sure that the equations balance when you're modifying individual parameters.
Alternatively, we could think about implementing some sort of max GPP that would cause GPP to saturate at a given number and help prevent parameterization insanity.
from ed2.
It looks like this problem happening in a site with short vegetation, or a young patch, right? I've had similar problems with grassland sites in the past. The canopy air space may become completely "detached" from the air aloft when conditions are stable and u* is weak, and ED tends to understimate u*. In grasslands the canopy air space is thin, so this issue could cause some unrealistic properties if it lasts very long.
Which settings have you been using for the surface layer model? If you are running with ISFCLYR=3 or 4, then you could try increasing the minimum wind (1m/s) and minimum u*(~0.1m/s), and reducing the maximum bulk Richardson number (about 0.5). It's not the best solution, but it should avoid these super-stable layers.
from ed2.
If the Vm0 is "extremely high", then I think you should adjust it rather than tweaking the respiration parameters since changing Vm0 affects other parameters (e.g. C/N ratio in leaf).
from ed2.
Hi @apourmok and @crollinson isn't "Vm0" the base rate (i.e. 15 degC) value of photosynthetic capacity (i.e. Vcmax)? Changing that value should change C:N ratio, for example. Also, the old way to calculate Rd0 in ED2 was:
Rd0 = dark_respiration_factor*Vm0
However, I thought we removed that old logic and went to a more appropriate approach of prescribing the Rd0 (i.e. leaf base respiration rate)? In other words, using the XML to prescribe both Vm0 and Rd0.
from ed2.
@mpaiao: You're right about the age thing: it usually happens in my youngest or next-youngest stand and had previously been tied to canopy roughness and decoupling of the canopy from the atmosphere.
@serbinsh: Based on my understanding of the mainline version of ED2, Vm0 = Vcmax and Rd0 is calculated as you have stated and not prescribed in the xml. You can specify D0 (which I believe is different) and dark_respiration_factor, but not Rd0 itself. For the most part when I bumped up vm0 from PFT 8, I left respiration alone since I didn't have any estimates for that PFT from Pecan. This could be part of my problem.
It looks like each of these solutions (decrease vm0, increase dark_rep, increase wind thresholds) fixes the problem without altering the code, so now I get to figure out which combination of factors makes the most sense. Thanks for everyone's help!
from ed2.
Related Issues (20)
- base flow goes to negative HOT 5
- Grass cohorts causing Floating-point exception at log() operations, several places HOT 11
- little mistake about the unit description in the output file HOT 5
- how to modify the settings to increase the NPLANT of pft 3 ? ? HOT 2
- Issue with HPC build HOT 6
- soil water off-track HOT 1
- Problem writing the T-file HOT 12
- NPLANT
- little mistake in adjust_veg_properties of rk4_misc HOT 3
- question about unit in reproduction.f90 HOT 3
- Typo in the events.f90 file HOT 3
- Docker image not being updated by GH actions
- `HYDAULIC_MORTALITY_SCHEME` missing from ED2IN HOT 4
- ED2 Installation error
- Patch Bfusion budget failed (CO2 Summary: FINE: F)
- hdf5 problem when running HISTORY mode HOT 4
- question about MPI HOT 1
- Crop yield estimation
- psi_max or psi_min in leaf phenology HOT 1
- typo in calc_plant_water_flux HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from ed2.