Re: Dear NCL users

From: Mateus Teixeira <mateus.teixeira_at_nyahnyahspammersnyahnyah>
Date: Sat Nov 21 2009 - 20:45:26 MST

Dear Si Dong,

Did you verify the indexes obtained in lines 24-27?
It seems that they are beyond of dimension limits.
Print these values to see they are inside you data dimension limits.

Best regards,

Mateus

2009/11/22 <sidong7041@hotmail.com>

> Dear NCL users!
>
> Is there something wrong with my ncl file?
> My data is from "NOAA Optimum Interpolation (OI) SST V2".
>
> Si Dong
> I'm using NCL 5.1.1
> -------------------------------------------------------------------
> Error message is :
>
> fatal:Subscript out of range, error in subscript #2
> fatal:An error occurred reading sst_ts
> fatal:Execute: Error occurred at or near line 32 in file sst_jja.ncl
> -------------------------------------------------------------------
> 1
> 2;*********************************************
> 3load "/cygdrive/E/cygwin/usr/local/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> 4load "/cygdrive/E/cygwin/usr/local/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> 5load
> "/cygdrive/E/cygwin/usr/local/lib/ncarg/nclscripts/csm/contributed.ncl"
> 6load "/cygdrive/E/cygwin/usr/local/lib/ncarg/nclscripts/csm/shea_util.ncl"
> 7
> 8begin
> 9;********************************************************
> 10 f = addfile("/cygdrive/G/data/NOAA_sst/OISST/sst.mnmean.nc","r")
> 11;**************************************************************
> 12 date=f->time
> 13 sst=short2flt(f->sst)
> 14
> 15 sst@_FillValue = 32767
> 16
> 17; compute SEASONAL means
> 18 sst_ts = sst(lat|:,lon|:,time|:) ; (lat,lon,time)
> 19 sst_ts = runave(sst_ts, 3, 0) ; (lat,lon,time) 3-mo run ave
> 20
> 21;********************************
> 22; indices for special dates
> 23;*********************************
> 24 ind8901 = ind(date.eq.198201) ; Jan 1982
> 25 ind8912 = ind(date.eq.199912) ; Dec 1999
> 26 ind0001 = ind(date.eq.200001) ; Jan 2000
> 27 ind0012 = ind(date.eq.200812) ; Dec 2008
> 28;********************************
> 29; compute [2-mon] climatologies
> 30; over specified periods
> 31;********************************
> 32 sstAve8299 = clmMonLLT ( sst_ts(:,:,ind8901:ind8912) )
> 33 sstStd8299 = stdMonLLT ( sst_ts(:,:,ind8901:ind8912) )
> 34
> 35 sstAve0008 = clmMonLLT ( sst_ts(:,:,ind0001:ind0012) )
> 36 sstStd0008 = stdMonLLT ( sst_ts(:,:,ind0001:ind0012) )
> ;********************************
> ; compute probabilities for means difference
> ;********************************
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>

-- 
Mateus da Silva Teixeira
Meteorologista
Instituto de Pesquisas Meteorológicas - IPMet
Universidade Estadual Paulista - UNESP
Av. Luis Edmundo Carrijo Coube, 14-01 - Cx. Postal 281 - CEP 17033-360
Bauru - SP - Brasil
Fone: +55 14 3103-6030 / fax: 3203-3649
Registered Linux User #466740 (http://counter.li.org/)

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat Nov 21 20:45:32 2009

This archive was generated by hypermail 2.1.8 : Mon Nov 23 2009 - 12:24:01 MST