display_volgrid6d display_volgrid6dv export_to_gridinfo export_to_gridinfov export_to_gridinfovv import_from_gridinfo import_from_gridinfovv v7d_v7d_transform v7d_volgrid6d_transform vg6d_c2a vg6d_reduce vg6d_rounding vg6d_wind_rot vg6d_wind_unrot vg6dv_rounding volgrid6d_alloc volgrid6d_alloc_vol volgrid6d_delete volgrid6d_export_to_file volgrid6d_import_from_file volgrid6d_init volgrid6d_read_from_file volgrid6d_transform volgrid6d_v7d_transform volgrid6d_v7d_transform_compute volgrid6d_write_on_file volgrid6dv_delete volgrid6dv_transform volgrid6dv_v7d_transform volgrid_get_vol_2d volgrid_get_vol_3d volgrid_set_vol_2d volgrid_set_vol_3d

volgrid6d_v7d_transform_compute() subroutine volgrid6d_class::volgrid6d_v7d_transform_compute ( type(grid_transform), intent(in)  this, type(volgrid6d), intent(in)  volgrid6d_in, type(vol7d), intent(inout)  vol7d_out, character(len=*), intent(in), optional  networkname, logical, intent(in), optional  noconvert  ) private Parametri [in] noconvert do not try to match variable and convert values during transform Definizione alla linea 2188 del file volgrid6d_class.F90. if (associated(volgrid6d_in%timerange)) ntimerange=size(volgrid6d_in%timerange) if (associated(volgrid6d_in%time)) then ntime=size(volgrid6d_in%time) if (time_definition /= volgrid6d_in%time_definition) then ! converto reference in validity allocate (validitytime(ntime,ntimerange),stat=stallo) if (stallo /=0)then call l4f_category_log(volgrid6d_in%category,l4f_fatal,"allocating memory") call raise_fatal_error() end if do itime=1,ntime do itimerange=1,ntimerange if (time_definition > volgrid6d_in%time_definition) then validitytime(itime,itimerange) = & volgrid6d_in%time(itime) + timedelta_new(sec=volgrid6d_in%timerange(itimerange)%p1) else validitytime(itime,itimerange) = & volgrid6d_in%time(itime) - timedelta_new(sec=volgrid6d_in%timerange(itimerange)%p1) end if end do end do ntime = count_distinct(reshape(validitytime,(/ntime*ntimerange/)), back=.true.) deallocate (validitytime) end if end if if (associated(volgrid6d_in%level)) nlevel=size(volgrid6d_in%level) if (associated(volgrid6d_in%var)) nvar=size(volgrid6d_in%var) CALL init(grid_trans, this, volgrid6d_in%griddim, v7d_locana, & maskgrid=maskgrid, maskbounds=maskbounds, find_index=find_index, & categoryappend=categoryappend) CALL init (vol7d_out,time_definition=time_definition) IF (c_e(grid_trans)) THEN nana=SIZE(v7d_locana%ana) CALL vol7d_alloc(vol7d_out, nana=nana, ntime=ntime, nlevel=nlevel, & ntimerange=ntimerange, ndativarr=nvar, nnetwork=nnetwork) vol7d_out%ana = v7d_locana%ana CALL get_val(grid_trans, output_point_index=point_index) IF (ALLOCATED(point_index)) THEN ! check that size(point_index) == nana? CALL vol7d_alloc(vol7d_out, nanavari=1) CALL init(vol7d_out%anavar%i(1), 'B01192') ENDIF CALL vol7d_alloc_vol(vol7d_out) IF (ALLOCATED(point_index)) THEN DO inetwork = 1, nnetwork vol7d_out%volanai(:,1,inetwork) = point_index(:) ENDDO ENDIF CALL compute(grid_trans, volgrid6d_in, vol7d_out, networkname, noconvert) ELSE CALL l4f_log(l4f_error, 'vg6d_v7d_transform: transformation not valid') CALL raise_error() ENDIF CALL delete(grid_trans) #ifdef HAVE_DBALLE ! set variables to a conformal state CALL vol7d_dballe_set_var_du(vol7d_out) #endif CALL delete(v7d_locana) END SUBROUTINE volgrid6d_v7d_transform SUBROUTINE volgrid6dv_v7d_transform(this, volgrid6d_in, vol7d_out, v7d, & maskgrid, maskbounds, networkname, noconvert, find_index, categoryappend) TYPE(transform_def),INTENT(in) :: this TYPE(volgrid6d),INTENT(inout) :: volgrid6d_in(:) TYPE(vol7d),INTENT(out) :: vol7d_out TYPE(vol7d),INTENT(in),OPTIONAL :: v7d REAL,INTENT(in),OPTIONAL :: maskgrid(:,:) REAL,INTENT(in),OPTIONAL :: maskbounds(:) CHARACTER(len=*),OPTIONAL,INTENT(in) :: networkname LOGICAL,OPTIONAL,INTENT(in) :: noconvert PROCEDURE(basic_find_index),POINTER,OPTIONAL :: find_index CHARACTER(len=*),INTENT(in),OPTIONAL :: categoryappend integer :: i TYPE(vol7d) :: v7dtmp CALL init(v7dtmp) CALL init(vol7d_out) DO i=1,SIZE(volgrid6d_in) CALL transform(this, volgrid6d_in(i), v7dtmp, v7d=v7d, & maskgrid=maskgrid, maskbounds=maskbounds, & networkname=networkname, noconvert=noconvert, find_index=find_index, & categoryappend=categoryappend) CALL vol7d_append(vol7d_out, v7dtmp) ENDDO END SUBROUTINE volgrid6dv_v7d_transform ! Internal method for performing sparse point to grid computations SUBROUTINE v7d_volgrid6d_transform_compute(this, vol7d_in, volgrid6d_out, networkname, gaid_template) TYPE(grid_transform),INTENT(in) :: this ! object specifying the specific transformation type(vol7d), INTENT(in) :: vol7d_in ! object to be transformed type(volgrid6d), INTENT(inout) :: volgrid6d_out ! transformed object CHARACTER(len=*),OPTIONAL,INTENT(in) :: networkname ! select the network to be processed from the \a vol7d input object, default the first network TYPE(grid_id),OPTIONAL,INTENT(in) :: gaid_template ! the template (typically grib_api) to be associated with output data, it also helps in improving variable conversion integer :: nana, ntime, ntimerange, nlevel, nvar INTEGER :: ilevel, itime, itimerange, ivar, inetwork REAL,POINTER :: voldatiout(:,:,:) type(vol7d_network) :: network TYPE(conv_func), pointer :: c_func(:) !TODO category sarebbe da prendere da vol7d #ifdef DEBUG CALL l4f_category_log(volgrid6d_out%category, l4f_debug, & 'start v7d_volgrid6d_transform_compute') #endif ntime=0 ntimerange=0 nlevel=0 nvar=0 IF (PRESENT(networkname)) THEN