diff --git a/code/m_parameters.f90 b/code/m_parameters.f90 index 18c6d50..f4d024f 100644 --- a/code/m_parameters.f90 +++ b/code/m_parameters.f90 @@ -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 diff --git a/code/post.f90 b/code/post.f90 index 8064bf9..706c745 100644 --- a/code/post.f90 +++ b/code/post.f90 @@ -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'