Changeset 1565 for LMDZ5/trunk/libf/dyn3d

Timestamp:

Aug 31, 2011, 2:53:29 PM (13 years ago)

Author:

jghattas

Message:

Added interface with chemestry model REPROBUS :

• Compile LMDZ together with Reprobus code (dependecies in both directions) and cpp key REPROBUS :

./makelmdz_fcm -ext_src my_path_to_reprobus -cpp REPROBUS ...

• For running, add type_trac=repr in run.def.

/Marion Marchand, JG

File:

• LMDZ5/trunk/libf/dyn3d/infotrac.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/dyn3d/infotrac.F90

@@ -32 +32 @@
   SUBROUTINE infotrac_init
     USE control_mod
-
+#ifdef REPROBUS
+    USE CHEM_REP, ONLY : Init_chem_rep_trac
+#endif
     IMPLICIT NONE
 !=======================================================================
@@ -94 +96 @@
           CALL abort_gcm('infotrac_init','Incoherence between type_trac and config_inca',1)
        END IF
+    ELSE IF (type_trac=='repr') THEN
+       WRITE(lunout,*) 'You have choosen to couple with REPROBUS chemestry model : type_trac=', type_trac
     ELSE IF (type_trac == 'lmdz') THEN
        WRITE(lunout,*) 'Tracers are treated in LMDZ only : type_trac=', type_trac
@@ -115 +119 @@
 !
 !-----------------------------------------------------------------------
-    IF (type_trac == 'lmdz') THEN
+    IF (type_trac == 'lmdz' .OR. type_trac == 'repr') THEN
        OPEN(90,file='traceur.def',form='formatted',status='old', iostat=ierr)
        IF(ierr.EQ.0) THEN
@@ -147 +151 @@
     END IF
     
+! Transfert number of tracers to Reprobus
+    IF (type_trac == 'repr') THEN
+#ifdef REPROBUS
+       CALL Init_chem_rep_trac(nbtr)
+#endif
+    END IF
        
 !
@@ -182 +192 @@
 !    Get choice of advection schema from file tracer.def or from INCA
 !---------------------------------------------------------------------
-    IF (type_trac == 'lmdz') THEN
+    IF (type_trac == 'lmdz' .OR. type_trac == 'repr') THEN
        IF(ierr.EQ.0) THEN
           ! Continue to read tracer.def