Closed rosiealice closed 6 years ago
jkshuman
commented
7 years ago Attaching an updated version of fates parameter file. Updated on April 12 with new values to spitfire specific parameters, and appears to fix issue #190 CWD_Stock_Col stabilizes just above 200 gC per m2 while biomass_stock_col stabilizes around 2600 gC per m2 for a single point run at Caxiuana, Brazil.
tropical_params_decomp_18Apr2017.nc.gz
This code is functional with same version referenced above.
ckoven
commented
7 years ago all- I've been playing around with parameter files again, trying to get (a) at a minimum, a healthy canopy with biomass in the 10-20 kg C / m2 range for a 1x1amazon default run, and (b) to see what I can do parametrically to help out the understory. I've got something that is doing a lot better, in that the biomass is up where it should be after spinup, and the understory isn't dying quite as fast as before. Not sure hwo else is playing around with parameters, but figured it may be of some interest.
I did this, starting from Jackie's file, by:
I also made it so that only a single PFT is active right now (via initd and seed_rain, the latter of which by the way we should turn off), and I had also played around with fates_flnr before I remembered that that doesn't do anything on the code right now (since vcmax25top is set directly via the fnitr parameter, which we ought to make more transparent)
And lastly, there were a few fairly recently added variables (fates_displar, fates_z0mr, fates_dbh2bl_dbh2carea_expnt_diff) and one recently-deleted variable (fates_profile_tol) that either are or aren't there now. For sanity's sake, I'll put the diff of the cdl files here:
1c1
< netcdf tropical_params_decomp_18Apr2017 {
---
> netcdf tropical_params_ckmod_01May2017 {
203a204,206
> double fates_dbh2bl_dbh2carea_expnt_diff(fates_pft) ;
> fates_dbh2bl_dbh2carea_expnt_diff:long_name = "difference between allometric DBH:bleaf and DBH:crown area exponents" ;
> fates_dbh2bl_dbh2carea_expnt_diff:units = "unitless" ;
331c334,341
< double fates_profile_tol(fates_scalar) ;
---
> double fates_displar(fates_pft) ;
> fates_displar:long_name = "Ratio of displacement height to canopy top height" ;
> fates_displar:units = "unitless" ;
> fates_displar:coordinates = "pftname" ;
> double fates_z0mr(fates_pft) ;
> fates_z0mr:long_name = "Ratio of momentum roughness length to canopy top height" ;
> fates_z0mr:units = "unitless" ;
> fates_z0mr:coordinates = "pftname" ;
361,367c371,375
< fates_flnr = 0, 0.0509, 0.0466, 0.0546, 0.0461, 0.0515, 0.0716, 0.1007,
< 0.1007, 0.0517, 0.0943, 0.0943, 0.1365, 0.1365, 0.09, 0.1758, 0.1758,
< 0.293, 0.293, 0.4102, 0.4102, 0.4102, 0.4102, 0.4102, 0.4102, 0.4102,
< 0.4102, 0.4102, 0.4102, 0.4102, 0.4102, 0.4102, 0.4102, 0, 0, 0, 0, 0, 0,
< 0, 0, 0.4102, 0.4102, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
< 0, 0.4102, 0.4102, 0, 0, 0, 0, 0.293, 0.293, 0, 0, 0, 0, 0, 0, 0.293,
< 0.293, 0.4102, 0.4102 ;
---
> fates_flnr = 0, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
> 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
> 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
> 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
> 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 ;
424c432
< fates_leafcn = 1, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30,
---
> fates_leafcn = 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30,
620,624c628,633
< fates_max_dbh = 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68,
< 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68,
< 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68,
< 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68,
< 68, 68, 68, 68, 68, 68, 68, 68, 68, 68 ;
---
> fates_max_dbh = 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150,
> 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150,
> 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150,
> 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150,
> 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150,
> 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150 ;
626,627c635,636
< fates_cushion = 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2,
< 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, _, _, _,
---
> fates_cushion = 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0,
> 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, _, _, _,
680,685c689,692
< fates_initd = 0.08, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2,
< 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2,
< 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2,
< 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2,
< 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2,
< 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2 ;
---
> fates_initd = 0, 0.2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ;
693,695c700,701
< fates_seed_rain = 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28,
< 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28,
< 0.28, 0.28, 0.28, 0.28, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
---
> fates_seed_rain = 0, 0.28, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
697c703
< _, _, _, _, _, _, _, _, _, _, _, _, _, _ ;
---
> _, _, _, _, _, _, _, _, _, _, _, _, _ ;
731c737
< fates_dbh2bl_a = 0.0419, 0.0419, 0.0419, 0.0419, _, _, _, _, _, _, _, _, _,
---
> fates_dbh2bl_a = .021, .021, .021, .021, _, _, _, _, _, _, _, _, _,
736c742
< fates_dbh2bl_b = 1.56, 1.56, 1.56, 1.56, _, _, _, _, _, _, _, _, _, _, _, _,
---
> fates_dbh2bl_b = 1.3, 1.3, 1.3, 1.3, _, _, _, _, _, _, _, _, _, _, _, _,
750a757,761
> fates_dbh2bl_dbh2carea_expnt_diff = 0, 0, 0, 0, _, _, _, _, _, _, _, _, _,
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _ ;
>
931c942,956
< fates_profile_tol = 0.1 ;
---
> fates_displar = 0, 0.67, 0.67, 0.67, 0.67, 0.67, 0.67, 0.67, 0.67, 0.68,
> 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68,
> 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68,
> 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68,
> 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68,
> 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68,
> 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68 ;
>
> fates_z0mr = 0, 0.055, 0.055, 0.055, 0.075, 0.075, 0.055, 0.055, 0.055,
> 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12,
> 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12,
> 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12,
> 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12,
> 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12,
> 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12 ;Thanks, Charlie. That is great news regarding understory!
On Mon, May 1, 2017 at 2:18 PM, Charlie Koven notifications@github.com wrote:
all- I've been playing around with parameter files again, trying to get (a) at a minimum, a healthy canopy with biomass in the 10-20 kg C / m2 range for a 1x1amazon default run, and (b) to see what I can do parametrically to help out the understory. I've got something that is doing a lot better, in that the biomass is up where it should be after spinup, and the understory isn't dying quite as fast as before. Not sure hwo else is playing around with parameters, but figured it may be of some interest.
I did this, starting from Jackie's file, by:
- increasing the cushion all the way to 2, so that the canopy trees have enough buffer to last through a few bad years before they die.
- decreasing the canopy exponent fates_dbh2bl_b to 1.3 so that it is closer to observed and also so that we can get more large trees into the canopy ( #206 https://github.com/NGEET/fates/issues/206 )
- decreasing the bleaf coefficient so that the canopy lai isn't so thick, and therefore the understory trees get a bit more life. I reduced it all the way to fates_dbh2bl_a = .021 here, which seems to be a point at which understory survival starts to be slightly less depressing.
I also made it so that only a single PFT is active right now (via initd and seed_rain, the latter of which by the way we should turn off), and I had also played around with fates_flnr before I remembered that that doesn't do anything on the code right now (since vcmax25top is set directly via the fnitr parameter, which we ought to make more transparent)
And lastly, there were a few fairly recently added variables (fates_displar, fates_z0mr, fates_dbh2bl_dbh2carea_expnt_diff) and one recently-deleted variable (fates_profile_tol) that either are or aren't there now. For sanity's sake, I'll put the diff of the cdl files here:
1c1 < netcdf tropical_params_decomp_18Apr2017 {
netcdf tropical_params_ckmod_01May2017 { 203a204,206 double fates_dbh2bl_dbh2carea_expnt_diff(fates_pft) ; fates_dbh2bl_dbh2carea_expnt_diff:long_name = "difference between allometric DBH:bleaf and DBH:crown area exponents" ; fates_dbh2bl_dbh2carea_expnt_diff:units = "unitless" ; 331c334,341 < double fates_profile_tol(fates_scalar) ;
double fates_displar(fates_pft) ; fates_displar:long_name = "Ratio of displacement height to canopy top height" ; fates_displar:units = "unitless" ; fates_displar:coordinates = "pftname" ; double fates_z0mr(fates_pft) ; fates_z0mr:long_name = "Ratio of momentum roughness length to canopy top height" ; fates_z0mr:units = "unitless" ; fates_z0mr:coordinates = "pftname" ;361,367c371,375 < fates_flnr = 0, 0.0509, 0.0466, 0.0546, 0.0461, 0.0515, 0.0716, 0.1007, < 0.1007, 0.0517, 0.0943, 0.0943, 0.1365, 0.1365, 0.09, 0.1758, 0.1758, < 0.293, 0.293, 0.4102, 0.4102, 0.4102, 0.4102, 0.4102, 0.4102, 0.4102, < 0.4102, 0.4102, 0.4102, 0.4102, 0.4102, 0.4102, 0.4102, 0, 0, 0, 0, 0, 0, < 0, 0, 0.4102, 0.4102, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, < 0, 0.4102, 0.4102, 0, 0, 0, 0, 0.293, 0.293, 0, 0, 0, 0, 0, 0, 0.293, < 0.293, 0.4102, 0.4102 ;
fates_flnr = 0, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 ; 424c432 < fates_leafcn = 1, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30,
fates_leafcn = 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 620,624c628,633 < fates_max_dbh = 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, < 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, < 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, < 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, < 68, 68, 68, 68, 68, 68, 68, 68, 68, 68 ;
fates_max_dbh = 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150, 150 ; 626,627c635,636 < fatescushion = 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, < 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, , , ,
fatescushion = 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, , , , 680,685c689,692 < fates_initd = 0.08, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, < 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, < 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, < 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, < 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, < 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2 ;
fates_initd = 0, 0.2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ; 693,695c700,701 < fates_seedrain = 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, < 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, < 0.28, 0.28, 0.28, 0.28, , , , , , , , , , , , , , , , _,
fates_seedrain = 0, 0.28, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, , , , , , , , , , , , , , , , , , 697c703 < , , , , , , , , , , , , , ;
_, _, _, _, _, _, _, _, _, _, _, _, _ ;731c737 < fates_dbh2bla = 0.0419, 0.0419, 0.0419, 0.0419, , , , , , , , , ,
fates_dbh2bla = .021, .021, .021, .021, , , , , , , , , , 736c742 < fates_dbh2blb = 1.56, 1.56, 1.56, 1.56, , , , , , , , , , , , _,
fates_dbh2blb = 1.3, 1.3, 1.3, 1.3, , , , , , , , , , , , _, 750a757,761 fates_dbh2bl_dbh2carea_expntdiff = 0, 0, 0, 0, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , ;
931c942,956 < fates_profile_tol = 0.1 ;
fates_displar = 0, 0.67, 0.67, 0.67, 0.67, 0.67, 0.67, 0.67, 0.67, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68, 0.68 ;
fates_z0mr = 0, 0.055, 0.055, 0.055, 0.075, 0.075, 0.055, 0.055, 0.055, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12, 0.12 ;
tropical_params_ckmod_01May2017.nc.zip https://github.com/NGEET/fates/files/968700/tropical_params_ckmod_01May2017.nc.zip
— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/NGEET/fates/issues/209#issuecomment-298421431, or mute the thread https://github.com/notifications/unsubscribe-auth/AVFDhp67pyAgeQdQZXC7hkKJoHykYw0_ks5r1j4KgaJpZM4M2Cks .
Dr. Jacquelyn Shuman
Terrestrial Sciences Section Climate and Global Dynamics National Center for Atmospheric Research 1850 Table Mesa Drive Boulder, Colorado, 80305 USA.
jkshuman@ucar.edu office: +1-303-497-1787
ckoven
commented
7 years ago All, just wanted to report that I've made some progress getting early/late successional coexistence using the strategy of specifying higher growth rates via vcmax in the early and forcing higher mortality rates in the early both via the background mortality as well as higher carbon starvation in the understory via increased respiration rates. file attached below, as well as diff of the file with respect to the May 01 file, and a couple figures that demonstrate the early/late coexistence seems stable. Note that this file (and possibly early/late coexistence in general) requires https://github.com/NGEET/fates/pull/245 for reasons that I haven't fully figured out yet.
I'll be out for a couple weeks but if others want to play around with these feel free...
[cdkoven@n0001 paramdata]$ diff tropical_params_ckmod_01May2017.cdl tropical_params_ckmod_2pfts_v14_strictppa_18jul2017.cdl
1c1
< netcdf tropical_params_ckmod_01May2017 {
---
> netcdf tropical_params_ckmod_2pfts_v9_13jul2017 {
378c378
< fates_fnitr = 0, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50,
---
> fates_fnitr = 0, 80, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50,
426c426
< fates_leaf_long = 0, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5,
---
> fates_leaf_long = 0, 0.75, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5,
433,437c433,437
< fates_leafcn = 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30,
< 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30,
< 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30,
< 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30,
< 30, 30, 30, 30, 30, 30, 30, 30, 30, 30 ;
---
> fates_leafcn = 40, 20, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40,
> 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40,
> 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40,
> 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40,
> 40, 40, 40, 40, 40, 40, 40, 40, 40, 40 ;
607,608c607,608
< fates_wood_density = 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7,
< 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, _,
---
> fates_wood_density = 0.9, 0.5, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9,
> 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, _,
689c689
< fates_initd = 0, 0.2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
---
> fates_initd = 0, 0.2, 0.2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
692c692
< 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ;
---
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ;
700,701c700,701
< fates_seed_rain = 0, 0.28, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
< 0, 0, 0, 0, 0, 0, 0, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
---
> fates_seed_rain = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
703c703
< _, _, _, _, _, _, _, _, _, _, _, _, _ ;
---
> _, _, _, _, _, _, _, _, _, _, _, _ ;
705c705
< fates_root_long = 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
---
> fates_root_long = 1, 0.5, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
737c737
< fates_dbh2bl_a = 0.028, 0.028, 0.028, 0.028, _, _, _, _, _, _, _, _, _, _,
---
> fates_dbh2bl_a = 0.021, 0.021, 0.021, 0.021, _, _, _, _, _, _, _, _, _, _,
787c787
< fates_bmort = 0.014, 0.014, 0.014, 0.014, _, _, _, _, _, _, _, _, _, _, _,
---
> fates_bmort = 0.014, 0.07, 0.014, 0.014, _, _, _, _, _, _, _, _, _, _, _, _,
790c790
< _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _ ;
---
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _ ;
854c854
< fates_comp_excln = 0.1 ;
---
> fates_comp_excln = -1.;
rgknox
commented
6 years ago I think I will close this, unless folks have points to still make?
ckoven
commented
6 years ago this problem (how to track default parameter file versions) is now being tracked via the git-controlled cdl files. closing as a thread.
This is the version I am using as of 6 April. It has the SAME VALUES for the first four PFTs which are a first attempt at a tropical-evergreen-broadleaf-tree specific file.
The file was built from fates_params.c170331.nc which was appended to issue #190 . The matlab script that made this file is also in issue #190.
The values are broadly just the values for TrEBrT from the CLM defaults.
This file also includes updated CWD decomposition parameters from @jkshuman that appear to fix issue #190 (to be confirmed once equilibrium is reached) as well as a couple of changes to the fire moisture parameters (script to follow)
tropical_params_decomp_apr06_2017.nc.gz
The code this works with is: $ git log commit bc571f9962db7254ee11830ab05ae82bf60de638 Merge: 8f2c36e 1b8ea3b (unsure which of these numbers is the most relevant. Maybe @bandre-ucar or @rgknox could advise).