v7d_v7d_transform v7d_volgrid6d_transform volgrid6d_transform volgrid6d_v7d_transform volgrid6dv_transform volgrid6dv_v7d_transform

volgrid6d_v7d_transform() subroutine 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 ) 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 2334 del file volgrid6d_class.F90. inetwork = index(vol7d_in%network,network) IF (inetwork <= 0) THEN CALL l4f_category_log(volgrid6d_out%category,l4f_warn, & 'network '//trim(networkname)//' not found, first network will be transformed') inetwork = 1 ENDIF ELSE inetwork = 1 ENDIF ! 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 !> Performs the specified abstract transformation on the data provided. !! The abstract transformation is specified by \a this parameter; the !! corresponding specifical transformation (\a grid_transform object) !! is created and destroyed internally. The output transformed object !! is created internally and it does not require preliminary !! initialisation. SUBROUTINE v7d_volgrid6d_transform(this, griddim, vol7d_in, volgrid6d_out, & networkname, gaid_template, categoryappend) TYPE(transform_def),INTENT(in) :: this !< object specifying the abstract transformation TYPE(griddim_def),INTENT(in),OPTIONAL :: griddim !< griddim specifying the output grid (required by most transformation types) ! TODO ripristinare intent(in) dopo le opportune modifiche in grid_class.F90 TYPE(vol7d),INTENT(inout) :: vol7d_in !< object to be transformed, it is not modified, despite the INTENT(inout) TYPE(volgrid6d),INTENT(out) :: volgrid6d_out !< transformed object, it does not require initialisation 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 CHARACTER(len=*),INTENT(in),OPTIONAL :: categoryappend !< append this suffix to log4fortran namespace category 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()