libsim  Versione 7.2.6

◆ volgrid_get_vol_3d()

subroutine, public volgrid6d_class::volgrid_get_vol_3d ( type(volgrid6d), intent(in)  this,
integer, intent(in)  itime,
integer, intent(in)  itimerange,
integer, intent(in)  ivar,
real, dimension(:,:,:), pointer  voldati 
)

Return a 3-d pointer to a x-y-z slice of a volume.

This method works both with volumes having allocated and non-allocated thisvoldati array, and it returns a pointer to a 3-d slice either from the allocated thisvoldati array or from the grid_id object on file or in memory. In the second case the pointer should be either ALLOCATE'd to the expected size or NULLIFY'ed, and if NULLIFY'ed, it is allocated within the method, thus it will have to be deallocated by the caller when not in use anymore. 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 pointer 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]thisobject from which the slice has to be retrieved
[in]itimeindex of time level of the slice
[in]itimerangeindex of timerange of the slice
[in]ivarindex of physical variable of the slice
voldatipointer to the data, if thisvoldati is already allocated, it will just point to the requested slice, otherwise it will be allocated if and only if it is nullified on entry

Definizione alla linea 662 del file volgrid6d_class.F90.

663 if (index(arg,'/',back=.true.) > 0) lfilename=lfilename(index(arg,'/',back=.true.)+1 : )
664 
665 if (present(filename))then
666  if (filename /= "")then
667  lfilename=filename
668  end if
669 end if
670 
671 if (present(filename_auto))filename_auto=lfilename
672 
673 
674 inquire(unit=lunit,opened=opened)
675 if (.not. opened) then
676  inquire(file=lfilename,exist=exist)
677  if (exist) CALL raise_error('file exist; cannot open new file')
678  if (.not.exist) open (unit=lunit,file=lfilename,form="UNFORMATTED")
679  !print *, "opened: ",lfilename
680 end if
681 
682 if (associated(this%time)) ntime=size(this%time)
683 if (associated(this%timerange)) ntimerange=size(this%timerange)
684 if (associated(this%level)) nlevel=size(this%level)
685 if (associated(this%var)) nvar=size(this%var)
686 
687 
688 write(unit=lunit)ldescription
689 write(unit=lunit)tarray
Index method.

Generated with Doxygen.