[gradsusr] Concatenation on Windows OS

ugo okoro uknac23 at yahoo.com
Tue Oct 3 05:04:39 EDT 2017


Dear All,I want to concatenate a number of daily_data.nc on Windows OS. Please can someone identify a script that works on OpenGrADS 2.0.
Thanks. Dr Ugochukwu K. Okoro
 

    On Monday, October 2, 2017 11:55 PM, "gradsusr-request at gradsusr.org" <gradsusr-request at gradsusr.org> wrote:
 

 Send gradsusr mailing list submissions to
    gradsusr at gradsusr.org

To subscribe or unsubscribe via the World Wide Web, visit
    http://gradsusr.org/mailman/listinfo/gradsusr
or, via email, send a message with subject or body 'help' to
    gradsusr-request at gradsusr.org

You can reach the person managing the list at
    gradsusr-owner at gradsusr.org

When replying, please edit your Subject line so it is more specific
than "Re: Contents of gradsusr digest..."


Today's Topics:

  1. Re: Maskout two different fields (Andrew Friedman)
  2. Re: Maskout two different fields (Daniel Nairn)
  3. Re: Maskout two different fields (Andrew Friedman)
  4. Re: Maskout two different fields (Daniel Nairn)
  5. Re: Maskout two different fields (Jeff Duda)


----------------------------------------------------------------------

Message: 1
Date: Sun, 1 Oct 2017 22:55:38 +0200
From: Andrew Friedman <andfried at gmail.com>
Subject: Re: [gradsusr] Maskout two different fields
To: GrADS Users Forum <gradsusr at gradsusr.org>
Message-ID: <140EC942-0F86-4C2E-8BCE-52CCFC174ED2 at gmail.com>
Content-Type: text/plain; charset=us-ascii

Hi Daniel,

What if you simplify it by just looking at the inner loop: maskout(cape,myhel-100)? Does the output look correct? Then you can build on it by adding the successive const loops.

Andrew

> On Oct 1, 2017, at 12:59 PM, Daniel Nairn <dannairn at gmail.com> wrote:
> 
> Hello
> 
> I am trying to plot a chart (from one ensemble member) that displays SBCAPE where 0-1km helicity is greater than 100. This is how my current script looks
> 
> 'define cape = capesfc(e=2)'
> 'define myhel = HLCY0_1000m(e=2)'
> 'define capehel = const(const(maskout(cape,myhel-100),1),0.0,-u)'
> 'd capehel'
> 
> I'm not convinced this is correct as my final image output does not show any CAPE values where expected. Can anybody help at all?
> 
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr




------------------------------

Message: 2
Date: Mon, 2 Oct 2017 08:00:59 +0100
From: Daniel Nairn <dannairn at gmail.com>
Subject: Re: [gradsusr] Maskout two different fields
To: GrADS Users Forum <gradsusr at gradsusr.org>
Message-ID:
    <CAPo+bHo9ut7EOP0hUDdW+0dBHfBUULKjkTzLpyDcCKKb9tyf3Q at mail.gmail.com>
Content-Type: text/plain; charset="utf-8"

Hello Andrew

Yes I simplified the script as suggested at that does look much better,
thank you. Not quite sure what I've done wrong with the const loops though.

On 1 October 2017 at 21:55, Andrew Friedman <andfried at gmail.com> wrote:

> Hi Daniel,
>
> What if you simplify it by just looking at the inner loop:
> maskout(cape,myhel-100)? Does the output look correct? Then you can build
> on it by adding the successive const loops.
>
> Andrew
>
> > On Oct 1, 2017, at 12:59 PM, Daniel Nairn <dannairn at gmail.com> wrote:
> >
> > Hello
> >
> > I am trying to plot a chart (from one ensemble member) that displays
> SBCAPE where 0-1km helicity is greater than 100. This is how my current
> script looks
> >
> > 'define cape = capesfc(e=2)'
> > 'define myhel = HLCY0_1000m(e=2)'
> > 'define capehel = const(const(maskout(cape,myhel-100),1),0.0,-u)'
> > 'd capehel'
> >
> > I'm not convinced this is correct as my final image output does not show
> any CAPE values where expected. Can anybody help at all?
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20171002/c1993263/attachment-0001.html 

------------------------------

Message: 3
Date: Mon, 2 Oct 2017 12:52:54 +0200
From: Andrew Friedman <andfried at gmail.com>
Subject: Re: [gradsusr] Maskout two different fields
To: GrADS Users Forum <gradsusr at gradsusr.org>
Message-ID: <877F230E-6AB9-43AF-8199-72820E086481 at gmail.com>
Content-Type: text/plain; charset=us-ascii

Hi Daniel,
Could it be an issue with contour shading? If you pick a point where you know there is a value, is the output correct?
Andrew

> On Oct 2, 2017, at 9:00 AM, Daniel Nairn <dannairn at gmail.com> wrote:
> 
> Hello Andrew
> 
> Yes I simplified the script as suggested at that does look much better, thank you. Not quite sure what I've done wrong with the const loops though.
> 
> On 1 October 2017 at 21:55, Andrew Friedman <andfried at gmail.com> wrote:
> Hi Daniel,
> 
> What if you simplify it by just looking at the inner loop: maskout(cape,myhel-100)? Does the output look correct? Then you can build on it by adding the successive const loops.
> 
> Andrew
> 
> > On Oct 1, 2017, at 12:59 PM, Daniel Nairn <dannairn at gmail.com> wrote:
> >
> > Hello
> >
> > I am trying to plot a chart (from one ensemble member) that displays SBCAPE where 0-1km helicity is greater than 100. This is how my current script looks
> >
> > 'define cape = capesfc(e=2)'
> > 'define myhel = HLCY0_1000m(e=2)'
> > 'define capehel = const(const(maskout(cape,myhel-100),1),0.0,-u)'
> > 'd capehel'
> >
> > I'm not convinced this is correct as my final image output does not show any CAPE values where expected. Can anybody help at all?
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> 
> 
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
> 
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr




------------------------------

Message: 4
Date: Mon, 2 Oct 2017 12:12:58 +0100
From: Daniel Nairn <dannairn at gmail.com>
Subject: Re: [gradsusr] Maskout two different fields
To: GrADS Users Forum <gradsusr at gradsusr.org>
Message-ID:
    <CAPo+bHpL0gBngdi7B1G1Mkx560L5sUSS=0VHCyCyU9RKXRgvCQ at mail.gmail.com>
Content-Type: text/plain; charset="utf-8"

No I don't think so, what I've done is break down the script (as below)
again and this now actually looks okay

'define cape = capesfc'
'define myhel = HLCY0_1000m'
'define myval = maskout(cape,myhel-100)'
'define myval2 = const(myval, 0, -u)'
'd myval2'

I'm still unsure as to what I'm doing wrong when its all combined into one
line, but at least this appears to be working. Thank you.



On 2 October 2017 at 11:52, Andrew Friedman <andfried at gmail.com> wrote:

> Hi Daniel,
> Could it be an issue with contour shading? If you pick a point where you
> know there is a value, is the output correct?
> Andrew
>
> > On Oct 2, 2017, at 9:00 AM, Daniel Nairn <dannairn at gmail.com> wrote:
> >
> > Hello Andrew
> >
> > Yes I simplified the script as suggested at that does look much better,
> thank you. Not quite sure what I've done wrong with the const loops though.
> >
> > On 1 October 2017 at 21:55, Andrew Friedman <andfried at gmail.com> wrote:
> > Hi Daniel,
> >
> > What if you simplify it by just looking at the inner loop:
> maskout(cape,myhel-100)? Does the output look correct? Then you can build
> on it by adding the successive const loops.
> >
> > Andrew
> >
> > > On Oct 1, 2017, at 12:59 PM, Daniel Nairn <dannairn at gmail.com> wrote:
> > >
> > > Hello
> > >
> > > I am trying to plot a chart (from one ensemble member) that displays
> SBCAPE where 0-1km helicity is greater than 100. This is how my current
> script looks
> > >
> > > 'define cape = capesfc(e=2)'
> > > 'define myhel = HLCY0_1000m(e=2)'
> > > 'define capehel = const(const(maskout(cape,myhel-100),1),0.0,-u)'
> > > 'd capehel'
> > >
> > > I'm not convinced this is correct as my final image output does not
> show any CAPE values where expected. Can anybody help at all?
> > >
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20171002/ab4fffdd/attachment-0001.html 

------------------------------

Message: 5
Date: Mon, 2 Oct 2017 10:50:04 -0500
From: Jeff Duda <jeffduda319 at gmail.com>
Subject: Re: [gradsusr] Maskout two different fields
To: GrADS Users Forum <gradsusr at gradsusr.org>
Message-ID:
    <CAAig09CAmXct31OiHAfwNM56jHMsA=9U0fKLD-b-B4+TLapdUA at mail.gmail.com>
Content-Type: text/plain; charset="utf-8"

Daniel,
Why are you applying the const() function in this problem? From your code,
myval2 should be a field of all zeros except for in the masked area. I
would think you would want to plot myval only. If myval looks wrong then
either your cape or helicity fields are full of bad data. I would display
each separately to confirm if that is a problem.

Jeff Duda

On Mon, Oct 2, 2017 at 6:12 AM, Daniel Nairn <dannairn at gmail.com> wrote:

> No I don't think so, what I've done is break down the script (as below)
> again and this now actually looks okay
>
> 'define cape = capesfc'
> 'define myhel = HLCY0_1000m'
> 'define myval = maskout(cape,myhel-100)'
> 'define myval2 = const(myval, 0, -u)'
> 'd myval2'
>
> I'm still unsure as to what I'm doing wrong when its all combined into one
> line, but at least this appears to be working. Thank you.
>
>
>
> On 2 October 2017 at 11:52, Andrew Friedman <andfried at gmail.com> wrote:
>
>> Hi Daniel,
>> Could it be an issue with contour shading? If you pick a point where you
>> know there is a value, is the output correct?
>> Andrew
>>
>> > On Oct 2, 2017, at 9:00 AM, Daniel Nairn <dannairn at gmail.com> wrote:
>> >
>> > Hello Andrew
>> >
>> > Yes I simplified the script as suggested at that does look much better,
>> thank you. Not quite sure what I've done wrong with the const loops though.
>> >
>> > On 1 October 2017 at 21:55, Andrew Friedman <andfried at gmail.com> wrote:
>> > Hi Daniel,
>> >
>> > What if you simplify it by just looking at the inner loop:
>> maskout(cape,myhel-100)? Does the output look correct? Then you can build
>> on it by adding the successive const loops.
>> >
>> > Andrew
>> >
>> > > On Oct 1, 2017, at 12:59 PM, Daniel Nairn <dannairn at gmail.com> wrote:
>> > >
>> > > Hello
>> > >
>> > > I am trying to plot a chart (from one ensemble member) that displays
>> SBCAPE where 0-1km helicity is greater than 100. This is how my current
>> script looks
>> > >
>> > > 'define cape = capesfc(e=2)'
>> > > 'define myhel = HLCY0_1000m(e=2)'
>> > > 'define capehel = const(const(maskout(cape,myhel-100),1),0.0,-u)'
>> > > 'd capehel'
>> > >
>> > > I'm not convinced this is correct as my final image output does not
>> show any CAPE values where expected. Can anybody help at all?
>> > >
>> > > _______________________________________________
>> > > gradsusr mailing list
>> > > gradsusr at gradsusr.org
>> > > http://gradsusr.org/mailman/listinfo/gradsusr
>> >
>> >
>> > _______________________________________________
>> > gradsusr mailing list
>> > gradsusr at gradsusr.org
>> > http://gradsusr.org/mailman/listinfo/gradsusr
>> >
>> > _______________________________________________
>> > gradsusr mailing list
>> > gradsusr at gradsusr.org
>> > http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>


-- 
Jeff Duda
Post-doctoral research fellow
University of Oklahoma School of Meteorology
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20171002/abd94d37/attachment.html 

------------------------------

_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org
http://gradsusr.org/mailman/listinfo/gradsusr


End of gradsusr Digest, Vol 92, Issue 2
***************************************


   
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20171003/d0a576d5/attachment-0001.html 


More information about the gradsusr mailing list