--- trunk/libf/dyn3d/comvert.f90	2011/07/01 15:00:48	47
+++ trunk/libf/dyn3d/comvert.f90	2012/01/30 12:54:02	57
@@ -27,7 +27,8 @@
     ! "pa" should be defined before this procedure is called.
 
     use nr_util, only: pi
-    use l_util, only: new_unit
+    use jumble, only: new_unit
+    use unit_nml_m, only: unit_nml
 
     REAL s(llm+1)
     ! "s(l)" is the atmospheric hybrid sigma-pressure coordinate at
@@ -37,7 +38,7 @@
     ! "ds(l)" : épaisseur de la couche "l" dans la coordonnée "s"
 
     INTEGER l, unit
-    REAL alpha, x(llm), trash
+    REAL alpha, x(llm)
 
     character(len=7):: s_sampling = "LMD5"
     ! (other allowed values are "param", "strato1", "strato2" and "read")
@@ -66,7 +67,7 @@
 
     print *, "Enter namelist 'disvert_nml'."
     read(unit=*, nml=disvert_nml)
-    write(unit=*, nml=disvert_nml)
+    write(unit_nml, nml=disvert_nml)
 
     select case (s_sampling)
     case ("param")
@@ -120,7 +121,7 @@
             position="rewind")
        read(unit, fmt=*) ! skip title line
        do l = 1, llm + 1
-          read(unit, fmt=*) trash, s(l)
+          read(unit, fmt=*) s(l)
        end do
        close(unit)
        ! Quick check: