diff --git a/code/m_parameters.f90 b/code/m_parameters.f90
index a9afe53..758e375 100644
--- a/code/m_parameters.f90
+++ b/code/m_parameters.f90
@@ -84,60 +84,31 @@ contains
        END IF
 
        OPEN(otape,FILE='otape')
-        !READ(itape,*) cdum,nx
-        !WRITE(otape,*) cdum,nx
         CALL READ_PARAMETER (nx)
-
-        READ(itape,*) cdum,l_0
-        WRITE(otape,*) cdum,l_0
-        READ(itape,*) cdum,int_pr
-        WRITE(otape,*) cdum,int_pr
-        READ(itape,*) cdum,tar_lo
-        WRITE(otape,*) cdum,tar_lo
-        READ(itape,*) cdum,dt
-        WRITE(otape,*) cdum,dt
-        READ(itape,*) cdum,sc
-        WRITE(otape,*) cdum,sc
-        READ(itape,*) cdum,vis
-        WRITE(otape,*) cdum,vis
-        READ(itape,*) cdum,pre
-        WRITE(otape,*) cdum,pre
-        READ(itape,*) cdum,ac
-        WRITE(otape,*) cdum,ac
-        READ(itape,*) cdum,bc
-        WRITE(otape,*) cdum,bc
-        READ(itape,*) cdum,u0
-        WRITE(otape,*) cdum,u0
-        READ(itape,*) cdum,tf
-        WRITE(otape,*) cdum,tf
-        READ(itape,*) cdum,dt_uf
-        WRITE(otape,*) cdum,dt_uf
-        READ(itape,*) cdum,ctmp
-        WRITE(otape,*) cdum,ctmp
-        READ(itape,*) cdum,c_cut
-        WRITE(otape,*) cdum,c_cut
-!        READ(itape,*) cdum,cs
-!        WRITE(otape,*) cdum,cs
-        READ(itape,*) cdum,c_ref
-        WRITE(otape,*) cdum,c_ref
-        READ(itape,*) cdum,min_wr
-        WRITE(otape,*) cdum,min_wr
-        READ(itape,*) cdum,prof_wr
-        WRITE(otape,*) cdum,prof_wr
-        READ(itape,*) cdum,lewis
-        WRITE(otape,*) cdum,lewis
-        READ(itape,*) cdum,lambda1
-        WRITE(otape,*) cdum,lambda1
-        READ(itape,*) cdum,lambda2
-        WRITE(otape,*) cdum,lambda2
-        READ(itape,*) cdum,beta1
-        WRITE(otape,*) cdum,beta1
-        READ(itape,*) cdum,hrp
-        WRITE(otape,*) cdum,hrp
-        READ(itape,*) cdum,le_a
-        WRITE(otape,*) cdum,le_a
-        READ(itape,*) cdum,le_x
-        WRITE(otape,*) cdum,le_x
+        CALL READ_PARAMETER (l_0)
+        CALL READ_PARAMETER (int_pr)
+        CALL READ_PARAMETER (tar_lo)
+        CALL READ_PARAMETER (dt)
+        CALL READ_PARAMETER (sc)
+        CALL READ_PARAMETER (vis)
+        CALL READ_PARAMETER (pre)
+        CALL READ_PARAMETER (ac)
+        CALL READ_PARAMETER (bc)
+        CALL READ_PARAMETER (u0)
+        CALL READ_PARAMETER (tf)
+        CALL READ_PARAMETER (dt_uf)
+        CALL READ_PARAMETER (ctmp)
+        CALL READ_PARAMETER (c_cut)
+        CALL READ_PARAMETER (c_ref)
+        CALL READ_PARAMETER (min_wr)
+        CALL READ_PARAMETER (prof_wr)
+        CALL READ_PARAMETER (lewis)
+        CALL READ_PARAMETER (lambda1)
+        CALL READ_PARAMETER (lambda2)
+        CALL READ_PARAMETER (beta1)
+        CALL READ_PARAMETER (hrp)
+        CALL READ_PARAMETER (le_a)
+        CALL READ_PARAMETER (le_x)
 
         ! READ(itape,*) cdum,table_size
         ! WRITE(otape,*) cdum,table_size