[gradsusr] GFS hourly clouds output (Ivan Toman)

Matt Rydzik mrydzik at commoditywx.com
Fri Sep 30 10:21:03 EDT 2016


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 at 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
****************************************



More information about the gradsusr mailing list