Changes between Version 5 and Version 6 of CleanOrchidee

Timestamp:

2011-10-07T16:36:02+02:00 (13 years ago)

Author:

dsolyga

Comment:

--

CleanOrchidee

@@ -85 +85 @@
     STOP 'Ou est-ce un bug?'
 }}}
- 
+- in '''hydrol.f90''' : drain_rest_cste is declared but not used
+{{{
+  REAL(r_std), PARAMETER                             :: drain_rest_cste = 15.0  !! time constant in days to return to free drainage after return flow
+}}}
+- in '''stomate_lpj.f90''' : the harvest subroutine is bad coded for vectorization 
+
+{{{
+    DO i = 1, npts
+       DO j=1, nvm
+          IF (j > 11) THEN
+             above_old = turnover_daily(i,j,ileaf) + turnover_daily(i,j,isapabove) + &
+                  &       turnover_daily(i,j,iheartabove) + turnover_daily(i,j,ifruit) + &
+                  &       turnover_daily(i,j,icarbres) + turnover_daily(i,j,isapbelow) + &
+                  &       turnover_daily(i,j,iheartbelow) + turnover_daily(i,j,iroot)
+
+             turnover_daily(i,j,ileaf) = turnover_daily(i,j,ileaf)*0.55
+             turnover_daily(i,j,isapabove) = turnover_daily(i,j,isapabove)*0.55
+             turnover_daily(i,j,isapbelow) = turnover_daily(i,j,isapbelow)*0.55
+             turnover_daily(i,j,iheartabove) = turnover_daily(i,j,iheartabove)*0.55
+             turnover_daily(i,j,iheartbelow) = turnover_daily(i,j,iheartbelow)*0.55
+             turnover_daily(i,j,iroot) = turnover_daily(i,j,iroot)*0.55
+             turnover_daily(i,j,ifruit) = turnover_daily(i,j,ifruit)*0.55
+             turnover_daily(i,j,icarbres) = turnover_daily(i,j,icarbres)*0.55
+             harvest_above(i)  = harvest_above(i) + veget_max(i,j) * above_old *0.45
+
+          ENDIF
+       ENDDO
+    ENDDO
+}}}
+