Re: What is wrong with this "ind" loop for a 2D array?

Thanks, Brian

I tried it, but it still does not work. Same error shows

fatal:Subscript out of range, error in subscript #1
fatal:An error occurred reading fog

I am wondering if the "ind" function only can be used with 1D array.

Thanks a lot.

Shawn




Wen.J.Qu
2012-09-26



·¢¼þÈË£º Brian Medeiros
·¢ËÍʱ¼ä£º 2012-09-26 17:02:11
ÊÕ¼þÈË£º Wen.J.Qu@gmail.com
³­ËÍ£º ncl-talk
Ö÷Ì⣺ Re: [ncl-talk] What is wrong with this "ind" loop for a 2D array?

Shawn,


I would bet that fog never equals 1, so the ind function returns "missing"


Maybe try
if .not. all(ismissing(ii)) then
  FOG = fog( ind(.not.ismissing(ii)), ind(.not.ismissing(jj)) ) ;; no reason to use ii and jj if they are the same.
else
  FOG = -9999 ;; or appropriate missing value
end if


.brian




On Sep 26, 2012, at 3:54 PM, Wen.J.Qu <wen.j.qu@gmail.com> wrote:


Hi, dear NCL users,

Would you please give me some hint of what is wrong with the "ind" loop for a 2D array as shown below?

Thanks a lot.

Shawn


################################################################

When I run it, it shows:
fatal:Subscript out of range, error in subscript #0
fatal:An error occurred reading fog
fatal:Execute: Error occurred at or near line "FOG = fog(ii,jj)" in the file

################################################################

;--fog is a nrows*ncols 2D array.
;--Get indices where fog value is 1, that is fog happens.
    ii = new(nrows,integer)
    jj = new(ncols,integer)
    do i=0,nrows-1
       do j=0,ncols-1
          ii(i) = ind(fog(i,j).eq.1)
          jj(j) = ind(fog(i,j).eq.1)
       end do
    end do

 FOG = fog(ii,jj)
 ;Debug print
 print(FOG)

#################################################################





Wen.J.Qu
2012-09-26
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Sep 27 00:35:41 2012