Hi,
Does anyone have a suggestion how to do this without using a do loop and if
statement, so it's faster?
T1=new ((/24, 848, 70,150/), double)
T1 (:,:,:,:) =-9999.
T1@_FillValue = -9999. ; sets _FillValue to -9999
m= -1
do i=0,15
do j=0,52
m=m+1
do k=0,69
do l=0,149
if(lat(k).ge.lat1(i).and.lat(k).le.lat1(i+1)) then
if(lon(l).ge.lon1(j).and.lon(l).le.lon1(j+1)) then
T1(:,m,k,l)= TLand(:,k,l)
end if
end if
end do
end do
print(m)
end do
end do
printVarSummary(T1)
Thanks
Debasish
----------------------------------------------------------------------------
-----
Debasish Pai Mazumder
Graduate Research Assistant
Geophysical Institute,University of Alaska Fairbanks,
903 Koyukuk Drive
Fairbanks, AK 99775
USA
Phone : +1 907 474 7618
Email : <mailto:debasish_at_gi.alaska.edu> debasish_at_gi.alaska.edu
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jun 06 2007 - 14:30:44 MDT
This archive was generated by hypermail 2.2.0 : Fri Jun 15 2007 - 12:34:56 MDT