example_vg6d_7.f90

Programma trasformazione da vol7d a volgrid6d

! Copyright (C) 2010  ARPA-SIM <urpsim@smr.arpa.emr.it>
! authors:
! Davide Cesari <dcesari@arpa.emr.it>
! Paolo Patruno <ppatruno@arpa.emr.it>

! This program is free software; you can redistribute it and/or
! modify it under the terms of the GNU General Public License as
! published by the Free Software Foundation; either version 2 of 
! the License, or (at your option) any later version.

! This program is distributed in the hope that it will be useful,
! but WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
! GNU General Public License for more details.

! You should have received a copy of the GNU General Public License
! along with this program.  If not, see <http://www.gnu.org/licenses/>.
program demo7

use log4fortran
use volgrid6d_class
use grid_class
use grid_transform_class
USE vol7d_dballe_class
USE vol7d_class
USE grid_id_class

implicit none

integer :: ier
type(l4f_handle) :: category
character(len=512):: a_name,filename="out.bufr"
type(volgrid6d),pointer  :: volgrid(:)
type(transform_def) :: trans
TYPE(vol7d_dballe) :: v7d_import
type(griddim_def) :: griddim_out

integer :: nx=40,ny=40,component_flag=0
type(grid_id) :: gaid_template
doubleprecision :: xmin=0., xmax=30., ymin=30., ymax=60.
doubleprecision :: latitude_south_pole=-32.5,longitude_south_pole=10.,angle_rotation=0.
character(len=80) :: type='regular_ll',trans_type='inter',sub_type='linear'

!questa chiamata prende dal launcher il nome univoco
call l4f_launcher(a_name,a_name_force="demo7")

!init di log4fortran
ier=l4f_init()

!imposta a_name
category=l4f_category_get_handle(a_name//".main")

call l4f_category_log(category,l4f_info,"inizio")

allocate (volgrid(1))

call init(griddim_out,&
 proj_type=type,nx=nx,ny=ny, &
 xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax, component_flag=component_flag, &
 latitude_south_pole=latitude_south_pole,longitude_south_pole=longitude_south_pole,angle_rotation=angle_rotation, &
 categoryappend=type)

call griddim_unproj(griddim_out)

print*,'grid di interpolazione >>>>>>>>>>>>>>>>>>>>'
call display(griddim_out)

gaid_template = grid_id_new(grib_api_template="regular_ll_sfc_grib1")

!trasformation object
call init(trans, trans_type=trans_type,sub_type=sub_type, categoryappend="trasformation")

! Chiamo il costruttore della classe vol7d_dballe per il mio oggetto in import
CALL init(v7d_import,file=.true.,write=.false.,filename=filename,&
 categoryappend="importBUFR",format="BUFR")

call import (v7d_import,var=(/"B12101"/),varkind=(/"r"/))

call display(v7d_import%vol7d)

call l4f_category_log(category,l4f_info,"trasformato")
call transform(trans,griddim_out, vol7d_in=v7d_import%vol7d, &
 volgrid6d_out=volgrid(1), gaid_template=gaid_template, &
 categoryappend="trasform->")

call l4f_category_log(category,l4f_info,"export to GRIB")
CALL export(volgrid, 'examp[le_v7d.grb', gaid_template=gaid_template,&
 categoryappend="volume scritto")

if (associated(volgrid)) call delete(volgrid)

call l4f_category_log(category,l4f_info,"terminato")

call delete (v7d_import)

!chiudo il logger
call l4f_category_delete(category)
ier=l4f_fini()

end program demo7