55 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
56 |
#include "PARAMS.h" |
#include "PARAMS.h" |
57 |
|
|
58 |
integer mythid,im,jm,maxtyp,nchp,nchptot,nchpland,Nsx,Nsy,Nxg,Nyg |
integer mythid,im,jm,maxtyp,nchp,Nsx,Nsy,Nxg,Nyg |
59 |
|
integer nchptot(Nsx,Nsy), nchpland(Nsx,Nsy) |
60 |
integer surftype(im,jm,maxtyp,Nsx,Nsy) |
integer surftype(im,jm,maxtyp,Nsx,Nsy) |
61 |
integer igrd(nchp,Nsx,Nsy),ityp(nchp,Nsx,Nsy) |
integer igrd(nchp,Nsx,Nsy),ityp(nchp,Nsx,Nsy) |
62 |
_RL tilefrac(im,jm,maxtyp,Nsx,Nsy) |
_RL tilefrac(im,jm,maxtyp,Nsx,Nsy) |
171 |
|
|
172 |
c land points |
c land points |
173 |
c ----------- |
c ----------- |
174 |
nchpland = 0 |
nchpland(bi,bj) = 0 |
175 |
do k=1,maxtyp |
do k=1,maxtyp |
176 |
do j=1,jm |
do j=1,jm |
177 |
do i=1,im |
do i=1,im |
178 |
if(surftype(i,j,k,bi,bj).lt.100 .and. |
if(surftype(i,j,k,bi,bj).lt.100 .and. |
179 |
. tilefrac(i,j,k,bi,bj).gt.0.) then |
. tilefrac(i,j,k,bi,bj).gt.0.) then |
180 |
nchpland = nchpland + 1 |
nchpland(bi,bj) = nchpland(bi,bj) + 1 |
181 |
igrd (nchpland,bi,bj) = i + (j-1)*im |
igrd (nchpland(bi,bj),bi,bj) = i + (j-1)*im |
182 |
ityp (nchpland,bi,bj) = surftype(i,j,k,bi,bj) |
ityp (nchpland(bi,bj),bi,bj) = surftype(i,j,k,bi,bj) |
183 |
chfr (nchpland,bi,bj) = tilefrac(i,j,k,bi,bj) |
chfr (nchpland(bi,bj),bi,bj) = tilefrac(i,j,k,bi,bj) |
184 |
chlon(nchpland,bi,bj) = lons(i,j,bi,bj) |
chlon(nchpland(bi,bj),bi,bj) = lons(i,j,bi,bj) |
185 |
chlt (nchpland,bi,bj) = lats(i,j,bi,bj) |
chlt (nchpland(bi,bj),bi,bj) = lats(i,j,bi,bj) |
186 |
endif |
endif |
187 |
enddo |
enddo |
188 |
enddo |
enddo |
190 |
|
|
191 |
c ocean points |
c ocean points |
192 |
c ------------ |
c ------------ |
193 |
nchptot = nchpland |
nchptot(bi,bj) = nchpland(bi,bj) |
194 |
|
|
195 |
do k=1,maxtyp |
do k=1,maxtyp |
196 |
do j=1,jm |
do j=1,jm |
197 |
do i=1,im |
do i=1,im |
198 |
if(surftype(i,j,k,bi,bj).ge.100 .and. |
if(surftype(i,j,k,bi,bj).ge.100 .and. |
199 |
. tilefrac(i,j,k,bi,bj).gt.0.) then |
. tilefrac(i,j,k,bi,bj).gt.0.) then |
200 |
nchptot = nchptot + 1 |
nchptot(bi,bj) = nchptot(bi,bj) + 1 |
201 |
igrd (nchptot,bi,bj) = i + (j-1)*im |
igrd (nchptot(bi,bj),bi,bj) = i + (j-1)*im |
202 |
ityp (nchptot,bi,bj) = surftype(i,j,k,bi,bj) |
ityp (nchptot(bi,bj),bi,bj) = surftype(i,j,k,bi,bj) |
203 |
chfr (nchptot,bi,bj) = tilefrac(i,j,k,bi,bj) |
chfr (nchptot(bi,bj),bi,bj) = tilefrac(i,j,k,bi,bj) |
204 |
chlon(nchptot,bi,bj) = lons(i,j,bi,bj) |
chlon(nchptot(bi,bj),bi,bj) = lons(i,j,bi,bj) |
205 |
chlt (nchptot,bi,bj) = lats(i,j,bi,bj) |
chlt (nchptot(bi,bj),bi,bj) = lats(i,j,bi,bj) |
206 |
endif |
endif |
207 |
enddo |
enddo |
208 |
enddo |
enddo |
209 |
enddo |
enddo |
210 |
|
|
211 |
if(bi.eq.1.and.bj.eq.1)then |
print *,'Number of Total Tiles for bi=',bi,': ',nchptot(bi,bj) |
212 |
print *, 'Number of Total Tiles: ',nchptot |
print *,'Number of Land Tiles for bi=',bi,': ',nchpland(bi,bj) |
213 |
print *, 'Number of Land Tiles: ',nchpland |
|
|
print * |
|
|
endif |
|
|
|
|
214 |
ENDDO |
ENDDO |
215 |
ENDDO |
ENDDO |
216 |
|
|