Re: Script to tally days below a given temperature - DO LOOP issue

Hi Gus

Yes, after I changed my do loop to

do iensemble = 0,8,1 ; ENSEMBLE LOOP MODELS 0,1...7,8

it then picked up the itime v it issue you mentioned so I changed
t2_uk = wrf_user_getvar(a,"ukmo",it) ; T2 in Kelvin to the following;

t2_uk = wrf_user_getvar(a[:],"ukmo",itime) ; T2 in Kelvin

but then I received the following error message.

fatal:_NclBuildArray: Arrays of files are not yet supported
fatal:Execute: Error occurred at or near line 72 in file plt_min2m_allmodels_allensembles_all_days.ncl

reference the ensemble member counter-is that related to the position on the end do statement, and if so should I have placed it before the tcrit line?

Kind regards

Dominique

________________________________________
From: ncl-talk-bounces@ucar.edu [ncl-talk-bounces@ucar.edu] on behalf of Gus Correa [gus@ldeo.columbia.edu]
Sent: 04 June 2012 20:43
Cc: NCL-talk
Subject: Re: Script to tally days below a given temperature - DO LOOP issue

Hi

You declared

ensemble = a[:]->number

which is not a scalar, then tried to use it as a loop counter:

    do ensemble = 0,9,1 ; ENSEMBLE LOOP MODELS 0,1...7,8

Have you tried to name the loop variable differently,
say:

    do emb = 0,9,1 ; ENSEMBLE LOOP MODELS 0,1...7,8

In addition, there may be more loose threads:
loop counter 'itime' vs. 'it' which appears inside the loops,
and the ensemble member counter doesn't seem to appear
inside the loop nest.

Gus Correa

On 06/04/2012 02:56 PM, Cannavina, Dominique wrote:
>
> Hi
>
> I am writing a script to count the number of occurances for given number
> of consecutive days that a temperature falls below a given temperature.
> The data I am am using is from the DEMETER project and each file
> contains 3 different models each running 9 ensembles to give daily
> minimum temperatures. The domain covering Europe.
>
> WhenI run the attached script I get the following error in the DO LOOP
>
> fatal:Loop variable must be scalar, can't execute loop
>
> I am not sure how to work around this as I need to loop through the days
> and ensemble runs for each model?
>
> Any help would be much appreciated.
>
> Kind regards
>
> Dominique
>
>
>
> _______________________________________________
> 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
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jun 4 14:04:46 2012