Index: trunk/LMDZ.MARS/libf/phymars/physiq_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 2907)
+++ trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 2908)
@@ -1033,5 +1033,5 @@
            CALL callradite(icount,ngrid,nlayer,nq,zday,zls,pq,
      &     albedo_meshavg,emis_meshavg,
-     &     mu0,zplev,zplay,pt,tsurf(:,1),fract,dist_sol,igout,
+     &     mu0,zplev,zplay,pt,tsurf_meshavg,fract,dist_sol,igout,
      &     zdtlw,zdtsw,fluxsurf_lw(:,iflat),fluxsurf_dn_sw(:,:,iflat),
      &     fluxsurf_up_sw,
@@ -1057,5 +1057,5 @@
                CALL callradite(icount,ngrid,nlayer,nq,zday,zls,pq,
      &              albedo_meshavg,emis_meshavg,mu0,zplev,zplay,pt,
-     &              tsurf(:,1),fract,
+     &              tsurf_meshavg,fract,
      &              dist_sol,igout,zdtlwclf,zdtswclf,
      &              fluxsurf_lwclf,fluxsurf_dn_swclf,fluxsurf_up_swclf,
@@ -1467,5 +1467,5 @@
 
           DO ig=1, ngrid
-             IF (zh(ig,1) .lt. tsurf(ig,1)) THEN
+             IF (zh(ig,1) .lt. tsurf_meshavg(ig)) THEN
                wstar(ig)=1.
                hfmax_th(ig)=0.2
@@ -2146,5 +2146,5 @@
       if (callthermos) then
         call thermosphere(ngrid,nlayer,nq,zplev,zplay,dist_sol,
-     $     mu0,ptimestep,ptime,zday,tsurf(:,iflat),zzlev,zzlay,
+     $     mu0,ptimestep,ptime,zday,tsurf_meshavg,zzlev,zzlay,
      &     pt,pq,pu,pv,pdt,pdq,
      $     zdteuv,zdtconduc,zdumolvis,zdvmolvis,zdqmoldiff,