Dimension mismatch on latitude variable ERA data extraction

Hi NCL-Talk

I am extracting data from an ERA data set and I get the following error messaage;

fatal:Dimension size mismatch on subscript #1, left-hand and right-hand side dimensions do not match

fatal:["Execute.c":7556]:Execute: Error occurred at or near line 185 in file /ECMWF_INTERIM/era_interim_spell_tot_0000hrs_run.ncl

The line it is refering to is;

spell_tot_unique_mn(nr,:,:) = dim_sum_n(cold_mn1({spell(nr)},:,:),0)

which is part of a do loop-see below (full script is attached)

do nr=0,nspell-1

spell_tot_unique_mn(nr,:,:) = dim_sum_n(cold_mn1({spell(nr)},:,:),0)

end do

I attached units and names to each part of both of the variables and for the dimension#1 i.e. the latitude they appear to be the same, so I am not sure what I am doing wrong in the script. The two variables are listed below (a summary of all the variables is also attached)

Any help on what I am doind wrong would be much appreciated.

Variable: cold_mn1

Type: integer

Total Size: 67145680 bytes

16786420 values

Number of Dimensions: 3

Dimensions and sizes: [run_length | 12235] x [lat | 28] x [lon | 49]

Coordinates:

run_length: [1..12235]

lat: [73.5..33]

lon: [-27..45]

Number Of Attributes: 2

cold_crit_values : 0

long_name : unique cold runs

$B!!(B

Variable: spell_tot_unique_mn

Type: float

Total Size: 548800 bytes

137200 values

Number of Dimensions: 3

Dimensions and sizes: [spell | 100] x [lat | 28] x [lon | 49]

Coordinates:

spell: [1..100]

lat: [73.5..33]

lon: [-27..45]

Number Of Attributes: 2

spell : 1

_FillValue : 9.96921e+36

(0) 100

(0)

(0) min=9.96921e+36 max=9.96921e+36

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