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() subroutine volgrid6d_class::volgrid6d_v7d_transform ( 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, dimension(:,:), intent(in), optional  maskgrid, real, dimension(:), intent(in), optional  maskbounds, character(len=*), intent(in), optional  networkname, logical, intent(in), optional  noconvert, procedure(basic_find_index), optional, pointer  find_index, character(len=*), intent(in), optional  categoryappend  ) private Performs the specified abstract transformation on the data provided. The abstract transformation is specified by this parameter; the corresponding specifical transformation (grid_transform object) is created and destroyed internally. The output transformed object is created internally and it does not require preliminary initialisation. Parametri [in] this object specifying the abstract transformation [in,out] volgrid6d_in object to be transformed, it is not modified, despite the INTENT(inout) [out] vol7d_out transformed object, it does not requires initialisation [in] v7d object containing a list of points over which transformation has to be done (required by some transformation types) [in] maskgrid 2D field to be used for defining subareas according to its values, it must have the same shape as the field to be interpolated (for transformation type 'maskinter') [in] maskbounds array of boundary values for defining subareas from the values of maskgrid, the number of subareas is SIZE(maskbounds) - 1, if not provided a default based on extreme values of makgrid is used [in] networkname set the output network name in vol7d_out (default='generic') [in] noconvert do not try to match variable and convert values during transform [in] categoryappend append this suffix to log4fortran namespace category Definizione alla linea 2346 del file volgrid6d_class.F90. ! no time_definition conversion implemented, output will be the same as input if (associated(vol7d_in%time))then ntime=size(vol7d_in%time) volgrid6d_out%time=vol7d_in%time end if if (associated(vol7d_in%timerange))then ntimerange=size(vol7d_in%timerange) volgrid6d_out%timerange=vol7d_in%timerange end if if (associated(vol7d_in%level))then nlevel=size(vol7d_in%level) volgrid6d_out%level=vol7d_in%level end if if (associated(vol7d_in%dativar%r))then nvar=size(vol7d_in%dativar%r) CALL varbufr2vargrib(vol7d_in%dativar%r, volgrid6d_out%var, c_func, gaid_template) end if nana=SIZE(vol7d_in%voldatir, 1) ! allocate once for speed IF (.NOT.ASSOCIATED(volgrid6d_out%voldati)) THEN ALLOCATE(voldatiout(volgrid6d_out%griddim%dim%nx, volgrid6d_out%griddim%dim%ny, & nlevel)) ENDIF DO ivar=1,nvar DO itimerange=1,ntimerange DO itime=1,ntime ! clone the gaid template where I have data IF (PRESENT(gaid_template)) THEN DO ilevel = 1, nlevel IF (any(c_e(vol7d_in%voldatir(:,itime,ilevel,itimerange,ivar,inetwork)))) THEN CALL copy(gaid_template, volgrid6d_out%gaid(ilevel,itime,itimerange,ivar)) ELSE volgrid6d_out%gaid(ilevel,itime,itimerange,ivar) = grid_id_new() ENDIF ENDDO ENDIF ! get data IF (ASSOCIATED(volgrid6d_out%voldati)) & ! improve!!!! CALL volgrid_get_vol_3d(volgrid6d_out, itime, itimerange, ivar, & voldatiout) ! do the interpolation CALL compute(this, & vol7d_in%voldatir(:,itime,:,itimerange,ivar,inetwork), voldatiout, & vol7d_in%dativar%r(ivar)) ! rescale valid data according to variable conversion table IF (ASSOCIATED(c_func)) THEN CALL compute(c_func(ivar), voldatiout(:,:,:)) ENDIF ! put data CALL volgrid_set_vol_3d(volgrid6d_out, itime, itimerange, ivar, & voldatiout) ENDDO ENDDO ENDDO IF (.NOT.ASSOCIATED(volgrid6d_out%voldati)) THEN DEALLOCATE(voldatiout) ENDIF IF (ASSOCIATED(c_func)) THEN DEALLOCATE(c_func) ENDIF END SUBROUTINE v7d_volgrid6d_transform_compute SUBROUTINE v7d_volgrid6d_transform(this, griddim, vol7d_in, volgrid6d_out, & networkname, gaid_template, categoryappend) TYPE(transform_def),INTENT(in) :: this TYPE(griddim_def),INTENT(in),OPTIONAL :: griddim ! TODO ripristinare intent(in) dopo le opportune modifiche in grid_class.F90 TYPE(vol7d),INTENT(inout) :: vol7d_in TYPE(volgrid6d),INTENT(out) :: volgrid6d_out CHARACTER(len=*),OPTIONAL,INTENT(in) :: networkname TYPE(grid_id),OPTIONAL,INTENT(in) :: gaid_template CHARACTER(len=*),INTENT(in),OPTIONAL :: categoryappend type(grid_transform) :: grid_trans integer :: ntime, ntimerange, nlevel, nvar !TODO la category sarebbe da prendere da vol7d !call l4f_category_log(vol7d_out%category,L4F_DEBUG,"start volgrid6d_transform") CALL vol7d_alloc_vol(vol7d_in) ! be safe ntime=SIZE(vol7d_in%time) ntimerange=SIZE(vol7d_in%timerange) nlevel=SIZE(vol7d_in%level) nvar=0 if (associated(vol7d_in%dativar%r)) nvar=size(vol7d_in%dativar%r) IF (nvar <= 0) THEN ! use vol7d category once it will be implemented CALL l4f_log(l4f_error, & "trying to transform a vol7d object incomplete or without real variables") CALL init(volgrid6d_out) ! initialize to empty CALL raise_error() RETURN ENDIF CALL init(grid_trans, this, vol7d_in, griddim, categoryappend=categoryappend) CALL init(volgrid6d_out, griddim, time_definition=vol7d_in%time_definition, & categoryappend=categoryappend) IF (c_e(grid_trans)) THEN CALL volgrid6d_alloc(volgrid6d_out, griddim%dim, ntime=ntime, nlevel=nlevel, &