--- trunk/libf/dyn3d/comvert.f90	2011/10/07 13:11:58	53
+++ trunk/libf/dyn3d/comvert.f90	2012/01/30 12:54:02	57
@@ -28,6 +28,7 @@
 
     use nr_util, only: pi
     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
@@ -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")