apply read parameter subroutine

This commit is contained in:
ignis 2019-09-15 03:08:21 +09:00
parent d6e2ad74e5
commit 2ca691010f

View file

@ -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