example_vg6d_6.f90

Programma trasformazione da volgrid6d a vol7d

! 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 demo6

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

implicit none

integer :: ier,i,nana
type(l4f_handle) :: category
character(len=512):: a_name,filename="out.bufr"
type(volgrid6d),pointer  :: volgrid(:),volgrid_out(:)
type(transform_def) :: trans
type(vol7d) :: v7d
type(vol7d) :: vol7d_out
TYPE(vol7d_dballe) :: v7d_exp


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

!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")

!!$ nana=5
!!$ allocate(ana(nana))
!!$ !target points 
!!$ call init(ana(1),lat=45.D0,lon=11.D0)
!!$ call init(ana(2),lat=45.6D0,lon=11.8D0)
!!$ call init(ana(3),lat=46.6D0,lon=12.8D0)
!!$ call init(ana(4),lat=40.6D0,lon=11.8D0)
!!$ call init(ana(5),lat=40.0D0,lon=10.0D0)

call init(v7d)

call import (v7d,filename="ana.v7d")

call display(v7d)

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

call import (volgrid,filename="in.grb",categoryappend="volume letto")

call transform(trans, volgrid6d_in=volgrid, vol7d_out=vol7d_out, v7d=v7d, categoryappend="trasforma")

call l4f_category_log(category,l4f_info,"trasformato")

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

call l4f_category_log(category,l4f_info,"export to BUFR")


! Chiamo il costruttore della classe vol7d_dballe per il mio oggetto in export
CALL init(v7d_exp,file=.true.,write=.true.,wipe=.true.,filename=filename,&
categoryappend="exportBUFR",format="BUFR",template="synop")

CALL display(vol7d_out)
v7d_exp%vol7d = vol7d_out
CALL export(v7d_exp)

CALL l4f_category_log(category,l4f_info,"terminato")

CALL delete (v7d_exp)

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

end program demo6