[ncl-talk] Changing forecast time iteration for GFS

[Barry Lynn]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190219/01553992/attachment.html>