Changeset 257 for branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/lpj_fire.f90

Timestamp:

2011-06-17T14:02:17+02:00 (13 years ago)

Author:

didier.solyga

Message:

Externalized version merged with the trunk

File:

• branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/lpj_fire.f90 (modified) (13 diffs)

branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/lpj_fire.f90

@@ -99 +99 @@
     !MM Shilong ??
 !!$    REAL(r_std), PARAMETER                                           :: tau_fire = 365.  ! GKtest
-
     ! fire perturbation
     REAL(r_std), DIMENSION(npts)                                     :: fire_disturb
@@ -273 +272 @@
        IF(.NOT.disable_fire.AND.natural(j))THEN
           WHERE ( aff(:) .GT. 0.1 )
-             firefrac(:,j) = 1. - ( 1. - aff(:) ) ** (dt/one_year)
+             firefrac(:,j) = un - ( un - aff(:) ) ** (dt/one_year)
           ELSEWHERE
              firefrac(:,j) = aff(:) * dt/one_year
@@ -315 +314 @@
           ! 4.2.1 Trees: always disturbed
 
-          fire_disturb(:) = ( 1. - resist(j) ) * firefrac(:,j)
+          fire_disturb(:) = ( un - resist(j) ) * firefrac(:,j)
 
        ELSE
@@ -323 +322 @@
           WHERE ( biomass(:,j,ileaf) .GT. min_stomate )
 
-             fire_disturb(:) = ( 1. - resist(j) ) * firefrac(:,j)
+             fire_disturb(:) = ( un - resist(j) ) * firefrac(:,j)
 
           ELSEWHERE
@@ -353 +352 @@
              ! 4.3.2 Determine the residue, in gC/m**2 of ground.
 
-             residue(:) = biomass(:,j,k) * fire_disturb(:) * ( 1. - co2frac(k) )
-             !MM in SZ ???        residue(:) = fire_disturb(:) * ( 1. - co2frac(k) )
+             residue(:) = biomass(:,j,k) * fire_disturb(:) * ( un - co2frac(k) )
+             !MM in SZ ???        residue(:) = fire_disturb(:) * ( un - co2frac(k) )
 
              ! 4.3.2.1 determine fraction of black carbon. Only for plant parts above the
@@ -400 +399 @@
           IF ( .NOT. ( ( .NOT. tree(j) ) .AND. ( ( k.EQ.iroot ) .OR. ( k.EQ.icarbres) ) ) ) THEN
 
-             biomass(:,j,k) = ( 1. - fire_disturb(:) ) * biomass(:,j,k)
+             biomass(:,j,k) = ( un - fire_disturb(:) ) * biomass(:,j,k)
 
           ENDIF
@@ -409 +408 @@
        !       individuals.
 
-       IF ( control%ok_dgvm .AND. tree(j) ) THEN
+       IF ( (control%ok_dgvm .OR. .NOT.lpj_gap_const_mort) .AND. tree(j) ) THEN
 
           ! fraction of plants that dies each day.
@@ -415 +414 @@
           firedeath(:,j) = fire_disturb(:) / dt
 
-          ind(:,j) = ( 1. - fire_disturb(:) ) * ind(:,j)
+          ind(:,j) = ( un - fire_disturb(:) ) * ind(:,j)
 
        ENDIF
@@ -440 +439 @@
 
        litter(:,imetabolic,j,iabove) = litter(:,imetabolic,j,iabove) * &
-            ( 1. - firefrac(:,j) )
+            ( un - firefrac(:,j) )
 
        !
@@ -455 +454 @@
        co2_fire(:,j) = co2_fire(:,j) + &
             litter(:,istructural,j,iabove) * firefrac(:,j) * &
-            ( 1. - struc_residual(:) )/ dt
+            ( un - struc_residual(:) )/ dt
 
        ! 5.2.3 determine residue (litter that undergoes fire, but is not transformed
        !       into CO2)
 
-       residue(:) = litter(:,istructural,j,iabove) * firefrac(:,j) * &
-            struc_residual(:)
+!NV,MM : We add this test to keep coherence with CMIP5 computations without DGVM.
+!        It has to be removed in trunk version after CMIP5.
+       IF (control%ok_dgvm .OR. .NOT.lpj_gap_const_mort) THEN
+          residue(:) = firefrac(:,j) * struc_residual(:)
+       ELSE
+          residue(:) = litter(:,istructural,j,iabove) * firefrac(:,j) * &
+               struc_residual(:)
+       ENDIF
+
+!       residue(:) = litter(:,istructural,j,iabove) * firefrac(:,j) * &
+!            struc_residual(:)
        !MM in SZ        residue(:) = firefrac(:,j) * struc_residual(:)
 
@@ -482 +490 @@
 
        litter(:,istructural,j,iabove) = &
-            litter(:,istructural,j,iabove) * ( 1. - firefrac(:,j) ) + &
-            residue(:) * ( 1. - bcfrac(:) )
-       !MM in SZ            residue(:) * ( 1. - bcfrac(:) ) * litter(:,iwoody,j,iabove)
+            litter(:,istructural,j,iabove) * ( un - firefrac(:,j) ) + &
+            residue(:) * ( un - bcfrac(:) )
+       !MM in SZ            residue(:) * ( un - bcfrac(:) ) * litter(:,iwoody,j,iabove)
 
     ENDDO  !  ground
@@ -496 +504 @@
 
        DO k = 1, nlitt
-          dead_leaves(:,j,k) = dead_leaves(:,j,k) * ( 1. - firefrac(:,j) )
+          dead_leaves(:,j,k) = dead_leaves(:,j,k) * ( un - firefrac(:,j) )
        ENDDO
 
@@ -543 +551 @@
 
     firefrac_result(:) = &
-!         x(:) * EXP( xm1(:) / ( -.13*xm1(:)*xm1(:)*xm1(:) + .6*xm1(:)*xm1(:) + .8*xm1(:) + .45 ) )
          x(:) * EXP( xm1(:) / ( -firefrac_coeff(4)*xm1(:)*xm1(:)*xm1(:) + firefrac_coeff(3)*xm1(:)*xm1(:) + firefrac_coeff(2)*xm1(:) + firefrac_coeff(1) ) )
 
-
   END FUNCTION firefrac_func