Changeset 3361

Timestamp:

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

Author:

oboucher

Message:

Embryon of changes for interactive CO2
type_trac=co2i

Location:

LMDZ6/trunk/libf

Files:

• dyn3d_common/infotrac.F90 (modified) (7 diffs)
• phylmd/phytrac_mod.F90 (modified) (4 diffs)
• phylmd/tracco2i_mod.F90 (added)

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

LMDZ6/trunk/libf/phylmd/phytrac_mod.F90

@@ -100 +100 @@
     USE print_control_mod, ONLY: lunout
     USE aero_mod, ONLY : naero_grp
+
+    USE tracco2i_mod
 
 #ifdef CPP_StratAer
@@ -474 +476 @@
        CASE('repr')
           source(:,:)=0.
+       CASE('co2i')
+          source(:,:)=0.
 #ifdef CPP_StratAer
        CASE('coag')
@@ -537 +541 @@
                 flag_cvltr(it)=.false.
 
+             CASE('co2i')
+                !--co2 tracers are not scavenged
+                flag_cvltr(it)=.false.
+
 #ifdef CPP_StratAer
              CASE('coag')
@@ -624 +632 @@
             t_seri, pplay, paprs, sh , &
             tr_seri)
+
+    CASE('co2i')
+       !   -- CO2 interactif --
+       !   -- source is updated with FF and BB emissions 
+       !   -- OB => PC need to add net flux from ocean and orchidee 
+       !   -- sign convention : positive into the atmosphere
+       CALL tracco2i(pdtphys, debutphy, &
+            xlat, xlon, pphis, pphi, &
+            t_seri, pplay, paprs, tr_seri, source)
 
 #ifdef CPP_StratAer