Changes between Version 2 and Version 3 of DevelopmentActivities/Bugs

Timestamp:

2015-12-02T14:38:51+01:00 (9 years ago)

Author:

ajornet

Comment:

--

DevelopmentActivities/Bugs

@@ -3 +3 @@
 List of bugs that need to be fixed. They are order by priority from the top to the bottom.
 
-== List ==
 
-=== PFT maps and Analytic spinup with DGVM ===
+== PFT maps and Analytic spinup with DGVM ==
+
+No trees with DGVM (without spitfire)
 
 You have a list in pft_parameters:
@@ -26 +27 @@
 This is independent from Chao's problem. I'm looking at DGVM without SPITFIRE. In dynamic mode, we have first crops but natural PFTs appear later: grasses and after a few years trees. At least that's how it behaves with the TRUNK version.
 
-=== Spitfire ===
+== Spitfire ==
 
-* There are regative values at fuel_XXhr variables
-* Multiple floating point exceptions at reat_of_spread
+=== There are regative values at fuel_XXhr variables ===
 
-=== VSD merge ===
+
+
+=== Multiple floating point exceptions at reat_of_spread ===
+
+
+There is a floating overflow exception at exp(XXXX) from lpj_spitfire.f90 -> rate_of_spread
+
+{{{
+! reaction intensity
+a(:)=8.9033*(sigma(:)**(-0.7913))
+WHERE (a(:).le.0.00001 .OR. bet(:).le.0.00001)
+   dummy(:)=0.0
+ELSEWHERE
+   dummy(:)=exp(a(:)*(1.0-bet(:)))
+ENDWHERE
+}}}
+
+Proposed solution: set sigma values to 10000. FAILED
+{{{
+     !calculate mass weighted surface-area-to-volume ratio by fuel types.
+     WHERE (dead_fuel(:).gt.min_stomate)
+       sigma(:)=(fuel_1hr_total(:) * sigma_1hr + &
+            fuel_10hr_total(:) * sigma_10hr + &
+            fuel_100hr_total(:) * sigma_100hr) / dead_fuel(:)
+     ELSEWHERE
+       sigma(:)=0.00001 -> 10000
+     ENDWHERE
+}}}
+
+
+Proposed solution: check sigma instead of a. OK (Accepted by Chao and Fabienne)
+
+{{{
+! reaction intensity
+a(:)=8.9033*(sigma(:)**(-0.7913))
+WHERE (sigma(:).le.0.00001 .OR. bet(:).le.0.00001)
+   dummy(:)=0.0
+ELSEWHERE
+   dummy(:)=exp(a(:)*(1.0-bet(:)))
+ENDWHERE
+}}}
+
+
+
+== VSD merge ==
 
 How to deal with already designed MICT soil layers ?