Silverfrost Forums

Welcome to our forums

Save Attributes

6 Aug 2010 2:41 #6735

6.8.2010

Dear Sirs,

when I compile the program Axreib the options /List and /Dump (see compall.txt), I find in the list Axreib.lis 12 Save attributes for local variables.

Axreib uses he module Global_Data. I have used the latest FTN95 version 5.50 as you can see from the files *.lis.

Since there is nowhere any Save statement I wonder where the Save attributes come from.

Best regards,

Klaus Lassmann

Compall.bat:

del comp.lis
del *.obj

ftn95  Global_Data.f95  /Checkmate /Full_Debug  >> comp.lis
ftn95  Axreib.f95       /Checkmate /Full_Debug /list /dump >> comp.lis

Axreib.f95:

! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
! ++++++++++++++                                        ++++++++++++++++
! ++++++++++++++        Subroutine   A X R E I B        ++++++++++++++++
! ++++++++++++++        Subroutine   A X R E I B        ++++++++++++++++
! ++++++++++++++        Subroutine   A X R E I B        ++++++++++++++++
! ++++++++++++++                                        ++++++++++++++++
! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 
 
!     #################
      Subroutine Axreib 
!     #################
 
 
!     Subroutine AXREIB calculates the local friction force FJ
!     and the total friction force PJ
 
!     ------------------------------------------------------------------
 
      Use  Global_Data 

      Implicit None
 
      Double Precision , Dimension ( m31, m31) :: a

      Double Precision , Dimension ( m31     ) :: b, x

      Double Precision                         :: summal,  &
                                                  sueps0,  &
                                                  summaf,  &
                                                  al    ,  &
                                                  summfj

      Integer                                  :: ierr, j, k, l
 
!     ------------------------------------------------------------------
 
      k = 1
 
!     *************
      Do j = 2, m31
!     *************
 
!       =============================
        If ( irakon (j) .eq. 2 ) Then
!       =============================
 
!         --- irakon (j) = 2
 
          summal = 0.d+00
          sueps0 = 0.d+00
 
!         *************
          Do l = k, j-1
!         *************
 
            summal = summal + hhref (l) / aequer (l,1)   &
                            + hbref (l) / aequer (l,2)
 
            sueps0 = sueps0 + ( Axrelo (l) + eps0b (l) ) * hbref (l)  &
                                           - eps0h (l)   * hhref (l)
 
!         ******
          End Do
!         ******

          a (j,j:m31) = summal
 
          b (j)       = sueps0 + zbref (j) - zhref (j)   &
                               - zbref (k) + zhref (k)   &
                               - wrel  (j) + wrel  (k)
 
          summaf = 0.d+00
 
!         *************
          Do l = k, j-2
!         *************
 
            al     = + hhref (l) / aequer (l,1)          &
                     + hbref (l) / aequer (l,2)
 
            summfj = sum ( Fj (l+1:j-1) )
 
            summaf = summaf + al * summfj 
 
!         ******
          End Do
!         ******
 
          b (j) = b (j) - summaf
 
          k = j
 
!       ====
        Else
!       ====

!         ---  irakon (j) = 0, 1 or -1
 
!         +++++++++++++++++++++++++++++
          If ( irakon (j) .eq. 0 ) Then
!         +++++++++++++++++++++++++++++
 
            Fj (j) = 0.
            b  (j) = 0.d+00

!         ++++
          Else
!         ++++
 
            Fj (j) = Real ( irakon (j) ) * Reib (j)

            b  (j) = Fj (j)
 
!         ++++++
          End If
!         ++++++
 
          a (j,j+1:m31) = 0.d+00
 
          a (j,j      ) = 1.d+00
 
!       ======
        End If
!       ======
 
!     ******
      End Do
!     ******


!     ###########
      Call Fgausl 
!     ###########
 
 
!     --- Store x to FJ
 
      Fj (2:m31) = x (2:m31)
 
!     =======================
      If ( ierr .eq. 1 ) Then
!     =======================
 
        istop = 1
 
!     ======
      End If
!     ======


!     ########
      Contains
!     ########

 
!       #################
        Subroutine FGAUSL 
!       #################
 
 
!         ******************************************************************
!         *                                                                *
!         * Test programs see  Tusubsys\\Fgausl                             *
!         *                                                                *
!         ******************************************************************
 
!         Backsubstitution of upper triangular matrix
 
!         Friction forces are calculated for section 2, 3, ..., m3+1
 
!         ------------------------------------------------------------------

          Implicit None

          Double Precision             :: sum
          Double Precision , Parameter :: Small_Value = Tiny (sum)
          Integer                      :: i
 
!         ------------------------------------------------------------------
 
            ierr = 0
 
!         ================================================
          If ( Abs ( a (m31,m31) ) .lt. Small_Value ) Then
!         ================================================


            Call Fgausl_Error
!           #################
 
 
            ierr = 1

            Return
!           ++++++

!         ======
          End If
!         ======
 
          x (m31) = b (m31) / a (m31,m31)
 
!         ******************
          Do i  = m3, 2, -1
!         ******************
 
            sum = b (i) - dot_product ( a (i,i+1:m31), x (i+1:m31) )
 
!           ============================================
            If ( Abs ( a (i,i) ) .lt. Small_Value ) Then
!           ============================================
        
        
              Call Fgausl_Error
!             #################
 
 
              ierr = 1
          
              Return
!             ++++++
          
!           ======
            End If
!           ======

            x (i) = sum / a (i,i)
 
 
!         ******
          End Do
!         ******
   
 
!       #####################
        End Subroutine Fgausl 
!       #####################
 
 
!       #######################
        Subroutine Fgausl_Error
!       #######################

 
          Implicit None

          Integer :: i, j

!         --- System of equations is singular

          write   ( nwrite,  1000  )
          write   ( nwrite, '(//)' )

!         *************
          Do i = 2, m31
!         *************

            write ( nwrite, 2010   )    i, b (i)

!         ******
          End Do
!         ******

          write   ( nwrite, '(//)' )

!         *************
          Do i = 2, m31
!         *************

            write ( nwrite, 2020 )      i, i
            write ( nwrite, 2030 ) ( a (i, j), j = i,m31 )

            write ( nwrite, '(/)' )

!         ******
          End Do
!         ******
 
!         ------------------------------------------------------------------
 
   1000   Format (////1x,130('=') /                                         &
                      1x,130('=') /                                         &
                      1x,130('=') ///                                       &
                      1x,'Programme stop in subroutine    F G A U S L   ',  &
                      1x,'called by subroutine AXREIB because'/             &
                      1x,'the system of equations to determine the local',  &
                      1x,'friction forces FJ is singular'///                &
                      1x,130('=') /                                         &
                      1x,130('=') /                                         &
                      1x,130('=') / )
 
   2010   Format ( ' B (',i5,' ) = ',   d16.7 )
 
   2020   Format ( ' A (',i5,', j =',i5,', m31 ) : ' )
 
   2030   Format ( 1x, 10d16.7 )
   
 
!       ###########################
        End Subroutine Fgausl_Error
!       ###########################
 
  
!     #####################
      End Subroutine Axreib
!     #####################

Modul Global_Data.f95:

! ##################
  Module Global_Data
! ##################
 
! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
! ++++++++++++++++                                  ++++++++++++++++++++
! ++++++++++++++++     TRANSURANUS Global Data      ++++++++++++++++++++
! ++++++++++++++++     TRANSURANUS Global Data      ++++++++++++++++++++
! ++++++++++++++++     TRANSURANUS Global Data      ++++++++++++++++++++
! ++++++++++++++++                                  ++++++++++++++++++++
! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 
!     ==================================================================
!                      Definition of symbolic constants
 
!     --- Maximum number of radial coarse zones
 
      Integer , Parameter :: igrmax = 20
 
!     --- Maximum number of mesh points in a coarse zone
 
      Integer , Parameter :: ifemax = 20
 
!     --- Maximum number of axial slices
 
      Integer , Parameter :: m3max  = 40
      Integer , Parameter :: m31max = m3max+1
      Integer , Parameter :: m32max = m3max+2
 
!     --- Maximum number of equations for the radial thermal analysis
 
      Integer , Parameter :: ieqnta = igrmax*ifemax+2+10
 
!     --- Maximum number of equations for the pore migration
!         ( see puredi )
!         and densification ( see densi5 and densi9 )
 
!CC   Parameter ( ieqnmi = igrmax * (ifemax-1) + 1 )
 
!     --- Maximum number of equations for mechanics
 
      Integer , Parameter :: ieqnme = 2*igrmax+2
 
!     --- Maximum number of points in fuel rod  ( see temp, melt )
 
!CC   Parameter (ieqnfr = igrmax*ifemax)
!     ==================================================================
!                     Global Data B00
 
!     --- Input and output units
 
      Integer      :: UnitNr_restart, kwrite,                             &
                      UnitNr_pli, UnitNr_mic, UnitNr_mac, UnitNr_sta,     &
                      nread0, nread, nwrite, UnitNr_sta_scratch,          &
                      nerror, UnitNr_sta_log,                             &
                      UnitNr_itlog1, UnitNr_itlog2
 
 
!     ==================================================================
!                     Global Data B01

!     --- Control variables
 
      Integer      :: IALPHA,IANZA,IAXKON(M3MAX),IBLOC,                   &
                      iBroyden_s,                                         &
                      iBroyden_u,                                         &
                      IBMECH,IBZTEI(30),IBZZEI(30),IDEGAS(M31MAX,10),     &
                      IDENSI,IEND1,IFAALT(M31MAX),IFALLL(M31MAX),IFORM,   &
                      IFREI(IGRMAX,IFEMAX,3),IFRONT,                      &
                      IGASZA(M31MAX,3),IGASZI(M31MAX,3),IHGAP,IGRNSZ,     &
                      IKUEHL, ikueka, INSTA,INTRUP,INUMB,IPLOT,           &
                      IRAKON(M31MAX),IRAND,IRELOC,                        &
                      ISCHR,ISCHRL(M31MAX),ISTATI,ISTOP,                  &
                      ISTREC(IGRMAX,IFEMAX,4),ISTRUK,ITEALT,ITEMTE,       &
                      ITERAX,ITERAG,                                      &
                      ITERGE,ITERGS,ITERS(M31MAX),ITERSG,                 &
                      ITERT,ITERTG,ITHEOC,                                &
                      IVAR1,IXMODE,inumbm,RecNr_mic,RecNr_mac,IZENKA,     &
                      KFRONT,KPLOT,KRIMAX(IGRMAX,IFEMAX),                 &
                      Init_one (20), Init_fuel (20), Init_cool (20),      &
                      Init_axi (M31MAX,20),LSCHNI,MEND,M1,                &
                      M1B,M1H,M1STR,M2(IGRMAX),M3,M31,M32,NFRONT,NKOMM,   &
                      NRSUDT,NRSUIT, NTIMAN (40), NTIMSU (40),            &
                      Number_Nodes_fuel (m31max),                         &
                      Number_Nodes_clad (m31max)
 
 
!     ==================================================================
!                     Global Data B01a
 
      Integer      :: iclfail, ICORRO, icrkpi, ifba  , ihbs  , iHe     ,  &
                      ihydd  , iloca , IMPLIC, intaxl, ioxide, ioxire  ,  &
                      iphasezr,iplnum, IPLOT1,                            &
                      IPLOT2 , IPOINT, IPURE ,                            &
                      ISCHR1 , ISCHR2, ISLICE, ISTAT , ISTZNE, isurfb  ,  &
                      ITERM  , ITERMG, itlog ,                            &
                      ITRNPT (M31MAX),                                    &
                      kokoko,                                             &
                      MaxRecNr_mac, MaxRecNr_mic, MaxRecNr_sta,           &
                      MSTEP , NRSUIS, NUMBER(5), RecNr_sta
 
!     ==================================================================
!                     Global Data B02 
 
!     --- Logical variables
 
      Logical      :: AXKONV,                                             &
                      CLADRP(M31MAX), ISCHMB(M31MAX), ISCHMH(M31MAX),     &
                      LXKONV(M31MAX), MKONV,          SKONV (M31MAX),     &
                      THANLY(M31MAX), TKONV, TRANSA,                      &
                      MacroTime        ,                                  &   ! indicates 'macro time'
                      MacroTime_zero   ,                                  &
                      MacroTime_end    ,                                  &
                      MacroTime_restart,                                  &
                      Normal_end,                                         &
                      Cladrp_glob,                                        &
                      Change_of_Geometry
 
 
!     ==================================================================
!                     Global Data B03 
 
!     --- Geometrical quantities
 
      Real         :: DAL(M31MAX),HB(M3MAX),HBREF(M3MAX),HH(M31MAX),      &
                      HHREF(M31MAX),PiNumb,RBA(M31MAX),RBI(M31MAX),       &
                      REF(IGRMAX,IFEMAX,2),                               &
                      REFBA(M31MAX),REFBA0(M31MAX),                       &
                      REFBI(M31MAX),REFBI0(M31MAX),                       &
                      REFHA(M31MAX),REFHA0(M31MAX),                       &
                      REFHI(M31MAX),REFHI0(M31MAX),                       &
                      RHA(M31MAX),RHI(M31MAX),                            &
                      HBSwidth (m31max,2),                                &
                      RISSL(M31MAX),                                      &
                      RMELT(M31MAX,2,2,2),RRVEK(IGRMAX,IFEMAX,2),         &
                      RSTRU(10,2),SDLTCL(M31MAX),                         &
                      SOXIDE(M31MAX,2),UDOTFI(M31MAX),                    &
                      UDOTFO(M31MAX),WB(M31MAX),WH(M32MAX),               &
                      WREL(M31MAX),ZB(M31MAX),ZBREF(M31MAX),              &
                      ZH(M32MAX),ZHREF(M32MAX)
 
 
!     ==================================================================
!                     Global Data B04 
 
!     --- Mechanics
 
      Real         :: AEQUER(M3MAX,2),AREACL(M31MAX),AREAFU(M31MAX),      &
                      AXRELO(M3MAX),C1(IGRMAX),C2(IGRMAX),C3(IGRMAX),     &
                      C3JB(M31MAX),C3JH(M31MAX),C3REFB(M31MAX),           &
                      C3REFH(M31MAX), cvdivt (m31max),                    &
                      EMOD  (IGRMAX), EMODCL (M31MAX),                    &
                      EMODFU(M31MAX),EPS0B(M3MAX),EPS0H(M3MAX),ERESGR,    &
                      ETA(IGRMAX,IFEMAX,3,12),ETAALT(M31MAX),ETACRP,      &
                      ETA10(IGRMAX,IFEMAX,3),                             &
                      ETA11(IGRMAX,IFEMAX,3),ETA12(IGRMAX,IFEMAX,3),      &
                      ETA3(IGRMAX,IFEMAX,3),ETA5(IGRMAX,IFEMAX,3),        &
                      ETA6(IGRMAX,IFEMAX,3),ETA7(IGRMAX,IFEMAX,3),        &
                      ETA8(IGRMAX,IFEMAX,3),ETA9(IGRMAX,IFEMAX,3),        &
                      FMUEF,FMUEH,FNUE(IGRMAX),FNUECL(M31MAX),            &
                      FNUEFU(M31MAX),G(18,IGRMAX),INT1(IGRMAX,IFEMAX),    &
                      INT2(IGRMAX,IFEMAX),INT3(IGRMAX,IFEMAX),            &
                      INT4(IGRMAX),PLAST1(IGRMAX,IFEMAX),RELO(M31MAX),    &
                      SIGAVE(2,4,M31MAX),SIGMA(IGRMAX,IFEMAX,4,2),        &
                      SIGS(IGRMAX,IFEMAX), RadDef_u (igrmax,ifemax),      &
                      VISCOK(IGRMAX,M31MAX),ZKMIN
 
 
!     ==================================================================
!                     Global Data B04a
 
!     --- Mechanics: average strains in cladding
 
!                     EpsEffCl_crp = average effective creep strain
!                     S_depsCl_crp = sum of the effective
!                                    creep strain increments
!                     EpsEffCl_pla = average effective plastic strain
!                     S_depsCl_pla = sum of the effective
!                                    plastic strain increments
!                     EpsTanCl_tot = average permanent tangential strain
 
      Real         :: EpsEffCl_crp (m31max)  ,                            &
                      S_depsCl_crp (m31max,2),                            &
                      EpsEffCl_pla (m31max)  ,                            &
                      S_depsCl_pla (m31max,2),                            &
                      EpsTanCl_tot (m31max)
 
 
!     ==================================================================
!                     Global Data B04b
 
!     --- Mechanics: time integral of effective stress
!                    applied for the primary creep of cladding
 
      Real         :: sigmaefftime (IGRMAX, IFEMAX, 2),                   &
                      ETAstat      (IGRMAX, IFEMAX, 2)
 
 
!     ==================================================================
!                     Global Data B06
 
!     --- Forces and weights
 
      Real         :: FJ(M31MAX),FSPRNG,GEWB(M31MAX),GEWH(M31MAX),        &
                      HAFT   (M31MAX),                                    &
                      PJ     (M31MAX),                                    &
                      PAXIAL (M31MAX),                                    &
                      REIB   (M31MAX)
 
 
! ######################
  End Module Global_Data
! ######################

Axreib.lis:

Silverfrost FTN95/WIN32 Ver 5.50.0  Axreib.F95  Fri Aug  6 16:16:29 2010

Compiler Options in Effect:
    CHECKMATE COLOUR DELETE_OBJ_ON_ERROR DUMP FULL_DEBUG LIST MINIMISE_REBUILD  NO_QUIT 
    NON_STANDARD SINGLE_THREADED 

   0001                                                                                          AT 0
   0002   ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++               AT 0
   0003   ! ++++++++++++++                                        ++++++++++++++++               AT 0
   0004   ! ++++++++++++++        Subroutine   A X R E I B        ++++++++++++++++               AT 0
   0005   ! ++++++++++++++        Subroutine   A X R E I B        ++++++++++++++++               AT 0
   0006   ! ++++++++++++++        Subroutine   A X R E I B        ++++++++++++++++               AT 0
   0007   ! ++++++++++++++                                        ++++++++++++++++               AT 0
   0008   ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++               AT 0
   0009                                                                                          AT 0
   0010                                                                                          AT 0
   0011   !     #################                                                                AT 0
   0012         Subroutine Axreib                                                                AT 0
   0013   !     #################                                                                AT 1d
   0014                                                                                          AT 1d
   0015                                                                                          AT 1d
   0016   !     Subroutine AXREIB calculates the local friction force FJ                         AT 1d
   0017   !     and the total friction force PJ                                                  AT 1d
   0018                                                                                          AT 1d
   0019   !     ------------------------------------------------------------------               AT 1d
   0020                                                                                          AT 1d
   0021         Use  Global_Data                                                                 AT 1d
   0022                                                                                          AT 1d
   0023         Implicit None                                                                    AT 1d
   0024                                                                                          AT 1d
   0025         Double Precision , Dimension ( m31, m31) :: a                                    AT 1d
   0026                                                                                          AT 6d
   0027         Double Precision , Dimension ( m31     ) :: b, x                                 AT 6d
   0028                                                                                          AT 176
   0029         Double Precision                         :: summal,  &                           AT 176
   0030                                                     sueps0,  &                           AT 176
   0031                                                     summaf,  &                           AT 176
   0032                                                     al    ,  &                           AT 176
   0033                                                     summfj                               AT 176
   0034                                                                                          AT 176
   0035         Integer                                  :: ierr, j, k, l                        AT 176
   0036                                                                                          AT 176
   0037   !     ------------------------------------------------------------------               AT 176
   0038                                                                                          AT 176
   0039         k = 1                                                                            AT 176
   0040                                                                                          AT 181
   0041   !     *************                                                                    AT 181
   0042         Do j = 2, m31                                                                    AT 181
   0043   !     *************                                                                    AT 1ca
   0044                                                                                          AT 1ca
   0045   !       =============================                                                  AT 1ca
   0046           If ( irakon (j) .eq. 2 ) Then                                                  AT 1ca
   0047   !       =============================                                                  AT 203
   0048                                                                                          AT 203
   0049   !         --- irakon (j) = 2                                                           AT 203
   0050                                                                                          AT 203
   0051             summal = 0.d+00                                                              AT 203
   0052             sueps0 = 0.d+00                                                              AT 218
   0053                                                                                          AT 22d
   0054   !         *************                                                                AT 22d
   0055             Do l = k, j-1                                                                AT 22d
   0056   !         *************                                                                AT 27a
   0057                                                                                          AT 27a
   0058               summal = summal + hhref (l) / aequer (l,1)   &                             AT 27a
   0059                               + hbref (l) / aequer (l,2)                                 AT 27a
   0060                                                                                          AT 342
   0061               sueps0 = sueps0 + ( Axrelo (l) + eps0b (l) ) * hbref (l)  &                AT 342
   0062                                              - eps0h (l)   * hhref (l)                   AT 342
   0063                                                                                          AT 434
   0064   !         ******                                                                       AT 434
   0065             End Do                                                                       AT 434
   0066   !         ******                                                                       AT 454
   0067                                                                                          AT 454
   0068             a (j,j:m31) = summal                                                         AT 454
   0069                                                                                          AT 4e1
   0070             b (j)       = sueps0 + zbref (j) - zhref (j)   &                             AT 4e1
   0071                                  - zbref (k) + zhref (k)   &                             AT 4e1
   0072                                  - wrel  (j) + wrel  (k)                                 AT 4e1
   0073                                                                                          AT 64c
   0074             summaf = 0.d+00                                                              AT 64c
   0075                                                                                          AT 661
   0076   !         *************                                                                AT 661
   0077             Do l = k, j-2                                                                AT 661
   0078   !         *************                                                                AT 6b0
   0079                                                                                          AT 6b0
   0080               al     = + hhref (l) / aequer (l,1)          &                             AT 6b0
   0081                        + hbref (l) / aequer (l,2)                                        AT 6b0
   0082                                                                                          AT 76a
   0083               summfj = sum ( Fj (l+1:j-1) )                                              AT 76a
   0084                                                                                          AT 86d
   0085               summaf = summaf + al * summfj                                              AT 86d
   0086                                                                                          AT 89d
   0087   !         ******                                                                       AT 89d
   0088             End Do                                                                       AT 89d
   0089   !         ******                                                                       AT 8bd
   0090                                                                                          AT 8bd
   0091             b (j) = b (j) - summaf                                                       AT 8bd
   0092                                                                                          AT 923
   0093             k = j                                                                        AT 923
   0094                                                                                          AT 935
   0095   !       ====                                                                           AT 935
   0096           Else                                                                           AT 935
   0097   !       ====                                                                           AT 935
   0098                                                                                          AT 935
   0099   !         ---  irakon (j) = 0, 1 or -1                                                 AT 935
   0100                                                                                          AT 935
   0101   !         +++++++++++++++++++++++++++++                                                AT 935
   0102             If ( irakon (j) .eq. 0 ) Then                                                AT 935
   0103   !         +++++++++++++++++++++++++++++                                                AT 96e
   0104                                                                                          AT 96e
   0105               Fj (j) = 0.                                                                AT 96e
   0106               b  (j) = 0.d+00                                                            AT 995
   0107                                                                                          AT 9cf
   0108   !         ++++                                                                         AT 9cf
   0109             Else                                                                         AT 9cf
   0110   !         ++++                                                                         AT 9cf
   0111                                                                                          AT 9cf
   0112               Fj (j) = Real ( irakon (j) ) * Reib (j)                                    AT 9cf
   0113                                                                                          AT a52
   0114               b  (j) = Fj (j)                                                            AT a52
   0115                                                                                          AT aa4
   0116   !         ++++++                                                                       AT aa4
   0117             End If                                                                       AT aa4
   0118   !         ++++++                                                                       AT aa4
   0119                                                                                          AT aa4
   0120             a (j,j+1:m31) = 0.d+00                                                       AT aa4
   0121                                                                                          AT b30
   0122             a (j,j      ) = 1.d+00                                                       AT b30
   0123                                                                                          AT b92
   0124   !       ======                                                                         AT b92
   0125           End If                                                                         AT b92
   0126   !       ======                                                                         AT b92
   0127                                                                                          AT b92
   0128   !     ******                                                                           AT b92
   0129         End Do                                                                           AT b92
   0130   !     ******                                                                           AT bb2
   0131                                                                                          AT bb2
   0132                                                                                          AT bb2
   0133   !     ###########                                                                      AT bb2
   0134         Call Fgausl                                                                      AT bb2
   0135   !     ###########                                                                      AT bc0
   0136                                                                                          AT bc0
   0137                                                                                          AT bc0
   0138   !     --- Store x to FJ                                                                AT bc0
   0139                                                                                          AT bc0
   0140         Fj (2:m31) = x (2:m31)                                                           AT bc0
   0141                                                                                          AT c1d
   0142   !     =======================                                                          AT c1d
   0143         If ( ierr .eq. 1 ) Then                                                          AT c1d
   0144   !     =======================                                                          AT c3a
   0145                                                                                          AT c3a
   0146           istop = 1                                                                      AT c3a


  Dump of all variables in AXREIB

  Local Variables
SAVE            INTEGER(KIND=3)      L                   @ 68  (at line 38)
SAVE            INTEGER(KIND=3)      K                   @ 72  (at line 38)
SAVE            INTEGER(KIND=3)      J                   @ 76  (at line 38)
SAVE            INTEGER(KIND=3)      IERR                @ 80  (at line 38)
SAVE            REAL(KIND=2)         SUMMFJ              @ 0   (at line 34)
SAVE            REAL(KIND=2)         AL                  @ 8   (at line 34)
SAVE            REAL(KIND=2)         SUMMAF              @ 16  (at line 34)
SAVE            REAL(KIND=2)         SUEPS0              @ 24  (at line 34)
SAVE            REAL(KIND=2)         SUMMAL              @ 32  (at line 34)
SAVE            REAL(KIND=2)         X(:)                @ 88  (at line 28)
SAVE            REAL(KIND=2)         B(:)                @ 96  (at line 28)
SAVE            REAL(KIND=2)         A(:, :)             @ 108 (at line 26)
                INTEGER(KIND=3)      UNITNR_RESTART      @ 15860 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      KWRITE              @ 15856 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      UNITNR_PLI          @ 15852 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      UNITNR_MIC          @ 15848 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      UNITNR_MAC          @ 15844 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      UNITNR_STA          @ 15840 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      NREAD0              @ 15836 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      NREAD               @ 15832 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      NWRITE              @ 15828 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      UNITNR_STA_SCRATCH  @ 15824 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      NERROR              @ 15820 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      UNITNR_STA_LOG      @ 15816 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      UNITNR_ITLOG1       @ 15812 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      UNITNR_ITLOG2       @ 15808 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IALPHA              @ 15804 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IANZA               @ 15800 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IAXKON(40)          @ 14080 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IBLOC               @ 15796 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IBROYDEN_S          @ 15792 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IBROYDEN_U          @ 15788 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IBMECH              @ 15784 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IBZTEI(30)          @ 15664 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IBZZEI(30)          @ 15544 in MODULE GLOBAL_DATA
/@@GLOBAL_DATA!IDEGAS/ INTEGER(KIND=3)      IDEGAS(41, 10)      @ 0   (at line 22)
                INTEGER(KIND=3)      IDENSI              @ 15540 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IEND1               @ 15536 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IFAALT(41)          @ 13916 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IFALLL(41)          @ 13752 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IFORM               @ 15532 in MODULE GLOBAL_DATA
/@@GLOBAL_DATA!IFREI/ INTEGER(KIND=3)      IFREI(20, 20, 3)    @ 0   (at line 22)
                INTEGER(KIND=3)      IFRONT              @ 15528 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IGASZA(41, 3)       @ 13260 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IGASZI(41, 3)       @ 12768 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IHGAP               @ 15524 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IGRNSZ              @ 15520 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IKUEHL              @ 15516 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IKUEKA              @ 15512 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      INSTA               @ 15508 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      INTRUP              @ 15504 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      INUMB               @ 15500 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IPLOT               @ 15496 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IRAKON(41)          @ 12604 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IRAND               @ 15492 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IRELOC              @ 15488 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ISCHR               @ 15484 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ISCHRL(41)          @ 12440 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ISTATI              @ 15480 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ISTOP               @ 15476 in MODULE GLOBAL_DATA
/@@GLOBAL_DATA!ISTREC/ INTEGER(KIND=3)      ISTREC(20, 20, 4)   @ 0   (at line 22)
                INTEGER(KIND=3)      ISTRUK              @ 15472 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ITEALT              @ 15468 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ITEMTE              @ 15464 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ITERAX              @ 15460 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ITERAG              @ 15456 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ITERGE              @ 15452 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ITERGS              @ 15448 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ITERS(41)           @ 12276 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ITERSG              @ 15444 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ITERT               @ 15440 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ITERTG              @ 15436 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ITHEOC              @ 15432 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IVAR1               @ 15428 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IXMODE              @ 15424 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      INUMBM              @ 15420 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      RECNR_MIC           @ 15416 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      RECNR_MAC           @ 15412 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IZENKA              @ 15408 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      KFRONT              @ 15404 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      KPLOT               @ 15400 in MODULE GLOBAL_DATA
/@@GLOBAL_DATA!KRIMAX/ INTEGER(KIND=3)      KRIMAX(20, 20)      @ 0   (at line 22)
                INTEGER(KIND=3)      INIT_ONE(20)        @ 15320 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      INIT_FUEL(20)       @ 15240 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      INIT_COOL(20)       @ 15160 in MODULE GLOBAL_DATA
/@@GLOBAL_DATA!INIT_AXI/ INTEGER(KIND=3)      INIT_AXI(41, 20)    @ 0   (at line 22)
                INTEGER(KIND=3)      LSCHNI              @ 15156 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      MEND                @ 15152 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      M1                  @ 15148 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      M1B                 @ 15144 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      M1H                 @ 15140 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      M1STR               @ 15136 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      M2(20)              @ 15056 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      M3                  @ 15052 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      M31                 @ 15048 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      M32                 @ 15044 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      NFRONT              @ 15040 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      NKOMM               @ 15036 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      NRSUDT              @ 15032 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      NRSUIT              @ 15028 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      NTIMAN(40)          @ 12116 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      NTIMSU(40)          @ 11956 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      NUMBER_NODES_FUEL(41) @ 11792 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      NUMBER_NODES_CLAD(41) @ 11628 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ICLFAIL             @ 15024 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ICORRO              @ 15020 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ICRKPI              @ 15016 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IFBA                @ 15012 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IHBS                @ 15008 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IHE                 @ 15004 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IHYDD               @ 15000 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ILOCA               @ 14996 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IMPLIC              @ 14992 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      INTAXL              @ 14988 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IOXIDE              @ 14984 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IOXIRE              @ 14980 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IPHASEZR            @ 14976 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IPLNUM              @ 14972 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IPLOT1              @ 14968 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IPLOT2              @ 14964 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IPOINT              @ 14960 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      IPURE               @ 14956 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ISCHR1              @ 14952 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ISCHR2              @ 14948 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ISLICE              @ 14944 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ISTAT               @ 14940 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ISTZNE              @ 14936 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ISURFB              @ 14932 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ITERM               @ 14928 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ITERMG              @ 14924 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ITLOG               @ 14920 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      ITRNPT(41)          @ 11464 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      KOKOKO              @ 14916 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      MAXRECNR_MAC        @ 14912 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      MAXRECNR_MIC        @ 14908 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      MAXRECNR_STA        @ 14904 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      MSTEP               @ 14900 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      NRSUIS              @ 14896 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      NUMBER(5)           @ 14876 in MODULE GLOBAL_DATA
                INTEGER(KIND=3)      RECNR_STA           @ 14872 in MODULE GLOBAL_DATA
                LOGICAL(KIND=3)      AXKONV              @ 14868 in MODULE GLOBAL_DATA
                LOGICAL(KIND=3)      CLADRP(41)          @ 11300 in MODULE GLOBAL_DATA
                LOGICAL(KIND=3)      ISCHMB(41)          @ 11136 in MODULE GLOBAL_DATA
                LOGICAL(KIND=3)      ISCHMH(41)          @ 10972 in MODULE GLOBAL_DATA
                LOGICAL(KIND=3)      LXKONV(41)          @ 10808 in MODULE GLOBAL_DATA
                LOGICAL(KIND=3)      MKONV               @ 14864 in MODULE GLOBAL_DATA
                LOGICAL(KIND=3)      SKONV(41)           @ 10644 in MODULE GLOBAL_DATA
                LOGICAL(KIND=3)      THANLY(41)          @ 10480 in MODULE GLOBAL_DATA
                LOGICAL(KIND=3)      TKONV               @ 14860 in MODULE GLOBAL_DATA
                LOGICAL(KIND=3)      TRANSA              @ 14856 in MODULE GLOBAL_DATA
                LOGICAL(KIND=3)      MACROTIME           @ 14852 in MODULE GLOBAL_DATA
                LOGICAL(KIND=3)      MACROTIME_ZERO      @ 14848 in MODULE GLOBAL_DATA
                LOGICAL(KIND=3)      MACROTIME_END       @ 14844 in MODULE GLOBAL_DATA
                LOGICAL(KIND=3)      MACROTIME_RESTART   @ 14840 in MODULE GLOBAL_DATA
                LOGICAL(KIND=3)      NORMAL_END          @ 14836 in MODULE GLOBAL_DATA
                LOGICAL(KIND=3)      CLADRP_GLOB         @ 14832 in MODULE GLOBAL_DATA
                LOGICAL(KIND=3)      CHANGE_OF_GEOMETRY  @ 14828 in MODULE GLOBAL_DATA
                REAL(KIND=1)         DAL(41)             @ 10316 in MODULE GLOBAL_DATA
                REAL(KIND=1)         HB(40)              @ 10156 in MODULE GLOBAL_DATA
                REAL(KIND=1)         HBREF(40)           @ 9996 in MODULE GLOBAL_DATA
                REAL(KIND=1)         HH(41)              @ 9832 in MODULE GLOBAL_DATA
                REAL(KIND=1)         HHREF(41)           @ 9668 in MODULE GLOBAL_DATA
                REAL(KIND=1)         PINUMB              @ 14824 in MODULE GLOBAL_DATA
                REAL(KIND=1)         RBA(41)             @ 9504 in MODULE GLOBAL_DATA
                REAL(KIND=1)         RBI(41)             @ 9340 in MODULE GLOBAL_DATA
/@@GLOBAL_DATA!REF/ REAL(KIND=1)         REF(20, 20, 2)      @ 0   (at line 22)
                REAL(KIND=1)         REFBA(41)           @ 9176 in MODULE GLOBAL_DATA
                REAL(KIND=1)         REFBA0(41)          @ 9012 in MODULE GLOBAL_DATA
                REAL(KIND=1)         REFBI(41)           @ 8848 in MODULE GLOBAL_DATA
                REAL(KIND=1)         REFBI0(41)          @ 8684 in MODULE GLOBAL_DATA
                REAL(KIND=1)         REFHA(41)           @ 8520 in MODULE GLOBAL_DATA
                REAL(KIND=1)         REFHA0(41)          @ 8356 in MODULE GLOBAL_DATA
                REAL(KIND=1)         REFHI(41)           @ 8192 in MODULE GLOBAL_DATA
                REAL(KIND=1)         REFHI0(41)          @ 8028 in MODULE GLOBAL_DATA
                REAL(KIND=1)         RHA(41)             @ 7864 in MODULE GLOBAL_DATA
                REAL(KIND=1)         RHI(41)             @ 7700 in MODULE GLOBAL_DATA
                REAL(KIND=1)         HBSWIDTH(41, 2)     @ 7372 in MODULE GLOBAL_DATA
                REAL(KIND=1)         RISSL(41)           @ 7208 in MODULE GLOBAL_DATA
/@@GLOBAL_DATA!RMELT/ REAL(KIND=1)         RMELT(41, 2, 2, 2)  @ 0   (at line 22)
/@@GLOBAL_DATA!RRVEK/ REAL(KIND=1)         RRVEK(20, 20, 2)    @ 0   (at line 22)
                REAL(KIND=1)         RSTRU(10, 2)        @ 14744 in MODULE GLOBAL_DATA
                REAL(KIND=1)         SDLTCL(41)          @ 7044 in MODULE GLOBAL_DATA
                REAL(KIND=1)         SOXIDE(41, 2)       @ 6716 in MODULE GLOBAL_DATA
                REAL(KIND=1)         UDOTFI(41)          @ 6552 in MODULE GLOBAL_DATA
                REAL(KIND=1)         UDOTFO(41)          @ 6388 in MODULE GLOBAL_DATA
                REAL(KIND=1)         WB(41)              @ 6224 in MODULE GLOBAL_DATA
                REAL(KIND=1)         WH(42)              @ 6056 in MODULE GLOBAL_DATA
                REAL(KIND=1)         WREL(41)            @ 5892 in MODULE GLOBAL_DATA
                REAL(KIND=1)         ZB(41)              @ 5728 in MODULE GLOBAL_DATA
                REAL(KIND=1)         ZBREF(41)           @ 5564 in MODULE GLOBAL_DATA
                REAL(KIND=1)         ZH(42)              @ 5396 in MODULE GLOBAL_DATA
                REAL(KIND=1)         ZHREF(42)           @ 5228 in MODULE GLOBAL_DATA
                REAL(KIND=1)         AEQUER(40, 2)       @ 4908 in MODULE GLOBAL_DATA
                REAL(KIND=1)         AREACL(41)          @ 4744 in MODULE GLOBAL_DATA
                REAL(KIND=1)         AREAFU(41)          @ 4580 in MODULE GLOBAL_DATA
                REAL(KIND=1)         AXRELO(40)          @ 4420 in MODULE GLOBAL_DATA
                REAL(KIND=1)         C1(20)              @ 14664 in MODULE GLOBAL_DATA
                REAL(KIND=1)         C2(20)              @ 14584 in MODULE GLOBAL_DATA
                REAL(KIND=1)         C3(20)              @ 14504 in MODULE GLOBAL_DATA
                REAL(KIND=1)         C3JB(41)            @ 4256 in MODULE GLOBAL_DATA
                REAL(KIND=1)         C3JH(41)            @ 4092 in MODULE GLOBAL_DATA
                REAL(KIND=1)         C3REFB(41)          @ 3928 in MODULE GLOBAL_DATA
                REAL(KIND=1)         C3REFH(41)          @ 3764 in MODULE GLOBAL_DATA
                REAL(KIND=1)         CVDIVT(41)          @ 3600 in MODULE GLOBAL_DATA
                REAL(KIND=1)         EMOD(20)            @ 14424 in MODULE GLOBAL_DATA
                REAL(KIND=1)         EMODCL(41)          @ 3436 in MODULE GLOBAL_DATA
                REAL(KIND=1)         EMODFU(41)          @ 3272 in MODULE GLOBAL_DATA
                REAL(KIND=1)         EPS0B(40)           @ 3112 in MODULE GLOBAL_DATA
                REAL(KIND=1)         EPS0H(40)           @ 2952 in MODULE GLOBAL_DATA
                REAL(KIND=1)         ERESGR              @ 14420 in MODULE GLOBAL_DATA
/@@GLOBAL_DATA!ETA/ REAL(KIND=1)         ETA(20, 20, 3, 12)  @ 0   (at line 22)
                REAL(KIND=1)         ETAALT(41)          @ 2788 in MODULE GLOBAL_DATA
                REAL(KIND=1)         ETACRP              @ 14416 in MODULE GLOBAL_DATA
/@@GLOBAL_DATA!ETA10/ REAL(KIND=1)         ETA10(20, 20, 3)    @ 0   (at line 22)
/@@GLOBAL_DATA!ETA11/ REAL(KIND=1)         ETA11(20, 20, 3)    @ 0   (at line 22)
/@@GLOBAL_DATA!ETA12/ REAL(KIND=1)         ETA12(20, 20, 3)    @ 0   (at line 22)
/@@GLOBAL_DATA!ETA3/ REAL(KIND=1)         ETA3(20, 20, 3)     @ 0   (at line 22)
/@@GLOBAL_DATA!ETA5/ REAL(KIND=1)         ETA5(20, 20, 3)     @ 0   (at line 22)
/@@GLOBAL_DATA!ETA6/ REAL(KIND=1)         ETA6(20, 20, 3)     @ 0   (at line 22)
/@@GLOBAL_DATA!ETA7/ REAL(KIND=1)         ETA7(20, 20, 3)     @ 0   (at line 22)
/@@GLOBAL_DATA!ETA8/ REAL(KIND=1)         ETA8(20, 20, 3)     @ 0   (at line 22)
/@@GLOBAL_DATA!ETA9/ REAL(KIND=1)         ETA9(20, 20, 3)     @ 0   (at line 22)
                REAL(KIND=1)         FMUEF               @ 14412 in MODULE GLOBAL_DATA
                REAL(KIND=1)         FMUEH               @ 14408 in MODULE GLOBAL_DATA
                REAL(KIND=1)         FNUE(20)            @ 14328 in MODULE GLOBAL_DATA
                REAL(KIND=1)         FNUECL(41)          @ 2624 in MODULE GLOBAL_DATA
                REAL(KIND=1)         FNUEFU(41)          @ 2460 in MODULE GLOBAL_DATA
/@@GLOBAL_DATA!G/ REAL(KIND=1)         G(18, 20)           @ 0   (at line 22)
/@@GLOBAL_DATA!INT1/ REAL(KIND=1)         INT1(20, 20)        @ 0   (at line 22)
/@@GLOBAL_DATA!INT2/ REAL(KIND=1)         INT2(20, 20)        @ 0   (at line 22)
/@@GLOBAL_DATA!INT3/ REAL(KIND=1)         INT3(20, 20)        @ 0   (at line 22)
                REAL(KIND=1)         INT4(20)            @ 14248 in MODULE GLOBAL_DATA
/@@GLOBAL_DATA!PLAST1/ REAL(KIND=1)         PLAST1(20, 20)      @ 0   (at line 22)
                REAL(KIND=1)         RELO(41)            @ 2296 in MODULE GLOBAL_DATA
/@@GLOBAL_DATA!SIGAVE/ REAL(KIND=1)         SIGAVE(2, 4, 41)    @ 0   (at line 22)
/@@GLOBAL_DATA!SIGMA/ REAL(KIND=1)         SIGMA(20, 20, 4, 2) @ 0   (at line 22)
/@@GLOBAL_DATA!SIGS/ REAL(KIND=1)         SIGS(20, 20)        @ 0   (at line 22)
/@@GLOBAL_DATA!RADDEF_U/ REAL(KIND=1)         RADDEF_U(20, 20)    @ 0   (at line 22)
/@@GLOBAL_DATA!VISCOK/ REAL(KIND=1)         VISCOK(20, 41)      @ 0   (at line 22)
                REAL(KIND=1)         ZKMIN               @ 14244 in MODULE GLOBAL_DATA
                REAL(KIND=1)         EPSEFFCL_CRP(41)    @ 2132 in MODULE GLOBAL_DATA
                REAL(KIND=1)         S_DEPSCL_CRP(41, 2) @ 1804 in MODULE GLOBAL_DATA
                REAL(KIND=1)         EPSEFFCL_PLA(41)    @ 1640 in MODULE GLOBAL_DATA
                REAL(KIND=1)         S_DEPSCL_PLA(41, 2) @ 1312 in MODULE GLOBAL_DATA
                REAL(KIND=1)         EPSTANCL_TOT(41)    @ 1148 in MODULE GLOBAL_DATA
/@@GLOBAL_DATA!SIGMAEFFTIME/ REAL(KIND=1)         SIGMAEFFTIME(20, 20, 2) @ 0   (at line 22)
/@@GLOBAL_DATA!ETASTAT/ REAL(KIND=1)         ETASTAT(20, 20, 2)  @ 0   (at line 22)
                REAL(KIND=1)         FJ(41)              @ 984 in MODULE GLOBAL_DATA
                REAL(KIND=1)         FSPRNG              @ 14240 in MODULE GLOBAL_DATA
                REAL(KIND=1)         GEWB(41)            @ 820 in MODULE GLOBAL_DATA
                REAL(KIND=1)         GEWH(41)            @ 656 in MODULE GLOBAL_DATA
                REAL(KIND=1)         HAFT(41)            @ 492 in MODULE GLOBAL_DATA
                REAL(KIND=1)         PJ(41)              @ 328 in MODULE GLOBAL_DATA
                REAL(KIND=1)         PAXIAL(41)          @ 164 in MODULE GLOBAL_DATA
                REAL(KIND=1)         REIB(41)            @ 0   in MODULE GLOBAL_DATA

  Total stack size: 96 bytes
  Total static size: 112 bytes


   0147                                                                                          AT 0
   0148   !     ======                                                                           AT 0
   0149         End If                                                                           AT 0
   0150   !     ======                                                                           AT 0
   0151                                                                                          AT 0
   0152                                                                                          AT 0
   0153   !     ########                                                                         AT 0
   0154         Contains                                                                         AT 0
   0155   !     ########                                                                         AT 0
   0156                                                                                          AT 0
   0157                                                                                          AT 0
   0158   !       #################                                                              AT 0
   0159           Subroutine FGAUSL                                                              AT 0
   0160   !       #################                                                              AT 23
   0161                                                                                          AT 23
   0162                                                                                          AT 23
   0163   !         ******************************************************************           AT 23
   0164   !         *                                                                *           AT 23
   0165   !         * Test programs see  Tusubsys\\Fgausl                             *           AT 23
   0166   !         *                                                                *           AT 23
   0167   !         ******************************************************************           AT 23
   0168                                                                                          AT 23
   0169   !         Backsubstitution of upper triangular matrix                                  AT 23
   0170                                                                                          AT 23
   0171   !         Friction forces are calculated for section 2, 3, ..., m3+1                   AT 23
   0172                                                                                          AT 23
   0173   !         ------------------------------------------------------------------           AT 23
   0174                                                                                          AT 23
   0175             Implicit None                                                                AT 23
   0176                                                                                          AT 23
   0177             Double Precision             :: sum                                          AT 23
   0178             Double Precision , Parameter :: Small_Value = Tiny (sum)                     AT 23
   0179             Integer                      :: i                                            AT 23
   0180                                                                                          AT 23
   0181   !         ------------------------------------------------------------------           AT 23
   0182                                                                                          AT 23
   0183               ierr = 0                                                                   AT 23
   0184                                                                                          AT 2e
   0185   !         ================================================                             AT 2e
   0186             If ( Abs ( a (m31,m31) ) .lt. Small_Value ) Then                             AT 2e
   0187   !         ================================================                             AT bc
   0188                                                                                          AT bc
   0189                                                                                          AT bc
   0190               Call Fgausl_Error                                                          AT bc
   0191   !           #################                                                          AT c8
   0192                                                                                          AT c8
   0193                                                                                          AT c8
   0194               ierr = 1                                                                   AT c8
   0195                                                                                          AT d3
   0196               Return                                                                     AT d3
   0197   !           ++++++                                                                     AT d9
   0198                                                                                          AT d9
   0199   !         ======                                                                       AT d9
   0200             End If                                                                       AT d9
   0201   !         ======                                                                       AT d9
   0202                                                                                          AT d9
   0203             x (m31) = b (m31) / a (m31,m31)                                              AT d9
   0204                                                                                          AT 1c2
   0205   !         ******************                                                           AT 1c2
   0206             Do i  = m3, 2, -1                                                            AT 1c2
   0207   !         ******************                                                           AT 1f0
   0208                                                                                          AT 1f0
   0209               sum = b (i) - dot_product ( a (i,i+1:m31), x (i+1:m31) )                   AT 1f0
   0210                                                                                          AT 2cd
   0211   !           ============================================                               AT 2cd
   0212               If ( Abs ( a (i,i) ) .lt. Small_Value ) Then                               AT 2cd
   0213   !           ============================================                               AT 335
   0214                                                                                          AT 335
   0215                                                                                          AT 335
   0216                 Call Fgausl_Error                                                        AT 335
   0217   !             #################                                                        AT 359
   0218                                                                                          AT 359
   0219                                                                                          AT 359
   0220                 ierr = 1                                                                 AT 359
   0221                                                                                          AT 364
   0222                 Return                                                                   AT 364
   0223   !             ++++++                                                                   AT 36a
   0224                                                                                          AT 36a
   0225   !           ======                                                                     AT 36a
   0226               End If                                                                     AT 36a
   0227   !           ======                                                                     AT 36a
   0228                                                                                          AT 36a
   0229               x (i) = sum / a (i,i)                                                      AT 36a
   0230                                                                                          AT 3f1
   0231                                                                                          AT 3f1
   0232   !         ******                                                                       AT 3f1
   0233             End Do                                                                       AT 3f1


  Dump of all variables in FGAUSL

  Local Variables
                INTEGER(KIND=3)      I                   @ -32 (at line 182)
                REAL(KIND=2)         SUM                 @ -24 (at line 177)

  Total stack size: 56 bytes
  Total static size: 124 bytes


   0234   !         ******                                                                       AT 0
   0235                                                                                          AT 0
   0236                                                                                          AT 0
   0237   !       #####################                                                          AT 0
   0238           End Subroutine Fgausl                                                          AT 0
   0239   !       #####################                                                          AT 0
   0240                                                                                          AT 0
   0241                                                                                          AT 0
   0242   !       #######################                                                        AT 0
   0243           Subroutine Fgausl_Error                                                        AT 0
   0244   !       #######################                                                        AT 23
   0245                                                                                          AT 23
   0246                                                                                          AT 23
   0247             Implicit None                                                                AT 23
   0248                                                                                          AT 23
   0249             Integer :: i, j                                                              AT 23
   0250                                                                                          AT 23
   0251   !         --- System of equations is singular                                          AT 23
   0252                                                                                          AT 23
   0253             write   ( nwrite,  1000  )                                                   AT 23
   0254             write   ( nwrite, '(//)' )                                                   AT 62
   0255                                                                                          AT a1
   0256   !         *************                                                                AT a1
   0257             Do i = 2, m31                                                                AT a1
   0258   !         *************                                                                AT d8
   0259                                                                                          AT d8
   0260               write ( nwrite, 2010   )    i, b (i)                                       AT d8
   0261                                                                                          AT 18a
   0262   !         ******                                                                       AT 18a
   0263             End Do                                                                       AT 18a
   0264   !         ******                                                                       AT 19e
   0265                                                                                          AT 19e
   0266             write   ( nwrite, '(//)' )                                                   AT 19e
   0267                                                                                          AT 1dd
   0268   !         *************                                                                AT 1dd
   0269             Do i = 2, m31                                                                AT 1dd
   0270   !         *************                                                                AT 214
   0271                                                                                          AT 214
   0272               write ( nwrite, 2020 )      i, i                                           AT 214
   0273               write ( nwrite, 2030 ) ( a (i, j), j = i,m31 )                             AT 297
   0274                                                                                          AT 39f
   0275               write ( nwrite, '(/)' )                                                    AT 39f
   0276                                                                                          AT 3f6
   0277   !         ******                                                                       AT 3f6
   0278             End Do                                                                       AT 3f6


  Dump of all variables in FGAUSL_ERROR

  Local Variables
                INTEGER(KIND=3)      J                   @ -60 (at line 252)
                INTEGER(KIND=3)      I                   @ -64 (at line 252)

  Total stack size: 96 bytes
  Total static size: 669 bytes


   0279   !         ******                                                                       AT 414
   0280                                                                                          AT 414
   0281   !         ------------------------------------------------------------------           AT 414
   0282                                                                                          AT 414
   0283      1000   Format (////1x,130('=') /                                         &          AT 414
   0284                         1x,130('=') /                                         &          AT 414
   0285                         1x,130('=') ///                                       &          AT 414
   0286                         1x,'Programme stop in subroutine    F G A U S L   ',  &          AT 414
   0287                         1x,'called by subroutine AXREIB because'/             &          AT 414
   0288                         1x,'the system of equations to determine the local',  &          AT 414
   0289                         1x,'friction forces FJ is singular'///                &          AT 414
   0290                         1x,130('=') /                                         &          AT 414
   0291                         1x,130('=') /                                         &          AT 414
   0292                         1x,130('=') / )                                                  AT 414
   0293                                                                                          AT 414
   0294      2010   Format ( ' B (',i5,' ) = ',   d16.7 )                                        AT 414
   0295                                                                                          AT 414
   0296      2020   Format ( ' A (',i5,', j =',i5,', m31 ) : ' )                                 AT 414
   0297                                                                                          AT 414
   0298      2030   Format ( 1x, 10d16.7 )                                                       AT 414
   0299                                                                                          AT 414
   0300                                                                                          AT 414
   0301   !       ###########################                                                    AT 414
   0302           End Subroutine Fgausl_Error                                                    AT 414
   0303   !       ###########################                                                    AT 414
   0304                                                                                          AT 414
   0305                                                                                          AT 414
   0306   !     #####################                                                            AT 414
   0307         End Subroutine Axreib                                                            AT 414
   0308   !     #####################                                                            AT 414
   0309                                                                                          AT 414




End of Compilation - Clocked 0.01 seconds

[/quote]

7 Aug 2010 6:53 #6737

Can you supply a short and complete program please?

7 Aug 2010 7:09 #6738

were the programs not transmitted completely?

7 Aug 2010 8:52 #6739

No. There appears to be limit to the number of lines that you can post.

8 Aug 2010 7:44 #6740

Sorry, I did not know about this restriction. However, I have of course tried to cut the program(s) down, but failed to reproduce this particular behaviour of FTN95. Since I have observed this behaviour many times, I'll continue to produce a shorter demo program.

Best regards,

Klaus

Please login to reply.