Changeset 1227 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/traclmdz_mod.F90

Timestamp:

Aug 10, 2009, 1:45:04 PM (15 years ago)

Author:

jghattas

Message:

• Inclusion d'un premier version du cycle de carbon dans LMDZ. Attention

!! Il s'agit d'un version ou les nouveaux cles cycle_carbon_tr et
cycle_carbon_cpl ne sont pas teste. Avec les ancinenes parametres le
modele donne les memes resultats qu'avant. L'interface avec ORCHIDEE n'a
pas encore etait modifie.

• physiq.F, phys_cal_mod.F90 : ajout d'un nouveau module qui contient qq parametres pour le calendrier et le pas de temps acutelle de la physiq. Ce module pourrait etre elargie plus tard / LF + JG

• infotrac.F90 : les noms du traceurs peut prendre un nom plus long (15 caracteres) dans traceur.def

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/traclmdz_mod.F90 (modified) (7 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/traclmdz_mod.F90

@@ -65 +65 @@
 
 
-  SUBROUTINE traclmdz_init(pctsrf, ftsol, aerosol, lessivage)
+  SUBROUTINE traclmdz_init(pctsrf, ftsol, tr_seri, aerosol, lessivage)
     ! This subroutine allocates and initialize module variables and control variables.
     USE dimphy
     USE infotrac
+    USE carbon_cycle_mod, ONLY : carbon_cycle_init, carbon_cycle_tr, carbon_cycle_cpl
 
     IMPLICIT NONE
@@ -75 +76 @@
 
 ! Input variables
-    REAL,DIMENSION(klon,nbsrf),INTENT(IN) :: pctsrf ! Pourcentage de sol f(nature du sol)
-    REAL,DIMENSION(klon,nbsrf),INTENT(IN) :: ftsol  ! Temperature du sol (surf)(Kelvin)
-    
+    REAL,DIMENSION(klon,nbsrf),INTENT(IN)     :: pctsrf ! Pourcentage de sol f(nature du sol)
+    REAL,DIMENSION(klon,nbsrf),INTENT(IN)     :: ftsol  ! Temperature du sol (surf)(Kelvin)
+    REAL,DIMENSION(klon,klev,nbtr),INTENT(IN) :: tr_seri! Concentration Traceur [U/KgA]  
+
 ! Output variables
     LOGICAL,DIMENSION(nbtr), INTENT(OUT) :: aerosol
@@ -150 +152 @@
     END IF
 
+!
+! Initialisation de module carbon_cycle_mod
+! ----------------------------------------------
+    IF (carbon_cycle_tr .OR. carbon_cycle_cpl) THEN
+       CALL carbon_cycle_init(tr_seri, aerosol, radio)
+    END IF
+
   END SUBROUTINE traclmdz_init
 
@@ -161 +170 @@
     USE dimphy
     USE infotrac
+    USE carbon_cycle_mod, ONLY : carbon_cycle, carbon_cycle_tr, carbon_cycle_cpl
     
     IMPLICIT NONE
@@ -199 +209 @@
     REAL,DIMENSION(klon,nbsrf),INTENT(IN) :: pctsrf ! Pourcentage de sol f(nature du sol)
 
+
 ! InOutput argument
     REAL,DIMENSION(klon,klev,nbtr),INTENT(INOUT)  :: tr_seri ! Concentration Traceur [U/KgA]  
@@ -206 +217 @@
     REAL,DIMENSION(klon,nbtr), INTENT(OUT)        :: source  ! a voir lorsque le flux de surface est prescrit 
     REAL,DIMENSION(klon,klev,nbtr), INTENT(OUT)   :: d_tr_cl ! Td couche limite/traceur
-
 
 !=======================================================================================
@@ -301 +311 @@
     END DO
 
+!======================================================================
+!   Calcul de cycle de carbon
+!======================================================================
+    IF (carbon_cycle_tr .OR. carbon_cycle_cpl) THEN
+       CALL carbon_cycle(nstep, pdtphys, pctsrf, tr_seri)
+    END IF
+
   END SUBROUTINE traclmdz