Index: LMDZ.3.3/trunk/libf/dyn3d/nudge.F
===================================================================
--- LMDZ.3.3/trunk/libf/dyn3d/nudge.F	(revision 346)
+++ LMDZ.3.3/trunk/libf/dyn3d/nudge.F	(revision 356)
@@ -26,5 +26,5 @@
 
 c   ...  Possibilite de choisir le shema de Van-leer pour l'advection de
-c         q  , en faisant iadv = 3  dans   traceur  (29/04/97) .
+c         q  , en faisant iadv = 10  dans   traceur  (29/04/97) .
 c
 c-----------------------------------------------------------------------
@@ -70,10 +70,14 @@
 c
       INTEGER itau,ij,l
+      integer ncidt,varidpl,nlev,status
+      integer rcod,rid 
       real ditau,tau,a
-
-      LOGICAL first
+      save nlev
+
+      LOGICAL first,ncep
       integer online
-      save first,online
+      save first,online,ncep
       data first,online/.true.,1/
+      data ncep/.false./
 
       save ucovrea1,vcovrea1,tetarea1,masserea1,psrea1
@@ -139,9 +143,17 @@
 
 c    itau_test    montre si l'importation a deja ete faite au rang itau
-
-
+c lecture d'un fichier netcdf pour determiner le nombre de niveaux
+         ncidt=NCOPN('T.nc',NCNOWRIT,rcod)
+         if (ncep) then
+          status=NF_INQ_DIMID(ncidt,'LEVEL',rid)
+         else
+          status=NF_INQ_DIMID(ncidt,'PRESSURE',rid)
+         endif
+          status=NF_INQ_DIMLEN(ncidt,rid,nlev)
+         print *,'nlev', nlev 
+          call ncclos(ncidt,rcod)
 c   Lecture du premier etat des reanalyses.
          call read_reanalyse(1
-     s   ,ucovrea2,vcovrea2,tetarea2,masserea2,ps,1)
+     s   ,ucovrea2,vcovrea2,tetarea2,masserea2,ps,1,nlev)
 
 c-----------------------------------------------------------------------
@@ -180,5 +192,5 @@
            itau_test=itau
            call read_reanalyse(step_rea
-     s     ,ucovrea2,vcovrea2,tetarea2,masserea2,ps,1)
+     s     ,ucovrea2,vcovrea2,tetarea2,masserea2,ps,1,nlev)
         endif
       else
Index: LMDZ.3.3/trunk/libf/dyn3d/read_reanalyse.F
===================================================================
--- LMDZ.3.3/trunk/libf/dyn3d/read_reanalyse.F	(revision 346)
+++ LMDZ.3.3/trunk/libf/dyn3d/read_reanalyse.F	(revision 356)
@@ -2,5 +2,5 @@
 c $Header$
 c
-      subroutine read_reanalyse(timestep,u,v,t,masse,ps,mode)
+      subroutine read_reanalyse(timestep,u,v,t,masse,ps,mode,nlevnc)
 
 c   mode=0 variables naturelles
@@ -15,6 +15,7 @@
 
       integer nlevnc
-      parameter (nlevnc=15)
-C pour annee 2000      parameter (nlevnc=21)
+c      parameter (nlevnc=15)
+C pour annee 2000      
+c      parameter (nlevnc=21)
       integer timestep,mode,l