Changeset 1808 for trunk/LMDZ.TITAN/libf/phytitan

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

File:

• trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90 (modified) (2 diffs)

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'