53 integer :: jjjjmmdd,ihmmss,i,ix,jy,j
54 real :: xp1,yp1,xp2,yp2
61 open(unitheader,file=path(2)(1:length(2))//
'header_nest', &
62 form=
'unformatted',err=998)
68 if (ldirect.eq.1)
then
69 write(unitheader) ibdate,ibtime,
'FLEXPART V8.2'
71 write(unitheader) iedate,ietime,
'FLEXPART V8.2'
77 write(unitheader) loutstep,loutaver,loutsample
82 write(unitheader) outlon0n,outlat0n,numxgridn,numygridn, &
84 write(unitheader) numzgrid,(outheight(i),i=1,numzgrid)
86 call
caldate(bdate,jjjjmmdd,ihmmss)
87 write(unitheader) jjjjmmdd,ihmmss
94 write(unitheader) 3*nspec,maxpointspec_act
96 write(unitheader) 1,
'WD_'//species(i)(1:7)
97 write(unitheader) 1,
'DD_'//species(i)(1:7)
98 write(unitheader) numzgrid,species(i)
105 write(unitheader) numpoint
107 write(unitheader) ireleasestart(i),ireleaseend(i),kindz(i)
108 xp1=xpoint1(i)*dx+xlon0
109 yp1=ypoint1(i)*dy+ylat0
110 xp2=xpoint2(i)*dx+xlon0
111 yp2=ypoint2(i)*dy+ylat0
112 write(unitheader) xp1,yp1,xp2,yp2,zpoint1(i),zpoint2(i)
113 write(unitheader) npart(i),1
114 if (numpoint.le.1000)
then
115 write(unitheader) compoint(i)
117 write(unitheader) compoint(1001)
120 write(unitheader) xmass(i,j)
121 write(unitheader) xmass(i,j)
122 write(unitheader) xmass(i,j)
129 write(unitheader) method,lsubgrid,lconvection, &
130 ind_source,ind_receptor
135 write(unitheader) nageclass,(lage(i),i=1,nageclass)
142 write(unitheader) (orooutn(ix,jy),jy=0,numygridn-1)
149 998
write(*,*)
' #### FLEXPART MODEL ERROR! THE FILE #### '
150 write(*,*)
' #### '//path(2)(1:length(2))//
'header'//
' #### '
151 write(*,*)
' #### CANNOT BE OPENED. IF A FILE WITH THIS #### '
152 write(*,*)
' #### NAME ALREADY EXISTS, DELETE IT AND START #### '
153 write(*,*)
' #### THE PROGRAM AGAIN. #### '
subroutine caldate(juldate, yyyymmdd, hhmiss)