Changeset 1808

Timestamp:

Oct 23, 2017, 10:08:56 PM (7 years ago)

Author:

jvatant

Message:

Minor modifs :
+ correct a bug for ifort compiling in muphy
+ added some outputs for chemistry
JVO

Location:

trunk/LMDZ.TITAN/libf

Files:

• muphytitan/strings.F90 (modified) (1 diff)
• phytitan/physiq_mod.F90 (modified) (2 diffs)

trunk/LMDZ.TITAN/libf/muphytitan/strings.F90

@@ -854 +854 @@
     !! If the size of the output string vector (i.e. the character length of the string elements within the
     !! vector) is too small, words can be truncated.
-    OBJECT(words), INTENT(in)                                :: this
+    OBJECT(words), INTENT(in)                                      :: this
       !! A words object reference
-    CHARACTER(len=*), INTENT(out), ALLOCATABLE, DIMENSION(:) :: ret
+    CHARACTER(len=st_slen), INTENT(out), ALLOCATABLE, DIMENSION(:) :: ret
       !! An allocatable vector of assumed length string with the words of __this__
     LOGICAL             :: ok

trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90

@@ -1460 +1460 @@
       ! Output tracers.
       if (tracer) then
-      
+
          if (callmufi) then ! For now we assume an given order for tracers !
             call writediagfi(ngrid,"mu_m0as","Spherical mode 0th order moment",'kg/kg',3,zq(:,:,1))
@@ -1468 +1468 @@
          endif ! end of 'callmufi'
 
+        if ( callchim .and. (.not. callmufi) ) then
+           call writediagfi(ngrid,"C2H2","C2H2",'kg/kg',3,zq(:,:,10))
+           call writediagfi(ngrid,"C2H4","C2H4",'kg/kg',3,zq(:,:,12))
+           call writediagfi(ngrid,"C2H6","C2H6",'kg/kg',3,zq(:,:,14))
+           call writediagfi(ngrid,"C4H2","C4H2",'kg/kg',3,zq(:,:,26))
+           call writediagfi(ngrid,"HCN","HCN",'kg/kg',3,zq(:,:,36))
+           call writediagfi(ngrid,"HC3N","HC3N",'kg/kg',3,zq(:,:,42))
+        else if ( callchim .and. callmufi ) then
+           call writediagfi(ngrid,"C2H2","C2H2",'kg/kg',3,zq(:,:,6))
+           call writediagfi(ngrid,"C2H4","C2H4",'kg/kg',3,zq(:,:,8))
+           call writediagfi(ngrid,"C2H6","C2H6",'kg/kg',3,zq(:,:,10))
+           call writediagfi(ngrid,"C4H2","C4H2",'kg/kg',3,zq(:,:,22))
+           call writediagfi(ngrid,"HCN","HCN",'kg/kg',3,zq(:,:,32))
+           call writediagfi(ngrid,"HC3N","HC3N",'kg/kg',3,zq(:,:,38))
+        endif
+
        endif ! end of 'tracer'