Changeset 3649 for LMDZ6/trunk/libf/phylmd

Timestamp:

Mar 17, 2020, 12:41:16 PM (5 years ago)

Author:

jghattas

Message:

Correction necessaire pour experience IPSLESM/CO2 pour eviter plantage aleatoire en mode debug : carbon_cpl_init doit etre appele avant premier appel a phys_output_write pour que tout les variables (fco2_land, etc..) soient allouees.

Location:

LMDZ6/trunk/libf/phylmd

Files:

• carbon_cycle_mod.F90 (modified) (1 diff)
• phytrac_mod.F90 (modified) (2 diffs)
• tracco2i_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/carbon_cycle_mod.F90

@@ -201 +201 @@
   
   SUBROUTINE carbon_cycle_init()
-! This subroutine is called from traclmdz_init, only at first timestep.
-! - Read controle parameters from .def input file
-! - Search for carbon tracers and set default values
-! - Allocate variables
-! - Test for compatibility
+    ! This subroutine is called from tracco2i_init, which is called from phytrac_init only at first timestep.
+    ! - Allocate variables. These variables must be allocated before first call to phys_output_write in physiq.
 
     USE dimphy

LMDZ6/trunk/libf/phylmd/phytrac_mod.F90

@@ -56 +56 @@
   SUBROUTINE phytrac_init()
     USE dimphy
-    USE infotrac_phy, ONLY: nbtr
+    USE infotrac_phy, ONLY: nbtr, type_trac
+    USE tracco2i_mod, ONLY: tracco2i_init
     IMPLICIT NONE
 
@@ -70 +71 @@
        ALLOCATE(d_tr_th(klon,klev,nbtr))
        ALLOCATE(d_tr_lessi_impa(klon,klev,nbtr),d_tr_lessi_nucl(klon,klev,nbtr))
+
+
+
+    !===============================================================================
+    !    -- Do specific treatment according to chemestry model or local LMDZ tracers
+    !      
+    !===============================================================================
+    SELECT CASE(type_trac)
+    CASE('co2i')
+       !   -- CO2 interactif --
+       CALL tracco2i_init()
+    END SELECT
+
 
   END SUBROUTINE phytrac_init

LMDZ6/trunk/libf/phylmd/tracco2i_mod.F90

@@ -14 +14 @@
 CONTAINS
 
+  SUBROUTINE tracco2i_init()
+    ! This subroutine calls carbon_cycle_init needed to be done before first call to phys_output_write in physiq.
+    USE carbon_cycle_mod, ONLY: carbon_cycle_init, carbon_cycle_cpl
+
+    ! Initialize carbon_cycle_mod
+    IF (carbon_cycle_cpl) THEN
+       CALL carbon_cycle_init()
+    ENDIF
+
+  END SUBROUTINE tracco2i_init
+
   SUBROUTINE tracco2i(pdtphys, debutphy, &
        xlat, xlon, pphis, pphi, &
@@ -21 +32 @@
     USE infotrac_phy
     USE geometry_mod, ONLY: cell_area
-    USE carbon_cycle_mod, ONLY: carbon_cycle_init
     USE carbon_cycle_mod, ONLY: id_CO2, nbcf_in, fields_in, cfname_in
     USE carbon_cycle_mod, ONLY: fco2_ocn_day, fco2_ff, fco2_bb, fco2_land, fco2_ocean
@@ -79 +89 @@
 !--convert 280 ppm into kg CO2 / kg air
     IF (debutphy) THEN
-
-! Initialisation de module carbon_cycle_mod
-      IF (carbon_cycle_cpl) THEN
-        CALL carbon_cycle_init()
-      ENDIF
 
 ! Initialisation de tr_seri(id_CO2) si pas initialise