56 real :: rluh(5,numclass),rgssh(5,numclass),rgsoh(5,numclass)
57 real :: rclsh(5,numclass),rcloh(5,numclass)
63 open(unitwesely,file=path(1)(1:length(1))//
'surfdepo.t', &
71 read(unitwesely,
'(8x,13f8.0)') (ri(i,j),j=1,numclass)
72 read(unitwesely,
'(8x,13f8.0)') (rluh(i,j),j=1,numclass)
73 read(unitwesely,
'(8x,13f8.0)') (rac(i,j),j=1,numclass)
74 read(unitwesely,
'(8x,13f8.0)') (rgssh(i,j),j=1,numclass)
75 read(unitwesely,
'(8x,13f8.0)') (rgsoh(i,j),j=1,numclass)
76 read(unitwesely,
'(8x,13f8.0)') (rclsh(i,j),j=1,numclass)
77 read(unitwesely,
'(8x,13f8.0)') (rcloh(i,j),j=1,numclass)
94 if (ri(i,j).eq.9999.) ri(i,j)=1.e25
95 if (rluh(i,j).eq.9999.) rluh(i,j)=1.e25
96 if (rac(i,j).eq.9999.) rac(i,j)=1.e25
97 if (rgssh(i,j).eq.9999.) rgssh(i,j)=1.e25
98 if (rgsoh(i,j).eq.9999.) rgsoh(i,j)=1.e25
99 if (rclsh(i,j).eq.9999.) rclsh(i,j)=1.e25
100 if (rcloh(i,j).eq.9999.) rcloh(i,j)=1.e25
108 ri(i,j)=max(ri(i,j),0.001)
109 rluh(i,j)=max(rluh(i,j),0.001)
110 rac(i,j)=max(rac(i,j),0.001)
111 rgssh(i,j)=max(rgssh(i,j),0.001)
112 rgsoh(i,j)=max(rgsoh(i,j),0.001)
113 rclsh(i,j)=max(rclsh(i,j),0.001)
114 rcloh(i,j)=max(rcloh(i,j),0.001)
124 if (reldiff(ic).gt.0.)
then
127 rlu(ic,i,j)=rluh(i,j)/(1.e-5*henry(ic)+f0(ic))
128 rgs(ic,i,j)=1./(henry(ic)/(10.e5*rgssh(i,j))+f0(ic)/ &
130 rcl(ic,i,j)=1./(henry(ic)/(10.e5*rclsh(i,j))+f0(ic)/ &
141 999
write(*,*)
'### FLEXPART ERROR! FILE ###'
142 write(*,*)
'### surfdepo.t DOES NOT EXIST. ###'