Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM	(revision 2020)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM	(revision 2021)
@@ -147,5 +147,38 @@
 ####MARS MARS MARS MARS MARS
 ####MARS MARS MARS MARS MARS
-
+###VENUS TRACER
+state  real  qco2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCO2"       "CO2 number mixing ratio"   "kg kg-1"
+state  real  qco       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCO"        "CO number mixing ratio"   "kg kg-1"
+state  real  qh2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2"        "H2 number mixing ratio"   "kg kg-1"
+state  real  qo1d      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO1D"        "O1d number mixing ratio"   "kg kg-1"
+state  real  qo        ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO"        "O number mixing ratio"   "kg kg-1"
+state  real  qo2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO2"        "O2 number mixing ratio"   "kg kg-1"
+state  real  qo2dg     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO2dg"        "O2dg number mixing ratio"   "kg kg-1"
+state  real  qo3       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO3"        "O3 number mixing ratio"   "kg kg-1"
+state  real  qh        ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH"        "H number mixing ratio"   "kg kg-1"
+state  real  qoh       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOH"        "OH number mixing ratio"   "kg kg-1"
+state  real  qho2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHO2"        "HO2 number mixing ratio"   "kg kg-1"
+state  real  qh2o2     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2O2"        "H2O2 number mixing ratio"   "kg kg-1"
+state  real  qcl       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCL"        "Cl number mixing ratio"   "kg kg-1"
+state  real  qclo      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLO"        "ClO number mixing ratio"   "kg kg-1"
+state  real  qcl2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCL2"        "Cl2 number mixing ratio"   "kg kg-1"
+state  real  qhcl      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHCL"        "HCl number mixing ratio"   "kg kg-1"
+state  real  qhocl     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHOCL"        "HOCl number mixing ratio"   "kg kg-1"
+state  real  qclco     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLCO"        "ClCO number mixing ratio"   "kg kg-1"
+state  real  qclco3    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLCO3"        "ClCO3 number mixing ratio"   "kg kg-1"
+state  real  qcocl2    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCOCL2"        "COCl2 number mixing ratio"   "kg kg-1"
+state  real  qss       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS"        "S number mixing ratio"   "kg kg-1"
+state  real  qso       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO"        "SO number mixing ratio"   "kg kg-1"
+state  real  qso2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO2"        "SO2 number mixing ratio"   "kg kg-1"
+state  real  qso3      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO3"        "SO3 number mixing ratio"   "kg kg-1"
+state  real  qs2o2     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS2O2"        "S2O2 number mixing ratio"   "kg kg-1"
+state  real  qocs      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOCS"        "OCS number mixing ratio"   "kg kg-1"
+state  real  qhso3     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHSO3"        "HSO3 number mixing ratio"   "kg kg-1"
+state  real  qh2so4    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2SO4"        "H2SO4 number mixing ratio"   "kg kg-1"
+state  real  qs2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS2"        "S2 number mixing ratio"   "kg kg-1"
+state  real  qclso2    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLSO2"        "ClSO2 number mixing ratio"   "kg kg-1"
+state  real  qoscl     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOSCL"        "OSCl number mixing ratio"   "kg kg-1"
+state  real  qh2oliq   ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2OLIQ"        "H2O liq number mixing ratio"   "kg kg-1"
+state  real  qh2so4liq ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2SO4LIQ"        "H2SO4 liq number mixing ratio"   "kg kg-1"
 
 
@@ -1517,4 +1550,5 @@
 package   radioac      mars==20                     -              scalar:qtrac1
 package   radioac2     mars==21                     -              scalar:upward,downward
+package   VenusPC      mars==34                     -              scalar:qco2,qco,qh2,qh2o,qo1d,qo,qo2,qo2dg,qo3,qh,qoh,qho2,qh2o2,qcl,qclo,qcl2,qhcl,qhocl,qclco,qclco3,qcocl2,qss,qso,qso2,qso3,qs2o2,qocs,qhso3,qh2so4,qs2,qclso2,qoscl,qh2oliq,qh2so4liq
 ##### MARS OPTIONS
 ##### MARS OPTIONS
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM.newphys
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM.newphys	(revision 2020)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM.newphys	(revision 2021)
@@ -144,4 +144,39 @@
 ####MARS MARS MARS MARS MARS
 ####MARS MARS MARS MARS MARS
+###VENUS TRACER
+state  real  qco2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCO2"       "CO2 number mixing ratio"   "kg kg-1"
+state  real  qco       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCO"        "CO number mixing ratio"   "kg kg-1"
+state  real  qh2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2"        "H2 number mixing ratio"   "kg kg-1"
+state  real  qo1d      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO1D"        "O1d number mixing ratio"   "kg kg-1"
+state  real  qo        ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO"        "O number mixing ratio"   "kg kg-1"
+state  real  qo2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO2"        "O2 number mixing ratio"   "kg kg-1"
+state  real  qo2dg     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO2dg"        "O2dg number mixing ratio"   "kg kg-1"
+state  real  qo3       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO3"        "O3 number mixing ratio"   "kg kg-1"
+state  real  qh        ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH"        "H number mixing ratio"   "kg kg-1"
+state  real  qoh       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOH"        "OH number mixing ratio"   "kg kg-1"
+state  real  qho2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHO2"        "HO2 number mixing ratio"   "kg kg-1"
+state  real  qh2o2     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2O2"        "H2O2 number mixing ratio"   "kg kg-1"
+state  real  qcl       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCL"        "Cl number mixing ratio"   "kg kg-1"
+state  real  qclo      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLO"        "ClO number mixing ratio"   "kg kg-1"
+state  real  qcl2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCL2"        "Cl2 number mixing ratio"   "kg kg-1"
+state  real  qhcl      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHCL"        "HCl number mixing ratio"   "kg kg-1"
+state  real  qhocl     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHOCL"        "HOCl number mixing ratio"   "kg kg-1"
+state  real  qclco     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLCO"        "ClCO number mixing ratio"   "kg kg-1"
+state  real  qclco3    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLCO3"        "ClCO3 number mixing ratio"   "kg kg-1"
+state  real  qcocl2    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCOCL2"        "COCl2 number mixing ratio"   "kg kg-1"
+state  real  qss       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS"        "S number mixing ratio"   "kg kg-1"
+state  real  qso       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO"        "SO number mixing ratio"   "kg kg-1"
+state  real  qso2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO2"        "SO2 number mixing ratio"   "kg kg-1"
+state  real  qso3      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO3"        "SO3 number mixing ratio"   "kg kg-1"
+state  real  qs2o2     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS2O2"        "S2O2 number mixing ratio"   "kg kg-1"
+state  real  qocs      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOCS"        "OCS number mixing ratio"   "kg kg-1"
+state  real  qhso3     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHSO3"        "HSO3 number mixing ratio"   "kg kg-1"
+state  real  qh2so4    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2SO4"        "H2SO4 number mixing ratio"   "kg kg-1"
+state  real  qs2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS2"        "S2 number mixing ratio"   "kg kg-1"
+state  real  qclso2    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLSO2"        "ClSO2 number mixing ratio"   "kg kg-1"
+state  real  qoscl     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOSCL"        "OSCl number mixing ratio"   "kg kg-1"
+state  real  qh2oliq   ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2OLIQ"        "H2O liq number mixing ratio"   "kg kg-1"
+state  real  qh2so4liq ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2SO4LIQ"        "H2SO4 liq number mixing ratio"   "kg kg-1"
+
 
 
@@ -1516,4 +1551,5 @@
 package   newwatersca  mars==12                     -              scalar:qh2o,qh2o_ice,qdust,qdustn,qccn,qccnn
 package   radioac      mars==20                     -              scalar:qtrac1
+package   VenusPC      mars==34                     -              scalar:qco2,qco,qh2,qh2o,qo1d,qo,qo2,qo2dg,qo3,qh,qoh,qho2,qh2o2,qcl,qclo,qcl2,qhcl,qhocl,qclco,qclco3,qcocl2,qss,qso,qso2,qso3,qs2o2,qocs,qhso3,qh2so4,qs2,qclso2,qoscl,qh2oliq,qh2so4liq
 ##### MARS OPTIONS
 ##### MARS OPTIONS
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_initialize_les.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_initialize_les.F	(revision 2020)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_initialize_les.F	(revision 2021)
@@ -87,9 +87,10 @@
    INTEGER, PARAMETER :: nl_max = 1000
    REAL, DIMENSION(nl_max) :: zk, p_in, theta, rho, u, v, qv, pd_in
+   REAL*8, DIMENSION(nl_max) :: pd_in8
    INTEGER :: nl_in
-
 
    INTEGER :: icm,jcm, ii, im1, jj, jm1, loop, error, fid, nxc, nyc
    REAL    :: u_mean,v_mean, f0, p_surf, p_level, qvf, z_at_v, z_at_u
+   REAL*8 :: p_level8
    REAL    :: xrad, yrad, zrad, rad, delt, cof1, cof2
 !   REAL, EXTERNAL :: interp_0
@@ -121,4 +122,5 @@
  REAL, DIMENSION(nl_max) :: venus_hrdyn
  REAL, DIMENSION(nl_max) :: altitude
+ REAL*8, DIMENSION(nl_max) :: trac
 !!MARS
 
@@ -969,4 +971,16 @@
     ENDIF
 
+    IF ((planet.eq."venus") .AND. ( config_flags%mars == 34 )) Then
+       pd_in8(:)=pd_in(:)
+       do i = 1,34
+         call read_tracer(trac,num_scalar,i,nl_in)
+         DO k=1,kte-1
+           p_level = grid%znu(k)*(pd_surf - grid%p_top) + grid%p_top
+           p_level8=p_level
+           scalar(its:ite,k,jts:jte,i+1) =  interp_0_log2( trac, pd_in8, p_level8, nl_in )
+         ENDDO
+       ENDDO
+       !close(unit=22)
+    ENDIF
 
     open(unit=21,file="altitude.txt",action="write")
@@ -1358,3 +1372,26 @@
       end subroutine read_lsf
 
+      subroutine read_tracer(trace,nq,qn,n)
+      implicit none
+      integer n,qn,nq ! qn : number of the tracer
+      real*8 tra(nq-1,n+1)
+      real*8 trace(n+1) !output
+      logical end_of_file
+
+      integer k,j
+
+! first element is the surface
+      open(unit=14,file='input_tracer',form='formatted',status='old')
+      rewind(14)
+      end_of_file = .false.
+        DO k=1,n
+          read(14,*) tra(:,k)
+          write(*,*) k,tra(qn,k)
+      ENDDO
+
+      close(14)
+      trace(:)=tra(qn,:)
+      end subroutine read_tracer
+
+
 END MODULE module_initialize_ideal
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/PREP_MARS/Venus/readmeteo.F90
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/PREP_MARS/Venus/readmeteo.F90	(revision 2020)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/PREP_MARS/Venus/readmeteo.F90	(revision 2021)
@@ -252,20 +252,41 @@
 allocate(aps(altlen),bps(altlen),levels(altlen))
 #ifdef PHOTOCHEM
-nchemtrac = 14
+nchemtrac = 34
 allocate(wtnom(nchemtrac))
-wtnom(1)  = "c_co2"
-wtnom(2)  = "c_co"
-wtnom(3)  = "c_o"
-wtnom(4)  = "c_o1d"
-wtnom(5)  = "c_o2"
-wtnom(6)  = "c_o3"
-wtnom(7)  = "c_h"
-wtnom(8)  = "c_h2"
-wtnom(9)  = "c_oh"
-wtnom(10) = "c_ho2"
-wtnom(11) = "c_h2o2"
-wtnom(12) = "c_ch4"
-wtnom(13) = "c_n2"
-wtnom(14) = "c_ar"
+print*,'PHOTOCHEM2.1'
+wtnom(1)  = "co2"
+wtnom(2)  = "co"
+wtnom(3)  = "h2"
+wtnom(4) = "h2o"
+wtnom(5)  = "o1d"
+wtnom(6)  = "o"
+wtnom(7)  = "o2"
+wtnom(8)  = "o2dg"
+wtnom(9)  = "o3"
+wtnom(10)  = "h"
+wtnom(11)  = "oh"
+wtnom(12)  = "ho2"
+wtnom(13) = "h2o2"
+wtnom(14)  = "cl"
+wtnom(15)  = "clo"
+wtnom(16)  = "cl2"
+wtnom(17)  = "hcl"
+wtnom(18)  = "hocl"
+wtnom(19)  = "clco"
+wtnom(20)  = "clco3"
+wtnom(21)  = "cocl2"
+wtnom(22)  = "s"
+wtnom(23)  = "so"
+wtnom(24)  = "so2"
+wtnom(25)  = "so3"
+wtnom(26)  = "s2o2"
+wtnom(27)  = "ocs"
+wtnom(28)  = "hso3"
+wtnom(29)  = "h2so4"
+wtnom(30)  = "s2"
+wtnom(31)  = "clso2"
+wtnom(32)  = "oscl"
+wtnom(33)  = "h2oliq"
+wtnom(34)  = "h2so4liq"
 allocate(chemtrac(lonlen,latlen,altlen,timelen,nchemtrac))
 chemtrac(:,:,:,:,:)=0
@@ -772,4 +793,7 @@
   vfile(j,i,:,:)=vfile(lonlen+1-j,latlen+1-i,:,:)
   tsoilfile(j,i,:,:)=tsoilfile(lonlen+1-j,latlen+1-i,:,:)
+#ifdef PHOTOCHEM
+ chemtrac(j,i,:,:,:)=chemtrac(lonlen+1-j,latlen+1-i,:,:,:)
+#endif
 ENDDO
 ENDDO  
@@ -1343,18 +1367,478 @@
 !------------------------!
 #ifdef PHOTOCHEM
-    DO i=1,nchemtrac
-       FIELD=wtnom(i)
-       UNITS='units'
-       DESC='desc'
-       XLVL=200100.
-       SLAB=chemtrac(:,:,1,time_out(l),i)
-       ! And now put everything in the destination file
-       ! ... Header
-       write(1) IFV
-       write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
-       write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
-       ! ... Data
-       write(1) SLAB
-    ENDDO
+        FIELD='CO2'
+        UNITS='kg/kg'
+        DESC='CO2 mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),1)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='CO'
+        UNITS='kg/kg'
+        DESC='CO mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),2)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='H2'
+        UNITS='kg/kg'
+        DESC='H2 mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),3)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='H2O'
+        UNITS='kg/kg'
+        DESC='H2O mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),4)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='O1D'
+        UNITS='kg/kg'
+        DESC='O1D mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),5)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='O'
+        UNITS='kg/kg'
+        DESC='O mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),6)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='O2'
+        UNITS='kg/kg'
+        DESC='O2 mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),7)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='O2DG'
+        UNITS='kg/kg'
+        DESC='O2DG mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),8)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='O3'
+        UNITS='kg/kg'
+        DESC='O3 mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),9)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='H'
+        UNITS='kg/kg'
+        DESC='H mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),10)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='OH'
+        UNITS='kg/kg'
+        DESC='OH mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),11)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='HO2'
+        UNITS='kg/kg'
+        DESC='hO2 mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),12)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='H2O2'
+        UNITS='kg/kg'
+        DESC='H2O2 mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),13)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='Cl'
+        UNITS='kg/kg'
+        DESC='Cl mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),14)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='ClO'
+        UNITS='kg/kg'
+        DESC='ClO mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),15)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='Cl2'
+        UNITS='kg/kg'
+        DESC='Cl2 mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),16)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='HCl'
+        UNITS='kg/kg'
+        DESC='HCl mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),17)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='HOCl'
+        UNITS='kg/kg'
+        DESC='HOCl mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),18)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='ClCO'
+        UNITS='kg/kg'
+        DESC='ClCO mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),19)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='ClCO3'
+        UNITS='kg/kg'
+        DESC='ClCO3 mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),20)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='COCl2'
+        UNITS='kg/kg'
+        DESC='COClC2 mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),21)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='S'
+        UNITS='kg/kg'
+        DESC='S mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),22)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='SO'
+        UNITS='kg/kg'
+        DESC='SO mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),23)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='SO2'
+        UNITS='kg/kg'
+        DESC='SO2 mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),24)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='SO3'
+        UNITS='kg/kg'
+        DESC='SO3 mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),25)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='S2O2'
+        UNITS='kg/kg'
+        DESC='S2O2 mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),26)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='OCS'
+        UNITS='kg/kg'
+        DESC='OCS mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),27)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='HSO3'
+        UNITS='kg/kg'
+        DESC='HSO3 mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),28)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='H2SO4'
+        UNITS='kg/kg'
+        DESC='H2SO4 mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),29)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='S2'
+        UNITS='kg/kg'
+        DESC='S2 mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),30)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='ClSO2'
+        UNITS='kg/kg'
+        DESC='ClSO2 mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),31)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='OSCl'
+        UNITS='kg/kg'
+        DESC='OSCl mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),32)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='H2Oliq'
+        UNITS='kg/kg'
+        DESC='H2Oliq mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),33)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
+
+
+        FIELD='H2SO4liq'
+        UNITS='kg/kg'
+        DESC='H2SO4liq mixing ratio'
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,1,time_out(l),34)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+        write(1) SLAB
 #endif
 
@@ -1686,20 +2170,513 @@
 !------------------------!
 #ifdef PHOTOCHEM
-    DO i=1,nchemtrac
-       FIELD=wtnom(i)
-       UNITS='units'
-       DESC='desc'
-       DO k = 1,altlen
-         XLVL=levels(k)
-         SLAB=chemtrac(:,:,k,time_out(l),i)
-         ! And now put everything in the destination file
-         ! ... Header
-         write(1) IFV
-         write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
-         write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
-         ! ... Data
-         write(1) SLAB
-       END DO
-    ENDDO
+FIELD='CO2'
+UNITS='kg/kg'
+DESC='CO2 mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),1)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='CO'
+UNITS='kg/kg'
+DESC='CO mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),2)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='H2'
+UNITS='kg/kg'
+DESC='H2 mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),3)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='H2O'
+UNITS='kg/kg'
+DESC='H2O mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),4)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='O1D'
+UNITS='kg/kg'
+DESC='O1D mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),5)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='O'
+UNITS='kg/kg'
+DESC='O mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),6)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='O2'
+UNITS='kg/kg'
+DESC='O2 mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),7)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='O2DG'
+UNITS='kg/kg'
+DESC='O2DG mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),8)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='O3'
+UNITS='kg/kg'
+DESC='O3 mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),9)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='H'
+UNITS='kg/kg'
+DESC='H mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),10)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='OH'
+UNITS='kg/kg'
+DESC='OH mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),11)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='HO2'
+UNITS='kg/kg'
+DESC='hO2 mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),12)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='H2O2'
+UNITS='kg/kg'
+DESC='H2O2 mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),13)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='Cl'
+UNITS='kg/kg'
+DESC='Cl mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),14)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='ClO'
+UNITS='kg/kg'
+DESC='ClO mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),15)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='Cl2'
+UNITS='kg/kg'
+DESC='Cl2 mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),16)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='HCl'
+UNITS='kg/kg'
+DESC='HCl mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),17)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='HOCl'
+UNITS='kg/kg'
+DESC='HOCl mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),18)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='ClCO'
+UNITS='kg/kg'
+DESC='ClCO mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),19)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='ClCO3'
+UNITS='kg/kg'
+DESC='ClCO3 mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),20)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='COCl2'
+UNITS='kg/kg'
+DESC='COClC2 mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),21)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='S'
+UNITS='kg/kg'
+DESC='S mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),22)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='SO'
+UNITS='kg/kg'
+DESC='SO mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),23)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='SO2'
+UNITS='kg/kg'
+DESC='SO2 mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),24)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='SO3'
+UNITS='kg/kg'
+DESC='SO3 mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),25)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='S2O2'
+UNITS='kg/kg'
+DESC='S2O2 mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),26)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='OCS'
+UNITS='kg/kg'
+DESC='OCS mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),27)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='HSO3'
+UNITS='kg/kg'
+DESC='HSO3 mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),28)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='H2SO4'
+UNITS='kg/kg'
+DESC='H2SO4 mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),29)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='S2'
+UNITS='kg/kg'
+DESC='S2 mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),30)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='ClSO2'
+UNITS='kg/kg'
+DESC='ClSO2 mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),31)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='OSCl'
+UNITS='kg/kg'
+DESC='OSCl mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),32)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='H2Oliq'
+UNITS='kg/kg'
+DESC='H2Oliq mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),33)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
+
+FIELD='H2SO4liq'
+UNITS='kg/kg'
+DESC='H2SO4liq mixing ratio'
+DO k = 1,altlen
+        XLVL=levels(k)
+        SLAB=chemtrac(:,:,k,time_out(l),34)
+                ! And now put everything in the destination file
+                ! ... Header
+        write(1) IFV
+        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
+        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
+                ! ... Data
+                write(1) SLAB
+END DO
 #endif
 
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM	(revision 2020)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM	(revision 2021)
@@ -171,4 +171,37 @@
 ####MARS MARS MARS MARS MARS
 ####MARS MARS MARS MARS MARS
+###VENUS TRACER
+state  real  qco       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCO"        "CO number mixing ratio"   "kg kg-1"
+state  real  qh2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2"        "H2 number mixing ratio"   "kg kg-1"
+state  real  qo1d      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO1D"        "O1d number mixing ratio"   "kg kg-1"
+state  real  qo        ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO"        "O number mixing ratio"   "kg kg-1"
+state  real  qo2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO2"        "O2 number mixing ratio"   "kg kg-1"
+state  real  qo2dg     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO2dg"        "O2dg number mixing ratio"   "kg kg-1"
+state  real  qo3       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO3"        "O3 number mixing ratio"   "kg kg-1"
+state  real  qh        ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH"        "H number mixing ratio"   "kg kg-1"
+state  real  qoh       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOH"        "OH number mixing ratio"   "kg kg-1"
+state  real  qho2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHO2"        "HO2 number mixing ratio"   "kg kg-1"
+state  real  qh2o2     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2O2"        "H2O2 number mixing ratio"   "kg kg-1"
+state  real  qcl       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCL"        "Cl number mixing ratio"   "kg kg-1"
+state  real  qclo      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLO"        "ClO number mixing ratio"   "kg kg-1"
+state  real  qcl2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCL2"        "Cl2 number mixing ratio"   "kg kg-1"
+state  real  qhcl      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHCL"        "HCl number mixing ratio"   "kg kg-1"
+state  real  qhocl     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHOCL"        "HOCl number mixing ratio"   "kg kg-1"
+state  real  qclco     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLCO"        "ClCO number mixing ratio"   "kg kg-1"
+state  real  qclco3    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLCO3"        "ClCO3 number mixing ratio"   "kg kg-1"
+state  real  qcocl2    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCOCL2"        "COCl2 number mixing ratio"   "kg kg-1"
+state  real  qss       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS"        "S number mixing ratio"   "kg kg-1"
+state  real  qso       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO"        "SO number mixing ratio"   "kg kg-1"
+state  real  qso2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO2"        "SO2 number mixing ratio"   "kg kg-1"
+state  real  qso3      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO3"        "SO3 number mixing ratio"   "kg kg-1"
+state  real  qs2o2     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS2O2"        "S2O2 number mixing ratio"   "kg kg-1"
+state  real  qocs      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOCS"        "OCS number mixing ratio"   "kg kg-1"
+state  real  qhso3     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHSO3"        "HSO3 number mixing ratio"   "kg kg-1"
+state  real  qh2so4    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2SO4"        "H2SO4 number mixing ratio"   "kg kg-1"
+state  real  qs2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS2"        "S2 number mixing ratio"   "kg kg-1"
+state  real  qclso2    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLSO2"        "ClSO2 number mixing ratio"   "kg kg-1"
+state  real  qoscl     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOSCL"        "OSCl number mixing ratio"   "kg kg-1"
+state  real  qh2oliq   ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2OLIQ"        "H2O liq number mixing ratio"   "kg kg-1"
+state  real  qh2so4liq ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2SO4LIQ"        "H2SO4 liq number mixing ratio"   "kg kg-1"
 
 
@@ -1392,4 +1425,5 @@
 package   newwater     mars==11                     -              scalar:qh2o,qh2o_ice,qdust,qdustn
 package   newwatersca  mars==12                     -              scalar:qh2o,qh2o_ice,qdust,qdustn,qccn,qccnn
+package   VenusPC      mars==34                     -              scalar:qco2,qco,qh2,qh2o,qo1d,qo,qo2,qo2dg,qo3,qh,qoh,qho2,qh2o2,qcl,qclo,qcl2,qhcl,qhocl,qclco,qclco3,qcocl2,qss,qso,qso2,qso3,qs2o2,qocs,qhso3,qh2so4,qs2,qclso2,qoscl,qh2oliq,qh2so4liq
 package   photochem    mars==42                     -              scalar:qco2,chem_co,chem_o,chem_o1d,chem_o2,chem_o3,chem_h,chem_h2,chem_oh,chem_ho2,chem_h2o2,chem_ch4,chem_n2,chem_ar,qh2o_ice,qh2o,qdust,qdustn
 ##### MARS OPTIONS
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM.newphys
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM.newphys	(revision 2020)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM.newphys	(revision 2021)
@@ -156,6 +156,37 @@
 ####MARS MARS MARS MARS MARS
 ####MARS MARS MARS MARS MARS
-
-
+###VENUS TRACER
+state  real  qco       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCO"        "CO number mixing ratio"   "kg kg-1"
+state  real  qh2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2"        "H2 number mixing ratio"   "kg kg-1"
+state  real  qo1d      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO1D"        "O1d number mixing ratio"   "kg kg-1"
+state  real  qo        ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO"        "O number mixing ratio"   "kg kg-1"
+state  real  qo2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO2"        "O2 number mixing ratio"   "kg kg-1"
+state  real  qo2dg     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO2dg"        "O2dg number mixing ratio"   "kg kg-1"
+state  real  qo3       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO3"        "O3 number mixing ratio"   "kg kg-1"
+state  real  qh        ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH"        "H number mixing ratio"   "kg kg-1"
+state  real  qoh       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOH"        "OH number mixing ratio"   "kg kg-1"
+state  real  qho2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHO2"        "HO2 number mixing ratio"   "kg kg-1"
+state  real  qh2o2     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2O2"        "H2O2 number mixing ratio"   "kg kg-1"
+state  real  qcl       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCL"        "Cl number mixing ratio"   "kg kg-1"
+state  real  qclo      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLO"        "ClO number mixing ratio"   "kg kg-1"
+state  real  qcl2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCL2"        "Cl2 number mixing ratio"   "kg kg-1"
+state  real  qhcl      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHCL"        "HCl number mixing ratio"   "kg kg-1"
+state  real  qhocl     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHOCL"        "HOCl number mixing ratio"   "kg kg-1"
+state  real  qclco     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLCO"        "ClCO number mixing ratio"   "kg kg-1"
+state  real  qclco3    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLCO3"        "ClCO3 number mixing ratio"   "kg kg-1"
+state  real  qcocl2    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCOCL2"        "COCl2 number mixing ratio"   "kg kg-1"
+state  real  qss       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS"        "S number mixing ratio"   "kg kg-1"
+state  real  qso       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO"        "SO number mixing ratio"   "kg kg-1"
+state  real  qso2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO2"        "SO2 number mixing ratio"   "kg kg-1"
+state  real  qso3      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO3"        "SO3 number mixing ratio"   "kg kg-1"
+state  real  qs2o2     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS2O2"        "S2O2 number mixing ratio"   "kg kg-1"
+state  real  qocs      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOCS"        "OCS number mixing ratio"   "kg kg-1"
+state  real  qhso3     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHSO3"        "HSO3 number mixing ratio"   "kg kg-1"
+state  real  qh2so4    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2SO4"        "H2SO4 number mixing ratio"   "kg kg-1"
+state  real  qs2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS2"        "S2 number mixing ratio"   "kg kg-1"
+state  real  qclso2    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLSO2"        "ClSO2 number mixing ratio"   "kg kg-1"
+state  real  qoscl     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOSCL"        "OSCl number mixing ratio"   "kg kg-1"
+state  real  qh2oliq   ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2OLIQ"        "H2O liq number mixing ratio"   "kg kg-1"
+state  real  qh2so4liq ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2SO4LIQ"        "H2SO4 liq number mixing ratio"   "kg kg-1"
 
 ################################################################################
@@ -1378,4 +1409,5 @@
 package   newwatersca  mars==12                     -              scalar:qh2o,qh2o_ice,qdust,qdustn,qccn,qccnn
 package   cloudco2     mars==32                     -              scalar:qh2o,qh2o_ice,qdust,qdustn,qccn,qccnn,qco2,qco2_ice,qccn_co2,qccnn_co2
+package   VenusPC      mars==34                     -              scalar:qco2,qco,qh2,qh2o,qo1d,qo,qo2,qo2dg,qo3,qh,qoh,qho2,qh2o2,qcl,qclo,qcl2,qhcl,qhocl,qclco,qclco3,qcocl2,qss,qso,qso2,qso3,qs2o2,qocs,qhso3,qh2so4,qs2,qclso2,qoscl,qh2oliq,qh2so4liq
 package   photochem    mars==42                     -              scalar:qco2,chem_co,chem_o,chem_o1d,chem_o2,chem_o3,chem_h,chem_h2,chem_oh,chem_ho2,chem_h2o2,chem_ch4,chem_n2,chem_ar,qh2o_ice,qh2o,qdust,qdustn
 ##### MARS OPTIONS
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_real.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_real.F	(revision 2020)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_real.F	(revision 2021)
@@ -1418,4 +1418,281 @@
 #endif
 
+!VENUS photochemistry
+if ( config_flags%mars == 34 ) then
+       print*,'grid%em_qco2_gc',grid%em_qco2_gc(0,:,0)
+       CALL vert_interp_old ( grid%em_qco2_gc , grid%em_pd_gc , scalar(:,:,:,2) , grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order , lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qco_gc , grid%em_pd_gc , scalar(:,:,:,3), grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order ,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qh2_gc , grid%em_pd_gc , scalar(:,:,:,4),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qh2o_gc , grid%em_pd_gc , scalar(:,:,:,5),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qo1d_gc , grid%em_pd_gc , scalar(:,:,:,6),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qo_gc , grid%em_pd_gc , scalar(:,:,:,7),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qo2_gc , grid%em_pd_gc , scalar(:,:,:,8),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qo2dg_gc , grid%em_pd_gc , scalar(:,:,:,9),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qo3_gc , grid%em_pd_gc , scalar(:,:,:,10),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qh_gc , grid%em_pd_gc , scalar(:,:,:,11),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qoh_gc , grid%em_pd_gc , scalar(:,:,:,12), grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qho2_gc , grid%em_pd_gc , scalar(:,:,:,13),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qh2o2_gc , grid%em_pd_gc , scalar(:,:,:,14),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qcl_gc , grid%em_pd_gc , scalar(:,:,:,15),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order ,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qclo_gc , grid%em_pd_gc , scalar(:,:,:,16),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qcl2_gc , grid%em_pd_gc , scalar(:,:,:,17),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qhcl_gc , grid%em_pd_gc , scalar(:,:,:,18),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qhocl_gc , grid%em_pd_gc , scalar(:,:,:,19),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qclco_gc , grid%em_pd_gc , scalar(:,:,:,20),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qclco3_gc , grid%em_pd_gc , scalar(:,:,:,21),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qcocl2_gc , grid%em_pd_gc , scalar(:,:,:,22),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qs_gc , grid%em_pd_gc , scalar(:,:,:,23), grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order ,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qso_gc , grid%em_pd_gc , scalar(:,:,:,24),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qso2_gc , grid%em_pd_gc , scalar(:,:,:,25), grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order ,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qo3_gc , grid%em_pd_gc , scalar(:,:,:,26),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order ,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qs2o2_gc , grid%em_pd_gc , scalar(:,:,:,27),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qocs_gc , grid%em_pd_gc , scalar(:,:,:,28),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qhso3_gc , grid%em_pd_gc , scalar(:,:,:,29),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qh2so4_gc , grid%em_pd_gc , scalar(:,:,:,30),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qs2_gc , grid%em_pd_gc , scalar(:,:,:,31),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qclso2_gc , grid%em_pd_gc , scalar(:,:,:,32),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qoscl_gc , grid%em_pd_gc , scalar(:,:,:,33),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qh2oliq_gc , grid%em_pd_gc , scalar(:,:,:,34),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+       CALL vert_interp_old ( grid%em_qh2so4liq_gc , grid%em_pd_gc , scalar(:,:,:,35),grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+
+endif
+
 !!! we want any scalar (i.e. tracer) to be positive
 !!! and because of interpolation it is possible that negative values occur...
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F	(revision 2020)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F	(revision 2021)
@@ -2751,4 +2751,5 @@
  .OR. (  config_flags%mars .eq. 11                    ) &    !! 11  pass EVERYTHING
  .OR. (  config_flags%mars .eq. 12                    ) &    !! 12  pass EVERYTHING 
+ .OR. (  config_flags%mars .eq. 34                    ) &
  .OR. config_flags%nested ) THEN                             !! *   pass all tracers if nested
          CALL relax_bdy_scalar ( scalar_tend(ims,kms,jms,is),            &
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/main/real_em.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/main/real_em.F	(revision 2020)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/main/real_em.F	(revision 2021)
@@ -568,4 +568,33 @@
    REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q5bdy3dtemp1
    REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q6bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q7bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q8bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q9bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q10bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q11bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q12bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q13bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q14bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q15bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q16bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q17bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q18bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q19bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q20bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q21bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q22bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q23bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q24bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q25bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q26bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q27bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q28bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q29bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q30bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q31bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q32bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q33bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q34bdy3dtemp1
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q35bdy3dtemp1
 !!!!***MARS <<
    REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: mbdy2dtemp1
@@ -577,4 +606,33 @@
    REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q5bdy3dtemp2
    REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q6bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q7bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q8bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q9bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q10bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q11bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q12bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q13bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q14bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q15bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q16bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q17bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q18bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q19bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q20bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q21bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q22bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q23bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q24bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q25bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q26bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q27bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q28bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q29bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q30bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q31bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q32bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q33bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q34bdy3dtemp2
+   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q35bdy3dtemp2
 !!!!***MARS <<
    REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: mbdy2dtemp2
@@ -631,4 +689,33 @@
       IF ( ALLOCATED ( q5bdy3dtemp1 ) ) DEALLOCATE ( q5bdy3dtemp1 )
       IF ( ALLOCATED ( q6bdy3dtemp1 ) ) DEALLOCATE ( q6bdy3dtemp1 )
+      IF ( ALLOCATED ( q7bdy3dtemp1 ) ) DEALLOCATE ( q7bdy3dtemp1 )
+      IF ( ALLOCATED ( q8bdy3dtemp1 ) ) DEALLOCATE ( q8bdy3dtemp1 )
+      IF ( ALLOCATED ( q9bdy3dtemp1 ) ) DEALLOCATE ( q9bdy3dtemp1 )
+      IF ( ALLOCATED ( q10bdy3dtemp1 ) ) DEALLOCATE ( q10bdy3dtemp1 )
+      IF ( ALLOCATED ( q11bdy3dtemp1 ) ) DEALLOCATE ( q11bdy3dtemp1 )
+      IF ( ALLOCATED ( q12bdy3dtemp1 ) ) DEALLOCATE ( q12bdy3dtemp1 )
+      IF ( ALLOCATED ( q13bdy3dtemp1 ) ) DEALLOCATE ( q13bdy3dtemp1 )
+      IF ( ALLOCATED ( q14bdy3dtemp1 ) ) DEALLOCATE ( q14bdy3dtemp1 )
+      IF ( ALLOCATED ( q15bdy3dtemp1 ) ) DEALLOCATE ( q15bdy3dtemp1 )
+      IF ( ALLOCATED ( q16bdy3dtemp1 ) ) DEALLOCATE ( q16bdy3dtemp1 )
+      IF ( ALLOCATED ( q17bdy3dtemp1 ) ) DEALLOCATE ( q17bdy3dtemp1 )
+      IF ( ALLOCATED ( q18bdy3dtemp1 ) ) DEALLOCATE ( q18bdy3dtemp1 )
+      IF ( ALLOCATED ( q19bdy3dtemp1 ) ) DEALLOCATE ( q19bdy3dtemp1 )
+      IF ( ALLOCATED ( q20bdy3dtemp1 ) ) DEALLOCATE ( q20bdy3dtemp1 )
+      IF ( ALLOCATED ( q21bdy3dtemp1 ) ) DEALLOCATE ( q21bdy3dtemp1 )
+      IF ( ALLOCATED ( q22bdy3dtemp1 ) ) DEALLOCATE ( q22bdy3dtemp1 )
+      IF ( ALLOCATED ( q23bdy3dtemp1 ) ) DEALLOCATE ( q23bdy3dtemp1 )
+      IF ( ALLOCATED ( q24bdy3dtemp1 ) ) DEALLOCATE ( q24bdy3dtemp1 )
+      IF ( ALLOCATED ( q25bdy3dtemp1 ) ) DEALLOCATE ( q25bdy3dtemp1 )
+      IF ( ALLOCATED ( q26bdy3dtemp1 ) ) DEALLOCATE ( q26bdy3dtemp1 )
+      IF ( ALLOCATED ( q27bdy3dtemp1 ) ) DEALLOCATE ( q27bdy3dtemp1 )
+      IF ( ALLOCATED ( q28bdy3dtemp1 ) ) DEALLOCATE ( q28bdy3dtemp1 )
+      IF ( ALLOCATED ( q29bdy3dtemp1 ) ) DEALLOCATE ( q29bdy3dtemp1 )
+      IF ( ALLOCATED ( q30bdy3dtemp1 ) ) DEALLOCATE ( q30bdy3dtemp1 )
+      IF ( ALLOCATED ( q31bdy3dtemp1 ) ) DEALLOCATE ( q31bdy3dtemp1 )
+      IF ( ALLOCATED ( q32bdy3dtemp1 ) ) DEALLOCATE ( q32bdy3dtemp1 )
+      IF ( ALLOCATED ( q33bdy3dtemp1 ) ) DEALLOCATE ( q33bdy3dtemp1 )
+      IF ( ALLOCATED ( q34bdy3dtemp1 ) ) DEALLOCATE ( q34bdy3dtemp1 )
+      IF ( ALLOCATED ( q35bdy3dtemp1 ) ) DEALLOCATE ( q35bdy3dtemp1 )
 !!!!***MARS <<
       IF ( ALLOCATED ( mbdy2dtemp1 ) ) DEALLOCATE ( mbdy2dtemp1 )
@@ -644,4 +731,33 @@
       IF ( ALLOCATED ( q5bdy3dtemp2 ) ) DEALLOCATE ( q5bdy3dtemp2 )
       IF ( ALLOCATED ( q6bdy3dtemp2 ) ) DEALLOCATE ( q6bdy3dtemp2 )
+      IF ( ALLOCATED ( q7bdy3dtemp2 ) ) DEALLOCATE ( q7bdy3dtemp2 )
+      IF ( ALLOCATED ( q8bdy3dtemp2 ) ) DEALLOCATE ( q8bdy3dtemp2 )
+      IF ( ALLOCATED ( q9bdy3dtemp2 ) ) DEALLOCATE ( q9bdy3dtemp2 )
+      IF ( ALLOCATED ( q10bdy3dtemp2 ) ) DEALLOCATE ( q10bdy3dtemp2 )
+      IF ( ALLOCATED ( q11bdy3dtemp2 ) ) DEALLOCATE ( q11bdy3dtemp2 )
+      IF ( ALLOCATED ( q12bdy3dtemp2 ) ) DEALLOCATE ( q12bdy3dtemp2 )
+      IF ( ALLOCATED ( q13bdy3dtemp2 ) ) DEALLOCATE ( q13bdy3dtemp2 )
+      IF ( ALLOCATED ( q14bdy3dtemp2 ) ) DEALLOCATE ( q14bdy3dtemp2 )
+      IF ( ALLOCATED ( q15bdy3dtemp2 ) ) DEALLOCATE ( q15bdy3dtemp2 )
+      IF ( ALLOCATED ( q16bdy3dtemp2 ) ) DEALLOCATE ( q16bdy3dtemp2 )
+      IF ( ALLOCATED ( q17bdy3dtemp2 ) ) DEALLOCATE ( q17bdy3dtemp2 )
+      IF ( ALLOCATED ( q18bdy3dtemp2 ) ) DEALLOCATE ( q18bdy3dtemp2 )
+      IF ( ALLOCATED ( q19bdy3dtemp2 ) ) DEALLOCATE ( q19bdy3dtemp2 )
+      IF ( ALLOCATED ( q20bdy3dtemp2 ) ) DEALLOCATE ( q20bdy3dtemp2 )
+      IF ( ALLOCATED ( q21bdy3dtemp2 ) ) DEALLOCATE ( q21bdy3dtemp2 )
+      IF ( ALLOCATED ( q22bdy3dtemp2 ) ) DEALLOCATE ( q22bdy3dtemp2 )
+      IF ( ALLOCATED ( q23bdy3dtemp2 ) ) DEALLOCATE ( q23bdy3dtemp2 )
+      IF ( ALLOCATED ( q24bdy3dtemp2 ) ) DEALLOCATE ( q24bdy3dtemp2 )
+      IF ( ALLOCATED ( q25bdy3dtemp2 ) ) DEALLOCATE ( q25bdy3dtemp2 )
+      IF ( ALLOCATED ( q26bdy3dtemp2 ) ) DEALLOCATE ( q26bdy3dtemp2 )
+      IF ( ALLOCATED ( q27bdy3dtemp2 ) ) DEALLOCATE ( q27bdy3dtemp2 )
+      IF ( ALLOCATED ( q28bdy3dtemp2 ) ) DEALLOCATE ( q28bdy3dtemp2 )
+      IF ( ALLOCATED ( q29bdy3dtemp2 ) ) DEALLOCATE ( q29bdy3dtemp2 )
+      IF ( ALLOCATED ( q30bdy3dtemp2 ) ) DEALLOCATE ( q30bdy3dtemp2 )
+      IF ( ALLOCATED ( q31bdy3dtemp2 ) ) DEALLOCATE ( q31bdy3dtemp2 )
+      IF ( ALLOCATED ( q32bdy3dtemp2 ) ) DEALLOCATE ( q32bdy3dtemp2 )
+      IF ( ALLOCATED ( q33bdy3dtemp2 ) ) DEALLOCATE ( q33bdy3dtemp2 )
+      IF ( ALLOCATED ( q34bdy3dtemp2 ) ) DEALLOCATE ( q34bdy3dtemp2 )
+      IF ( ALLOCATED ( q35bdy3dtemp2 ) ) DEALLOCATE ( q35bdy3dtemp2 )
 !!!!***MARS <<
       IF ( ALLOCATED ( mbdy2dtemp2 ) ) DEALLOCATE ( mbdy2dtemp2 )
@@ -658,4 +774,33 @@
       ALLOCATE ( q5bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
       ALLOCATE ( q6bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q7bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q8bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q9bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q10bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q11bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q12bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q13bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q14bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q15bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q16bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q17bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q18bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q19bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q20bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q21bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q22bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q23bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q24bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q25bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q26bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q27bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q28bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q29bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q30bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q31bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q32bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q33bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q34bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q35bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
 !!!!***MARS <<
       ALLOCATE ( mbdy2dtemp1(ims:ime,1:1,    jms:jme) )
@@ -671,4 +816,33 @@
       ALLOCATE ( q5bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
       ALLOCATE ( q6bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q7bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q8bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q9bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q10bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q11bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q12bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q13bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q14bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q15bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q16bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q17bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q18bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q19bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q20bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q21bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q22bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q23bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q24bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q25bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q26bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q27bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q28bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q29bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q30bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q31bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q32bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q33bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q34bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
+      ALLOCATE ( q35bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
 !!!!***MARS <<
       ALLOCATE ( mbdy2dtemp2(ims:ime,1:1,    jms:jme) )
@@ -741,4 +915,76 @@
       CALL couple ( grid%em_mu_2 , grid%em_mub , q6bdy3dtemp1 , grid%scalar(:,:,:,7) , 't' , grid%msft , &
                     ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms, kme, ips, ipe, jps, jpe, kps, kpe )
+ENDIF
+
+!VENUS
+IF (config_flags%mars .eq. 34) THEN
+      CALL couple ( grid%em_mu_2 , grid%em_mub , qbdy3dtemp1 ,grid%scalar(:,:,:,2) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms, kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q2bdy3dtemp1,grid%scalar(:,:,:,3) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q3bdy3dtemp1,grid%scalar(:,:,:,4) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q4bdy3dtemp1,grid%scalar(:,:,:,5) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q5bdy3dtemp1,grid%scalar(:,:,:,6) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q6bdy3dtemp1,grid%scalar(:,:,:,7) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q7bdy3dtemp1,grid%scalar(:,:,:,8) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q8bdy3dtemp1,grid%scalar(:,:,:,9) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q9bdy3dtemp1,grid%scalar(:,:,:,10) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q10bdy3dtemp1,grid%scalar(:,:,:,11) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q11bdy3dtemp1,grid%scalar(:,:,:,12) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q12bdy3dtemp1,grid%scalar(:,:,:,13) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q13bdy3dtemp1,grid%scalar(:,:,:,14) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q14bdy3dtemp1,grid%scalar(:,:,:,15) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q15bdy3dtemp1,grid%scalar(:,:,:,16) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q16bdy3dtemp1,grid%scalar(:,:,:,17) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q17bdy3dtemp1,grid%scalar(:,:,:,18) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q18bdy3dtemp1,grid%scalar(:,:,:,19) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q19bdy3dtemp1,grid%scalar(:,:,:,20) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q20bdy3dtemp1,grid%scalar(:,:,:,21) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q21bdy3dtemp1,grid%scalar(:,:,:,22) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q22bdy3dtemp1,grid%scalar(:,:,:,23) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q23bdy3dtemp1,grid%scalar(:,:,:,24) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q24bdy3dtemp1,grid%scalar(:,:,:,25) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q25bdy3dtemp1,grid%scalar(:,:,:,26) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q26bdy3dtemp1,grid%scalar(:,:,:,27) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q27bdy3dtemp1,grid%scalar(:,:,:,28) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q28bdy3dtemp1,grid%scalar(:,:,:,29) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q29bdy3dtemp1,grid%scalar(:,:,:,30) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q30bdy3dtemp1,grid%scalar(:,:,:,31) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q31bdy3dtemp1,grid%scalar(:,:,:,32) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q32bdy3dtemp1,grid%scalar(:,:,:,33) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q33bdy3dtemp1,grid%scalar(:,:,:,34) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q34bdy3dtemp1,grid%scalar(:,:,:,35) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
 ENDIF
 !!!!!MARS!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -833,4 +1079,144 @@
                                                                  ips , ipe , jps , jpe , kps , kpe )
 ENDIF
+
+!VENUS
+IF (config_flags%mars .eq. 34) THEN
+      CALL stuff_bdy     ( qbdy3dtemp1 , grid%scalar_b(:,:,:,:,2)   , 'T' , ijds, ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q2bdy3dtemp1 , grid%scalar_b(:,:,:,:,3)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q3bdy3dtemp1 , grid%scalar_b(:,:,:,:,4)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q4bdy3dtemp1 , grid%scalar_b(:,:,:,:,5)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q5bdy3dtemp1 , grid%scalar_b(:,:,:,:,6)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q6bdy3dtemp1 , grid%scalar_b(:,:,:,:,7)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q7bdy3dtemp1 , grid%scalar_b(:,:,:,:,8)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q8bdy3dtemp1 , grid%scalar_b(:,:,:,:,9)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q9bdy3dtemp1 , grid%scalar_b(:,:,:,:,10)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q10bdy3dtemp1 , grid%scalar_b(:,:,:,:,11)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q11bdy3dtemp1 , grid%scalar_b(:,:,:,:,12)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q12bdy3dtemp1 , grid%scalar_b(:,:,:,:,13)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q13bdy3dtemp1 , grid%scalar_b(:,:,:,:,14)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q14bdy3dtemp1 , grid%scalar_b(:,:,:,:,15)   , 'T' ,ijds , ijde , spec_bdy_width      , &                                                                 
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q15bdy3dtemp1 , grid%scalar_b(:,:,:,:,16)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q16bdy3dtemp1 , grid%scalar_b(:,:,:,:,17)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q17bdy3dtemp1 , grid%scalar_b(:,:,:,:,18)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q18bdy3dtemp1 , grid%scalar_b(:,:,:,:,19)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q19bdy3dtemp1 , grid%scalar_b(:,:,:,:,20)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q20bdy3dtemp1 , grid%scalar_b(:,:,:,:,21)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q21bdy3dtemp1 , grid%scalar_b(:,:,:,:,22)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q22bdy3dtemp1 , grid%scalar_b(:,:,:,:,23)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q23bdy3dtemp1 , grid%scalar_b(:,:,:,:,24)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q24bdy3dtemp1 , grid%scalar_b(:,:,:,:,25)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q25bdy3dtemp1 , grid%scalar_b(:,:,:,:,26)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q26bdy3dtemp1 , grid%scalar_b(:,:,:,:,27)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q27bdy3dtemp1 , grid%scalar_b(:,:,:,:,28)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q28bdy3dtemp1 , grid%scalar_b(:,:,:,:,29)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q29bdy3dtemp1 , grid%scalar_b(:,:,:,:,30)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q30bdy3dtemp1 , grid%scalar_b(:,:,:,:,31)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q31bdy3dtemp1 , grid%scalar_b(:,:,:,:,33)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q32bdy3dtemp1 , grid%scalar_b(:,:,:,:,33)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q33bdy3dtemp1 , grid%scalar_b(:,:,:,:,34)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q34bdy3dtemp1 , grid%scalar_b(:,:,:,:,35)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+ENDIF
 !!!!MARS <<
 
@@ -911,4 +1297,76 @@
                     ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms, kme, ips, ipe, jps, jpe, kps, kpe )
 ENDIF
+
+!VENUS
+IF (config_flags%mars .eq. 34) THEN
+      CALL couple ( grid%em_mu_2 , grid%em_mub , qbdy3dtemp2 ,grid%scalar(:,:,:,2) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms, kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q2bdy3dtemp2,grid%scalar(:,:,:,3) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q3bdy3dtemp2,grid%scalar(:,:,:,4) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q4bdy3dtemp2,grid%scalar(:,:,:,5) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q5bdy3dtemp2,grid%scalar(:,:,:,6) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q6bdy3dtemp2,grid%scalar(:,:,:,7) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q7bdy3dtemp2,grid%scalar(:,:,:,8) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q8bdy3dtemp2,grid%scalar(:,:,:,9) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q9bdy3dtemp2,grid%scalar(:,:,:,10) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q10bdy3dtemp2,grid%scalar(:,:,:,11) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q11bdy3dtemp2,grid%scalar(:,:,:,12) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q12bdy3dtemp2,grid%scalar(:,:,:,13) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q13bdy3dtemp2,grid%scalar(:,:,:,14) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q14bdy3dtemp2,grid%scalar(:,:,:,15) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q15bdy3dtemp2,grid%scalar(:,:,:,16) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q16bdy3dtemp2,grid%scalar(:,:,:,17) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q17bdy3dtemp2,grid%scalar(:,:,:,18) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q18bdy3dtemp2,grid%scalar(:,:,:,19) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q19bdy3dtemp2,grid%scalar(:,:,:,20) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q20bdy3dtemp2,grid%scalar(:,:,:,21) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q21bdy3dtemp2,grid%scalar(:,:,:,22) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q22bdy3dtemp2,grid%scalar(:,:,:,23) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q23bdy3dtemp2,grid%scalar(:,:,:,24) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q24bdy3dtemp2,grid%scalar(:,:,:,25) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q25bdy3dtemp2,grid%scalar(:,:,:,26) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q26bdy3dtemp2,grid%scalar(:,:,:,27) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q27bdy3dtemp2,grid%scalar(:,:,:,28) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q28bdy3dtemp2,grid%scalar(:,:,:,29) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q29bdy3dtemp2,grid%scalar(:,:,:,30) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q30bdy3dtemp2,grid%scalar(:,:,:,31) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q31bdy3dtemp2,grid%scalar(:,:,:,32) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q32bdy3dtemp2,grid%scalar(:,:,:,33) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q33bdy3dtemp2,grid%scalar(:,:,:,34) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+      CALL couple ( grid%em_mu_2 , grid%em_mub , q34bdy3dtemp2,grid%scalar(:,:,:,35) , 't' , grid%msft , &
+                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
+ENDIF
 !!!!!MARS <<
 
@@ -1010,4 +1468,179 @@
                                                             ims , ime , jms , jme , kms , kme , &
                                                             ips , ipe , jps , jpe , kps , kpe )
+ENDIF
+
+!VENUS
+IF (config_flags%mars .eq. 34) THEN
+      CALL stuff_bdytend ( qbdy3dtemp2 , qbdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,2) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q2bdy3dtemp2 , q2bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,3) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q3bdy3dtemp2 , q3bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,4) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q4bdy3dtemp2 , q4bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,5) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q5bdy3dtemp2 , q5bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,6) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q6bdy3dtemp2 , q6bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,7) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q7bdy3dtemp2 , q7bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,8) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q8bdy3dtemp2 , q8bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,9) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q9bdy3dtemp2 , q9bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,10) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q10bdy3dtemp2 , q10bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,11) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q11bdy3dtemp2 , q11bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,12) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q12bdy3dtemp2 , q12bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,13) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q13bdy3dtemp2 , q13bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,14) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q14bdy3dtemp2 , q14bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,15) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q15bdy3dtemp2 , q15bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,16) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q16bdy3dtemp2 , q16bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,17) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q17bdy3dtemp2 , q17bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,18) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q18bdy3dtemp2 , q18bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,19) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q19bdy3dtemp2 , q19bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,20) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q20bdy3dtemp2 , q20bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,21) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q21bdy3dtemp2 , q21bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,22) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q22bdy3dtemp2 , q22bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,23) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q23bdy3dtemp2 , q23bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,24) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q24bdy3dtemp2 , q24bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,25) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q25bdy3dtemp2 , q25bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,26) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q26bdy3dtemp2 , q26bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,27) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q27bdy3dtemp2 , q27bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,28) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q28bdy3dtemp2 , q28bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,29) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q29bdy3dtemp2 , q29bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,30) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q30bdy3dtemp2 , q30bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,31) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q31bdy3dtemp2 , q31bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,32) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q32bdy3dtemp2 , q32bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,33) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q33bdy3dtemp2 , q33bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,34) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+      CALL stuff_bdytend ( q34bdy3dtemp2 , q34bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,35) , 'T' , &
+                                                            ijds , ijde , spec_bdy_width      , &
+                                                            ids , ide , jds , jde , kds , kde , &
+                                                            ims , ime , jms , jme , kms , kme , &
+                                                            ips , ipe , jps , jpe , kps , kpe )
+
 ENDIF
 !!!!!MARS <<
@@ -1101,4 +1734,33 @@
                   q5bdy3dtemp1(i,k,j) = q5bdy3dtemp2(i,k,j)
                   q6bdy3dtemp1(i,k,j) = q6bdy3dtemp2(i,k,j)
+                  q7bdy3dtemp1(i,k,j) = q7bdy3dtemp2(i,k,j)
+                  q8bdy3dtemp1(i,k,j) = q8bdy3dtemp2(i,k,j)
+                  q9bdy3dtemp1(i,k,j) = q9bdy3dtemp2(i,k,j)
+                  q10bdy3dtemp1(i,k,j) = q10bdy3dtemp2(i,k,j)
+                  q11bdy3dtemp1(i,k,j) = q11bdy3dtemp2(i,k,j)
+                  q12bdy3dtemp1(i,k,j) = q12bdy3dtemp2(i,k,j)
+                  q13bdy3dtemp1(i,k,j) = q13bdy3dtemp2(i,k,j)
+                  q14bdy3dtemp1(i,k,j) = q14bdy3dtemp2(i,k,j)
+                  q15bdy3dtemp1(i,k,j) = q15bdy3dtemp2(i,k,j)
+                  q16bdy3dtemp1(i,k,j) = q16bdy3dtemp2(i,k,j)
+                  q17bdy3dtemp1(i,k,j) = q17bdy3dtemp2(i,k,j)
+                  q18bdy3dtemp1(i,k,j) = q18bdy3dtemp2(i,k,j)
+                  q19bdy3dtemp1(i,k,j) = q19bdy3dtemp2(i,k,j)
+                  q20bdy3dtemp1(i,k,j) = q20bdy3dtemp2(i,k,j)
+                  q21bdy3dtemp1(i,k,j) = q21bdy3dtemp2(i,k,j)
+                  q22bdy3dtemp1(i,k,j) = q22bdy3dtemp2(i,k,j)
+                  q23bdy3dtemp1(i,k,j) = q23bdy3dtemp2(i,k,j)
+                  q24bdy3dtemp1(i,k,j) = q24bdy3dtemp2(i,k,j)
+                  q25bdy3dtemp1(i,k,j) = q25bdy3dtemp2(i,k,j)
+                  q26bdy3dtemp1(i,k,j) = q26bdy3dtemp2(i,k,j)
+                  q27bdy3dtemp1(i,k,j) = q27bdy3dtemp2(i,k,j)
+                  q28bdy3dtemp1(i,k,j) = q28bdy3dtemp2(i,k,j)
+                  q29bdy3dtemp1(i,k,j) = q29bdy3dtemp2(i,k,j)
+                  q30bdy3dtemp1(i,k,j) = q30bdy3dtemp2(i,k,j)
+                  q31bdy3dtemp1(i,k,j) = q31bdy3dtemp2(i,k,j)
+                  q32bdy3dtemp1(i,k,j) = q32bdy3dtemp2(i,k,j)
+                  q33bdy3dtemp1(i,k,j) = q33bdy3dtemp2(i,k,j)
+                  q34bdy3dtemp1(i,k,j) = q34bdy3dtemp2(i,k,j)                  
+                  q35bdy3dtemp1(i,k,j) = q35bdy3dtemp2(i,k,j)
 !!!!MARS <<
                END DO
@@ -1192,4 +1854,144 @@
                                                                     ips , ipe , jps , jpe , kps , kpe )
 ENDIF
+
+!VENUS
+IF (config_flags%mars .eq. 34) THEN
+      CALL stuff_bdy     ( qbdy3dtemp1 , grid%scalar_b(:,:,:,:,2)   , 'T' , ijds, ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q2bdy3dtemp1 , grid%scalar_b(:,:,:,:,3)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q3bdy3dtemp1 , grid%scalar_b(:,:,:,:,4)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q4bdy3dtemp1 , grid%scalar_b(:,:,:,:,5)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q5bdy3dtemp1 , grid%scalar_b(:,:,:,:,6)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q6bdy3dtemp1 , grid%scalar_b(:,:,:,:,7)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q7bdy3dtemp1 , grid%scalar_b(:,:,:,:,8)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q8bdy3dtemp1 , grid%scalar_b(:,:,:,:,9)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q9bdy3dtemp1 , grid%scalar_b(:,:,:,:,10)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q10bdy3dtemp1 , grid%scalar_b(:,:,:,:,11)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q11bdy3dtemp1 , grid%scalar_b(:,:,:,:,12)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q12bdy3dtemp1 , grid%scalar_b(:,:,:,:,13)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q13bdy3dtemp1 , grid%scalar_b(:,:,:,:,14)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q14bdy3dtemp1 , grid%scalar_b(:,:,:,:,15)   , 'T' ,ijds , ijde , spec_bdy_width      , &                                                                 
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q15bdy3dtemp1 , grid%scalar_b(:,:,:,:,16)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q16bdy3dtemp1 , grid%scalar_b(:,:,:,:,17)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q17bdy3dtemp1 , grid%scalar_b(:,:,:,:,18)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q18bdy3dtemp1 , grid%scalar_b(:,:,:,:,19)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q19bdy3dtemp1 , grid%scalar_b(:,:,:,:,20)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q20bdy3dtemp1 , grid%scalar_b(:,:,:,:,21)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q21bdy3dtemp1 , grid%scalar_b(:,:,:,:,22)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q22bdy3dtemp1 , grid%scalar_b(:,:,:,:,23)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q23bdy3dtemp1 , grid%scalar_b(:,:,:,:,24)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q24bdy3dtemp1 , grid%scalar_b(:,:,:,:,25)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q25bdy3dtemp1 , grid%scalar_b(:,:,:,:,26)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q26bdy3dtemp1 , grid%scalar_b(:,:,:,:,27)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q27bdy3dtemp1 , grid%scalar_b(:,:,:,:,28)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q28bdy3dtemp1 , grid%scalar_b(:,:,:,:,29)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q29bdy3dtemp1 , grid%scalar_b(:,:,:,:,30)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q30bdy3dtemp1 , grid%scalar_b(:,:,:,:,31)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q31bdy3dtemp1 , grid%scalar_b(:,:,:,:,33)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q32bdy3dtemp1 , grid%scalar_b(:,:,:,:,33)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q33bdy3dtemp1 , grid%scalar_b(:,:,:,:,34)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+      CALL stuff_bdy     ( q34bdy3dtemp1 , grid%scalar_b(:,:,:,:,35)   , 'T' ,ijds , ijde , spec_bdy_width      , &
+                                                                 ids , ide , jds, jde , kds , kde , &
+                                                                 ims , ime , jms, jme , kms , kme , &
+                                                                 ips , ipe , jps, jpe , kps , kpe )
+ENDIF
 !!!!MARS <<
          CALL stuff_bdy     ( mbdy2dtemp1 , grid%em_mu_b    , 'M' , ijds , ijde , spec_bdy_width      , &
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_venus_lmd_new/update_inputs_physiq_mod.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_venus_lmd_new/update_inputs_physiq_mod.F	(revision 2020)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_venus_lmd_new/update_inputs_physiq_mod.F	(revision 2021)
@@ -70,13 +70,51 @@
 SUBROUTINE update_inputs_physiq_tracers(nq,MARS_MODE)
 
-  !use infotrac
-  
+  use infotrac, only: nqtot,tname
+  use logic_mod, only: iflag_trac
   implicit none
   INTEGER, INTENT(IN) :: nq,MARS_MODE
 
-  ALLOCATE(traceurs(nq))
-  traceurs(:)='zolbxs' 
-  !call infotrac_init ! Need traceur.def
-
+  IF (.not.ALLOCATED(tname)) ALLOCATE(tname(nq))
+  nqtot=nq
+
+  IF (MARS_MODE .EQ. 34) THEN
+    tname(1)="co2"
+    tname(2)="co"
+    tname(3)="h2"
+    tname(4)="h2o"
+    tname(5)="o1d"
+    tname(6)="o"
+    tname(7)="o2"
+    tname(8)="o2dg"
+    tname(9)="o3"
+    tname(10)="h"
+    tname(11)="oh"
+    tname(12)="ho2"
+    tname(13)="h2o2"
+    tname(14)="cl"
+    tname(15)="clo"
+    tname(16)="cl2"
+    tname(17)="hcl"
+    tname(18)="hocl"
+    tname(19)="clco"
+    tname(20)="clco3"
+    tname(21)="cocl2"
+    tname(22)="s"
+    tname(23)="so"
+    tname(24)="so2"
+    tname(25)="so3"
+    tname(26)="s2o2"
+    tname(27)="ocs"
+    tname(28)="hso3"
+    tname(29)="h2so4"
+    tname(30)="s2"
+    tname(31)="clso2"
+    tname(32)="oscl"
+    tname(33)="h2oliq"
+    tname(34)="h2so4liq"
+    iflag_trac=1
+  ELSE
+    tname(:)='zolbxs'
+  ENDIF
 
 END SUBROUTINE update_inputs_physiq_tracers
