alloc_subset_str.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"
#include "scr17_c.inc"
#include "tabsiz_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	alloc_subset_str (subset)

Function/Subroutine Documentation

alloc_subset_str()

subroutine alloc_subset_str

(

type (subset_), dimension(nsubs)

subset

)

Definition at line 34 of file alloc_subset_str.F.

C-----------------------------------------------
C   M o d u l e s
C----------------------------------------------- 
      USE message_mod
      USE groupdef_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
#include      "scr17_c.inc"
#include      "tabsiz_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      TYPE (SUBSET_) , DIMENSION(NSUBS) :: SUBSET
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J,ID,ISU,NCHILD,NBPART,NBTPART,ICHILD,IPART,IPARTTOT,
     .        PARENT,LEVEL,TH_PART,L_SUB,ERR,ITITLE(LTITR),THIAD,IVARTH,
     .        LEN
      CHARACTER(LEN=LTITR) ::  TITR,NEW_TITLE(NSUBS)
      INTEGER, ALLOCATABLE, DIMENSION (:)  ::  ISUB
!=======================================================================
      DO isu=1,nsubs
        CALL read_i_c(ititle,ltitr)
        CALL fretitl2(titr,ititle,ltitr)
        len = 0
        DO j=1,ltitr
          IF (titr(j:j) /= ' ') len = j
        ENDDO
        new_title(isu) = titr(1:len)
      ENDDO ! DO ISU=1,NSUBS
!!------  Allocate and read subset data from restart
      ALLOCATE (isub(lensubs) ,stat=err)
      CALL read_i_c(isub, lensubs)
!=======================================================================
!
      l_sub = 0
!
      DO isu=1,nsubs
        id       = isub(l_sub+1)
          l_sub = l_sub+1
        nchild   = isub(l_sub+1)
          l_sub = l_sub+1
        nbpart   = isub(l_sub+1)
          l_sub = l_sub+1
        nbtpart  = isub(l_sub+1)
          l_sub = l_sub+1
        parent   = isub(l_sub+1)
          l_sub = l_sub+1
        level   = isub(l_sub+1)
          l_sub = l_sub+1
        thiad   = isub(l_sub+1)
          l_sub = l_sub+1
!!        TH_FLAG = ISUB(L_SUB+1)
!!          L_SUB = L_SUB+1
!!        ICHAR(TITR)    = ISUB(L_SUB+1)
!!          L_SUB = L_SUB+1
!---
       subset(isu)%ID     = id
       subset(isu)%NCHILD = nchild
       subset(isu)%NPART  = nbpart
       subset(isu)%NTPART = nbtpart
       subset(isu)%PARENT = parent
       subset(isu)%LEVEL  = level
       subset(isu)%THIAD  = thiad
!!       SUBSET(ISU)%TH_FLAG = TH_FLAG
       subset(isu)%TITLE = new_title(isu)
       ALLOCATE(subset(isu)%CHILD(nchild))
       ALLOCATE(subset(isu)%PART(nbpart))
       ALLOCATE(subset(isu)%TPART(nbtpart))
       ALLOCATE(subset(isu)%NVARTH(10))  !  /iTH (9) + /TH (1) -> SUBSET(I)%NVARTH(ITHFLAG)
!---
! NOMBRE DE SUBSETS FILS
        DO j=1,nchild
          ichild = isub(l_sub+1)
            l_sub = l_sub+1
          subset(isu)%CHILD(j) = ichild
        ENDDO
! NOMBRE DE PARTS FILLES
        DO j=1,nbpart
          ipart = isub(l_sub+1)
            l_sub = l_sub+1
          subset(isu)%PART(j) = ipart
        ENDDO
! NOMBRE TOTAL DE PARTS (DESCENDANCE)
        DO j=1,nbtpart
          iparttot = isub(l_sub+1)
            l_sub = l_sub+1
          subset(isu)%TPART(j) = iparttot
        ENDDO
!---
!   - TIME HISTORY - /iTH (9) + /TH (1)
!---
        DO j=1,10
          ivarth =  isub(l_sub+1)
            l_sub = l_sub+1
          subset(isu)%NVARTH(j) = ivarth
        ENDDO
      ENDDO ! DO ISU=1,NSUBS
!---------
      DEALLOCATE (isub)
!---------
      RETURN