trying to extract a seasonal subset of data

[Katherin Kullgren]

Thank you everyone for all your input.  It helped me look at things
from a different angle.  More than one way to do it will allow me to
check my results!!!
Thanks again!
Katherin







On Aug 15, 2006, at 3:46 PM, Diane Stokes wrote:

> Hi, Katherin.
>
> If you change:
>
> 'set t ' imap pmap
>   to
> 'set t ' imap ' ' pmap
>
>   and, change:
>
> pmap = pmap + 14
>   to
> pmap = pmap + 12
>   or
> pmap = imap + 2
>
>
> the last iteration of your loop would output the fields for:
>    T = 600 to 602
>
> The resulting output file should have 3 fields per year for 50 years:
>   dec1950
>   jan1951
>   feb1951
>   dec1951
>   jan1952
>   feb1952
>     .
>     .
>     .
>   dec1999
>   jan2000
>   feb2000
>
> Is that what you were after?
>
> If you wanted seasonal /averages/ more changes are needed.  (And I see
> Chuck just responded on that topic).
>
>     Diane
>
>
> Katherin Kullgren wrote:
>> Hi,
>>
>> I am trying to extract seasonal subsets of data from my dataset.
>> I  want DJF MAM JJA SON
>>
>> So this is what I tried
>>
>> 'open anom.ctl'
>> 'set gxout fwrite'
>>
>> * This creates the DJF season from Dec1950/Feb1951 - Dec1999/Feb2000
>>
>> 'set fwrite djf.data'
>> imap = 12
>> pmap = 14
>> while (imap <= 600)
>> 'set t ' imap pmap
>> 'd anom'
>> imap = imap + 12
>> pmap = pmap + 14
>> endwhile
>>
>> 'disable fwrite'
>>
>> This does not work. I am not sure the correct way  define the
>> increments that I want.
>> I  want to have
>> set t 12 14
>> set t 24 26
>>  ..........
>> set t 603 605
>>
>> I keep getting the error
>>
>> Data Request Warning:  Request beyond file limits
>> Data Request Warning:  Request beyond file limits
>> Data Request Warning:  Request beyond file limits
>> Data Request Warning:  Request beyond file limits
>> Data Request Warning:  Request beyond file limits
>> Data Request Warning:  Request beyond file limits
>> Data Request Warning:  Request beyond file limits
>> and so on
>>
>> If someone could suggest a way that this would work.  I would
>> appreciate it.
>> BTW- its monthly gridded data.
>>
>> Thanks
>> Katherin
>>
>
> --
> Diane Stokes
> Environmental Modeling Center
> National Weather Service/NOAA

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20060815/b1c9fda7/attachment.html