MAKE_SDOF_OUTPUT_FILES.f90

Go to the documentation of this file.

!    Copyright (C) 2012 The SPEED FOUNDATION
!    Author: Ilario Mazzieri
!
!    This file is part of SPEED.
!
!    SPEED is free software; you can redistribute it and/or modify it
!    under the terms of the GNU Affero General Public License as
!    published by the Free Software Foundation, either version 3 of the
!    License, or (at your option) any later version.
!
!    SPEED 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
!    Affero General Public License for more details.
!
!    You should have received a copy of the GNU Affero General Public License
!    along with SPEED.  If not, see <http://www.gnu.org/licenses/>.
 
 
subroutine make_sdof_output_files
 
  use speed_sci
  use speed_timeloop
  use speed_par
 
  implicit none
  character*100 :: file
  integer*4 :: temp, sfs, ii, SDOFmon
 
  sdofmon = 701 + mpi_id
 
  if (n_bld.gt.0) then
 
    sfs = 0
    do ii = 1, n_bld
      temp = sys(ii)%SFS
      if(temp.eq.1) sfs = temp
    enddo
 
    if (sdofout(1).eq.1) then     ! displacement
 
      if(sfs.eq.0) then
 
        !!! structure relative displacement
 
        file="SDOF000000.DX"
 
        if(len_trim(monitor_file) .ne. 70) then
          sdofdisplx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          sdofdisplx = file
        endif
 
        open(sdofmon,file=sdofdisplx,status='replace')
        close(sdofmon)
 
        file="SDOF000000.DY"
 
        if(len_trim(monitor_file) .ne. 70) then
          sdofdisply = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          sdofdisply = file
        endif
 
        open(sdofmon,file=sdofdisply,status='replace')
        close(sdofmon)
 
        file="SDOF000000.DZ"
 
        if(len_trim(monitor_file) .ne. 70) then
          sdofdisplz = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          sdofdisplz = file
        endif
 
        open(sdofmon,file=sdofdisplz,status='replace')
        close(sdofmon)
 
        !!! soil displacement
 
        file="SDOF000000.GDX"
 
        if(len_trim(monitor_file) .ne. 70) then
          sdofgrdisplx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          sdofgrdisplx = file
        endif
 
        open(sdofmon,file=sdofgrdisplx,status='replace')
        close(sdofmon)
 
        file="SDOF000000.GDY"
 
        if(len_trim(monitor_file) .ne. 70) then
          sdofgrdisply = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          sdofgrdisply = file
        endif
 
        open(sdofmon,file=sdofgrdisply,status='replace')
        close(sdofmon)
 
        file="SDOF000000.GDZ"
 
        if(len_trim(monitor_file) .ne. 70) then
          sdofgrdisplz = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          sdofgrdisplz = file
        endif
 
        open(sdofmon,file=sdofgrdisplz,status='replace')
        close(sdofmon)
 
      elseif(sfs.eq.1) then
 
        !!! structure relative displacement
 
        file="STR000000.DX"
 
        if(len_trim(monitor_file) .ne. 70) then
          strdisplx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          strdisplx = file
        endif
 
        open(sdofmon,file=strdisplx,status='replace')
        close(sdofmon)
 
        file="STR000000.DY"
 
        if(len_trim(monitor_file) .ne. 70) then
          strdisply = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          strdisply = file
        endif
 
        open(sdofmon,file=strdisply,status='replace')
        close(sdofmon)
 
        !!! soil displacement
 
        file="GRD000000.DX"
 
        if(len_trim(monitor_file) .ne. 70) then
          grddisplx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          grddisplx = file
        endif
 
        open(sdofmon,file=grddisplx,status='replace')
        close(sdofmon)
 
        file="GRD000000.DY"
 
        if(len_trim(monitor_file) .ne. 70) then
          grddisply = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          grddisply = file
        endif
 
        open(sdofmon,file=grddisply,status='replace')
        close(sdofmon)
 
        file="GRD000000.DZ"
 
        if(len_trim(monitor_file) .ne. 70) then
          grddisplz = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          grddisplz = file
        endif
 
        open(sdofmon,file=grddisplz,status='replace')
        close(sdofmon)
 
        !!! foundation displacement
 
        file="FND000000.DX"
 
        if(len_trim(monitor_file) .ne. 70) then
          fnddisplx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          fnddisplx = file
        endif
 
        open(sdofmon,file=fnddisplx,status='replace')
        close(sdofmon)
 
        file="FND000000.DY"
 
        if(len_trim(monitor_file) .ne. 70) then
          fnddisply = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          fnddisply = file
        endif
 
        open(sdofmon,file=fnddisply,status='replace')
        close(sdofmon)
 
        file="FND000000.DRX"
 
        if(len_trim(monitor_file) .ne. 70) then
          fnddisplrx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          fnddisplrx = file
        endif
 
        open(sdofmon,file=fnddisplrx,status='replace')
        close(sdofmon)
 
        file="FND000000.DRY"
 
        if(len_trim(monitor_file) .ne. 70) then
          fnddisplry = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          fnddisplry = file
        endif
 
        open(sdofmon,file=fnddisplry,status='replace')
        close(sdofmon)
 
        file="FND000000.DZX"
 
        if(len_trim(monitor_file) .ne. 70) then
          fnddisplzx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          fnddisplzx = file
        endif
 
        open(sdofmon,file=fnddisplzx,status='replace')
        close(sdofmon)
 
        file="FND000000.DZY"
 
        if(len_trim(monitor_file) .ne. 70) then
          fnddisplzy = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          fnddisplzy = file
        endif
 
        open(sdofmon,file=fnddisplzy,status='replace')
        close(sdofmon)
      endif
    endif
 
    if (sdofout(2).eq.1) then       ! acceleration
 
      if(sfs.eq.0) then
 
        !!! structure total acceleration
 
        file="SDOF000000.AX"
 
        if(len_trim(monitor_file) .ne. 70) then
          sdofaccx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          sdofaccx = file
        endif
 
        open(sdofmon,file=sdofaccx,status='replace')
        close(sdofmon)
 
        file="SDOF000000.AY"
 
        if(len_trim(monitor_file) .ne. 70) then
          sdofaccy = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          sdofaccy = file
        endif
 
        open(sdofmon,file=sdofaccy,status='replace')
        close(sdofmon)
 
        file="SDOF000000.AZ"
 
        if(len_trim(monitor_file) .ne. 70) then
          sdofaccz = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          sdofaccz = file
        endif
 
        open(sdofmon,file=sdofaccz,status='replace')
        close(sdofmon)
 
        !!! soil acceleration
 
        file="SDOF000000.GAX"
 
        if(len_trim(monitor_file) .ne. 70) then
          sdofgraccx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          sdofgraccx = file
        endif
 
        open(sdofmon,file=sdofgraccx,status='replace')
        close(sdofmon)
 
        file="SDOF000000.GAY"
 
        if(len_trim(monitor_file) .ne. 70) then
          sdofgraccy = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          sdofgraccy = file
        endif
 
        open(sdofmon,file=sdofgraccy,status='replace')
        close(sdofmon)
 
        file="SDOF000000.GAZ"
 
        if(len_trim(monitor_file) .ne. 70) then
          sdofgraccz = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          sdofgraccz = file
        endif
 
        open(sdofmon,file=sdofgraccz,status='replace')
        close(sdofmon)
      elseif(sfs.eq.1) then
 
        !!! structure acceleration
 
        file="STR000000.AX"
 
        if(len_trim(monitor_file) .ne. 70) then
          straccx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          straccx = file
        endif
 
        open(sdofmon,file=straccx,status='replace')
        close(sdofmon)
 
        file="STR000000.AY"
 
        if(len_trim(monitor_file) .ne. 70) then
          straccy = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          straccy = file
        endif
 
        open(sdofmon,file=straccy,status='replace')
        close(sdofmon)
 
        !!! soil acceleration
 
        file="GRD000000.AX"
 
        if(len_trim(monitor_file) .ne. 70) then
          grdaccx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          grdaccx = file
        endif
 
        open(sdofmon,file=grdaccx,status='replace')
        close(sdofmon)
 
        file="GRD000000.AY"
 
        if(len_trim(monitor_file) .ne. 70) then
          grdaccy = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          grdaccy = file
        endif
 
        open(sdofmon,file=grdaccy,status='replace')
        close(sdofmon)
 
        file="GRD000000.AZ"
 
        if(len_trim(monitor_file) .ne. 70) then
          grdaccz = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          grdaccz = file
        endif
 
        open(sdofmon,file=grdaccz,status='replace')
        close(sdofmon)
 
        !!! foundation acceleration
 
        file="FND000000.AX"
 
        if(len_trim(monitor_file) .ne. 70) then
          fndaccx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          fndaccx = file
        endif
 
        open(sdofmon,file=fndaccx,status='replace')
        close(sdofmon)
 
        file="FND000000.AY"
 
        if(len_trim(monitor_file) .ne. 70) then
          fndaccy = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          fndaccy = file
        endif
 
        open(sdofmon,file=fndaccy,status='replace')
        close(sdofmon)
 
        file="FND000000.ARX"
 
        if(len_trim(monitor_file) .ne. 70) then
          fndaccrx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          fndaccrx = file
        endif
 
        open(sdofmon,file=fndaccrx,status='replace')
        close(sdofmon)
 
        file="FND000000.ARY"
 
        if(len_trim(monitor_file) .ne. 70) then
          fndaccry = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          fndaccry = file
        endif
 
        open(sdofmon,file=fndaccry,status='replace')
        close(sdofmon)
 
        file="FND000000.AZX"
 
        if(len_trim(monitor_file) .ne. 70) then
          fndacczx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          fndacczx = file
        endif
 
        open(sdofmon,file=fndacczx,status='replace')
        close(sdofmon)
 
        file="FND000000.AZY"
 
        if(len_trim(monitor_file) .ne. 70) then
          fndacczy = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          fndacczy = file
        endif
 
        open(sdofmon,file=fndacczy,status='replace')
        close(sdofmon)
      endif
    endif
 
    if (sdofout(3).eq.1) then       ! interaction force
 
      if (sfs.eq.0) then
 
        !!! structure shear force
 
        file="SDOF000000.FX"
 
        if(len_trim(monitor_file) .ne. 70) then
          sdoffx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          sdoffx = file
        endif
 
        open(sdofmon,file=sdoffx,status='replace')
        close(sdofmon)
 
        file="SDOF000000.FY"
 
        if(len_trim(monitor_file) .ne. 70) then
          sdoffy = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          sdoffy = file
        endif
 
        open(sdofmon,file=sdoffy,status='replace')
        close(sdofmon)
 
        file="SDOF000000.FZ"
 
        if(len_trim(monitor_file) .ne. 70) then
          sdoffz = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          sdoffz = file
        endif
 
        open(sdofmon,file=sdoffz,status='replace')
        close(sdofmon)
      elseif(sfs.eq.1) then
 
        !!! structure shear force
 
        file="STR000000.FX"
 
        if(len_trim(monitor_file) .ne. 70) then
          strfx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          strfx = file
        endif
 
        open(sdofmon,file=strfx,status='replace')
        close(sdofmon)
 
        file="STR000000.FY"
 
        if(len_trim(monitor_file) .ne. 70) then
          strfy = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          strfy = file
        endif
 
        open(sdofmon,file=strfy,status='replace')
        close(sdofmon)
 
        !!! foundation shear force
 
        file="FND000000.FX"
 
        if(len_trim(monitor_file) .ne. 70) then
          fndfx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          fndfx = file
        endif
 
        open(sdofmon,file=fndfx,status='replace')
        close(sdofmon)
 
        file="FND000000.FY"
 
        if(len_trim(monitor_file) .ne. 70) then
          fndfy = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          fndfy = file
        endif
 
        open(sdofmon,file=fndfy,status='replace')
        close(sdofmon)
 
        !!! interaction force
 
        file="INT000000.FX"
 
        if(len_trim(monitor_file) .ne. 70) then
          intfx = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          intfx = file
        endif
 
        open(sdofmon,file=intfx,status='replace')
        close(sdofmon)
 
        file="INT000000.FY"
 
        if(len_trim(monitor_file) .ne. 70) then
          intfy = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          intfy = file
        endif
 
        open(sdofmon,file=intfy,status='replace')
        close(sdofmon)
 
        file="INT000000.FZ"
 
        if(len_trim(monitor_file) .ne. 70) then
          intfz = monitor_file(1:len_trim(monitor_file)) // '/' // file
        else
          intfz = file
        endif
 
        open(sdofmon,file=intfz,status='replace')
        close(sdofmon)
      endif
    endif
  endif
 
  return
end subroutine make_sdof_output_files