[gradsusr] Printing values to a text file

Wed Mar 11 14:55:09 EDT 2020

ugrdprs and vgrdprs are assumed names of your u-wind and v-wind grids. In
operational NWP forecasts in the US, these are standard naming conventions.
But the names of the u- and v-wind grids may be very different in your
framework.

Just substitute the actual names representing those fields in your
particular data files where I put them in my examples.

Jeff

On Wed, Mar 11, 2020 at 12:33 PM Jan Ising <jising at aggies.ncat.edu> wrote:

> Hello Jeff,
>
>     Thank you for your reply, I apologize for my delayed response, I
> became very busy suddenly with unrelated tasks this week. I have been
> working multiple angles at this and will try what you had suggested next. I
> have discovered a couple more routes I can go after this but I want to try
> to avoid a manual CAPE calculation if possible, especially with GrADS.
> Should it come down to it I will reach out again with more details.
>     On a side note, could you please explain a little bit more on the
> commands "ugrdprs" and "vgrdprs"? I have not worked with these before.
>
>
>
> -Thanks,
> Jon
>
>  _______________________________________________
> Message: 2
> Date: Sat, 7 Mar 2020 12:18:59 -0700
> From: Jeff Duda <jeffduda319 at gmail.com>
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Subject: Re: [gradsusr] gradsusr Digest, Vol 121, Issue 6
> Message-ID:
>         <CAAig09CqJW7Kv=NZ1Y4VVTcE+W35n_HkfmkbDyFH7Yi_1=
> e+Bw at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Jon,
> Based on what you are telling me, I believe I did understand what you
> asked, and my original suggestion remains a viable option. In order to get
> the numerical values of a grid printed out in ASCII format, you need to
> display one of the fields in your data file. The code snippet you attached
> in your original email does not have that line anywhere. Because of that
> omission, you are getting output that differs from your expectations.
>
> But, given your additional information, this may turn out to be much easier
> than it first sounded:
>
> If you do have the fundamental components of BRN (CAPE and shear
> basically), and especially if you already have an array of CAPE values
> (thank goodness, because computing that manually from the T/qv grid would
> be no fun at all), you can calculate BRN directly in grads in one single
> command (or two or three, depending on how you want to break up the
> calculation):
>
> 'define shear = mag(ugrdprs(lev=*lev1*)-ugrdprs(lev=*lev2*),vgrdprs(lev=
> *lev1*)-vgrdprs(lev=*lev2*)'
> 'define BRN = CAPE / pow(shear,2)'
> 'display BRN'
>
> QED.
>
> Try that out and see what works. If you're missing CAPE in your data file,
> then we can talk more about your originally suggested method.
>
> Jeff Duda
>
> On Sat, Mar 7, 2020 at 11:56 AM Jan Ising <jising at aggies.ncat.edu> wrote:
>
> > Hello Jeff,
> >     Thank you for your response however, I believe there was some
> > confusion. Let me explain with more detail. I am trying to calculate the
> > Bulk Richardson Number (my data file (a .ctl file) does not have this
> > predefined), to do so I need to calculate each individual component of
> the
> > BRN from the different variables in my file (MUCAPE, Wind, etc...)
> > separately for each gridpoint, since GrADS does not allow for 2D arrays
> > that the user makes (say for instance in a script). This all works fine
> and
> > the values are outputted to my .txt. file, but only as 1 column with x
> rows
> > (where x is the number of grid points there are in my file). But I need
> it
> > to be in rows AND columns so that I can use the new grid, made only of
> > numbers (i.e the text file) and continue working with it in MATLAB (say
> for
> > example plotting it). Sorry if this was unclear before, I tried to
> simplify
> > it so my question is better understood, but it seems I may have
> > achieved the opposite.
> >
> >
> >> ----------------------------------------------------------------------
> >>
> >> Date: Fri, 6 Mar 2020 12:30:35 -0700
> >> From: Jeff Duda <jeffduda319 at gmail.com>
> >> To: GrADS Users Forum <gradsusr at gradsusr.org>
> >> Subject: Re: [gradsusr] Printing values to a text file
> >> Message-ID:1
> >>         <
> >> CAAig09DcTUY1j+k7hNDGMpHr1JhuS4D5k_3PEiVSG20aPV17cw at mail.gmail.com>
> >> Content-Type: text/plain; charset="utf-8"
> >>
> >> I think the problem is that you're not displaying a data field in your
> >> code. It looks like you're telling Grads just to spit out a sequence of
> >> integers. You don't need Grads for that. Also, the gxout print and
> prnopts
> >> probably aren't meant to be used for that. You could easily just do
> >>
> >> a = 1
> >> while (a < 400)
> >>  say a
> >>  a = a + 1
> >> endwhile
> >>
> >> Then copy your screen to a text file, or just use grads -H [*logfile*]
> to
> >> tell Grads to log commands + output to *logfile*
> >>
> >> Jeff Duda
> >>
> >> On Fri, Mar 6, 2020 at 11:47 AM Jan Ising <jising at aggies.ncat.edu>
> wrote:
> >>
> >> > Hello,
> >> >
> >> > How would I display multiple values per line in a .txt file from a
> >> > single-valued changing variable (see code)? Based on the code below it
> >> is
> >> > currently displaying 1 value per line instead of 200 values per line:
> >> >
> >> >  'set gxout print'
> >> > 'set prnopts %f 200 1'
> >> > a = 0
> >> >
> >> > while (a < 400)
> >> > a = a + 1
> >> > dummy = write(test_output, a, append)
> >> > endwhile
> >> >
> >> >
> >> > Result is:
> >> > 1
> >> > 2
> >> > 3
> >> > 4
> >> > 5
> >> > ...
> >> >
> >> > Instead of:
> >> > 1 2 3 4 5
> >> > 6 7 8 9 10
> >> >
> >> > any help would be greatly appreciated
> >> >
> >> > -Thanks,
> >> > Jon
> >> > _______________________________________________
> >> > gradsusr mailing list
> >> > gradsusr at gradsusr.org
> >> > http://gradsusr.org/mailman/listinfo/gradsusr
> >> >
> >>
> >>
> >> --
> >> Jeff Duda, Research Scientist
> >> University of Colorado Boulder
> >> Cooperative Institute for Research in Environmental Sciences
> >> NOAA/OAR/ESRL/Global Systems Division
> >> Boulder, CO
> >> -------------- next part --------------
> >> An HTML attachment was scrubbed...
> >> URL: <
> >>
> http://gradsusr.org/pipermail/gradsusr/attachments/20200306/31eab672/attachment-0001.html
> >> >
> >>
> >> ------------------------------
> >>
> >>
> >
> >
> > -Thanks,
> > Jon
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
>
>
> --
> Jeff Duda, Research Scientist
> University of Colorado Boulder
> Cooperative Institute for Research in Environmental Sciences
> NOAA/OAR/ESRL/Global Systems Division
> Boulder, CO
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <
> http://gradsusr.org/pipermail/gradsusr/attachments/20200307/f04f89ed/attachment-0001.html
> >
>
> ------------------------------
>
> Subject: Digest Footer
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> ------------------------------
>
> End of gradsusr Digest, Vol 121, Issue 7
> ****************************************
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>

-- 
Jeff Duda, Research Scientist
University of Colorado Boulder
Cooperative Institute for Research in Environmental Sciences
NOAA/OAR/ESRL/Global Systems Division
Boulder, CO
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://gradsusr.org/pipermail/gradsusr/attachments/20200311/31b3e303/attachment-0001.html>