Changeset 1405 for LMDZ4/branches/LMDZ4_AR5/libf/phylmd/phys_output_mod.F90

Timestamp:

Jul 1, 2010, 3:16:29 PM (14 years ago)

Author:

jghattas

Message:

• Added tracer "Age of stratospheric air", activated by putting Aga in tracer.def.
• The logical rnpb is now controled during execution. rnpb is true only if both tracers RN and PB existe in tracer.def. RN and PB can now be removed from tracer.def
• Modification in output. It is now possible to have 3 tracers. But still only tracer number 3 and tracer number 4(if existing) will be written.

F Lott, JG

File:

• LMDZ4/branches/LMDZ4_AR5/libf/phylmd/phys_output_mod.F90 (modified) (1 diff)

LMDZ4/branches/LMDZ4_AR5/libf/phylmd/phys_output_mod.F90

@@ -1224 +1224 @@
      ENDIF
 
-      if (nqtot>=3) THEN
-!Attention    DO iq=3,nqtot
-    DO iq=3,4  
-       iiq=niadv(iq)
-! CALL histdef3d (iff, o_trac%flag,'o_'//tnom(iq)%name,ttext(iiq), "-" )
-  CALL histdef3d (iff, o_trac(iq-2)%flag,o_trac(iq-2)%name,ttext(iiq), "-" )
-    ENDDO
-      endif
+
+!Attention : sorties uniquement pour traceurs 3 et 4
+     IF (nqtot>=3) THEN
+        IF (tname(3) == "Aga" .OR. tname(3)=="AGA") THEN
+           o_trac(1)%name= "Aga"
+           ttext(3)= "Age stratospheric air"
+        END IF
+        CALL histdef3d (iff, o_trac(1)%flag,o_trac(1)%name,ttext(3), "-" )
+     END IF
+     
+     IF (nqtot>=4) CALL histdef3d (iff, o_trac(2)%flag,o_trac(2)%name,ttext(4), "-" )
+
 
         CALL histend(nid_files(iff))