◆ volgrid_set_vol_3d()

subroutine, public volgrid_set_vol_3d	(	type(volgrid6d), intent(inout)	this,
		integer, intent(in)	itime,
		integer, intent(in)	itimerange,
		integer, intent(in)	ivar,
		real, dimension(:,:,:), intent(in)	voldati )

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 thisvoldati array, and it updates the requested slice. In case thisvoldati 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 this, so it should be called only when everything has been checked to be in good shape.

Parametri

[in,out]	this	object in which slice has to be updated
[in]	itime	index of time level of the slice
[in]	itimerange	index of timerange of the slice
[in]	ivar	index of physical variable of the slice
[in]	voldati	updated values of the slice

Definizione alla linea 727 del file volgrid6d_class.F90.

integer,intent(out),optional :: tarray(8) !< vettore come definito da "date_and_time" della data di scrittura del volume
 
integer :: ntime, ntimerange, nlevel, nvar
 
character(len=254) :: ldescription,lfilename,arg
integer :: ltarray(8),lunit
logical :: opened,exist
 
#ifdef DEBUG
call l4f_category_log(this%category,l4f_debug,"read from file")
#endif
 
call getarg(0,arg)
 
if (.not. present(unit))then
  lunit=getunit()
else
  if (unit==0)then
    lunit=getunit()
    unit=lunit
  else
    lunit=unit
  end if