|
◆ 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.
2337 if ( associated(vol7d_in%time)) then
2338 ntime= size(vol7d_in%time)
2339 volgrid6d_out%time=vol7d_in%time
2342 if ( associated(vol7d_in%timerange)) then
2343 ntimerange= size(vol7d_in%timerange)
2344 volgrid6d_out%timerange=vol7d_in%timerange
2347 if ( associated(vol7d_in%level)) then
2348 nlevel= size(vol7d_in%level)
2349 volgrid6d_out%level=vol7d_in%level
2352 if ( associated(vol7d_in%dativar%r)) then
2353 nvar= size(vol7d_in%dativar%r)
2354 CALL varbufr2vargrib(vol7d_in%dativar%r, volgrid6d_out%var, c_func, gaid_template)
2357 nana= SIZE(vol7d_in%voldatir, 1)
2359 IF (.NOT. ASSOCIATED(volgrid6d_out%voldati)) THEN
2360 ALLOCATE(voldatiout(volgrid6d_out%griddim%dim%nx, volgrid6d_out%griddim%dim%ny, &
2365 DO itimerange=1,ntimerange
2369 IF ( PRESENT(gaid_template)) THEN
2370 DO ilevel = 1, nlevel
2371 IF (any(c_e(vol7d_in%voldatir(:,itime,ilevel,itimerange,ivar,inetwork)))) THEN
2372 CALL copy(gaid_template, volgrid6d_out%gaid(ilevel,itime,itimerange,ivar))
2374 volgrid6d_out%gaid(ilevel,itime,itimerange,ivar) = grid_id_new()
2380 IF ( ASSOCIATED(volgrid6d_out%voldati)) &
2381 CALL volgrid_get_vol_3d(volgrid6d_out, itime, itimerange, ivar, &
2384 CALL compute(this, &
2385 vol7d_in%voldatir(:,itime,:,itimerange,ivar,inetwork), voldatiout, &
2386 vol7d_in%dativar%r(ivar))
2388 IF ( ASSOCIATED(c_func)) THEN
2389 CALL compute(c_func(ivar), voldatiout(:,:,:))
2392 CALL volgrid_set_vol_3d(volgrid6d_out, itime, itimerange, ivar, &
2399 IF (.NOT. ASSOCIATED(volgrid6d_out%voldati)) THEN
2400 DEALLOCATE(voldatiout)
2402 IF ( ASSOCIATED(c_func)) THEN
2406 END SUBROUTINE v7d_volgrid6d_transform_compute
2415 SUBROUTINE v7d_volgrid6d_transform(this, griddim, vol7d_in, volgrid6d_out, &
2416 networkname, gaid_template, categoryappend)
2417 TYPE(transform_def), INTENT(in) :: this
2418 TYPE(griddim_def), INTENT(in), OPTIONAL :: griddim
2420 TYPE(vol7d), INTENT(inout) :: vol7d_in
2421 TYPE(volgrid6d), INTENT(out) :: volgrid6d_out
2422 CHARACTER(len=*), OPTIONAL, INTENT(in) :: networkname
2423 TYPE(grid_id), OPTIONAL, INTENT(in) :: gaid_template
2424 CHARACTER(len=*), INTENT(in), OPTIONAL :: categoryappend
2426 type(grid_transform) :: grid_trans
2427 integer :: ntime, ntimerange, nlevel, nvar
2433 CALL vol7d_alloc_vol(vol7d_in)
2434 ntime= SIZE(vol7d_in%time)
2435 ntimerange= SIZE(vol7d_in%timerange)
2436 nlevel= SIZE(vol7d_in%level)
2438 if ( associated(vol7d_in%dativar%r)) nvar= size(vol7d_in%dativar%r)
2441 CALL l4f_log(l4f_error, &
2442 "trying to transform a vol7d object incomplete or without real variables")
2443 CALL init(volgrid6d_out)
2448 CALL init(grid_trans, this, vol7d_in, griddim, categoryappend=categoryappend)
2449 CALL init(volgrid6d_out, griddim, time_definition=vol7d_in%time_definition, &
2450 categoryappend=categoryappend)
2452 IF (c_e(grid_trans)) THEN
2454 CALL volgrid6d_alloc(volgrid6d_out, griddim%dim, ntime=ntime, nlevel=nlevel, &
|