!----------------------------------------------------------------- ! FPL (Fortran Parameter List) ! Copyright (c) 2015 Santiago Badia, Alberto F. Martín, ! Javier Principe and Víctor Sande. ! All rights reserved. ! ! This library is free software; you can redistribute it and/or ! modify it under the terms of the GNU Lesser General Public ! License as published by the Free Software Foundation; either ! version 3.0 of the License, or (at your option) any later version. ! ! This library is distributed in the hope that it will be useful, ! but WITHOUT ANY WARRANTY; without even the implied warranty of ! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU ! Lesser General Public License for more details. ! ! You should have received a copy of the GNU Lesser General Public ! License along with this library. !----------------------------------------------------------------- module DimensionsWrapper5D_R4P USE DimensionsWrapper5D USE PENF, only: I4P, R4P, str, byte_size USE ErrorMessages implicit none private type, extends(DimensionsWrapper5D_t) :: DimensionsWrapper5D_R4P_t real(R4P), allocatable :: Value(:,:,:,:,:) contains private procedure, public :: Set => DimensionsWrapper5D_R4P_Set procedure, public :: Get => DimensionsWrapper5D_R4P_Get procedure, public :: GetShape => DimensionsWrapper5D_R4P_GetShape procedure, public :: GetPointer => DimensionsWrapper5D_R4P_GetPointer procedure, public :: GetPolymorphic => DimensionsWrapper5D_R4P_GetPolymorphic procedure, public :: DataSizeInBytes=> DimensionsWrapper5D_R4P_DataSizeInBytes procedure, public :: isOfDataType => DimensionsWrapper5D_R4P_isOfDataType procedure, public :: toString => DimensionsWrapper5D_R4P_toString procedure, public :: Print => DimensionsWrapper5D_R4P_Print procedure, public :: Free => DimensionsWrapper5D_R4P_Free final :: DimensionsWrapper5D_R4P_Final end type public :: DimensionsWrapper5D_R4P_t contains subroutine DimensionsWrapper5D_R4P_Final(this) !----------------------------------------------------------------- !< Final procedure of DimensionsWrapper5D !----------------------------------------------------------------- type(DimensionsWrapper5D_R4P_t), intent(INOUT) :: this !----------------------------------------------------------------- call this%Free() end subroutine subroutine DimensionsWrapper5D_R4P_Set(this, Value) !----------------------------------------------------------------- !< Set R4P Wrapper Value !----------------------------------------------------------------- class(DimensionsWrapper5D_R4P_t), intent(INOUT) :: this class(*), intent(IN) :: Value(:,:,:,:,:) integer :: err !----------------------------------------------------------------- select type (Value) type is (real(R4P)) allocate(this%Value(size(Value,dim=1), & size(Value,dim=2), & size(Value,dim=3), & size(Value,dim=4), & size(Value,dim=5)), & stat=err) this%Value = Value if(err/=0) & call msg%Error( txt='Setting Value: Allocation error ('//& str(no_sign=.true.,n=err)//')', & file=__FILE__, line=__LINE__ ) class Default call msg%Warn( txt='Setting value: Expected data type (R4P)', & file=__FILE__, line=__LINE__ ) end select end subroutine subroutine DimensionsWrapper5D_R4P_Get(this, Value) !----------------------------------------------------------------- !< Get R4P Wrapper Value !----------------------------------------------------------------- class(DimensionsWrapper5D_R4P_t), intent(IN) :: this class(*), intent(OUT) :: Value(:,:,:,:,:) integer(I4P), allocatable :: ValueShape(:) !----------------------------------------------------------------- select type (Value) type is (real(R4P)) call this%GetShape(ValueShape) if(all(ValueShape == shape(Value))) then Value = this%Value else call msg%Warn(txt='Getting value: Wrong shape ('//& str(no_sign=.true.,n=ValueShape)//'/='//& str(no_sign=.true.,n=shape(Value))//')',& file=__FILE__, line=__LINE__ ) endif class Default call msg%Warn(txt='Getting value: Expected data type (R4P)',& file=__FILE__, line=__LINE__ ) end select end subroutine subroutine DimensionsWrapper5D_R4P_GetShape(this, ValueShape) !----------------------------------------------------------------- !< Get Wrapper Value Shape !----------------------------------------------------------------- class(DimensionsWrapper5D_R4P_t), intent(IN) :: this integer(I4P), allocatable, intent(INOUT) :: ValueShape(:) !----------------------------------------------------------------- if(allocated(ValueShape)) deallocate(ValueShape) allocate(ValueShape(this%GetDimensions())) ValueShape = shape(this%Value, kind=I4P) end subroutine function DimensionsWrapper5D_R4P_GetPointer(this) result(Value) !----------------------------------------------------------------- !< Get Unlimited Polymorphic pointer to Wrapper Value !----------------------------------------------------------------- class(DimensionsWrapper5D_R4P_t), target, intent(IN) :: this class(*), pointer :: Value(:,:,:,:,:) !----------------------------------------------------------------- Value => this%Value end function subroutine DimensionsWrapper5D_R4P_GetPolymorphic(this, Value) !----------------------------------------------------------------- !< Get Unlimited Polymorphic Wrapper Value !----------------------------------------------------------------- class(DimensionsWrapper5D_R4P_t), intent(IN) :: this class(*), allocatable, intent(OUT) :: Value(:,:,:,:,:) !----------------------------------------------------------------- allocate(Value(size(this%Value,dim=1), & size(this%Value,dim=2), & size(this%Value,dim=3), & size(this%Value,dim=4), & size(this%Value,dim=5)), & source=this%Value) end subroutine subroutine DimensionsWrapper5D_R4P_Free(this) !----------------------------------------------------------------- !< Free a DimensionsWrapper5D !----------------------------------------------------------------- class(DimensionsWrapper5D_R4P_t), intent(INOUT) :: this integer :: err !----------------------------------------------------------------- if(allocated(this%Value)) then deallocate(this%Value, stat=err) if(err/=0) call msg%Error(txt='Freeing Value: Deallocation error ('// & str(no_sign=.true.,n=err)//')', & file=__FILE__, line=__LINE__ ) endif end subroutine function DimensionsWrapper5D_R4P_DataSizeInBytes(this) result(DataSizeInBytes) !----------------------------------------------------------------- !< Return the size of the stored data in bytes !----------------------------------------------------------------- class(DimensionsWrapper5D_R4P_t), intent(IN) :: this !< Dimensions wrapper 5D integer(I4P) :: DataSizeInBytes !< Size of lthe stored data in bytes !----------------------------------------------------------------- DataSizeInBytes = byte_size(this%value(1,1,1,1,1))*size(this%value) end function DimensionsWrapper5D_R4P_DataSizeInBytes function DimensionsWrapper5D_R4P_isOfDataType(this, Mold) result(isOfDataType) !----------------------------------------------------------------- !< Check if Mold and Value are of the same datatype !----------------------------------------------------------------- class(DimensionsWrapper5D_R4P_t), intent(IN) :: this !< Dimensions wrapper 5D class(*), intent(IN) :: Mold !< Mold for data type comparison logical :: isOfDataType !< Boolean flag to check if Value is of the same data type as Mold !----------------------------------------------------------------- isOfDataType = .false. select type (Mold) type is (real(R4P)) isOfDataType = .true. end select end function DimensionsWrapper5D_R4P_isOfDataType function DimensionsWrapper5D_R4P_toString(this, Separator) result(String) !----------------------------------------------------------------- !< Return the wrapper value as a string !----------------------------------------------------------------- class(DimensionsWrapper5D_R4P_t), intent(IN) :: this character(len=1), optional, intent(IN) :: Separator character(len=:), allocatable :: String character(len=1) :: Sep integer(I4P) :: idx2,idx3,idx4,idx5 !----------------------------------------------------------------- String = '' Sep = ',' if(allocated(this%Value)) then if(present(Separator)) Sep = Separator do idx5=1, size(this%Value,5) do idx4=1, size(this%Value,4) do idx3=1, size(this%Value,3) do idx2=1, size(this%Value,2) String = String // trim(str(n=this%Value(:,idx2,idx3,idx4,idx5))) // Sep enddo enddo enddo enddo String = trim(adjustl(String(:len(String)-1))) endif end function subroutine DimensionsWrapper5D_R4P_Print(this, unit, prefix, iostat, iomsg) !----------------------------------------------------------------- !< Print Wrapper !----------------------------------------------------------------- class(DimensionsWrapper5D_R4P_t), intent(IN) :: this !< DimensionsWrapper integer(I4P), intent(IN) :: unit !< Logic unit. character(*), optional, intent(IN) :: prefix !< Prefixing string. integer(I4P), optional, intent(OUT) :: iostat !< IO error. character(*), optional, intent(OUT) :: iomsg !< IO error message. character(len=:), allocatable :: prefd !< Prefixing string. integer(I4P) :: iostatd !< IO error. character(500) :: iomsgd !< Temporary variable for IO error message. !----------------------------------------------------------------- prefd = '' ; if (present(prefix)) prefd = prefix write(unit=unit,fmt='(A,$)',iostat=iostatd,iomsg=iomsgd) prefd//' Data Type = R4P'//& ', Dimensions = '//trim(str(no_sign=.true., n=this%GetDimensions()))//& ', Bytes = '//trim(str(no_sign=.true., n=this%DataSizeInBytes()))//& ', Value = ' write(unit=unit,fmt=*,iostat=iostatd,iomsg=iomsgd) this%toString() if (present(iostat)) iostat = iostatd if (present(iomsg)) iomsg = iomsgd end subroutine DimensionsWrapper5D_R4P_Print end module DimensionsWrapper5D_R4P