libsim  Versione 7.2.6

◆ v7d_volgrid6d_transform()

subroutine v7d_volgrid6d_transform ( type(transform_def), intent(in)  this,
type(griddim_def), intent(in), optional  griddim,
type(vol7d), intent(inout)  vol7d_in,
type(volgrid6d), intent(out)  volgrid6d_out,
character(len=*), intent(in), optional  networkname,
type(grid_id), intent(in), optional  gaid_template,
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]thisobject specifying the abstract transformation
[in]griddimgriddim specifying the output grid (required by most transformation types)
[in,out]vol7d_inobject to be transformed, it is not modified, despite the INTENT(inout)
[out]volgrid6d_outtransformed object, it does not require initialisation
[in]networknameselect the network to be processed from the vol7d input object, default the first network
[in]gaid_templatethe template (typically grib_api) to be associated with output data, it also helps in improving variable conversion
[in]categoryappendappend this suffix to log4fortran namespace category

Definizione alla linea 2613 del file volgrid6d_class.F90.

2615  var_coord_in = -1
2616  IF (PRESENT(vol7d_coord_in)) THEN
2617  IF (ASSOCIATED(vol7d_coord_in%voldatir) .AND. &
2618  ASSOCIATED(vol7d_coord_in%dativar%r)) THEN
2619 
2620 ! strictly 1 time, 1 timerange and 1 network
2621  IF (SIZE(vol7d_coord_in%voldatir,2) /= 1 .OR. &
2622  SIZE(vol7d_coord_in%voldatir,4) /= 1 .OR. &
2623  SIZE(vol7d_coord_in%voldatir,6) /= 1) THEN
2624  CALL l4f_log(l4f_error, &
2625  'volume providing constant input vertical coordinate must have &
2626  &only 1 time, 1 timerange and 1 network')
2627  CALL raise_error()
2628  RETURN
2629  ENDIF
2630 
2631 ! search for variable providing vertical coordinate
2632  CALL get_val(this, output_levtype=output_levtype)
2633  vcoord_var = vol7d_var_new(vol7d_level_to_var(output_levtype))
2634  IF (.NOT.c_e(vcoord_var)) THEN
2635  CALL l4f_log(l4f_error, &
2636  'requested output level type '//t2c(output_levtype%level1)// &
2637  ' does not correspond to any known physical variable for &
2638  &providing vertical coordinate')
2639  CALL raise_error()
2640  RETURN
2641  ENDIF
2642 
2643  var_coord_in = index(vol7d_coord_in%dativar%r, vcoord_var)
2644 
2645  IF (var_coord_in <= 0) THEN
2646  CALL l4f_log(l4f_error, &
2647  'volume providing constant input vertical coordinate contains no &
2648  &real variables matching output level type '//t2c(output_levtype%level1))
2649  CALL raise_error()
2650  RETURN
2651  ENDIF
2652  CALL l4f_log(l4f_info, &
2653  'Coordinate for vertint found in coord volume at position '// &
2654  t2c(var_coord_in))
2655 
2656 ! check vertical coordinate system
2657  CALL get_val(this, input_levtype=input_levtype)
2658  mask_in = & ! implicit allocation
2659  (vol7d_coord_in%level(:)%level1 == input_levtype%level1) .AND. &
2660  (vol7d_coord_in%level(:)%level2 == input_levtype%level2)
2661  ulstart = firsttrue(mask_in)
2662  ulend = lasttrue(mask_in)
2663  IF (ulstart == 0 .OR. ulend == 0) THEN
2664  CALL l4f_log(l4f_error, &
2665  'coordinate file does not contain levels of type '// &
2666  t2c(input_levtype%level1)//'/'//t2c(input_levtype%level2)// &
2667  ' specified for input data')
Index method.

Generated with Doxygen.