readsdv.f90

Go to the documentation of this file.

! *********************************************************************
! * Rocstar Simulation Suite                                          *
! * Copyright@2015, Illinois Rocstar LLC. All rights reserved.        *
! *                                                                   *
! * Illinois Rocstar LLC                                              *
! * Champaign, IL                                                     *
! * www.illinoisrocstar.com                                           *
! * sales@illinoisrocstar.com                                         *
! *                                                                   *
! * License: See LICENSE file in top level of distribution package or *
! * http://opensource.org/licenses/NCSA                               *
! *********************************************************************
! *********************************************************************
! * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,   *
! * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES   *
! * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND          *
! * NONINFRINGEMENT.  IN NO EVENT SHALL THE CONTRIBUTORS OR           *
! * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER       *
! * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,   *
! * Arising FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE    *
! * USE OR OTHER DEALINGS WITH THE SOFTWARE.                          *
! *********************************************************************
PROGRAM readsdv

  IMPLICIT NONE

  include 'roccomf90.h'
  include 'mpif.h'

  INTEGER,PARAMETER :: pid = 1 ! will change to input late


  INTEGER :: hdl_read, hdl_obtain, hdl_all

  CHARACTER(*), PARAMETER :: overlaywin = "Overlay"

  INTEGER :: nsubn, nsubf, nn, nf

  INTEGER :: comm_self
  CHARACTER(*), PARAMETER :: prefix = "A"

  CHARACTER(LEN=5) :: sdv_material, sdv_wname
  CHARACTER(*), PARAMETER :: fname = 'A_101_sdv.hdf'

! obtain function handle ------------------------------------------------------

  print*,'Loading module'

  CALL com_init

  CALL rocin_load_module( "SDV_IN")

  print*,'before hdl_read '
  hdl_read   = com_get_function_handle( 'SDV_IN.read_windows')
  print*,'after hdl_read'
  hdl_obtain = com_get_function_handle( 'SDV_IN.obtain_attribute')

  print*,'after hdl_obtain'

! Define the base-window and sdv-window names


  sdv_material = prefix//'_sdv'
  sdv_wname = sdv_material

  CALL com_new_window( overlaywin )

!  // Read the pane from the given file. Note that the file contains both
!  // the parent and the subdivided windows. Read only the subdivided one.

  comm_self = mpi_comm_self
  CALL com_call_function( hdl_read, 4, fname, sdv_wname, &
              sdv_material, comm_self)
  hdl_all = com_get_attribute_handle( "A_sdv.all")
  CALL com_call_function( hdl_obtain, 3, hdl_all, hdl_all, pid) 

!  // Obtain number of sub-nodes, sub-faces, nodes, and faces
  
  CALL com_get_size( overlaywin//'.sn_parent_fcID', pid, nsubn)

  print*,nsubn

!!$  CALL COM_get_size( OverlayWin//'.:t3', pid, nsubf)
!!$
!!$  CALL COM_get_size( OverlayWin//'.sn_subID', pid, nn)
!!$  CALL COM_get_size( OverlayWin//'.sf_offset', pid, nf)


  CALL com_finalize

END PROGRAM readsdv