example_vg6d_7.f90

Programma trasformazione da vol7d a volgrid6d.

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