grid_transform_grid_vol7d_init grid_transform_init grid_transform_levtype_levtype_init grid_transform_vol7d_grid_init grid_transform_vol7d_vol7d_init transform_init

transform_init() subroutine transform_init ( type(transform_def), intent(out) this, character(len=*) trans_type, character(len=*) sub_type, integer, intent(in), optional ix, integer, intent(in), optional iy, integer, intent(in), optional fx, integer, intent(in), optional fy, doubleprecision, intent(in), optional ilon, doubleprecision, intent(in), optional ilat, doubleprecision, intent(in), optional flon, doubleprecision, intent(in), optional flat, integer, intent(in), optional npx, integer, intent(in), optional npy, doubleprecision, intent(in), optional boxdx, doubleprecision, intent(in), optional boxdy, doubleprecision, intent(in), optional radius, type(arrayof_georef_coord_array), optional poly, doubleprecision, intent(in), optional percentile, real, intent(in), optional interv_gt, real, intent(in), optional interv_ge, real, intent(in), optional interv_lt, real, intent(in), optional interv_le, logical, intent(in), optional extrap, integer, intent(in), optional time_definition, type(vol7d_level), intent(in), optional input_levtype, type(vol7d_var), intent(in), optional input_coordvar, type(vol7d_level), intent(in), optional output_levtype, character(len=*), intent(in), optional categoryappend ) Constructor for a transform_def object, defining an abstract transformation between gridded and/or sparse point data. The parameters trans_type and sub_type define the type of transformation, while all the other following parameters are optional, they have to be passed in keyword mode and those required by the transformation type and subtype chosen have to be present. Parametri [out] this transformation object trans_type type of transformation, can be 'zoom', 'boxregrid', 'interp', 'vertint' ... sub_type sub type of transformation, it depends on trans_type [in] ix index of initial point of new grid on x (for zoom) [in] iy index of initial point of new grid on y (for zoom) [in] fx index of final point of new grid on x (for zoom) [in] fy index of final point of new grid on y (for zoom) [in] ilon coordinate of initial point of new grid or of bounding box on x (for zoom and metamorphosis) [in] ilat coordinate of initial point of new grid or of bounding box on y (for zoom and metamorphosis) [in] flon coordinate of final point of new grid or of bounding box on x (for zoom and metamorphosis) [in] flat coordinate of final point of new grid or of bounding box on y (for zoom and metamorphosis) [in] npx number of points to average along x direction (for boxregrid) [in] npy number of points to average along y direction (for boxregrid) [in] boxdx longitudinal/x extension of the box for box interpolation, default the target x grid step (unimplemented !) [in] boxdy latitudinal/y extension of the box for box interpolation, default the target y grid step (unimplemented !) [in] radius radius of stencil in grid points (also fractionary values) for stencil interpolation poly array of polygons indicating areas over which to interpolate (for transformations 'polyinter' or 'metamorphosis:poly') [in] percentile percentile [0,100.] of the distribution of points in the box to use as interpolated value for 'percentile' subtype [in] interv_gt greater than condition for defining interval [in] interv_ge greater equal condition for defining interval [in] interv_lt less than condition for defining interval [in] interv_le less equal condition for defining interval [in] extrap activate extrapolation outside input domain (use with care!) [in] time_definition time definition for output vol7d object 0=time is reference time ; 1=time is validity time [in] input_levtype type of vertical level of input data to be vertically interpolated (only type of first and second surface are used, level values are ignored) [in] input_coordvar variable that defines the vertical coordinate in the input volume for vertical interpolation, if missing, the value of the vertical level defined with input_levtype is used [in] output_levtype type of vertical level to which data should be vertically interpolated (only type of first and second surface are used, level values are ignored) [in] categoryappend suffix to append to log4fortran namespace category Definizione alla linea 651 del file grid_transform_class.F90. ENDIF IF (this%trans_type == 'polyinter') THEN IF (this%poly%arraysize <= 0) THEN CALL l4f_category_log(this%category,l4f_error, & "polyinter: poly parameter missing or empty") CALL raise_fatal_error() ENDIF ENDIF IF (this%trans_type == 'stencilinter') THEN IF (.NOT.c_e(this%area_info%radius)) THEN CALL l4f_category_log(this%category,l4f_error, & "stencilinter: radius parameter missing") CALL raise_fatal_error() ENDIF ENDIF IF (this%sub_type == 'average' .OR. this%sub_type == 'stddev' & .OR. this%sub_type == 'stddevnm1') THEN this%stat_info%percentile = rmiss ELSE IF (this%sub_type == 'max') THEN this%stat_info%percentile = 101. ELSE IF (this%sub_type == 'min') THEN this%stat_info%percentile = -1. ELSE IF (this%sub_type == 'percentile') THEN IF (.NOT.c_e(this%stat_info%percentile)) THEN CALL l4f_category_log(this%category,l4f_error,trim(this%trans_type)// & ':percentile: percentile value not provided') CALL raise_fatal_error() ELSE IF (this%stat_info%percentile >= 100.) THEN this%sub_type = 'max' ELSE IF (this%stat_info%percentile <= 0.) THEN this%sub_type = 'min' ENDIF ELSE IF (this%sub_type == 'frequency') THEN IF (.NOT.c_e(this%interval_info%gt) .AND. .NOT.c_e(this%interval_info%gt)) THEN CALL l4f_category_log(this%category,l4f_error,trim(this%trans_type)// & ':frequency: lower and/or upper limit not provided') CALL raise_fatal_error() ENDIF ELSE CALL sub_type_error() RETURN ENDIF ELSE IF (this%trans_type == 'maskgen')THEN IF (this%sub_type == 'poly') THEN IF (this%poly%arraysize <= 0) THEN CALL l4f_category_log(this%category,l4f_error,"maskgen:poly poly parameter missing or empty") CALL raise_fatal_error() ENDIF ELSE IF (this%sub_type == 'grid') THEN ! nothing to do for now ELSE CALL sub_type_error() RETURN ENDIF ELSE IF (this%trans_type == 'vertint') THEN IF (this%vertint%input_levtype == vol7d_level_miss) THEN CALL l4f_category_log(this%category,l4f_error, & 'vertint parameter input_levtype not provided') CALL raise_fatal_error() ENDIF IF (this%vertint%output_levtype == vol7d_level_miss) THEN CALL l4f_category_log(this%category,l4f_error, & 'vertint parameter output_levtype not provided') CALL raise_fatal_error() ENDIF IF (this%sub_type == 'linear' .OR. this%sub_type == 'linearsparse') THEN ! nothing to do here ELSE CALL sub_type_error() RETURN ENDIF ELSE IF (this%trans_type == 'metamorphosis') THEN IF (this%sub_type == 'all') THEN ! nothing to do here ELSE IF (this%sub_type == 'coordbb')THEN IF (c_e(this%rect_coo%ilon) .AND. c_e(this%rect_coo%ilat) .AND. & c_e(this%rect_coo%flon) .AND. c_e(this%rect_coo%flat)) THEN ! coordinates given ELSE CALL l4f_category_log(this%category,l4f_error,"metamorphosis: coordbb parameters missing") CALL raise_fatal_error() ENDIF ELSE IF (this%sub_type == 'poly')THEN IF (this%poly%arraysize <= 0) THEN CALL l4f_category_log(this%category,l4f_error,"metamorphosis:poly: poly parameter missing or empty") CALL raise_fatal_error() ENDIF ELSE IF (this%sub_type == 'mask' .OR. this%sub_type == 'maskvalid' .OR. & this%sub_type == 'maskinvalid' .OR. this%sub_type == 'setinvalidto' .OR. & this%sub_type == 'settoinvalid' .OR. this%sub_type == 'lemaskinvalid' .OR. & this%sub_type == 'ltmaskinvalid' .OR. this%sub_type == 'gemaskinvalid' .OR. & this%sub_type == 'gtmaskinvalid' .OR. this%sub_type == 'index') THEN ! nothing to do here ELSE CALL sub_type_error() RETURN ENDIF ELSE CALL trans_type_error() RETURN ENDIF CONTAINS SUBROUTINE sub_type_error() CALL l4f_category_log(this%category, l4f_error, trim(this%trans_type) & //': sub_type '//trim(this%sub_type)//' is not defined') CALL raise_fatal_error() END SUBROUTINE sub_type_error SUBROUTINE trans_type_error() CALL l4f_category_log(this%category, l4f_error, 'trans_type '//this%trans_type & //' is not defined') CALL raise_fatal_error() END SUBROUTINE trans_type_error END SUBROUTINE transform_init !> Destructor of \a tranform_def object. !! It releases any memory and data associated to the \a transform_def !! object \a this, the logger category will be deleted too. SUBROUTINE transform_delete(this) TYPE(transform_def),INTENT(inout) :: this !< transformation object this%trans_type=cmiss this%sub_type=cmiss this%rect_ind%ix=imiss this%rect_ind%iy=imiss this%rect_ind%fx=imiss this%rect_ind%fy=imiss this%rect_coo%ilon=dmiss this%rect_coo%ilat=dmiss this%rect_coo%flon=dmiss this%rect_coo%flat=dmiss this%box_info%npx=imiss this%box_info%npy=imiss this%extrap=.false. !chiudo il logger CALL l4f_category_delete(this%category) END SUBROUTINE transform_delete !> Method for returning the contents of the object. SUBROUTINE transform_get_val(this, time_definition, trans_type, sub_type, & input_levtype, output_levtype) type(transform_def),intent(in) :: this !< object to examine INTEGER,INTENT(out),OPTIONAL :: time_definition !< 0=time is reference time, 1=time is validity time CHARACTER(len=*),INTENT(out),OPTIONAL :: trans_type !< type of transformation CHARACTER(len=*),INTENT(out),OPTIONAL :: sub_type !< subtype of transformation TYPE(vol7d_level),INTENT(out),OPTIONAL :: input_levtype !< type of vertical level of input data (only type of first and second surface are used, level values are ignored) TYPE(vol7d_level),INTENT(out),OPTIONAL :: output_levtype !< type of vertical level of output data (only type of first and second surface are used, level values are ignored) IF (PRESENT(time_definition)) time_definition=this%time_definition IF (PRESENT(trans_type)) trans_type = this%trans_type IF (PRESENT(sub_type)) sub_type = this%sub_type IF (PRESENT(input_levtype)) input_levtype = this%vertint%input_levtype IF (PRESENT(output_levtype)) output_levtype = this%vertint%output_levtype END SUBROUTINE transform_get_val !> Constructor for a \a grid_transform object, defining a particular !! vertical transformation. !! It defines an object describing a transformation involving !! operations on the vertical direction only, thus it applies in the !! same way to grid-to-grid and to sparse points-to-sparse points !! transformations; the abstract type of the transformation is !! described in the transformation object \a trans (type !! transform_def) which must have been properly initialised. The !! additional information required here is the list of input and !! output levels and an optional 3-d field indicating the vertical !! coordinates of the input dataset. !! !! The input level list \a lev_in is a 1-d array of \a vol7d_level !! type, describing the vertical coordinate system of the whole input !! dataset, only the vertical levels or layers matching the level type !! indicated when initialising the transformation object \a trans will !! be used, the others will be discarded in the vertical !! transformation. However the relevant vertical levels must be !! contiguous and sorted accordingly to the default \a vol7d_level !! sort order. The argument \a lev_out describes the vertical !! coordinate system of the output dataset. A particular case to be !! considered is when \a SIZE(lev_out)==0, this means that the output !! coordinates have to be computed automatically in the current !! subroutine, this is supported only for hybrid levels/layers. !! !! When the input and output level types are different, the \a !! coord_3d_in array must be provided, indicating the vertical !! coordinate of every input grid point expressed in terms of the !! output vertical coordinate system (e.g. height if interpolating to !! constant height levels or pressure if interpolating to isobaric !! levels). This array must contain, in the 3rd vertical dimension, !! only the those levels/layers of \a lev_in that are actually used !! for interpolation. The storage space of \a coord_3d_in is "stolen" !! by this method, so the array will appear as unallocated and !! unusable to the calling procedure after return from this !! subroutine. !! !! Layers in the grib2 sense (i.e. layers between two surfaces) can be !! handled by this class only when the upper and lower surfaces are of !! the same type; in these cases the coordinate assigned to every !! layer fro interpolation is the average (or log-average in case of !! isobaric surfaces) between the coordinates of the corresponding !! upper and lower surfaces. SUBROUTINE grid_transform_levtype_levtype_init(this, trans, lev_in, lev_out, & coord_3d_in, categoryappend) TYPE(grid_transform),INTENT(out) :: this !< grid_transformation object TYPE(transform_def),INTENT(in) :: trans !< transformation object TYPE(vol7d_level),INTENT(in) :: lev_in(:) !< vol7d_level from input object TYPE(vol7d_level),INTENT(in) :: lev_out(:) !< vol7d_level object defining target vertical grid REAL,INTENT(inout),OPTIONAL,ALLOCATABLE :: coord_3d_in(:,:,:) !< vertical coordinates of each input point in target reference system CHARACTER(len=*),INTENT(in),OPTIONAL :: categoryappend !< append this suffix to log4fortran namespace category DOUBLE PRECISION :: coord_in(SIZE(lev_in)) DOUBLE PRECISION,ALLOCATABLE :: coord_out(:) LOGICAL :: mask_in(SIZE(lev_in)) LOGICAL,ALLOCATABLE :: mask_out(:) LOGICAL :: dolog INTEGER :: i, j, icache, inused, istart, iend, ostart, oend CALL grid_transform_init_common(this, trans, categoryappend) #ifdef DEBUG CALL l4f_category_log(this%category, l4f_debug, "grid_transform vertint") #endif IF (this%trans%trans_type == 'vertint') THEN IF (c_e(trans%vertint%input_levtype%level2) .AND. & trans%vertint%input_levtype%level1 /= trans%vertint%input_levtype%level2) THEN CALL l4f_category_log(this%category, l4f_error, & 'vertint: input upper and lower surface must be of the same type, '// & t2c(trans%vertint%input_levtype%level1)//'/='// & t2c(trans%vertint%input_levtype%level2)) CALL raise_error() RETURN ENDIF IF (c_e(trans%vertint%output_levtype%level2) .AND. & trans%vertint%output_levtype%level1 /= trans%vertint%output_levtype%level2) THEN CALL l4f_category_log(this%category, l4f_error, & 'vertint: output upper and lower surface must be of the same type'// & t2c(trans%vertint%output_levtype%level1)//'/='// & t2c(trans%vertint%output_levtype%level2)) CALL raise_error() RETURN ENDIF mask_in(:) = (lev_in(:)%level1 == trans%vertint%input_levtype%level1) .AND. & (lev_in(:)%level2 == trans%vertint%input_levtype%level2) CALL make_vert_coord(lev_in, mask_in, coord_in, dolog) this%innz = SIZE(lev_in) istart = firsttrue(mask_in) iend = lasttrue(mask_in) inused = iend - istart + 1 IF (inused /= count(mask_in)) THEN CALL l4f_category_log(this%category, l4f_error, & 'grid_transform_levtype_levtype_init: input levels badly sorted '//& t2c(inused)//'/'//t2c(count(mask_in))) CALL raise_error() RETURN ENDIF this%levshift = istart-1 this%levused = inused IF (trans%vertint%input_levtype%level1 /= trans%vertint%output_levtype%level1) THEN #ifdef DEBUG CALL l4f_category_log(this%category, l4f_debug, & 'vertint: different input and output level types '// & t2c(trans%vertint%input_levtype%level1)//' '// & t2c(trans%vertint%output_levtype%level1)) #endif ALLOCATE(mask_out(SIZE(lev_out)), this%vcoord_out(SIZE(lev_out))) mask_out(:) = (lev_out(:)%level1 == trans%vertint%output_levtype%level1) .AND. & (lev_out(:)%level2 == trans%vertint%output_levtype%level2) CALL make_vert_coord(lev_out, mask_out, this%vcoord_out, dolog) this%outnz = SIZE(mask_out) DEALLOCATE(mask_out) IF (.NOT.PRESENT(coord_3d_in)) THEN CALL l4f_category_log(this%category, l4f_warn, & 'vertint: different input and output level types & &and no coord_3d_in, expecting vert. coord. in volume') this%dolog = dolog ! a little bit dirty, I must compute log later ELSE IF (SIZE(coord_3d_in,3) /= inused) THEN CALL l4f_category_log(this%category, l4f_error, & 'vertint: vertical size of coord_3d_in (vertical coordinate) & &different from number of input levels suitable for interpolation') CALL l4f_category_log(this%category, l4f_error, & 'coord_3d_in: '//t2c(SIZE(coord_3d_in,3))// & ', input levels for interpolation: '//t2c(inused)) CALL raise_error() RETURN