◆ volgrid6d_alloc_vol()

subroutine, public volgrid6d_alloc_vol	(	type(volgrid6d), intent(inout)	this,
		logical, intent(in), optional	ini,
		logical, intent(in), optional	inivol,
		logical, intent(in), optional	decode )

Allocate the data array of the volgrid6d object.

This method allocates the main 6-dimensional data array thisvoldati and the 4-dimensional grid_id array thisgaid with a shape dictated by the previous call(s) to vol7d_alloc(). if any descriptor (except horizontal grid) has not been allocated yet, it is allocated here with a size of 1. This method should be explicitly used only in rare cases, it is usually called implicitly through the import interface.

Parametri

[in,out]	this	object whose decriptors should be allocated
[in]	ini	if provided and .TRUE., for each dimension descriptor not yet allocated and allocated here the constructor is called without extra parameters, thus initializing the element as missing value
[in]	inivol	if provided and .FALSE., the allocated volumes will not be initialized to missing values
[in]	decode	if provided and .TRUE., the thisvoldati volume is allocated, otherwise only thisgaid will be allocated

Definizione alla linea 520 del file volgrid6d_class.F90.

ELSE
  CALL grid_id_encode_data(this%gaid(ilevel,itime,itimerange,ivar), voldati)
ENDIF
 
END SUBROUTINE volgrid_set_vol_2d
 
 
!> Reset a 3-d x-y-z slice of a volume after the data have been modified.
!! This method works both with volumes having allocated and
!! non-allocated this%voldati array, and it updates the requested
!! slice.  In case \a this%voldati is already allocated, this is a
!! no-operation while in the other case this method encodes the field
!! provided into the grid_id object on file or in memory. Since this
!! method may be called many times by a program, it is optimized for
!! speed and it does not make any check about the matching size of the
!! field and the array or about the allocation status of \a this, so
!! it should be called only when everything has been checked to be in
!! good shape.
SUBROUTINE volgrid_set_vol_3d(this, itime, itimerange, ivar, voldati)
TYPE(volgrid6d),INTENT(inout) :: this !< object in which slice has to be updated
INTEGER,INTENT(in) :: itime !< index of time level of the slice
INTEGER,INTENT(in) :: itimerange !< index of timerange of the slice
INTEGER,INTENT(in) :: ivar !< index of physical variable of the slice
REAL,INTENT(in) :: voldati(:,:,:) !< updated values of the slice
 
INTEGER :: ilevel
 
IF (ASSOCIATED(this%voldati)) THEN
  RETURN
ELSE
!$OMP PARALLEL DEFAULT(SHARED)
!$OMP MASTER
  DO ilevel = 1, SIZE(this%level)
!$OMP TASK FIRSTPRIVATE(ilevel)
    CALL grid_id_encode_data(this%gaid(ilevel,itime,itimerange,ivar), &
     voldati(:,:,ilevel))
!$OMP END TASK
  ENDDO
!$OMP END MASTER
!$OMP END PARALLEL
ENDIF
 
END SUBROUTINE volgrid_set_vol_3d
 
 
!> Destructor, it releases every information and memory buffer
!! associated with the object. It should be called also for objects
!! crated through the \a import interface.
SUBROUTINE volgrid6d_delete(this)
TYPE(volgrid6d),INTENT(inout) :: this
 
INTEGER :: i, ii, iii, iiii
 
#ifdef DEBUG
call l4f_category_log(this%category,l4f_debug,"delete")
#endif
 
if (associated(this%gaid))then
 
  DO i=1 ,SIZE(this%gaid,1)
    DO ii=1 ,SIZE(this%gaid,2)
      DO iii=1 ,SIZE(this%gaid,3)
        DO iiii=1 ,SIZE(this%gaid,4)
          CALL delete(this%gaid(i,ii,iii,iiii))
        ENDDO
      ENDDO
    ENDDO
  ENDDO
  DEALLOCATE(this%gaid)
 
end if
 
call delete(this%griddim)
 
!  call delete(this%time)
!  call delete(this%timerange)
!  call delete(this%level)
!  call delete(this%var)
 