Changeset 6319 for branches/ORCHIDEE_2_2/ORCHIDEE/src_stomate/stomate_io.f90

Timestamp:

2019-11-08T13:00:52+01:00 (5 years ago)

Author:

josefine.ghattas

Message:

Improvment for the ESM CO2 configuration:

• Separate variable fco2_lu into 3 parts: fco2_lu, fco2_wh and fco2_ha
• Move calculation of co2_flux from dt_sechiba time-step to daily time-step (in the part for do_slow)
• Removed co2_flux and fco2_lu from stomate_intialize argument list. These variables were never used in the intialization phase.
• Add co2_flux, and fco2_wh, fco2_ha to restart file
• Corrected output unit for nee to be consistent with LMDZ and stomate output variables. It is now in kgC/m2/s.
• Corrected output for znetco2
• Added fCO2_fWoodharvest and fCO2_fHarvest as new possible tracers in LMDZ (intersurf).
• Added diagnostic output for fCO2_fWoodharvest and fCO2_fHarvest

16. Cadule

File:

• branches/ORCHIDEE_2_2/ORCHIDEE/src_stomate/stomate_io.f90 (modified) (6 diffs)

branches/ORCHIDEE_2_2/ORCHIDEE/src_stomate/stomate_io.f90

@@ -80 +80 @@
        &  litterpart, litter, dead_leaves, &
        &  carbon, lignin_struc,turnover_time, &
-       &  fco2_lu, &
+       &  co2_flux, fco2_lu, fco2_wh, fco2_ha, &
        &  prod10,prod100,flux10, flux100, &
        &  convflux, cflux_prod10, cflux_prod100, &
@@ -327 +327 @@
     ! products remaining in the 10/100 year-turnover pool after the annual release for each compartment
     ! (10 or 100 + 1 : input from year of land cover change)
+    REAL(r_std),DIMENSION(npts, nvm),INTENT(out)                           :: co2_flux
     REAL(r_std),DIMENSION(npts),INTENT(out)                                :: fco2_lu
+    REAL(r_std),DIMENSION(npts),INTENT(out)                                :: fco2_wh
+    REAL(r_std),DIMENSION(npts),INTENT(out)                                :: fco2_ha
     REAL(r_std),DIMENSION(npts,0:10),INTENT(out)                           :: prod10
     REAL(r_std),DIMENSION(npts,0:100),INTENT(out)                          :: prod100
@@ -1034 +1037 @@
     !-
     ! Read from restart file or set to zero if the variables or restart file were not found
+
+    var_name = 'co2_flux'
+    CALL restget_p (rest_id_stomate, var_name, nbp_glo, nvm, 1, itime, &
+         .TRUE., co2_flux, 'gather', nbp_glo, index_g)
+    IF (ALL(co2_flux(:,:) == val_exp)) co2_flux(:,:) = zero
+
     var_name = 'fco2_lu'
     CALL restget_p (rest_id_stomate, var_name, nbp_glo, 1, 1, itime, &
          .TRUE., fco2_lu, 'gather', nbp_glo, index_g)
     IF (ALL(fco2_lu(:) == val_exp)) fco2_lu(:) = zero
+
+    var_name = 'fco2_wh'
+    CALL restget_p (rest_id_stomate, var_name, nbp_glo, 1, 1, itime, &
+         .TRUE., fco2_wh, 'gather', nbp_glo, index_g)
+    IF (ALL(fco2_wh(:) == val_exp)) fco2_wh(:) = zero
+
+    var_name = 'fco2_ha'
+    CALL restget_p (rest_id_stomate, var_name, nbp_glo, 1, 1, itime, &
+         .TRUE., fco2_ha, 'gather', nbp_glo, index_g)
+    IF (ALL(fco2_ha(:) == val_exp)) fco2_ha(:) = zero
+
        
     IF (vegetmap_reset) THEN
@@ -1332 +1352 @@
        &  litterpart, litter, dead_leaves, &
        &  carbon, lignin_struc, turnover_time, &
-       &  fco2_lu, &
+       &  co2_flux, fco2_lu, fco2_wh, fco2_ha,  &
        &  prod10,prod100 ,flux10, flux100, &
        &  convflux, cflux_prod10, cflux_prod100, & 
@@ -1570 +1590 @@
     ! products remaining in the 10/100 year-turnover pool after the annual release for each compartment
     ! (10 or 100 + 1 : input from year of land cover change)
+    REAL(r_std),DIMENSION(npts,nvm),INTENT(in)                            :: co2_flux
     REAL(r_std),DIMENSION(npts),INTENT(in)                                :: fco2_lu
+    REAL(r_std),DIMENSION(npts),INTENT(in)                                :: fco2_wh
+    REAL(r_std),DIMENSION(npts),INTENT(in)                                :: fco2_ha
     REAL(r_std),DIMENSION(npts,0:10),INTENT(in)                           :: prod10
     REAL(r_std),DIMENSION(npts,0:100),INTENT(in)                          :: prod100
@@ -2074 +2097 @@
     ! 18 land cover change
     !-
+    var_name = 'co2_flux'
+    CALL restput_p (rest_id_stomate, var_name, nbp_glo, nvm, 1, itime, &
+         &                co2_flux, 'scatter', nbp_glo, index_g)
     var_name = 'fco2_lu'
     CALL restput_p (rest_id_stomate, var_name, nbp_glo, 1, 1, itime, &
          &                fco2_lu, 'scatter', nbp_glo, index_g)
+    var_name = 'fco2_wh'
+    CALL restput_p (rest_id_stomate, var_name, nbp_glo, 1, 1, itime, &
+         &                fco2_wh, 'scatter', nbp_glo, index_g)
+    var_name = 'fco2_ha'
+    CALL restput_p (rest_id_stomate, var_name, nbp_glo, 1, 1, itime, &
+         &                fco2_ha, 'scatter', nbp_glo, index_g)
+
     var_name = 'prod10'
     CALL restput_p (rest_id_stomate, var_name, nbp_glo, 11, 1, itime, &