46 real :: xp1,yp1,xp2,yp2
52 open(unitouttraj,file=path(2)(1:length(2))//
'trajectories.txt', &
53 form=
'formatted',err=998)
55 if (ldirect.eq.1)
then
56 write(unitouttraj,
'(i8,1x,i6,1x,a)') ibdate,ibtime,
'FLEXPART V8.2'
58 write(unitouttraj,
'(i8,1x,i6,1x,a)') iedate,ietime,
'FLEXPART V8.2'
60 write(unitouttraj,*) method,lsubgrid,lconvection
61 write(unitouttraj,*) numpoint
63 xp1=xpoint1(i)*dx+xlon0
64 yp1=ypoint1(i)*dy+ylat0
65 xp2=xpoint2(i)*dx+xlon0
66 yp2=ypoint2(i)*dy+ylat0
67 write(unitouttraj,*) ireleasestart(i),ireleaseend(i), &
68 xp1,yp1,xp2,yp2,zpoint1(i),zpoint2(i),kindz(i),npart(i)
69 if (numpoint.le.1000)
then
70 write(unitouttraj,
'(a)') compoint(i)(1:40)
72 write(unitouttraj,
'(a)') compoint(1001)(1:40)
78 998
write(*,*)
' #### FLEXPART MODEL ERROR! THE FILE #### '
79 write(*,*)
' #### trajectories.txt #### '
80 write(*,*)
' #### CANNOT BE OPENED. IF A FILE WITH THIS #### '
81 write(*,*)
' #### NAME ALREADY EXISTS, DELETE IT AND START #### '
82 write(*,*)
' #### THE PROGRAM AGAIN. #### '