30 integer :: category,ier
31 character(len=512):: a_name,filename=
"out.bufr"
37 integer :: nx=40,ny=40,component_flag=0
39 doubleprecision :: xmin=0., xmax=30., ymin=30., ymax=60.
40 doubleprecision :: latitude_south_pole=-32.5,longitude_south_pole=10.,angle_rotation=0.
41 character(len=80) :: type=
'regular_ll',trans_type=
'inter',sub_type=
'linear'
44 call l4f_launcher(a_name,a_name_force=
"demo7")
50 category=l4f_category_get(a_name//
".main")
56 call
init(griddim_out,&
57 proj_type=type,nx=nx,ny=ny, &
58 xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax, component_flag=component_flag, &
59 latitude_south_pole=latitude_south_pole,longitude_south_pole=longitude_south_pole,angle_rotation=angle_rotation, &
62 call griddim_unproj(griddim_out)
64 print*,
'grid di interpolazione >>>>>>>>>>>>>>>>>>>>'
67 gaid_template = grid_id_new(grib_api_template=
"regular_ll_sfc_grib1")
70 call
init(trans, trans_type=trans_type,sub_type=sub_type, categoryappend=
"trasformation")
73 CALL
init(v7d_import,file=.true.,write=.false.,filename=filename,&
74 categoryappend=
"importBUFR",format=
"BUFR")
76 call
import(v7d_import,var=(/
"B12101"/),varkind=(/
"r"/))
81 call
transform(trans,griddim_out, vol7d_in=v7d_import%vol7d, &
82 volgrid6d_out=volgrid(1), gaid_template=gaid_template, &
83 categoryappend=
"trasform->")
86 CALL
export(volgrid,
'examp[le_v7d.grb', gaid_template=gaid_template,&
87 categoryappend=
"volume scritto")
89 if (
associated(volgrid)) call
delete(volgrid)
96 call l4f_category_delete(category)