[gradsusr] GFS hourly clouds output (Ivan Toman)

Bimal Gyawali bml.gya at gmail.com
Tue Oct 4 16:24:48 EDT 2016


Hello All,
I need help to calculate lag correlation. I have two monthly time series. I
want to plot lag-lead correlation lag-30 to lead-30. Would you please help
me to calculate lag-lead correlation. I want to plot similar to the plot
which I have attached here with.
.[image: Inline image 1]

Thank you.

*Bimal Gyawali*





* Graduate Research Assistant Department of Physical and Environment
Science Natural Resources Center Texas A&M University-Corpus Christi Corpus
Christi, TX 78412 Tel: 435-881-3603*

On Sat, Oct 1, 2016 at 5:48 PM, Ivan Toman <ivtoman at inet.hr> wrote:

> Hello,
>
> Thanks, I see. I suspected something like that from visual
> investigation. Thank you again for clarification.
>
> Best regards,
> Ivan
>
> On 09/30/2016 04:21 PM, Matt Rydzik wrote:
> > Hi Ivan,
> >
> > The forecast fields directly from NCEP are averaged.
> >
> > Forecast Hour 22:
> > 316:174518205:d=2016093006:TCDC:low cloud layer:18-22 hour ave fcst:
> > 317:174968393:d=2016093006:TCDC:middle cloud layer:18-22 hour ave fcst:
> > 318:175285665:d=2016093006:TCDC:high cloud layer:18-22 hour ave fcst:
> > 319:175676409:d=2016093006:TCDC:entire atmosphere:18-22 hour ave fcst:
> > http://nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/
> gfs.2016093006/gfs.t06z.pgrb2.0p25.f022.idx
> >
> > Forecast Hour 24
> > 316:175355590:d=2016093006:TCDC:low cloud layer:18-24 hour ave fcst:
> > 317:175812820:d=2016093006:TCDC:middle cloud layer:18-24 hour ave fcst:
> > 318:176137930:d=2016093006:TCDC:high cloud layer:18-24 hour ave fcst:
> > 319:176531519:d=2016093006:TCDC:entire atmosphere:18-24 hour ave fcst:
> > http://nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/
> gfs.2016093006/gfs.t06z.pgrb2.0p25.f024.idx
> >
> > Matthew Rydzik
> > Commodity Weather Group, LLC
> > Meteorological Applications Developer
> > 4933 Auburn Ave, 1st Floor, Bethesda, MD 20814
> > Work: (240) 207-2500 | Cell (540) 216-1940 | Fax: (301) 656-1425
> >
> > -----Original Message-----
> > From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces@
> gradsusr.org] On Behalf Of gradsusr-request at gradsusr.org
> > Sent: Friday, September 30, 2016 9:53 AM
> > To: gradsusr at gradsusr.org
> > Subject: gradsusr Digest, Vol 79, Issue 30
> >
> > 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: GFS hourly clouds output (James T. Potemra)
> >     2. Re: GFS hourly clouds output (Ivan Toman)
> >     3. Fog/Visibility (Chorley Weather)
> >     4. Re: gradsusr Digest, Vol 79, Issue 29 (ramakrishna surireddi)
> >     5. How to calculate lead lag correlation? (Bimal Gyawali)
> >
> >
> > ----------------------------------------------------------------------
> >
> > Message: 1
> > Date: Thu, 29 Sep 2016 08:25:06 -1000
> > From: "James T. Potemra" <jimp at hawaii.edu>
> > Subject: Re: [gradsusr] GFS hourly clouds output
> > To: GrADS Users Forum <gradsusr at gradsusr.org>
> > Message-ID: <344f8aa0-1ddd-46bd-dfb9-4ff26a70c95c at hawaii.edu>
> > Content-Type: text/plain; charset=windows-1252; format=flowed
> >
> > Hi Ivan,
> >
> > Are you getting these from a single forecast run?  If not, and you are
> looking at an aggregation of different forecasts, you might be seeing an
> impact from re-initialization of the model (I think the GFS is run at 00,
> 06, 12 and 18 daily).
> >
> > Jim
> >
> >
> > On 9/27/16 10:26 AM, Ivan Toman wrote:
> >> Hello always helpful GrADS community;
> >>
> >> I have problem with displaying GFS 0.25deg hourly clouds. It seems to
> >> me that something is broken. It is hard to explain, but best
> >> explanation is this 24h animation:
> >> http://www.meteoadriatic.net/pub/clouds_01.gif
> >>
> >> So, it seems that clouds somehow "jump" at every 6th hour, and
> >> inbetween are less and less defined, something like averaged more and
> >> more every hour, until they "jump" little bit and become defined again
> >> well. And then again in loop that repeats it's behaviour every 6
> hours...
> >>
> >> What is going on there?
> >>
> >> Many thanks,
> >>
> >> Ivan Toman
> >>
> >> _______________________________________________
> >> gradsusr mailing list
> >> gradsusr at gradsusr.org
> >> http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> > ------------------------------
> >
> > Message: 2
> > Date: Thu, 29 Sep 2016 22:19:41 +0200
> > From: Ivan Toman <ivtoman at inet.hr>
> > Subject: Re: [gradsusr] GFS hourly clouds output
> > To: gradsusr at gradsusr.org
> > Message-ID: <2bc73df9-f68e-d18f-d2fe-5dd64f47b441 at inet.hr>
> > Content-Type: text/plain; charset=windows-1252; format=flowed
> >
> > Hello,
> >
> > Thank you for comment;
> >
> > Yes of course it is single run.
> >
> > Just to explain briefly what I do, because there is no reply, I started
> to believe it is my mistake somewhere. But I do think I do everything OK
> because all other things I postprocess are ok, including various ARW, NMM,
> NMMB, "old" GFS, and so on.. and also all other fields from GFS is OK,
> exept cloud cover.
> >
> > So, I download series of grib2 GFS files from grib filter using time
> advancing for loop in bash. After gribs are in place, I concatenate them
> into single grib, by:
> > cat *pgrb2f* > $GRIBDIR/gfs.grib2
> >
> > After that I run g2ctl and gribmap on that grib:
> > g2ctl -verf gfs.grib2 > gfs.grib2.ctl
> > gribmap -i gfs.grib2.ctl
> >
> > Finally, I plot fields using GrADS, for cloud cover using TCDChcll,
> TCDCmcll, TCDClcll and TCDCccll variables, in time advancing loop inside
> GrADS over t:
> >
> > while ( i <= last )
> > 'set t 'i
> > ...
> > endwhile
> >
> > Do I miss something?
> >
> > Thanks again,
> > Ivan Toman
> >
> >
> > On 09/29/2016 08:25 PM, James T. Potemra wrote:
> >> Hi Ivan,
> >>
> >> Are you getting these from a single forecast run?  If not, and you are
> >> looking at an aggregation of different forecasts, you might be seeing
> >> an impact from re-initialization of the model (I think the GFS is run
> >> at 00, 06, 12 and 18 daily).
> >>
> >> Jim
> >>
> >>
> >> On 9/27/16 10:26 AM, Ivan Toman wrote:
> >>> Hello always helpful GrADS community;
> >>>
> >>> I have problem with displaying GFS 0.25deg hourly clouds. It seems to
> >>> me that something is broken. It is hard to explain, but best
> >>> explanation is this 24h animation:
> >>> http://www.meteoadriatic.net/pub/clouds_01.gif
> >>>
> >>> So, it seems that clouds somehow "jump" at every 6th hour, and
> >>> inbetween are less and less defined, something like averaged more and
> >>> more every hour, until they "jump" little bit and become defined
> >>> again well. And then again in loop that repeats it's behaviour every 6
> hours...
> >>>
> >>> What is going on there?
> >>>
> >>> Many thanks,
> >>>
> >>> Ivan Toman
> >>>
> >>> _______________________________________________
> >>> 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: 3
> > Date: Thu, 29 Sep 2016 21:51:33 +0100
> > From: Chorley Weather <weatherstu at chorleyweather.com>
> > Subject: [gradsusr] Fog/Visibility
> > To: GrADS Users Forum <gradsusr at gradsusr.org>
> > Message-ID: <333159b105b51ab519b2976821eedb13 at chorleyweather.com>
> > Content-Type: text/plain; charset="us-ascii"
> >
> >
> >
> > Hi All,
> >
> > I wonder if anybody can help. Does anybody have the calculation or
> script to chart fog/visibility using GFS OP25 data? Many thanks in
> advance.....
> > --
> > Regards,
> >
> > Stuart Markham
> >
> >
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160929/
> b650a6c5/attachment-0001.html
> >
> > ------------------------------
> >
> > Message: 4
> > Date: Fri, 30 Sep 2016 00:33:43 +0000 (UTC)
> > From: ramakrishna surireddi <ssvs_rk at yahoo.co.in>
> > Subject: Re: [gradsusr] gradsusr Digest, Vol 79, Issue 29
> > To: "gradsusr at gradsusr.org" <gradsusr at gradsusr.org>
> > Message-ID: <827342991.6009858.1475195623463 at mail.yahoo.com>
> > Content-Type: text/plain; charset="utf-8"
> >
> > Many thanks for the command.best Wishesramakrishna
> ?Dr.S.S.V.S.Ramakrishna Professor & Chairman Department of Meteorology &
> Oceanography Andhra University Visakhapatnam-530003 Andhra Pradesh India
> >
> >        From: "gradsusr-request at gradsusr.org" <
> gradsusr-request at gradsusr.org>
> >   To: gradsusr at gradsusr.org
> >   Sent: Friday, 30 September 2016 1:00 AM
> >   Subject: gradsusr Digest, Vol 79, Issue 29
> >
> > 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: GRB files for Grads (Boddapati Anil) ? 2. decimal point as
> comma int csv output (AMS EmsiWx Byt?a)
> >
> >
> > ----------------------------------------------------------------------
> >
> > Message: 1
> > Date: Thu, 29 Sep 2016 07:19:50 -0400
> > From: Boddapati Anil <boddapati.anil2 at gmail.com>
> > Subject: Re: [gradsusr] GRB files for Grads
> > To: GrADS Users Forum <gradsusr at gradsusr.org>
> > Message-ID:
> > ??? <CAJdYfFGXRX3oRTm9ojQeTSQa=J4phw-y0c7ea9KZ2Os2ZjmOaA at mail.gmail.com>
> > Content-Type: text/plain; charset="utf-8"
> >
> > Dear Sir,
> >
> > grib2ctl.pl grib_file >grib_file.ctl
> > $ gribmap -i grib_file.ctl -0? ? <-- this is a zero $ grads Landscape
> mode? (no for portrait):
> > ga-> open grib_file.ctl
> >
> >
> > try above commands for creating ctl file and then ctl file I think it
> may useful to you
> >
> > With Regards
> > Boddapati Anil
> > Project Fellow
> > CSIR Fourth paradigm Institute
> > Bengaluru-17
> >
> > On Thu, Sep 29, 2016 at 7:06 AM, Jennifer M Adams <jadams21 at gmu.edu>
> wrote:
> >
> >> Here is the documentation on how to use GRIB with GrADS:
> >> http://cola.gmu.edu/grads/gadoc/grib.html
> >> ?Jennifer
> >>
> >> On Sep 29, 2016, at 2:33 AM, ramakrishna surireddi
> >> <ssvs_rk at yahoo.co.in>
> >> wrote:
> >>
> >> Dear Grads Users,
> >> Can any body help me with using a GRB file in Grads.
> >> Thank you,
> >> With Best Regards
> >> Ramakrishna
> >>
> >> Dr.S.S.V.S.Ramakrishna Professor & Chairman Department of Meteorology
> >> & Oceanography Andhra University Visakhapatnam-530003 Andhra Pradesh
> >> India _______________________________________________
> >> gradsusr mailing list
> >> gradsusr at gradsusr.org
> >> http://gradsusr.org/mailman/listinfo/gradsusr
> >>
> >>
> >> --
> >> Jennifer Miletta Adams
> >> Center for Ocean-Land-Atmosphere Studies (COLA) George Mason
> >> University
> >>
> >>
> >>
> >>
> >> _______________________________________________
> >> 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/20160929/
> 7030b063/attachment-0001.html
> >
> > ------------------------------
> >
> > Message: 2
> > Date: Thu, 29 Sep 2016 16:38:24 +0200
> > From: AMS EmsiWx Byt?a <emsiwx at pocasie-bytca.sk>
> > Subject: [gradsusr] decimal point as comma int csv output
> > To: GrADS Users Forum <gradsusr at gradsusr.org>
> > Message-ID: <33c1c379-4c70-61aa-89b1-5cdfee2c1810 at pocasie-bytca.sk>
> > Content-Type: text/plain; charset="utf-8"
> >
> > Dear Grads users,
> >
> > I have been using (please see example below) to output data to csv:
> >
> > getseries(TMP2m,T2m,1000)
> > 'fprintf T2m-273.16 'name'_'date'.csv %.1f, 49 0'
> > getseries(UGRD10m,ugrd,1000)
> > 'fprintf ugrd 'name'_'date'.csv %.1f, 49 0'
> > getseries(VGRD10m,vgrd,1000)
> >
> > so I get a line of a single parameter for a period of time
> >
> > Now I would like to get something like this:
> >
> > Line 1:? ? ? ? ? ? ? dd.mm.yyyy; model run; HH; zip code (this can be
> done in command line as input parameters); parameter1 (e.g. T2m);
> > parameter2 (e.g. T850); parameter3 (etc.) Line 2:? ? ? ? ? ? dd.mm.yyyy;
> model run; HH (this means next forecast hour); zip code (this can be done
> in command line as input parameters);
> > parameter1 (e.g. T2m); parameter2 (e.g. T850); parameter3 (etc.)
> >
> > So I would like to ask:
> >
> > 1. How to output dd.mm.yyyy (probably using as input parameter, if you
> know any other way please comment) 2. How to output model run (probably
> using as input parameter, if you know any other way please comment)
> >
> > 3. How to output HH - means forecast hour
> >
> > 4. Since in Europe, we preferably use semi-colon in csv to separate the
> values, change comma to semi-colon as a separator
> >
> > 5. Output comma "," instead of decimal point? "." ..... meaning 2.5
> should be 2,5
> >
> >
> > Thank your very much for your time and help. I highly appreciate it.
> >
> >
> > Kindest regards,
> >
> >
> > Marian
> >
> >
> >
> >
> >
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160929/
> ff932459/attachment-0001.html
> >
> > ------------------------------
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> > End of gradsusr Digest, Vol 79, Issue 29
> > ****************************************
> >
> >
> >
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160930/
> b101ec3a/attachment-0001.html
> >
> > ------------------------------
> >
> > Message: 5
> > Date: Fri, 30 Sep 2016 08:52:42 -0500
> > From: Bimal Gyawali <bml.gya at gmail.com>
> > Subject: [gradsusr] How to calculate lead lag correlation?
> > To: GrADS Users Forum <gradsusr at gradsusr.org>
> > Message-ID:
> >       <CAHVD8=pM1cJXNyXwS+BrcMmzDrv7yvRQz=C4Ev=ZmczY_LO+
> +A at mail.gmail.com>
> > Content-Type: text/plain; charset="utf-8"
> >
> > Hello All,
> > I am trying to calculate lead lag correlation and plot it. I want to
> > calculate nlag=-30 to nlag=30. How can calculate and plot this? please
> help
> > me with this. I have two time series with n=730. I want to plot similar
> to
> > this.
> > Thank you.[image: Inline image 1]
> >
> > *Bimal Gyawali*
> >
> >
> >
> >
> >
> > * Graduate Research Assistant Department of Physical and Environment
> > Science Natural Resources Center Texas A&M University-Corpus Christi
> Corpus
> > Christi, TX 78412 Tel: 435-881-3603*
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160930/
> 1e34f00e/attachment.html
> > -------------- next part --------------
> > A non-text attachment was scrubbed...
> > Name: ....PNG
> > Type: image/png
> > Size: 32804 bytes
> > Desc: not available
> > Url : http://gradsusr.org/pipermail/gradsusr/attachments/20160930/
> 1e34f00e/attachment.png
> >
> > ------------------------------
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> > End of gradsusr Digest, Vol 79, Issue 30
> > ****************************************
> >
> > _______________________________________________
> > 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/20161004/a6cd41a7/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ....PNG
Type: image/png
Size: 32804 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20161004/a6cd41a7/attachment-0001.png 


More information about the gradsusr mailing list