research/incomp-flame-post
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

move subroutine read_intro to m_parameters

Browse source
This commit is contained in:
ignis 2019-04-12 14:08:44 +09:00
parent d20a384d06
commit 7a060505f3
2 changed files with 68 additions and 65 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

68
code/m_parameters.f90
View file

@ -41,4 +41,72 @@ logical function to_omit(num)
end function to_omit
subroutine read_intro
character(LEN=10) :: cdum
integer :: ierr,i
open (100, FILE='post-edge-cold-bc-hybrid-intro')
open (101, FILE='otape')
read (100,*) cdum,l_0
write(101,*) cdum,l_0
read (100,*) cdum,startnum
write(101,*) cdum,startnum
read (100,*) cdum,endnum
write(101,*) cdum,endnum
read (100,*) cdum,skipnum
write(101,*) cdum,skipnum
read (100,*) cdum,shiftnum
write(101,*) cdum,shiftnum
read (100,*) cdum,vis0p
write(101,*) cdum,vis0p
read (100,*) cdum,scp
write(101,*) cdum,scp
read (100,*) cdum,lep
write(101,*) cdum,lep
read (100,*) cdum,min_wr
write(101,*) cdum,min_wr
read (100,*) cdum,prof_wr
write(101,*) cdum,prof_wr
read (100,*) cdum,min_fsd
write(101,*) cdum,min_fsd
read (100,*) cdum,min_c
write(101,*) cdum,min_c
read (100,*) cdum,pre
write(101,*) cdum,pre
read (100,*) cdum,ac
write(101,*) cdum,ac
read (100,*) cdum,bc
write(101,*) cdum,bc
read (100,*) cdum,c_cut
write(101,*) cdum,c_cut
read (100,*) cdum,c_ref
write(101,*) cdum,c_ref
read (100,*) cdum,syp
write(101,*) cdum,syp
read (100,*) cdum,eyp
write(101,*) cdum,eyp
read (100,*) cdum,SL_u
write(101,*) cdum,SL_u
read (100,*) cdum,omitnum
write(101,*) cdum,omitnum
if (omitnum.gt.0) THEN
allocate (omit_t(omitnum,2), stat=ierr)
do i=1,omitnum
read (100,*) cdum,omit_t(i,1) ! begining fullsave number
write(101,*) cdum,omit_t(i,1)
read (100,*) cdum,omit_t(i,2) ! ending fullsave number
write(101,*) cdum,omit_t(i,2)
enddo
endif
prp=scp/lep
close (100)
close (101)
end subroutine read_intro
end module m_parameters

65
code/post.f90
View file

@ -118,71 +118,6 @@
! End of main routine
!========================================================================================
SUBROUTINE READ_INTRO
CHARACTER(LEN=10) :: cdum
INTEGER :: ierr,i
OPEN (100,FILE='post-edge-cold-bc-hybrid-intro')
OPEN (101,FILE='otape')
READ (100,*) cdum,l_0
WRITE(101,*) cdum,l_0
READ (100,*) cdum,startnum
WRITE(101,*) cdum,startnum
READ (100,*) cdum,endnum
WRITE(101,*) cdum,endnum
READ (100,*) cdum,skipnum
WRITE(101,*) cdum,skipnum
READ (100,*) cdum,shiftnum
WRITE(101,*) cdum,shiftnum
READ (100,*) cdum,vis0p
WRITE(101,*) cdum,vis0p
READ (100,*) cdum,scp
WRITE(101,*) cdum,scp
READ (100,*) cdum,lep
WRITE(101,*) cdum,lep
READ (100,*) cdum,min_wr
WRITE(101,*) cdum,min_wr
READ (100,*) cdum,prof_wr
WRITE(101,*) cdum,prof_wr
READ (100,*) cdum,min_fsd
WRITE(101,*) cdum,min_fsd
READ (100,*) cdum,min_c
WRITE(101,*) cdum,min_c
READ (100,*) cdum,pre
WRITE(101,*) cdum,pre
READ (100,*) cdum,ac
WRITE(101,*) cdum,ac
READ (100,*) cdum,bc
WRITE(101,*) cdum,bc
READ (100,*) cdum,c_cut
WRITE(101,*) cdum,c_cut
READ (100,*) cdum,c_ref
WRITE(101,*) cdum,c_ref
READ (100,*) cdum,syp
WRITE(101,*) cdum,syp
READ (100,*) cdum,eyp
WRITE(101,*) cdum,eyp
READ (100,*) cdum,SL_u
WRITE(101,*) cdum,SL_u
READ (100,*) cdum,omitnum
WRITE(101,*) cdum,omitnum
IF (omitnum.gt.0) THEN
ALLOCATE(omit_t(omitnum,2),STAT=ierr)
DO i=1,omitnum
READ (100,*) cdum,omit_t(i,1) ! begining fullsave number
WRITE(101,*) cdum,omit_t(i,1)
READ (100,*) cdum,omit_t(i,2) ! ending fullsave number
WRITE(101,*) cdum,omit_t(i,2)
ENDDO
ENDIF
prp=scp/lep
CLOSE (100)
CLOSE (101)
END SUBROUTINE READ_INTRO
SUBROUTINE PRINT_BANNER
WRITE(*,*) ' This program, x-edge-cold-bc-5-hybrid, is written by D. Kim, 2018'
WRITE(*,*) ' It is to study the statistics of the flame parameters at the leading edge'