Changeset 3549 for LMDZ6/trunk/libf/phylmd

Timestamp:

Aug 5, 2019, 9:40:19 PM (5 years ago)

Author:

oboucher

Message:

Allocate CO2 flux variables at first timestep

File:

• LMDZ6/trunk/libf/phylmd/tracco2i_mod.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/tracco2i_mod.F90

@@ -65 +65 @@
 !--convert 280 ppm into kg CO2 / kg air
     IF (debutphy) THEN
+
       IF (MAXVAL(tr_seri(:,:,id_CO2)).LT.1.e-15) THEN 
         !!tr_seri(:,:,id_CO2)=280.e-6/RMD*RMCO2
         tr_seri(:,:,id_CO2)=400.e-6/RMD*RMCO2 !--initialised to 400 ppm for a test
       ENDIF 
+
+      ALLOCATE(fco2_ff(klon))
+      ALLOCATE(fco2_bb(klon))
+      ALLOCATE(fco2_land(klon))
+      ALLOCATE(fco2_ocean(klon))
+
     ENDIF
 
@@ -87 +94 @@
     fco2_ocean(:)=0.0
     DO nb=1, nbcf_in
+      print *,'nb tracco2=', nb, cfname_in(nb)
 !--fCO2_nep comes in unit of kg C m-2 s-1
 !--converting to kg CO2 m-2 s-1
-      IF (cfname_in(nb) == "fCO2_nep" )   fco2_land(:)=fields_in(:,nb)*RMCO2/RMC*pctsrf(:,is_ter)
+      IF (cfname_in(nb) == "fCO2_nbp" )   fco2_land(:)=fields_in(:,nb)*RMCO2/RMC*pctsrf(:,is_ter)
 !--fCO2_fgco2 comes in unit of mol C02 m-2 s-1
 !--converting to kg CO2 m-2 s-1 + change sign
@@ -106 +114 @@
       CALL gather(m_air,m_air_glo)
 !$OMP MASTER
+
 !--compute a global mean CO2 value and print its value in ppm
        IF (is_mpi_root) THEN
@@ -138 +147 @@
     USE netcdf, ONLY: nf90_get_var, nf90_noerr, nf90_nowrite
 
-    USE carbon_cycle_mod, ONLY : fco2_ff, fco2_bb
+    USE carbon_cycle_mod, ONLY : fco2_ff, fco2_bb, fco2_land, fco2_ocean
 
     IMPLICIT NONE
@@ -259 +268 @@
     PRINT *,'probleme avec le mois dans co2_ini =', mth_cur
   ENDIF
-  IF (.NOT.ALLOCATED(fco2_ff)) ALLOCATE(fco2_ff(klon))
-  IF (.NOT.ALLOCATED(fco2_bb)) ALLOCATE(fco2_bb(klon))
+
   fco2_ff(:) = flx_co2ff(:,mth_cur)
   fco2_bb(:) = flx_co2bb(:,mth_cur)