dunno

stineb · Feb 1, 2021 · 9715c2a · 9715c2a
2 parents 6ae9228 + 9701069
commit 9715c2a
Show file tree

Hide file tree

Showing 5 changed files with 496 additions and 75 deletions.
diff --git a/params_std/params_gpp_pmodel.dat b/params_std/params_gpp_pmodel.dat
@@ -7,13 +7,13 @@
 # beta   244.033
 beta                                    146.000000
 rd_to_vcmax                             0.01400000
-soilm_par_a                             0.00000000
-soilm_par_b                             0.78880321
+soilm_par_a                             0.33349283
+soilm_par_b                             1.45602286
 tau_acclim                              10.0000000
 
-kphio_Gr3                               0.08573142
+kphio_Gr3                               0.09423773
 kphio_GN3                               0.08573142
-kphio_Gr4                               0.08573142
+kphio_Gr4                               1.00000000
 kphio_TND                               0.08573142
 kphio_TNE                               0.08573142
 kphio_TrD                               0.08573142

diff --git a/params_std/params_nimpl.dat b/params_std/params_nimpl.dat
@@ -25,14 +25,31 @@ vpd_alnpp                               -1.0192
 intersect_alnpp                         -9.2375
 
 #leaf C/N model
-lma_leafcn                              0.4470
-vcmax25_leafcn                          -0.1068
-intersect_leafcn                        1.5878
+lma_leafcn                              0.0161
+vcmax25_leafcn                          0.0041
+cmass_leafcn                            0.5000
+
+#nre model
+tg_nre                                  -0.0679
+vpd_nre                                 0.4217
+intersect_nre                           1.4541
 
 #constant ratio - derived from TRY database (N=2365 for root; N=13 for wood)
 root_cn                                 42
 wood_cn                                 97
 
+#Grassland models
+#BP/GPP model
+cue_grass                               0.4046
+
+#ANPP/NPP model
+tg_anpp_grass                           0.0385
+alpha_anpp_grass                        1.8720
+intersect_anpp_grass                    -2.5642
+
+#Constant of root C/N
+root_cn_grass                           46
+
 #                                      |
 # write value *behind* this position:  |->
 ## EOF 
diff --git a/src/gpp_pmodel.mod.f90 b/src/gpp_pmodel.mod.f90
@@ -771,7 +771,11 @@ function pmodel( kphio, ppfd, co2, tc, vpd, patm, c4, method_optci, method_jmaxl
       jmax = 0.0
       jmax25 = 0.0
     else
-      fact_jmaxlim = vcmax * (ci + 2.0 * gammastar) / (kphio * ppfd * (ci + kmm))
+      if (kphio==0.0) then
+        fact_jmaxlim = 0.0
+      else
+        fact_jmaxlim = vcmax * (ci + 2.0 * gammastar) / (kphio * ppfd * (ci + kmm))
+      end if
       !print*,'fact_jmaxlim       ', fact_jmaxlim
       if (fact_jmaxlim >= 1 .or. fact_jmaxlim <= 0) then
         jmax = dummy
@@ -1740,6 +1744,11 @@ function calc_ftemp_kphio( dtemp, c4 ) result( ftemp )
 
     if (c4) then
       ftemp = (-0.008 + 0.00375 * dtemp - 0.58e-4 * dtemp**2) * 8.0 ! Based on calibrated values by Shirley
+      if (ftemp < 0.0) then
+        ftemp = 0.0
+      else
+        ftemp = ftemp
+      end if
     else
       ftemp = 0.352 + 0.022 * dtemp - 3.4e-4 * dtemp**2  ! Based on Bernacchi et al., 2003
     end if