[ncl-talk] Changing forecast time iteration for GFS

Tue Feb 19 13:01:26 MST 2019

Hi:

Could you please show us what happens when you run the program.  What the
output is.

Are you able to print("fname_in " + fname_in)

Are  you able to print("files = " + files)

This might help in the debugging.

Thank you

On Tue, Feb 19, 2019 at 7:55 PM tschonholz <tschonholz at smn.gov.ar> wrote:

> On 2019-02-19 17:15, Barry Lynn wrote:
>
> Hi:
>
> Please copy ncl-talk on these e-mails.
>
> I am not sure that you are using the sprinti command correctly.  I think
> you need a comma before the gfs name.  However, it would be best if you
> provided context for what you are trying to do.  For instance, a complete
> sample script that tries to read the filename.  Also, you could include a
> file for reference.
>
> https://www.ncl.ucar.edu/Document/Functions/Built-in/sprinti.shtml
>
> Barry
>
>
> On Tue, Feb 19, 2019 at 5:39 PM tschonholz <tschonholz at smn.gov.ar> wrote:
>
>> On 2019-02-19 15:35, Barry Lynn wrote:
>>
>> So, you want to read a string of files on your computer? In this case, it
>> doesn't matter the actual time step, just that you want to read named
>> files.
>>
>> Barry
>>
>> On Tue, Feb 19, 2019 at 5:29 PM tschonholz <tschonholz at smn.gov.ar> wrote:
>>
>>> On 2019-02-19 15:14, Barry Lynn wrote:
>>>
>>> Do you mean that you want to download hourly data, or just read it?
>>>
>>> Barry
>>>
>>> On Tue, Feb 19, 2019 at 7:15 AM Barry Lynn <barry.h.lynn at gmail.com>
>>> wrote:
>>>
>>>> Hello:
>>>>
>>>> Could you please explain in greater detail what you would like to do.
>>>>
>>>> What do you mean by changing the forecast time iteration?
>>>>
>>>> Barry
>>>>
>>>> On Mon, Feb 18, 2019 at 10:15 PM tschonholz <tschonholz at smn.gov.ar>
>>>> wrote:
>>>>
>>>>> Hi Everybody, i have been looking for help but couldn´t find.
>>>>>
>>>>> My gfs grib2 have forecast every 6 hs, so the name of the files are
>>>>> (for example): gep19.t00z.pgrb2f384.grb2 in which 384 is the time of
>>>>> forecast which comes by 6 hs.
>>>>>
>>>>> I have this part of my lib_ensambleGFS.ncl:
>>>>>
>>>>>
>>>>>
>>>>> function get_filelist_forecast(BASEDIR[1]:string, \
>>>>> nens:numeric, \
>>>>> y:numeric, \
>>>>> m:numeric, \
>>>>> d:numeric, \
>>>>> h:numeric, \
>>>>> hf:numeric)
>>>>> local files,path_in,fname_in,dd,fdate, \
>>>>> y_start,m_start,d_start,h_start, \
>>>>> y_end,m_end,d_end,h_end
>>>>> begin
>>>>> hf at units = sprinti("hours since %0.4i-", y) + \
>>>>> sprinti("%0.2i-", m) + \
>>>>> sprinti("%0.2i ", d) + \
>>>>> sprinti("%0.2i:00:00", h)
>>>>> dd = ispan(1,nens,1)
>>>>> if(hf.lt.0) then
>>>>> ;---Guess Case
>>>>> hf = -1
>>>>> fdate = ut_calendar(hf, 0)
>>>>> y_start = tointeger(fdate(:,0))
>>>>> m_start = tointeger(fdate(:,1))
>>>>> d_start = tointeger(fdate(:,2))
>>>>> h_start = tointeger(fdate(:,3))
>>>>> y_end = y
>>>>> m_end = m
>>>>> d_end = d
>>>>> h_end = h
>>>>> else
>>>>> ;---Forecast Case
>>>>> y_start = y
>>>>> m_start = m
>>>>> d_start = d
>>>>> h_start = h
>>>>> fdate = ut_calendar(hf, 0)
>>>>> y_end = tointeger(fdate(:,0))
>>>>> m_end = tointeger(fdate(:,1))
>>>>> d_end = tointeger(fdate(:,2))
>>>>> h_end = tointeger(fdate(:,3))
>>>>> end if
>>>>> path_in = sprinti("%0.4i", y_start) + \
>>>>> sprinti("%0.2i", m_start) + \
>>>>> sprinti("%0.2i", d_start) + \
>>>>> sprinti("%0.2i/", h_start)
>>>>> fname_in = "gep" + \
>>>>> sprinti("%0.2i.t00z.pgrb2f", dd) + \
>>>>> sprinti("%0.2i.grb2",h_end )
>>>>> files = BASEDIR + "/gfs/" + \
>>>>> path_in + \
>>>>> fname_in
>>>>> return files
>>>>> end
>>>>>
>>>>>
>>>>>
>>>>> -------------
>>>>>
>>>>> Can anyone please help me to change the way it reads? I have
>>>>> f06,12,18..and so on.
>>>>>
>>>>> Thanks a lot! I am new to the program and want to improve.
>>>>>
>>>>>
>>>>>
>>>>> Tamara.-
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> ncl-talk at ucar.edu
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>>
>>>>
>>>> --
>>>> Barry H. Lynn, Ph.D
>>>> Senior Associate Scientist, Lecturer,
>>>> The Institute of the Earth Science,
>>>> The Hebrew University of Jerusalem,
>>>> Givat Ram, Jerusalem 91904, Israel
>>>> Tel: 972 547 231 170
>>>> Fax: (972)-25662581
>>>>
>>>> C.E.O, Weather It Is, LTD
>>>> Weather and Climate Focus
>>>> http://weather-it-is.com
>>>> Jerusalem, Israel
>>>> Local: 02 930 9525
>>>> Cell: 054 7 231 170
>>>> Int-IS: x972 2 930 9525
>>>>
>>>>
>>>
>>> --
>>> Barry H. Lynn, Ph.D
>>> Senior Associate Scientist, Lecturer,
>>> The Institute of the Earth Science,
>>> The Hebrew University of Jerusalem,
>>> Givat Ram, Jerusalem 91904, Israel
>>> Tel: 972 547 231 170
>>> Fax: (972)-25662581
>>>
>>> C.E.O, Weather It Is, LTD
>>> Weather and Climate Focus
>>> http://weather-it-is.com
>>> Jerusalem, Israel
>>> Local: 02 930 9525
>>> Cell: 054 7 231 170
>>> Int-IS: x972 2 930 9525
>>>
>>> just read, i have it on my server
>>> --
>>>
>>>
>>
>> --
>> Barry H. Lynn, Ph.D
>> Senior Associate Scientist, Lecturer,
>> The Institute of the Earth Science,
>> The Hebrew University of Jerusalem,
>> Givat Ram, Jerusalem 91904, Israel
>> Tel: 972 547 231 170
>> Fax: (972)-25662581
>>
>> C.E.O, Weather It Is, LTD
>> Weather and Climate Focus
>> http://weather-it-is.com
>> Jerusalem, Israel
>> Local: 02 930 9525
>> Cell: 054 7 231 170
>> Int-IS: x972 2 930 9525
>>
>> the problem is that my script reads the name of the file, which contains
>> the forecast time, which I read in my script
>>
>> sprinti("%0.2i.t00z.pgrb2f", dd) + \
>>
>>
>>
>> it is dd
>> --
>>
>>
>
> --
> Barry H. Lynn, Ph.D
> Senior Associate Scientist, Lecturer,
> The Institute of the Earth Science,
> The Hebrew University of Jerusalem,
> Givat Ram, Jerusalem 91904, Israel
> Tel: 972 547 231 170
> Fax: (972)-25662581
>
> C.E.O, Weather It Is, LTD
> Weather and Climate Focus
> http://weather-it-is.com
> Jerusalem, Israel
> Local: 02 930 9525
> Cell: 054 7 231 170
> Int-IS: x972 2 930 9525
>
>  This is my file. I need the file to read my GFS ensemble forecast (every
> 6 hours) to be read and the program is set so that it reads every hour, so
> it won´t do anything at all (it keeps telling me there are missing files).
>
> Can you please help?
> --
>
>

-- 
Barry H. Lynn, Ph.D
Senior Associate Scientist, Lecturer,
The Institute of the Earth Science,
The Hebrew University of Jerusalem,
Givat Ram, Jerusalem 91904, Israel
Tel: 972 547 231 170
Fax: (972)-25662581

C.E.O, Weather It Is, LTD
Weather and Climate Focus
http://weather-it-is.com
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190219/260e4f60/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 35ae030d.jpeg
Type: image/jpeg
Size: 13273 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190219/260e4f60/attachment.jpeg>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 35ae030d.jpeg
Type: image/jpeg
Size: 13273 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190219/260e4f60/attachment-0001.jpeg>