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_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 2176 del file volgrid6d_class.F90. call get_val(this,time_definition=time_definition) if (.not. c_e(time_definition)) then time_definition=1 ! default to validity time endif IF (PRESENT(v7d)) THEN CALL vol7d_copy(v7d, v7d_locana) ELSE CALL init(v7d_locana) ENDIF 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 !> Performs the specified abstract transformation on the arrays of !! 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. The transformation performed on !! each element of the input \a volgrid6d array object is merged into !! a single \a vol7d output object. SUBROUTINE volgrid6dv_v7d_transform(this, volgrid6d_in, vol7d_out, v7d, & maskgrid, maskbounds, networkname, noconvert, find_index, categoryappend) TYPE(transform_def),INTENT(in) :: this !< object specifying the abstract transformation TYPE(volgrid6d),INTENT(inout) :: volgrid6d_in(:) !< object to be transformed, it is an array of volgrid6d objects, each of which will be transformed, it is not modified, despite the INTENT(inout) TYPE(vol7d),INTENT(out) :: vol7d_out !< transformed object, it does not require initialisation TYPE(vol7d),INTENT(in),OPTIONAL :: v7d !< object containing a list of points over which transformation has to be done (required by some transformation types) REAL,INTENT(in),OPTIONAL :: 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') REAL,INTENT(in),OPTIONAL :: maskbounds(:) !< array of boundary values for defining subareas from the values of \a maskgrid, the number of subareas is SIZE(maskbounds) - 1, if not provided a default based on extreme values of \a makgrid is used CHARACTER(len=*),OPTIONAL,INTENT(in) :: networkname !< set the output network name in vol7d_out (default='generic') LOGICAL,OPTIONAL,INTENT(in) :: noconvert !< do not try to match variable and convert values during transform PROCEDURE(basic_find_index),POINTER,OPTIONAL :: find_index CHARACTER(len=*),INTENT(in),OPTIONAL :: categoryappend !< append this suffix to log4fortran namespace category 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