51 character(len=16) :: receptor
57 #if defined WITH_CTBTO_PATCHES
63 if (ldirect.lt.0)
then
73 open(unitreceptor,file=path(1)(1:length(1))//
'RECEPTORS', &
76 call
skplin(5,unitreceptor)
84 read(unitreceptor,*,end=99)
85 read(unitreceptor,*,end=99)
86 read(unitreceptor,*,end=99)
87 read(unitreceptor,
'(4x,a16)',end=99) receptor
88 call
skplin(3,unitreceptor)
89 read(unitreceptor,
'(4x,f11.4)',end=99) x
90 call
skplin(3,unitreceptor)
91 read(unitreceptor,
'(4x,f11.4)',end=99) y
92 if ((x.eq.0.).and.(y.eq.0.).and. &
93 (receptor.eq.
' '))
then
97 if (j.gt.maxreceptor)
then
98 write(*,*)
' #### FLEXPART MODEL ERROR! TOO MANY RECEPTOR #### '
99 write(*,*)
' #### POINTS ARE GIVEN. #### '
100 write(*,*)
' #### MAXIMUM NUMBER IS ',maxreceptor,
' #### '
101 write(*,*)
' #### PLEASE MAKE CHANGES IN FILE RECEPTORS #### '
103 receptorname(j)=receptor
104 xreceptor(j)=(x-xlon0)/dx
105 yreceptor(j)=(y-ylat0)/dy
106 xm=r_earth*cos(y*pi/180.)*dx/180.*pi
107 ym=r_earth*dy/180.*pi
108 receptorarea(j)=xm*ym
116 999
write(*,*)
' #### FLEXPART MODEL ERROR! FILE "RECEPTORS" #### '
117 write(*,*)
' #### CANNOT BE OPENED IN THE DIRECTORY #### '
118 write(*,
'(a)') path(1)(1:length(1))
subroutine skplin(nlines, iunit)