wrtdes.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    wrtdes         ../engine/source/output/th/wrtdes.F
!||--- called by ------------------------------------------------------
!||    hist1          ../engine/source/output/th/hist1.F
!||    hist13         ../engine/source/output/th/hist13.F
!||    hist2          ../engine/source/output/th/hist2.F
!||    lecnoise       ../engine/source/general_controls/computation/lecnoise.F
!||    noise          ../engine/source/general_controls/computation/noise.f
!||    thcluster      ../engine/source/output/th/thcluster.F
!||    thkin          ../engine/source/output/th/thkin.F
!||    thmonv         ../engine/source/output/th/thmonv.F
!||    thsechecksum   ../engine/source/output/th/thchecksum.F90
!||    thsens         ../engine/source/output/th/thsens.F
!||    thsurf         ../engine/source/output/th/thsurf.F
!||    write_th       ../engine/source/output/th/write_th.F
!||    wrtdes0        ../engine/source/output/th/wrtdes0.F
!||--- calls      -----------------------------------------------------
!||    eor_c          ../common_source/tools/input_output/write_routtines.c
!||    write_i_c      ../common_source/tools/input_output/write_routtines.c
!||    write_r_c      ../common_source/tools/input_output/write_routtines.c
!||====================================================================
      SUBROUTINE wrtdes(A,IA,L,IFORM,IR)
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER L, IFORM, IR
      INTEGER IA(*)
      my_real a(*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "scr13_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, L1, L2, II
      REAL R4
      CHARACTER EOR*8, CH8(1000)*8
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
      CHARACTER STRR*8, STRI*8
C-----------------------------------------------
      EXTERNAL strr,stri
      DATA eor/'ZZZZZEOR'/
      IF(l==0)RETURN
 
      IF (iform==0)THEN
       IF (ir==1)THEN
        WRITE(iunit) (a(i),i=1,l)
       ELSE
        WRITE(iunit) (ia(i),i=1,l)
       ENDIF
       
      ELSEIF (iform==1)THEN
       l1=1
       l2=1000
       DO
         l2=min(l,l2)
         ii=0
         IF(ir==1)THEN
          DO i=l1,l2
           ii=ii+1
           ch8(ii)=strr(a(i))
          ENDDO
         ELSE
          DO i=l1,l2
            ii=ii+1
            ch8(ii)=stri(ia(i))
          ENDDO
         ENDIF
         IF(l2<l)THEN
          WRITE(iunit,'(10A8)') (ch8(i),i=1,ii)
          l1=l1+1000
          l2=l2+1000
          cycle
         ELSEIF(ii<1000) THEN
          WRITE(iunit,'(10A8)') (ch8(i),i=1,ii),eor
         ELSE
          WRITE(iunit,'(10A8)') (ch8(i),i=1,ii)
          WRITE(iunit,'(10A8)') eor
         ENDIF
       ENDDO
 
      ELSEIF (iform==2)THEN
       IF(ir==1)THEN
        WRITE(iunit,'(A,I5,A)')eor,l,'R'
        WRITE(iunit,'((5(1X,1PE15.8)))')(a(i),i=1,l)
       ELSE
        WRITE(iunit,'(A,I5,A)')eor,l,'I'
        WRITE(iunit,'(16I10)')          (ia(i),i=1,l)
       ENDIF
 
      ELSEIF (iform==3)THEN
       CALL eor_c(4*l)
       IF(ir==1)THEN
        DO i=1,l
          r4 = a(i)
          CALL write_r_c(r4,1)
        ENDDO
       ELSE
        CALL write_i_c(ia,l)
       ENDIF
       CALL eor_c(4*l)
       
      ENDIF
 
      RETURN
      END