121 |
tmpij(i,j) = 0. |
tmpij(i,j) = 0. |
122 |
enddo |
enddo |
123 |
enddo |
enddo |
124 |
do i = 1,nchptot(bi,bj) |
do i = 1,nchpland(bi,bj) |
125 |
tmpchp(i) = tcanopy(i,bi,bj) |
tmpchp(i) = tcanopy(i,bi,bj) |
126 |
enddo |
enddo |
127 |
call msc2grd(igrd(1,bi,bj),chfr(1,bi,bj),tmpchp, |
call msc2grd(igrd(1,bi,bj),chfr(1,bi,bj),tmpchp, |
128 |
. nchp,nchptot(bi,bj),fracl,tmpij,im2,jm2) |
. nchp,nchpland(bi,bj),fracl,tmpij,im2,jm2) |
129 |
do j = jm1,jm2 |
do j = jm1,jm2 |
130 |
do i = im1,im2 |
do i = im1,im2 |
131 |
tgz(i,j,bi,bj) = tmpij(i,j) |
tgz(i,j,bi,bj) = tmpij(i,j) |