v7d_v7d_transform v7d_volgrid6d_transform volgrid6d_transform volgrid6d_v7d_transform volgrid6dv_transform volgrid6dv_v7d_transform

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  ) 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 2786 del file volgrid6d_class.F90. 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 subroutine vg6d_wind_unrot(this) type(volgrid6d) :: this 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 subroutine vg6d_wind_rot(this) type(volgrid6d) :: this 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 subroutine vg6d_c2a (this) TYPE(volgrid6d),INTENT(inout) :: this(:) integer :: ngrid,igrid,jgrid,ugrid,vgrid,tgrid doubleprecision :: xmin, xmax, ymin, ymax doubleprecision :: xmin_t, xmax_t, ymin_t, ymax_t doubleprecision :: step_lon_t,step_lat_t character(len=80) :: type_t,type TYPE(griddim_def):: griddim_t