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
90 open(unitland,file=path(1)(1:length(1)) &
91 //
'IGBP_int1.dat',status=
'old', &
93 form=
'UNFORMATTED', err=998, convert=
'little_endian')
95 read (unitland) (ilr_buffer(i),i=1,2160000)
106 ilr=ilr_buffer(irecread)
125 landinvent(ix,jy,k)=lu_cat
126 landinvent(ix,jy,k+3)=lu_perc
135 open(unitsurfdata,file=path(1)(1:length(1))//
'surfdata.t', &
136 status=
'old',err=999)
142 read(unitsurfdata,
'(45x,f15.3)') z0(i)
151 998
write(*,*)
' #### FLEXPART ERROR! FILE CONTAINING ####'
152 write(*,*)
' #### LANDUSE INVENTORY DOES NOT EXIST ####'
155 999
write(*,*)
' #### FLEXPART ERROR! FILE CONTAINING ####'
156 write(*,*)
' #### RELATION LANDUSE,z0 DOES NOT EXIST ####'