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

transform_init() subroutine grid_transform_class::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  ) private 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 667 del file grid_transform_class.F90. 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 SUBROUTINE transform_delete(this) TYPE(transform_def),INTENT(inout) :: this 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 SUBROUTINE transform_get_val(this, time_definition, trans_type, sub_type, & input_levtype, output_levtype) type(transform_def),intent(in) :: this INTEGER,INTENT(out),OPTIONAL :: time_definition CHARACTER(len=*),INTENT(out),OPTIONAL :: trans_type CHARACTER(len=*),INTENT(out),OPTIONAL :: sub_type TYPE(vol7d_level),INTENT(out),OPTIONAL :: input_levtype TYPE(vol7d_level),INTENT(out),OPTIONAL :: output_levtype 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 SUBROUTINE grid_transform_levtype_levtype_init(this, trans, lev_in, lev_out, & coord_3d_in, categoryappend) TYPE(grid_transform),INTENT(out) :: this TYPE(transform_def),INTENT(in) :: trans TYPE(vol7d_level),INTENT(in) :: lev_in(:) TYPE(vol7d_level),INTENT(in) :: lev_out(:) REAL,INTENT(inout),OPTIONAL,ALLOCATABLE :: coord_3d_in(:,:,:) CHARACTER(len=*),INTENT(in),OPTIONAL :: categoryappend 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 ENDIF CALL move_alloc(coord_3d_in, this%coord_3d_in) ! steal allocation IF (dolog) THEN WHERE(c_e(this%coord_3d_in) .AND. this%coord_3d_in > 0.0) this%coord_3d_in = log(this%coord_3d_in) ELSE WHERE this%coord_3d_in = rmiss END WHERE ENDIF