alloc_surf_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_surf_str (igrsurf)

Function/Subroutine Documentation

alloc_surf_str()

subroutine alloc_surf_str

(

type(surf_), dimension(nsurf)

igrsurf

)

Definition at line 35 of file alloc_surf_str.F.

C-----------------------------------------------
C   M o d u l e s
C----------------------------------------------- 
      USE message_mod
      USE groupdef_mod
      USE names_and_titles_mod, ONLY: nchartitle
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(SURF_)   ,DIMENSION(NSURF) :: IGRSURF
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER ISU,I,J,K,L_SURF,ID,NSEG,TYPE,ID_MADYMO,IAD_BUFR,
     .        NB_MADYMO,TYPE_MADYMO,LEVEL,TH_SURF,ISH4N3N,NSEG_R2R_ALL,
     .        NSEG_R2R_SHARE,NODE,ELTYP,ELEM,ERR,ITITLE(LTITR),LEN
      CHARACTER(LEN=NCHARTITLE) :: TITR
      INTEGER, ALLOCATABLE, DIMENSION (:)  ::  ISURFI
!
!=======================================================================
      DO isu=1,nsurf
        titr(1:nchartitle) = ' '
        CALL read_i_c(ititle,ltitr)
        CALL fretitl2(titr,ititle,ltitr)
        len = 0
        DO j=1,ltitr
          IF (titr(j:j) /= ' ') len = j
        ENDDO
        igrsurf(isu)%TITLE = titr
      ENDDO ! DO ISU=1,NSURF
!!------  Allocate and read subset data from restart
!
      ALLOCATE (isurfi(lenisurf) ,stat=err)
      CALL read_i_c(isurfi, lenisurf)
!=======================================================================
      l_surf = 0
!
      DO isu=1,nsurf
        id             = isurfi(l_surf+1)
          l_surf = l_surf+1
        nseg           = isurfi(l_surf+1)
          l_surf = l_surf+1
        TYPE           = isurfi(l_surf+1)
          l_surf = l_surf+1
        id_madymo      = isurfi(l_surf+1)
          l_surf = l_surf+1
        iad_bufr       = isurfi(l_surf+1)
          l_surf = l_surf+1
        nb_madymo      = isurfi(l_surf+1)
          l_surf = l_surf+1
        type_madymo    = isurfi(l_surf+1)
          l_surf = l_surf+1
        level          = isurfi(l_surf+1)
          l_surf = l_surf+1
        th_surf        = isurfi(l_surf+1)
          l_surf = l_surf+1
        ish4n3n        = isurfi(l_surf+1)
          l_surf = l_surf+1
        nseg_r2r_all   = isurfi(l_surf+1)
          l_surf = l_surf+1
        nseg_r2r_share = isurfi(l_surf+1)
          l_surf = l_surf+1
!
!        ICHAR(TITR)    = ISURFI(L_SURF+1)
!          L_SURF = L_SURF+1
!---
       igrsurf(isu)%ID             = id
       igrsurf(isu)%NSEG           = nseg
       igrsurf(isu)%TYPE           = TYPE
       igrsurf(isu)%ID_MADYMO      = id_madymo
       igrsurf(isu)%IAD_BUFR       = iad_bufr
       igrsurf(isu)%NB_MADYMO      = nb_madymo
       igrsurf(isu)%TYPE_MADYMO    = type_madymo
       igrsurf(isu)%LEVEL          = level
       igrsurf(isu)%TH_SURF        = th_surf
       igrsurf(isu)%ISH4N3N        = ish4n3n
       igrsurf(isu)%NSEG_R2R_ALL   = nseg_r2r_all
       igrsurf(isu)%NSEG_R2R_SHARE = nseg_r2r_share
!
! SURF ENTITIES (NODES, ELTYP, ELEM)
!
       ALLOCATE(igrsurf(isu)%NODES(nseg,4))
       ALLOCATE(igrsurf(isu)%ELTYP(nseg))
       ALLOCATE(igrsurf(isu)%ELEM(nseg))
!
        DO j=1,nseg
          DO k=1,4
            node = isurfi(l_surf+1)
              l_surf = l_surf+1
            igrsurf(isu)%NODES(j,k) = node
          ENDDO
          eltyp = isurfi(l_surf+1)
            l_surf = l_surf+1
          igrsurf(isu)%ELTYP(j) = eltyp
          elem = isurfi(l_surf+1)
            l_surf = l_surf+1
          igrsurf(isu)%ELEM(j) = elem
        ENDDO
      ENDDO ! DO ISU=1,NSURF
!---------
      DEALLOCATE (isurfi)
!---------
      RETURN