example_alchimiavg6d.f03

Sample program to demostrate the alchimia module with volgrid6d.

Sample program to demostrate the alchimia module with volgrid6d.

program alchimiavg6d
 
USE alchimia
USE termo
USE volgrid6d_class
USE volgrid6d_alchimia_class
use log4fortran
 
IMPLICIT NONE
 
type(fndsv) :: vfn, vfnoracle
character(len=10), allocatable:: mybout(:)
type(volgrid6d),pointer :: myin(:),myout(:)
 
character(len=255) :: filenamein,filenameout
 
TYPE(l4f_handle) :: category
INTEGER :: ier
CHARACTER(len=512):: a_name
 
!questa chiamata prende dal launcher il nome univoco
call l4f_launcher(a_name)
 
!init di log4fortran
ier=l4f_init()
 
!imposta a_name
category=l4f_category_get_handle(a_name)
 
call l4f_category_log(category,l4f_info,"Start")
 
mybout = [character(len=10) :: "B12192"]
filenamein="../data/t_p.grb"
filenameout="../data/tp.grb"
 
call register_termo(vfn)
 
CALL import(myin,filename=filenamein,decode=.true., time_definition=0, categoryappend="input")
 
call display(myin)
 
if (alchemy(myin,vfn,mybout,myout,copy=.true.,vfnoracle=vfnoracle) /= 0) then
  print*, "I cannot make ",mybout
  
  if (.not. shoppinglist(mybout,vfn,vfnoracle)) then
    print*, " error shoppinglist"
    stop 2
  else
        call display(compile_sl(vfnoracle))
    stop 3
  end if
end if
 
call display(vfnoracle)
 
call export(myout,filenameout)
call display(myout)
 
call delete(myout)
call delete(myin)
 
!chiudo il logger
call l4f_category_delete(category)
ier=l4f_fini()
 
end program alchimiavg6d