Changeset 7326 for branches/ORCHIDEE_2_2/ORCHIDEE/src_stomate/stomate_lpj.f90

Timestamp:

2021-10-20T18:39:22+02:00 (3 years ago)

Author:

josefine.ghattas

Message:

Corrected bug on carbon balance closure. See ticket #785
Integration in branch 2_2 done by P. Cadule

File:

• branches/ORCHIDEE_2_2/ORCHIDEE/src_stomate/stomate_lpj.f90 (modified) (13 diffs)

branches/ORCHIDEE_2_2/ORCHIDEE/src_stomate/stomate_lpj.f90

@@ -875 +875 @@
     sum_cVegTree = zero
 
-    DO j=2,nvm
+    DO j=1,nvm
 
        tot_litter_carb(:,j) = tot_litter_carb(:,j) + (litter(:,istructural,j,iabove,icarbon) + &
@@ -958 +958 @@
     CALL xios_orchidee_send_field("LAI",lai)
     CALL xios_orchidee_send_field("VEGET_COV_MAX",veget_cov_max)
-    CALL xios_orchidee_send_field("NPP_STOMATE",npp_daily)
-    CALL xios_orchidee_send_field("GPP",gpp_daily)
+    CALL xios_orchidee_send_field("NPP_STOMATE",npp_daily+co2_to_bm)
+    CALL xios_orchidee_send_field("GPP",gpp_daily+co2_to_bm)
     CALL xios_orchidee_send_field("IND",ind)
     CALL xios_orchidee_send_field("CN_IND",cn_ind)
@@ -1053 +1053 @@
     
     ! Carbon fluxes transformed from gC/m2/d into kgC/m2/s
-    CALL xios_orchidee_send_field("gpp_ipcc",SUM(gpp_daily*veget_cov_max,dim=2)/1e3/one_day)
+    CALL xios_orchidee_send_field("gpp_ipcc",SUM((gpp_daily+co2_to_bm)*veget_cov_max,dim=2)/1e3/one_day)
     CALL xios_orchidee_send_field("ra",SUM((resp_maint+resp_growth)*veget_cov_max,dim=2)/1e3/one_day)
     vartmp(:)=zero
@@ -1079 +1079 @@
     CALL xios_orchidee_send_field("raTree",vartmp/1e3/one_day)
 
-    CALL xios_orchidee_send_field("npp_ipcc",SUM(npp_daily*veget_cov_max,dim=2)/1e3/one_day)
+    CALL xios_orchidee_send_field("npp_ipcc",SUM((npp_daily+co2_to_bm)*veget_cov_max,dim=2)/1e3/one_day)
+
     vartmp(:)=zero
     DO j = 2, nvm
        IF ( .NOT. is_tree(j) .AND. natural(j) ) THEN
-          vartmp(:) = vartmp(:) + npp_daily(:,j)*veget_cov_max(:,j)
+          vartmp(:) = vartmp(:) +( npp_daily(:,j)+co2_to_bm(:,j))*veget_cov_max(:,j)
        ENDIF
     ENDDO
@@ -1090 +1091 @@
     DO j = 2, nvm
        IF ( (.NOT. is_tree(j)) .AND. (.NOT. natural(j)) ) THEN
-          vartmp(:) = vartmp(:) + npp_daily(:,j)*veget_cov_max(:,j)
+          vartmp(:) = vartmp(:) + (npp_daily(:,j)+co2_to_bm(:,j))*veget_cov_max(:,j)
        ENDIF
     ENDDO
@@ -1098 +1099 @@
     DO j = 2, nvm
        IF ( is_tree(j) ) THEN
-          vartmp(:) = vartmp(:) + npp_daily(:,j)*veget_cov_max(:,j)
+          vartmp(:) = vartmp(:) + (npp_daily(:,j)+co2_to_bm(:,j))*veget_cov_max(:,j)
        ENDIF
     ENDDO
@@ -1196 +1197 @@
     CALL xios_orchidee_send_field("flulccatmlut",flulccatmlut)
 
-   ! co2_to_bm is not added as it is already included in gpp
-    CALL xios_orchidee_send_field("nbp",(SUM((gpp_daily-(resp_maint+resp_growth+resp_hetero)-co2_fire) * &
+    CALL xios_orchidee_send_field("nbp",(SUM((gpp_daily+co2_to_bm-(resp_maint+resp_growth+resp_hetero)-co2_fire) * &
           veget_cov_max,dim=2)-cflux_prod_total-cflux_prod_harvest_total-harvest_above)/1e3/one_day)
     CALL xios_orchidee_send_field("fVegLitter",SUM((tot_bm_to_litter(:,:,icarbon) + tot_turnover(:,:,icarbon))*&
@@ -1267 +1267 @@
                resp_hetero(:,j)*veget_cov_max(:,j)/1e3/one_day
           npplut(:,id_psl) = npplut(:,id_psl) + &
-               npp_daily(:,j)*veget_cov_max(:,j)/1e3/one_day
+               (npp_daily(:,j)+co2_to_bm(:,j))*veget_cov_max(:,j)/1e3/one_day
        ELSE
           clitterlut(:,id_crp) = clitterlut(:,id_crp) + tot_litter_carb(:,j)*veget_cov_max(:,j)/1e3
@@ -1278 +1278 @@
                resp_hetero(:,j)*veget_cov_max(:,j)/1e3/one_day
           npplut(:,id_crp) = npplut(:,id_crp) + &
-               npp_daily(:,j)*veget_cov_max(:,j)/1e3/one_day
+               (npp_daily(:,j)+co2_to_bm(:,j))*veget_cov_max(:,j)/1e3/one_day
        END IF
     END DO
@@ -1422 +1422 @@
          veget_cov_max, npts*nvm, horipft_index)
     CALL histwrite_p (hist_id_stomate, 'NPP', itime, &
-         npp_daily, npts*nvm, horipft_index)
+         npp_daily+co2_to_bm, npts*nvm, horipft_index)
     CALL histwrite_p (hist_id_stomate, 'GPP', itime, &
-         gpp_daily, npts*nvm, horipft_index)
+         gpp_daily+co2_to_bm, npts*nvm, horipft_index)
     CALL histwrite_p (hist_id_stomate, 'IND', itime, &
          ind, npts*nvm, horipft_index)
@@ -1528 +1528 @@
        CALL histwrite_p (hist_id_stomate_IPCC, "lai", itime, &
             vartmp, npts, hori_index)
-       vartmp(:)=SUM(gpp_daily*veget_cov_max,dim=2)/1e3/one_day
+       vartmp(:)=SUM((gpp_daily+co2_to_bm)*veget_cov_max,dim=2)/1e3/one_day
        CALL histwrite_p (hist_id_stomate_IPCC, "gpp", itime, &
             vartmp, npts, hori_index)
@@ -1534 +1534 @@
        CALL histwrite_p (hist_id_stomate_IPCC, "ra", itime, &
             vartmp, npts, hori_index)
-       vartmp(:)=SUM(npp_daily*veget_cov_max,dim=2)/1e3/one_day
+       vartmp(:)=SUM((npp_daily+co2_to_bm)*veget_cov_max,dim=2)/1e3/one_day
        CALL histwrite_p (hist_id_stomate_IPCC, "npp", itime, &
             vartmp, npts, hori_index)
@@ -1553 +1553 @@
             vartmp, npts, hori_index)
        ! co2_to_bm is not added as it is already included in gpp
-       vartmp(:)=(SUM((gpp_daily-(resp_maint+resp_growth+resp_hetero)-co2_fire) &
+       vartmp(:)=(SUM((gpp_daily+co2_to_bm-(resp_maint+resp_growth+resp_hetero)-co2_fire) &
             &        *veget_cov_max,dim=2)-cflux_prod_total-cflux_prod_harvest_total-harvest_above)/1e3/one_day
        CALL histwrite_p (hist_id_stomate_IPCC, "nbp", itime, &