masking an area of china, there is a error or a bug.

Hello all,
  
    I try to mask an area of china with the shp file, but something goes wrong, it showed as:
  
     "fatal:Subscript out of range, error in subscript #0
     fatal:An error occurred reading lat1d
     fatal:Execute: Error occurred at or near line 73 in file mask_9_cn.ncl"
  
     David Brown told me that,
     " the problem starts with the line: ilt_mn = ind(min_mrb_lat.gt.lat1d)
        The variable ilt_mn gets set to a missing value because there are no values of the array lat1d where the value of min_mrb_lat is greater. Therefore lit_mn gets set to a missing value and therefore causes an out-of-range subscript in the line that reports the error:
                  if(gc_inout(lat1d(ilt),lon1d(iln),mrb_lat,mrb_lon)) then "
  
     I guest that the problem is from the shp file, Maybe there are two or more areas other than one area. so the function gc_inout() connot do it. I am trying to make a new shp file now, but I'm not sure this will work then.

  
     Is there anyone would like to help me?
     
     Thanks!
     
     Dazhen
  
  ------------------
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 Dazhen HUANG, Mr. #(»Æ´óÕê)
 Department of Atmospheric Sciences,
 Sun Yat-sen University,
 No.135, Xingang West Road. Guangzhou, 510275, P.R.China
 Email: huangdzh@foxmail.com
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk