grid_transform_c_e grid_transform_compute grid_transform_delete grid_transform_get_val grid_transform_grid_vol7d_init grid_transform_init grid_transform_levtype_levtype_init grid_transform_v7d_grid_compute grid_transform_vol7d_grid_init grid_transform_vol7d_vol7d_init interval_info_new transform_delete transform_get_val transform_init

grid_transform_vol7d_grid_init() subroutine grid_transform_vol7d_grid_init ( type(grid_transform), intent(out) this, type(transform_def), intent(in) trans, type(vol7d), intent(in) v7d_in, type(griddim_def), intent(in) out, character(len=*), intent(in), optional categoryappend ) private Constructor for a grid_transform object, defining a particular sparse points-to-grid transformation. It defines an object describing a transformation from a set of sparse points to a rectangular grid; the abstract type of transformation is described in the transformation object trans (type transform_def) which must have been properly initialised. The additional information required here is the list of the input sparse points in the form of a vol7d object (parameter v7d_in), which can be the same volume that will be successively used for interpolation, or a volume with just the same coordinate data, and the description of the output grid griddim (a griddim_def object). The generated grid_transform object is specific to the sparse point list and grid provided. The function c_e can be used in order to check whether the object has been successfully initialised, if the result is .FALSE., it should not be used further on. Parametri [out] this grid transformation object [in] trans transformation object [in] v7d_in vol7d object with the coordinates of the sparse point to be used as input (only information about coordinates is used) [in] out griddim object defining target grid [in] categoryappend append this suffix to log4fortran namespace category Definizione alla linea 2626 del file grid_transform_class.F90. point = georef_coord_new(x=lon1, y=lat1) DO n = 1, this%trans%poly%arraysize IF (inside(point, this%trans%poly%array(n))) THEN ! stop at the first matching polygon this%inter_index_x(i,1) = n EXIT ENDIF ENDDO ENDDO this%outnx=this%trans%poly%arraysize this%outny=1 CALL vol7d_alloc(v7d_out, nana=this%outnx) ! setup output point list, equal to average of polygon points ! warning, in case of concave areas points may coincide! CALL poly_to_coordinates(this%trans%poly, v7d_out) this%valid = .true. ! warning, no check of subtype ELSE IF (this%trans%trans_type == 'metamorphosis') THEN ! common to all metamorphosis subtypes this%innx = SIZE(v7d_in%ana) this%inny = 1 ! allocate index array ALLOCATE(this%point_index(this%innx,this%inny)) this%point_index(:,:) = imiss IF (this%trans%sub_type == 'all' ) THEN CALL metamorphosis_all_setup() ELSE IF (this%trans%sub_type == 'coordbb' ) THEN ALLOCATE(lon(this%innx),lat(this%innx)) ! count and mark points falling into requested bounding-box this%outnx = 0 this%outny = 1 CALL getval(v7d_in%ana(:)%coord,lon=lon,lat=lat) DO i = 1, this%innx ! IF (geo_coord_inside_rectang() IF (lon(i) > this%trans%rect_coo%ilon .AND. & lon(i) < this%trans%rect_coo%flon .AND. & lat(i) > this%trans%rect_coo%ilat .AND. & lat(i) < this%trans%rect_coo%flat) THEN ! improve! this%outnx = this%outnx + 1 this%point_index(i,1) = this%outnx ENDIF ENDDO IF (this%outnx <= 0) THEN CALL l4f_category_log(this%category,l4f_warn, & "metamorphosis:coordbb: no points inside bounding box "//& trim(to_char(this%trans%rect_coo%ilon))//","// & trim(to_char(this%trans%rect_coo%flon))//","// & trim(to_char(this%trans%rect_coo%ilat))//","// & trim(to_char(this%trans%rect_coo%flat))) ENDIF CALL vol7d_alloc(v7d_out, nana=this%outnx) ! collect coordinates of points falling into requested bounding-box n = 0 DO i = 1, this%innx IF (c_e(this%point_index(i,1))) THEN n = n + 1 CALL init(v7d_out%ana(n),lon=lon(i),lat=lat(i)) ENDIF ENDDO DEALLOCATE(lon, lat) this%valid = .true. ELSE IF (this%trans%sub_type == 'poly' ) THEN ! count and mark points falling into requested polygon this%outnx = 0 this%outny = 1 DO i = 1, this%innx ! temporary, improve!!!! CALL getval(v7d_in%ana(i)%coord,lon=lon1,lat=lat1) point = georef_coord_new(x=lon1, y=lat1) DO n = 1, this%trans%poly%arraysize IF (inside(point, this%trans%poly%array(n))) THEN ! stop at the first matching polygon this%outnx = this%outnx + 1 this%point_index(i,1) = n EXIT