Changes between Version 51 and Version 52 of Branches/MergeOCN/Goll


Ignore:
Timestamp:
2015-06-15T13:13:45+02:00 (9 years ago)
Author:
dgoll
Comment:

--

Legend:

Unmodified
Added
Removed
Modified

  • Branches/MergeOCN/Goll

    v51 v52  
    567567             !DSGsimplify 
    568568             DO ielem = 1,nelements 
    569                  ! leaves 
     569                 ! Stand level turnover (gC m-2) for 
     570 
     571                 ! - leaves, ... 
    570572                 dturnover(:) = biomass(:,ivm,ileaf,ielem) * leaf_frac(:,ivm,ilage) * turnover_rate(:) 
    571573                 biomass(:,ivm,ileaf,ielem)  = biomass(:,ivm,ileaf,ielem) - dturnover(:) 
    … …  
    584586                     STOP   
    585587                 END IF 
    586                   
    587                  ! roots   
     588      
     589                 ! - roots, ...   
    588590                 dturnover(:) = biomass(:,ivm,iroot,ielem) * leaf_frac(:,ivm,ilage) * turnover_rate(:) 
    589591                 biomass(:,ivm,iroot,ielem)  = biomass(:,ivm,iroot,ielem) - dturnover(:) 
    … …  
    603605                 END IF 
    604606 
    605                  ! roots   
    606                  dturnover(:) = biomass(:,ivm,iroot,ielem) * leaf_frac(:,ivm,ilage) * turnover_rate(:) 
    607                  biomass(:,ivm,iroot,ielem)  = biomass(:,ivm,iroot,ielem) - dturnover(:) 
    608                   
    609                  IF (ielem == icarbon) THEN ! no recycling 
    610                      turnover(:,ivm,iroot,ielem) = turnover(:,ivm,iroot,ielem) + dturnover(:) 
    611                      ! save root mass change; carbon only 
    612                      delta_lm(:,ilage) = - dturnover(:) 
    613                  ELSEIF (ielem == initrogen) THEN ! recycle 
    614                      biomass(:,ivm,ilabile,ielem) = biomass(:,ivm,ilabile,ielem) + recycle_root(ivm) * dturnover(:) 
    615                      turnover(:,ivm,iroot, ielem) = turnover(:,ivm,iroot, ielem) + ( un - recycle_root(ivm) ) * dturnover(:) 
    616                  ELSEIF (ielem == iphosphorus) THEN ! recycle 
    617                      biomass(:,ivm,ilabile,ielem) = biomass(:,ivm,ilabile,ielem) + p_recycle_root(ivm) * dturnover(:) 
    618                      turnover(:,ivm,iroot, ielem) = turnover(:,ivm,iroot, ielem) + ( un - p_recycle_root(ivm) ) * dturnover(:) 
    619                  ELSE 
    620                      STOP 
    621                  END IF 
    622  
    623                  ! Fruits 
     607                 ! ... and fruits. 
    624608                 dturnover(:) = biomass(:,ivm,ifruit,ielem) * leaf_frac(:,ivm,ilage) * turnover_rate(:) 
    625609                 biomass(:,ivm,ifruit,ielem) = biomass(:,ivm,ifruit,ielem) - dturnover(:)