Changeset 3361 for LMDZ6/trunk/libf/dyn3d_common

Timestamp:

Jul 6, 2018, 9:01:31 AM (6 years ago)

Author:

oboucher

Message:

Embryon of changes for interactive CO2
type_trac=co2i

File:

• LMDZ6/trunk/libf/dyn3d_common/infotrac.F90 (modified) (7 diffs)

LMDZ6/trunk/libf/dyn3d_common/infotrac.F90

@@ -152 +152 @@
        CALL abort_gcm('infotrac_init','You must compile with cpp key REPROBUS',1)
 #endif
+    ELSE IF (type_trac == 'co2i') THEN
+       WRITE(lunout,*) 'You have chosen to run with CO2 cycle: type_trac=', type_trac
     ELSE IF (type_trac == 'coag') THEN
        WRITE(lunout,*) 'Tracers are treated for COAGULATION tests : type_trac=', type_trac
@@ -177 +179 @@
 !
 !-----------------------------------------------------------------------
-    IF (type_trac == 'lmdz' .OR. type_trac == 'repr' .OR. type_trac == 'coag') THEN
+    IF (type_trac == 'lmdz' .OR. type_trac == 'repr' .OR. type_trac == 'coag' .OR. type_trac == 'co2i') THEN
        OPEN(90,file='traceur.def',form='formatted',status='old', iostat=ierr)
        IF(ierr.EQ.0) THEN
@@ -281 +283 @@
 !    Get choice of advection schema from file tracer.def or from INCA
 !---------------------------------------------------------------------
-    IF (type_trac == 'lmdz' .OR. type_trac == 'repr' .OR. type_trac == 'coag') THEN
+    IF (type_trac == 'lmdz' .OR. type_trac == 'repr' .OR. type_trac == 'coag' .OR. type_trac == 'co2i') THEN
        IF(ierr.EQ.0) THEN
           ! Continue to read tracer.def
@@ -673 +675 @@
             if (ipere.gt.nqtot) then
                 WRITE(lunout,*) 'Le traceur',iq,'appele ', &
-      &          trim(tnom_0(iq)),', est orpelin.'
+      &          trim(tnom_0(iq)),', est orphelin.'
                 CALL abort_gcm('infotrac_init','Un traceur est orphelin',1)
             endif !if (ipere.gt.nqtot) then
@@ -684 +686 @@
     WRITE(lunout,*) 'iqpere=',iqpere
     WRITE(lunout,*) 'iqfils=',iqfils
+    WRITE(lunout,*) 'infotrac'
 
 ! Calculer le nombre de descendants à partir de iqfils et de nbfils
@@ -702 +705 @@
         endif !if (ipere.gt.0) then
       enddo !do while (continu)    
+      WRITE(lunout,*) iq
       WRITE(lunout,*) 'Le traceur ',iq,', appele ',trim(tnom_0(iq)),' est un traceur de generation: ',generation
     enddo !DO iq=1,nqtot
@@ -726 +730 @@
     enddo !do iq=1,nqtot
 
+    WRITE(lunout,*) 'infotrac init fin'
 
 ! detecter quels sont les traceurs isotopiques parmi des traceurs