display_volgrid6d display_volgrid6dv export_to_gridinfo export_to_gridinfov export_to_gridinfovv import_from_gridinfo import_from_gridinfovv v7d_v7d_transform v7d_volgrid6d_transform vg6d_c2a vg6d_reduce vg6d_rounding vg6d_wind_rot vg6d_wind_unrot vg6dv_rounding volgrid6d_alloc volgrid6d_alloc_vol volgrid6d_delete volgrid6d_export_to_file volgrid6d_import_from_file volgrid6d_init volgrid6d_read_from_file volgrid6d_transform volgrid6d_v7d_transform volgrid6d_v7d_transform_compute volgrid6d_write_on_file volgrid6dv_delete volgrid6dv_transform volgrid6dv_v7d_transform volgrid_get_vol_2d volgrid_get_vol_3d volgrid_set_vol_2d volgrid_set_vol_3d

v7d_v7d_transform() subroutine v7d_v7d_transform ( type(transform_def), intent(in) this, type(vol7d), intent(inout) vol7d_in, type(vol7d), intent(out) vol7d_out, type(vol7d), intent(in), optional v7d, real, dimension(:), intent(in), optional maskbounds, type(vol7d_level), dimension(:), intent(in), optional, target lev_out, type(vol7d), intent(in), optional vol7d_coord_in, character(len=*), intent(in), optional categoryappend ) private 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. The success of the transformation can be checked with the c_e method: c_e(vol7d_out). Parametri [in] this object specifying the abstract transformation [in,out] vol7d_in object to be transformed, it is not modified, despite the INTENT(inout) [out] vol7d_out transformed object, it does not require initialisation [in] v7d object containing a list of points over which transformation has to be done (required by some transformation types) [in] maskbounds array of boundary values for defining a subset of valid points where the values of maskgrid are within the first and last value of maskbounds (for transformation type 'metamorphosis:maskfill') [in] lev_out vol7d_level object defining target vertical grid, for vertical interpolations [in] vol7d_coord_in object providing time constant input vertical coordinate for some kind of vertical interpolations [in] categoryappend append this suffix to log4fortran namespace category Definizione alla linea 2774 del file volgrid6d_class.F90. CALL get_val(grid_trans, point_mask=point_mask, output_point_index=point_index) IF (ALLOCATED(point_index)) THEN CALL vol7d_alloc(vol7d_out, nanavari=1) CALL init(vol7d_out%anavar%i(1), 'B01192') ENDIF CALL vol7d_alloc_vol(vol7d_out) IF (ALLOCATED(point_index)) THEN DO inetwork = 1, SIZE(vol7d_in%network) vol7d_out%volanai(:,1,inetwork) = point_index(:) ENDDO ENDIF CALL compute(grid_trans, vol7d_in, vol7d_out) IF (ALLOCATED(point_mask)) THEN ! keep full ana IF (SIZE(point_mask) /= SIZE(vol7d_in%ana)) THEN CALL l4f_log(l4f_warn, & 'v7d_v7d_transform: inconsistency in point size: '//t2c(SIZE(point_mask)) & //':'//t2c(SIZE(vol7d_in%ana))) ELSE #ifdef DEBUG CALL l4f_log(l4f_debug, 'v7d_v7d_transform: merging ana from in to out') #endif CALL vol7d_copy(vol7d_in, vol7d_tmpana, & lana=point_mask, lnetwork=(/.true./), & ltime=(/.false./), ltimerange=(/.false./), llevel=(/.false./)) CALL vol7d_append(vol7d_out, vol7d_tmpana) ENDIF ENDIF ELSE CALL l4f_log(l4f_error, 'v7d_v7d_transform: transformation not valid') CALL raise_error() ENDIF ENDIF CALL delete (grid_trans) IF (.NOT. PRESENT(v7d)) CALL delete(v7d_locana) END SUBROUTINE v7d_v7d_transform !> Unrotate the wind components. !! It converts u and v components of vector quantities relative to the !! defined grid in the direction of increasing x and y coordinates to !! u and v components relative to easterly and notherly direction. The !! original fields are overwritten. !! \todo Check and correct wind component flag (to be moved in !! griddim_def?) subroutine vg6d_wind_unrot(this) type(volgrid6d) :: this !< object containing wind to be unrotated integer :: component_flag call get_val(this%griddim,component_flag=component_flag) if (component_flag == 1) then call l4f_category_log(this%category,l4f_info, & "unrotating vector components") call vg6d_wind__un_rot(this,.false.) call set_val(this%griddim,component_flag=0) else call l4f_category_log(this%category,l4f_info, & "no need to unrotate vector components") end if end subroutine vg6d_wind_unrot !> Rotate the wind components. !! It converts u and v components of vector quantities !! relative to easterly and notherly direction to !! defined grid in the direction of increasing x and y coordinates. !! The original fields are overwritten. subroutine vg6d_wind_rot(this) type(volgrid6d) :: this !< object containing wind to be rotated integer :: component_flag call get_val(this%griddim,component_flag=component_flag) if (component_flag == 0) then call l4f_category_log(this%category,l4f_info, & "rotating vector components") call vg6d_wind__un_rot(this,.true.) call set_val(this%griddim,component_flag=1) else call l4f_category_log(this%category,l4f_info, & "no need to rotate vector components") end if end subroutine vg6d_wind_rot ! Generic UnRotate the wind components. SUBROUTINE vg6d_wind__un_rot(this,rot) TYPE(volgrid6d) :: this ! object containing wind to be (un)rotated LOGICAL :: rot ! if .true. rotate else unrotate INTEGER :: i, j, k, l, a11, a12, a21, a22, stallo double precision,pointer :: rot_mat(:,:,:) real,allocatable :: tmp_arr(:,:) REAL,POINTER :: voldatiu(:,:), voldativ(:,:) INTEGER,POINTER :: iu(:), iv(:) IF (.NOT.ASSOCIATED(this%var)) THEN CALL l4f_category_log(this%category, l4f_error, & "trying to unrotate an incomplete volgrid6d object") CALL raise_fatal_error() ! RETURN ENDIF CALL volgrid6d_var_hor_comp_index(this%var, iu, iv) IF (.NOT.ASSOCIATED(iu)) THEN CALL l4f_category_log(this%category,l4f_error, & "unrotation impossible") CALL raise_fatal_error() ! RETURN ENDIF ! Temporary workspace ALLOCATE(tmp_arr(this%griddim%dim%nx, this%griddim%dim%ny),stat=stallo) IF (stallo /= 0) THEN CALL l4f_category_log(this%category, l4f_fatal, "allocating memory") CALL raise_fatal_error() ENDIF ! allocate once for speed IF (.NOT.ASSOCIATED(this%voldati)) THEN ALLOCATE(voldatiu(this%griddim%dim%nx, this%griddim%dim%ny), & voldativ(this%griddim%dim%nx, this%griddim%dim%ny)) ENDIF CALL griddim_unproj(this%griddim) CALL wind_unrot(this%griddim, rot_mat) a11=1 if (rot)then a12=2 a21=3 else a12=3 a21=2 end if a22=4 DO l = 1, SIZE(iu) DO k = 1, SIZE(this%timerange) DO j = 1, SIZE(this%time) DO i = 1, SIZE(this%level) ! get data CALL volgrid_get_vol_2d(this, i, j, k, iu(l), voldatiu) CALL volgrid_get_vol_2d(this, i, j, k, iv(l), voldativ) ! convert units forward ! CALL compute(conv_fwd(iu(l)), voldatiu) ! CALL compute(conv_fwd(iv(l)), voldativ) ! multiply wind components by rotation matrix WHERE(voldatiu /= rmiss .AND. voldativ /= rmiss) tmp_arr(:,:) = real(voldatiu(:,:)*rot_mat(:,:,a11) + & voldativ(:,:)*rot_mat(:,:,a12)) voldativ(:,:) = real(voldatiu(:,:)*rot_mat(:,:,a21) + & voldativ(:,:)*rot_mat(:,:,a22)) voldatiu(:,:) = tmp_arr(:,:) END WHERE ! convert units backward ! CALL uncompute(conv_fwd(iu(l)), voldatiu) ! CALL uncompute(conv_fwd(iv(l)), voldativ) ! put data CALL volgrid_set_vol_2d(this, i, j, k, iu(l), voldatiu) CALL volgrid_set_vol_2d(this, i, j, k, iv(l), voldativ) ENDDO ENDDO ENDDO ENDDO IF (.NOT.ASSOCIATED(this%voldati)) THEN DEALLOCATE(voldatiu, voldativ) ENDIF DEALLOCATE(rot_mat, tmp_arr, iu, iv) END SUBROUTINE vg6d_wind__un_rot !!$ try to understand the problem: !!$ !!$ case: !!$ !!$ 1) we have only one volume: we have to provide the direction of shift !!$ compute H and traslate on it !!$ 2) we have two volumes: !!$ 1) volume U and volume V: compute H and traslate on it !!$ 2) volume U/V and volume H : translate U/V on H !!$ 3) we have tree volumes: translate U and V on H !!$ !!$ strange cases: !!$ 1) do not have U in volume U !!$ 2) do not have V in volume V !!$ 3) we have others variables more than U and V in volumes U e V !!$ !!$ !!$ so the steps are: !!$ 1) find the volumes !!$ 2) define or compute H grid !!$ 3) trasform the volumes in H !!$ N.B. !!$ case 1) for only one vg6d (U or V) is not managed, but !!$ the not pubblic subroutines will work but you have to know what you want to do !> Convert grids type C to type A. !! Use this to interpolate data from grid type C to grid type A !! Grids type are defined by Arakawa. !! !! We need to find these types of area in a vg6d array !! T area of points with temterature etc. !! U area of points with u components of winds !! V area of points with v components of winds !! !! this method works if it finds !! two volumes: !! 1) volume U and volume V: compute H and traslate on it !! 2) volume U/V and volume H : translate U/V on H !! three volumes: translate U and V on H !! !! try to work well on more datasets at once subroutine vg6d_c2a (this)