Programma esempio semplice per la lettura di file grib. Programma che legge i grib contenuti in un file e li organizza in un vettore di oggetti gridinfo
27 integer :: category,ier
28 character(len=512):: a_name
29 type(arrayof_gridinfo) :: gridinfo
31 TYPE(grid_file_id) :: ifile
36 call l4f_launcher(a_name,a_name_force=
"demo2")
39 category=l4f_category_get(a_name//
".main")
47 ifile = grid_file_id_new(
'../data/in.grb',
'r')
50 gaid = grid_id_new(ifile)
51 IF (.NOT.
c_e(gaid))
EXIT 60 "Numero totale di grib: "//
to_char(ngrib))
63 CALL insert(gridinfo, nelem=ngrib)
67 ifile = grid_file_id_new(
'../data/in.grb',
'r')
70 gaid = grid_id_new(ifile)
71 IF (.NOT.
c_e(gaid))
EXIT 75 CALL init (gridinfo%array(ngrib), gaid=gaid)
76 CALL import(gridinfo%array(ngrib))
87 call l4f_category_delete(category)