Index: trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F
===================================================================
--- trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F	(revision 1801)
+++ trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F	(revision 1802)
@@ -6,4 +6,5 @@
       use mod_grid_phy_lmdz, only : regular_lonlat
       use infotrac, only: nqtot, tname
+      use tracer_h, only: noms
       use surfdat_h, only: albedodat, phisfi, dryness, watercaptag,
      &                     zmea, zstd, zsig, zgam, zthe,
@@ -23,5 +24,5 @@
      &                            nday, iphysiq
       use callkeys_mod, only: tracer,check_cpp_match,rings_shadow,
-     &	                      specOLR,water,pceil,ok_slab_ocean
+     &	                      specOLR,water,pceil,ok_slab_ocean,photochem
       USE comvert_mod, ONLY: ap,bp,aps,bps,pa,preff, sig,
      &                       presnivs,pseudoalt,scaleheight
@@ -138,4 +139,5 @@
 !     added by AS to avoid the use of adv trac common
       character*20,allocatable :: nametrac(:)   ! name of the tracer (no need for adv trac common)
+      character*20,allocatable :: nametmp(:)   !
 
       real :: latitude(1), longitude(1), cell_area(1)
@@ -254,4 +256,5 @@
           if (nq>0) then
             allocate(tname(nq))
+            allocate(noms(nq))
             allocate(q(llm,nq))
             allocate(qsurf(nq))
@@ -267,4 +270,5 @@
           ! minimal version, just read in the tracer names, 1 per line
             read(90,*,iostat=ierr) tname(iq)
+            noms(iq)=tname(iq)
             if (ierr.ne.0) then
               write(*,*) 'rcm1d: error reading tracer names...'
@@ -292,4 +296,5 @@
         endif
         close(90)
+
 
       else ! of if (tracer)
@@ -816,8 +821,8 @@
 ! -----------------
       do isoil=0,nsoil-1
-        mlayer(isoil)=2.e-4*(2.**(isoil-0.5)) ! mid-layer depth
+        mlayer(isoil)=3.e-2*(2.**(isoil-0.5)) ! mid-layer depth
       enddo
       do isoil=1,nsoil
-        layer(isoil)=2.e-4*(2.**(isoil-1)) ! layer depth
+        layer(isoil)=3.e-2*(2.**(isoil-1)) ! layer depth
       enddo
 
@@ -838,4 +843,15 @@
       endif
 
+
+! Initialize chemical species
+! -----------------
+      if(tracer.and.photochem) then
+           call initracer(1,nq,tname)
+           allocate(nametmp(nq))
+           nametmp(1:nq)=tname(1:nq)
+           call inichim_1D(nq, q, qsurf, psurf, 0, 0)
+           tname(1:nq)=nametmp(1:nq)
+           noms(1:nq)=nametmp(1:nq)
+      endif ! tracer and photochem