Commit 00decbe2 authored by dingwell's avatar dingwell
Browse files

Added debug information

parent 7eeb484a
......@@ -83,6 +83,10 @@
write(*,*) ' NOTICE: RELEASE POINT OUT OF DOMAIN DETECTED.'
write(*,*) ' IT IS REMOVED NOW ... '
write(*,*) ' DEBUG INFO:'
write(*,*) ' Point number:',i
write(*,*) ' xpoint1,ypoint1:',xpoint1,ypoint1
write(*,*) ' xpoint2,ypoint2:',xpoint2,ypoint2
! write(*,*) ' COMMENT: ',compoint(i)
if (i.lt.numpoint) then
......
......@@ -92,7 +92,7 @@ subroutine readinput
! print*,'path2'
! print*,'path',unitpath
! open(unitpath,file='flexwrf.input',status='old',err=801)
write(*,*) "Opening ",inputname," for reading"
write(*,*) "Opening '",inputname,"' for reading"
open(unitpath,file=inputname,status='old',err=801)
! jdf start read pathnames
......@@ -234,6 +234,8 @@ subroutine readinput
read(unitcommand,*) iouttype
read(unitcommand,*) ncnumrec
read(unitcommand,*) option_verbose
write(*,*) "option_verbose=",option_verbose
option_verbose=10
ifine=max(ifine,1)
......@@ -735,6 +737,7 @@ subroutine readinput
! within the modelling period.
!*********************************************************************
write(*,*) 'Opening file: ',path(3)(1:length(3)),' for reading'
open(unitavailab,file=path(3)(1:length(3)),status='old', err=804)
do i=1,3
......@@ -769,6 +772,7 @@ subroutine readinput
!*********************************************************************
do k=1,numbnests
write(*,*) 'Opening file: ',path(numpath+2*k)(1:length(numpath+2*k)),' for reading'
open(unitavailab,file=path(numpath+2*(k-1)+2) &
(1:length(numpath+2*(k-1)+2)),status='old',err=803)
......@@ -923,7 +927,9 @@ subroutine readinput
memtime(i)=999999999
enddo
if(option_verbose.ge.10) write(*,*) "readinput() calling gridcheck()"
call gridcheck()
if(option_verbose.ge.10) write(*,*) "readinput() calling gridcheck_nests()"
call gridcheck_nests()
! jdf end read available
......@@ -965,6 +971,7 @@ subroutine readinput
! Open the OUTGRID file and read output grid specifications
!**********************************************************
write(*,*) "Reading OUTGRID"
call skplin(1,unitpath)
......@@ -1152,6 +1159,7 @@ subroutine readinput
youtshift=ymet0-out_ym0
if(nested_output .eq. 1) then
write(*,*) "Reading OUTGRID nest"
call skplin(1,unitpath)
!read(unitpath,'(f15.8)') xtmp1
!read(unitpath,'(f15.8)') ytmp1
......@@ -1288,6 +1296,7 @@ subroutine readinput
! *
!*******************************************************************************
write(*,*) "Reading RECEPTORS"
call skplin(1,unitpath)
read(unitpath,*) numreceptor
! For backward runs, do not allow receptor output. Thus, set number of receptors to zero
......@@ -1374,6 +1383,7 @@ subroutine readinput
! Open the SPECIES file and read species names and properties
!************************************************************
write(*,*) 'Reading SPECIES'
read(unitpath,*)
read(unitpath,*) numtable
read(unitpath,*)
......@@ -1407,6 +1417,7 @@ subroutine readinput
endif
endif
write(*,*) "option_verbose:",option_verbose
if (partsig(i).eq.1.) partsig(i)=1.0001 ! avoid realing exception
if (partsig(i).eq.0.) partsig(i)=1.0001 ! avoid realing exception
......@@ -1481,6 +1492,7 @@ subroutine readinput
! or using a formatted mask)
! Use of formatted mask is assumed if line 10 contains the word 'DIRECTION'
!**************************************************************************
write(*,*) 'Reading RELEASES'
call skplin(1,unitpath)
! Read the number of species and the link to the species information table
......@@ -1699,6 +1711,7 @@ subroutine readinput
! if (numpoint.gt.maxpoint) goto 997
if(option_verbose.ge.1) write(*,*) 'numpoint =',numpoint
do j=1,numpoint
read(unitpath,*) id1,it1
read(unitpath,*) id2,it2
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment