71 integer :: ix,jy,i,k,lu_cat,lu_perc
72 integer(kind=1) :: ilr
73 integer(kind=1) :: ilr_buffer(2160000)
74 integer :: il,irecread
91 write (*,*)
'reading: ',path(1)(1:length(1))
92 open(unitland,file=path(1)(1:length(1)) &
93 //
'IGBP_int1.dat',status=
'old', &
94 form=
'UNFORMATTED', err=998)
95 read (unitland) (ilr_buffer(i),i=1,2160000)
97 write (*,*)
'reading: '
107 ilr=ilr_buffer(irecread)
125 landinvent(ix,jy,k)=lu_cat
126 landinvent(ix,jy,k+3)=lu_perc
127 if ((jy.lt.10).and.(ix.lt.10))
then
128 write(*,*)
'reading: ', ix, jy, lu_cat, lu_perc
137 open(unitsurfdata,file=path(1)(1:length(1))//
'surfdata.t', &
138 status=
'old',err=999)
144 read(unitsurfdata,
'(45x,f15.3)') z0(i)
153 998
write(*,*)
' #### FLEXPART ERROR! FILE CONTAINING ####'
154 write(*,*)
' #### LANDUSE INVENTORY DOES NOT EXIST ####'
157 999
write(*,*)
' #### FLEXPART ERROR! FILE CONTAINING ####'
158 write(*,*)
' #### RELATION LANDUSE,z0 DOES NOT EXIST ####'