Index: trunk/LMDZ.MARS/libf/phymars/callradite.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/callradite.F	(revision 342)
+++ trunk/LMDZ.MARS/libf/phymars/callradite.F	(revision 353)
@@ -2,5 +2,5 @@
      $     emis,mu0,pplev,pplay,pt,tsurf,fract,dist_sol,igout,
      $     dtlw,dtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw,fluxtop_sw,
-     &     tauref,tau,aerosol,ccn,rdust,rice,nuice)
+     &     tauref,tau,aerosol,ccn,rdust,rice,nuice,co2ice)
 
        IMPLICIT NONE
@@ -178,4 +178,6 @@
       REAL rice(ngridmx,nlayermx)   ! Ice geometric mean radius (m)
       REAL nuice(ngridmx,nlayermx)  ! Estimated effective variance
+      REAL co2ice(ngridmx)           ! co2 ice surface layer (kg.m-2)
+
 c
 c    Local variables :
@@ -491,5 +493,6 @@
      .        ,fluxsurf_lw(ig0+1),fluxtop_lw(ig0+1)
      .        ,znetrad
-     &        ,zQIRsQREF3d,zomegaIR3d,zgIR3d)
+     &        ,zQIRsQREF3d,zomegaIR3d,zgIR3d
+     &        ,co2ice(ig0+1))
 
 c       Shortwave ("sw") radiative transfer (= solar radiation)
Index: trunk/LMDZ.MARS/libf/phymars/lwmain.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/lwmain.F	(revision 342)
+++ trunk/LMDZ.MARS/libf/phymars/lwmain.F	(revision 353)
@@ -4,5 +4,6 @@
      .                   ,fluxground,fluxtop
      .                   ,netrad
-     &                   ,QIRsQREF3d,omegaIR3d,gIR3d)
+     &                   ,QIRsQREF3d,omegaIR3d,gIR3d
+     &                   ,co2ice)
 
 c----------------------------------------------------------------------
@@ -76,4 +77,7 @@
       real newcoolrate(ndlon,nflev) ! cooling rate (K/s) / with implicite scheme
 
+      REAL co2ice(ndlo2)           ! co2 ice surface layer (kg.m-2)
+      REAL emis_gaz(ndlo2)         ! emissivity for gaz computations
+
       integer jk,jkk,ja,jl
 
@@ -103,4 +107,12 @@
       endif
 
+      DO jl=1 , kdlon
+         IF(co2ice(jl) .GT. 20.e-3) THEN
+             emis_gaz(jl)=1. 
+         ELSE
+             emis_gaz(jl)=emis(jl)
+         ENDIF
+      ENDDO
+
 c----------------------------------------------------------------------
 c         1.0   planck function
@@ -129,5 +141,5 @@
 
 c     print*, 'CALL of DISTANTS'
-      call lwxd ( ig0, kdlon, kflev, emis
+      call lwxd ( ig0, kdlon, kflev, emis_gaz
      .          , aer_t, co2_u, co2_up)
 
@@ -148,5 +160,5 @@
 
 c     print*, 'CALL of BOUNDARIES'
-      call lwxb ( ig0, kdlon, kflev, emis
+      call lwxb ( ig0, kdlon, kflev, emis_gaz
      .          , aer_t, co2_u, co2_up)
 
Index: trunk/LMDZ.MARS/libf/phymars/physiq.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/physiq.F	(revision 342)
+++ trunk/LMDZ.MARS/libf/phymars/physiq.F	(revision 353)
@@ -548,5 +548,5 @@
      $     emis,mu0,pplev,pplay,pt,tsurf,fract,dist_sol,igout,
      $     zdtlw,zdtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw,fluxtop_sw,
-     &     tauref,tau,aerosol,ccn,rdust,rice,nuice)
+     &     tauref,tau,aerosol,ccn,rdust,rice,nuice,co2ice)
 
 c          Outputs for basic check (middle of domain)
@@ -1692,5 +1692,5 @@
 
 
-         z_out=1. !0.
+         z_out=0.
          if (calltherm .and. (z_out .gt. 0.)) then
          call surflayer_interpol(ngrid,nlayer,z0,g,zzlay,zu,zv,wmax_th