Index: /LMDZ5/trunk/libf/phy1d/1DUTILS.h
===================================================================
--- /LMDZ5/trunk/libf/phy1d/1DUTILS.h	(revision 1952)
+++ /LMDZ5/trunk/libf/phy1d/1DUTILS.h	(revision 1953)
@@ -3017,6 +3017,4 @@
 
        open(ilesfile,file='trac.inp.001',status='old',iostat=ierr)
-       nt1=ntrac
-       nt2=ntrac
         if (ierr /= 0) then
             print*,'WARNING : trac.inp does not exist'
Index: /LMDZ5/trunk/libf/phy1d/1D_read_forc_cases.h
===================================================================
--- /LMDZ5/trunk/libf/phy1d/1D_read_forc_cases.h	(revision 1952)
+++ /LMDZ5/trunk/libf/phy1d/1D_read_forc_cases.h	(revision 1953)
@@ -3,4 +3,8 @@
 ! forcing_radconv = .T. : Pure radiative-convective equilibrium:
 !----------------------------------------------------------------------
+
+
+      nq1=0
+      nq2=0
 
       if (forcing_les .or. forcing_radconv
@@ -62,5 +66,5 @@
         ug(l)   = ugprof(kmax)-frac*( ugprof(kmax)- ugprof(kmax-1))
         vg(l)   = vgprof(kmax)-frac*( vgprof(kmax)- vgprof(kmax-1))
-        q(l,nq1:nq2)=qprof(kmax,nq1:nq2)
+        IF (nq2>0) q(l,nq1:nq2)=qprof(kmax,nq1:nq2)
      s               -frac*(qprof(kmax,nq1:nq2)-qprof(kmax-1,nq1:nq2))
         omega(l)=   wfls(kmax)-frac*(   wfls(kmax)-   wfls(kmax-1))
@@ -87,5 +91,5 @@
             ug(l)   = ugprof(k)-frac*( ugprof(k)- ugprof(k-1))
             vg(l)   = vgprof(k)-frac*( vgprof(k)- vgprof(k-1))
-            q(l,nq1:nq2)=qprof(k,nq1:nq2)
+            IF (nq2>0) q(l,nq1:nq2)=qprof(k,nq1:nq2)
      s                   -frac*(qprof(k,nq1:nq2)-qprof(k-1,nq1:nq2))
             omega(l)=   wfls(k)-frac*(   wfls(k)-   wfls(k-1))
@@ -108,5 +112,5 @@
             vg(l)   = vgprof(1)
             omega(l)=   wfls(1)
-            q(l,nq1:nq2)=qprof(1,nq1:nq2)
+            IF (nq2>0) q(l,nq1:nq2)=qprof(1,nq1:nq2)
             dq_dyn(l,1)  =dqtdtls(1)
             dt_cooling(l)=thlpcar(1)