53 real,
parameter :: eps=1.e-4
60 open(unitoutgrid,file=path(1)(1:length(1))//
'OUTGRID_NEST', &
71 read(unitoutgrid,
'(4x,f11.4)') outlon0n
73 read(unitoutgrid,
'(4x,f11.4)') outlat0n
75 read(unitoutgrid,
'(4x,i5)') numxgridn
77 read(unitoutgrid,
'(4x,i5)') numygridn
79 read(unitoutgrid,
'(4x,f12.5)') dxoutn
81 read(unitoutgrid,
'(4x,f12.5)') dyoutn
84 allocate(orooutn(0:numxgridn-1,0:numygridn-1) &
86 if (stat.ne.0)
write(*,*)
'ERROR: could not allocate outh'
87 allocate(arean(0:numxgridn-1,0:numygridn-1) &
89 if (stat.ne.0)
write(*,*)
'ERROR: could not allocate outh'
90 allocate(volumen(0:numxgridn-1,0:numygridn-1,numzgrid) &
92 if (stat.ne.0)
write(*,*)
'ERROR: could not allocate outh'
97 xr=outlon0n+
real(numxgridn)*dxoutn
98 yr=outlat0n+
real(numygridn)*dyoutn
99 xr1=xlon0+
real(nxmin1)*dx
100 yr1=ylat0+
real(nymin1)*dy
101 if ((outlon0n+eps.lt.xlon0).or.(outlat0n+eps.lt.ylat0) &
102 .or.(xr.gt.xr1+eps).or.(yr.gt.yr1+eps))
then
103 write(*,*)
' #### FLEXPART MODEL ERROR! PART OF OUTPUT ####'
104 write(*,*)
' #### NEST IS OUTSIDE MODEL DOMAIN. CHANGE ####'
105 write(*,*)
' #### FILE OUTGRID IN DIRECTORY ####'
106 write(*,
'(a)') path(1)(1:length(1))
110 xoutshiftn=xlon0-outlon0n
111 youtshiftn=ylat0-outlat0n
116 999
write(*,*)
' #### FLEXPART MODEL ERROR! FILE OUTGRID_NEST #### '
117 write(*,*)
' #### CANNOT BE OPENED IN THE DIRECTORY #### '
118 write(*,*)
' #### xxx/flexpart/options #### '
subroutine skplin(nlines, iunit)
subroutine readoutgrid_nest