iniconstant.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/.
!||====================================================================
!||    iniconstant   ../starter/source/system/iniconstant.F
!||--- called by ------------------------------------------------------
!||    starter0      ../starter/source/starter/starter0.F
!||====================================================================
      SUBROUTINE iniconstant
C      
#include "implicit_f.inc"   
c-------------------------------------------------------------  
c++++ les constantes suivantes ne sont pas dans le common
c     mais peuvent y etre ajoutees si necessaire
c-------------------------------------------------------------  
C--------------------------------------
C      Entiers
C--------------------------------------
       infinity       = 1e20
       zero           = 0.
       one            = 1.
       two            = 2.
       three          = 3.
       four           = 4.
       five           = 5.
       six            = 6.
       seven          = 7.
       eight          = 8.
       nine           = 9.
       ten            = 10.
       eleven         = 11.
       twelve         = 12.
       fourteen       = 14.
       fifteen        = 15.
       sixteen        = 16.
       twenty         = 20.
       twenty4        = 24.
       twenty5        = 25.
       twenty7        = 27.
       thirty         = 30.
       thirty2        = 32.
       thirty6        = 36.
       fourty         = 40.
       fourty5        = 45.
       fourty8        = 48.
       fifty          = 50.
       sixty4         = 64.
       eighty         = 80.
       ninty          = 90.
       eighty16       = 96.
       eighty19       = 99.
       hundred        = 100.
       hundred28      = 128.
       hundred50      = 150.
       hundred80      = 180.
       three100       = 300.
       twohundred25   = 225.
       twohundred95   = 295.
       twohundred98   = 298.
       twohundred93   = 293.
       threehundred24 = 324.
       fivehundred12  = 512.
       three1000      = 3000.
C--------------------------------------
C      10^n
C--------------------------------------
       ep02  = hundred     
       ep03  = ep02 * ten    
       ep04  = ep03 * ten   
       ep05  = ep04 * ten    
       ep06  = ep05 * ten     
       ep07  = ep06 * ten              
       ep08  = ep07 * ten              
       ep09  = ep08 * ten              
       ep10  = ep09 * ten              
       ep11  = ep10 * ten           
       ep12  = ep11 * ten           
       ep13  = ep12 * ten           
       ep14  = ep13 * ten           
       ep15  = ep14 * ten           
       ep16  = ep15 * ten             
       ep17  = ep16 * ten             
       ep18  = ep17 * ten             
       ep19  = ep18 * ten             
       ep20  = ep19 * ten             
       ep29  = ep20 * ep09           
       ep30  = ep20 * ep10            
       ep31  = ep30 * ten             
       ep35  = ep30 * ep05             
       ep36  = ep30 * ep06             
       ep37  = ep30 * ep07               
C--------------------------------------
C      10^-n
C--------------------------------------
       em01 = one / ten  
       em02 = one / ep02 
       em03 = one / ep03 
       em04 = one / ep04 
       em05 = one / ep05 
       em06 = one / ep06 
       em07 = one / ep07 
       em08 = one / ep08 
       em09 = one / ep09 
       em10 = one / ep10 
       em11 = one / ep11 
       em12 = one / ep12 
       em13 = one / ep13 
       em14 = one / ep14 
       em15 = one / ep15 
       em16 = one / ep16 
       em17 = one / ep17 
       em18 = one / ep18 
       em19 = one / ep19 
       em20 = one / ep20 
       em30 = one / ep30 
       em35 = one / ep35 
C--------------------------------------
C      m10^-n
C--------------------------------------
       zep1     = em01 
       zep2     = two / ten 
       zep3     = three / ten 
       zep4     = four / ten 
       zep5     = five / ten 
       zep6     = one - zep4 
       zep7     = one - zep3 
       zep8     = one - zep2 
       zep9     = one - zep1 
C
       zep01    = em02
       zep02    = two / ep02
       zep03    = three / ep02
       zep04    = four / ep02
       zep05    = five / ep02
       zep06    = zep1 - zep04
       zep07    = zep1 - zep03
       zep08    = zep1 - zep02
       zep09    = zep1 - zep01   
C
       twoem3   = two / ep03
       threeem3 = three / ep03
       fourem3  = four / ep03
       fiveem3  = five / ep03
       sixem3   = zep01 - fourem3
       sevenem3 = zep01 - threeem3 
       eightem3 = zep01 - twoem3
       nineem3  = zep01 - em03      
C
       twoem4   = two / ep04
       threeem4 = three / ep04
       fourem4  = four / ep04
       fiveem4  = five / ep04
       sixem4   = em03 - fourem4 
       sevenem4 = em03 - threeem4 
       eightem4 = em03 - twoem4
       nineem4  = em03 - em04             
C
       twoem5   = two / ep05
       threeem5 = three / ep05
       fourem5  = four / ep05
       fiveem5  = five / ep05
       sixem5   = em04 - fourem5
       sevenem5 = em04 - threeem5 
       eightem5 = em04 - twoem5
       nineem5  = em04 - em05             
C
       twoem6   = two / ep06
       threeem6 = three / ep06
       fourem6  = four / ep06
       fiveem6  = five / ep06
       sixem6   = em05 - fourem6
       sevenem6 = em05 - threeem6
       eightem6 = em05 - twoem6
       nineem6  = em05 - em06        
C--------------------------------------
C      m10^n
C--------------------------------------
       twoep30   = two * ep30
       threeep30 = three * ep30
       fourep30  = four * ep30
       zep9ep30  = nine  *ep29
C--------------------------------------
C      Fractions
C--------------------------------------
       half          = one/two         
       third         = one/three        
       fourth        = one/four       
       one_fifth     = one/five         
       one_over_6    = one/six          
       one_over_8    = one/eight         
       one_over_9    = one/nine         
       one_over_12   = one/twelve        
       one_over_16   = one/sixteen        
       one_over_64   = one/sixty4    
       one_over_512  = one/fivehundred12 
       two_third     = two/three      
       three_half    = three/two      
       three_over_4  = three/four    
       trhee_over_14 = three/fourteen  
       four_over_3   = four/three    
       four_over_5   = four/five     
       five_over_3   = five/three      
       five_over_4   = five/four     
       five_over_6   = five/six        
       six_over_5    = six/five        
       seven_half    = seven/two       
       seven_over_5  = seven/five       
       seven_over_9  = seven/nine       
       nine_over_2   = nine/two       
       nine_over_20  = nine/twenty      
C--------------------------------------
C      0.xxxx
C--------------------------------------
       zep11     = zep1 + zep01
       zep12     = zep1 + zep02
       zep13     = zep1 + zep03
       zep14     = zep1 + zep04
       zep15     = zep1 + zep05
       zep16     = zep1 + zep06
       zep17     = zep1 + zep07
       zep18     = zep1 + zep08
       zep19     = zep1 + zep09
C
       zep21     = zep2  + zep01
       zep23     = zep2  + zep03
       zep25     = zep2  + zep05
       zep26     = zep2  + zep06
       zep213    = zep21 + threeem3
       zep29     = zep2  + zep09
C
       zep33     = zep3  + zep03
       zep34     = zep3  + zep04
       zep35     = zep3  + zep05
       zep37     = zep3  + zep07
       zep38     = zep3  + zep08
       zep333    = zep33 + threeem3
       zep373    = zep37 + threeem3         
C
       zep41     = zep4  + zep01
       zep4187   = zep41 + eightem3 + sevenem4
       zep43     = zep4  + zep03
       zep45     = zep4  + zep05
       zep46     = zep4  + zep06
       zep47     = zep4  + zep07
       zep475    = zep4  + zep07 + fiveem3
       zep499    = zep4  + zep09 + nineem3     
C
       zep62     = one   - zep38
       zep65     = one   - zep35
       zep66     = one   - zep34 
       zep67     = one   - zep33
c
       zep71     = one   - zep29
       zep74     = one   - zep26         
       zep707    = one   - zep3 + sevenem3
C
       zep81     = one   - zep19
       zep82     = one   - zep18
       zep84     = one   - zep16         
       zep85     = one   - zep15
       zep87     = one   - zep13
       zep89     = one   - zep11
C
       zep92     = one   - zep08
       zep93     = one   - zep07       
       zep94     = one   - zep06
       zep95     = one   - zep05 
       zep96     = one   - zep04
       zep97     = one   - zep03
       zep98     = one   - zep02
       zep99     = one   - zep01
       zep995    = one   - fiveem3
       zep997    = one   - threeem3
       zep999    = one   - em03
       zep9999   = one   - em04
       zep99999  = one   - em05
       zep931    = zep93 + em03
c
       zep015    = zep15 / ten
       zep015625 = one / sixty4
       zep025    = one / fourty
       zep037    = zep03 + sevenem3                                          
       zep072169 = zep07 + twoem3 + em04 + sixem5 + nineem6
       zep075    = three / fourty
C
       onep01    = one + zep01
       onep001   = one + em03
       onep0001  = one + em04        
       onep00001 = one + em05
       onep02    = one + zep02
       onep03    = one + zep03
       onep04    = one + zep04
       onep05    = one + zep05
       onep06    = one + zep06
       onep08    = one + zep08
       onep09    = one + zep09
C
       onep1   = one    + zep1
       onep11  = onep1  + zep01 
       onep12  = onep1  + zep02
       onep13  = onep1  + zep03
       onep14  = onep1  + zep04
       onep15  = onep1  + zep05 
       onep16  = onep1  + zep06
       onep17  = onep1  + zep07
       onep18  = onep1  + zep08
       onep19  = onep1  + zep09
       onep136 = onep13 + sixem3
       onep138 = onep13 + eightem3
C
       onep2   = one + zep2 
       onep21  = onep2  + zep01
       onep23  = onep2  + zep03
       onep24  = onep2  + zep04
       onep25  = five / four      
       onep26  = onep2  + zep06
       onep27  = onep2  + zep07
       onep28  = onep2  + zep08        
       onep251 = onep25 + em03
       onep252 = onep25 + twoem3
C
       onep3   = one + zep3
       onep31  = onep3  + zep01  
       onep32  = onep3  + zep02
       onep33  = onep3  + zep03
       onep34  = onep3  + zep04
       onep35  = onep3  + zep05
       onep36  = onep3  + zep06
       onep37  = onep3  + zep07
       onep38  = onep3  + zep08
       onep39  = onep3  + zep09
       onep333 = onep33 + threeem3
       onep3348= onep33 + fourem3 + eightem4
       onep369 = onep36 + nineem3
C
       onep4   = one + zep4
       onep41  = onep4  + zep01
       onep42  = onep4  + zep02
       onep44  = onep4  + zep04         
       onep45  = onep4  + zep05
       onep46  = onep4  + zep06
       onep48  = onep4  + zep08
       onep49  = onep4  + zep09
       onep414 = onep41 + fourem3
       onep484 = onep48 + fourem3
C
       onep5   = three_half
       onep51  = onep5 + zep01
       onep52  = onep5 + zep02
       onep57  = onep5 + zep07
       onep59  = onep5 + zep09
C
       onep6   = one + zep6
       onep62  = onep6 + zep02
       onep65  = onep6 + zep05
       onep66  = onep6 + zep06
       onep67  = onep6 + zep07
       onep68  = onep6 + zep08 
C
       onep7   = one + zep7
       onep72  = onep7 + zep02
       onep73  = onep7 + zep03        
       onep75  = onep7 + zep05                
       onep78  = onep7 + zep08
       onep732 = onep73 + twoem3
C
       onep8    = one + zep8
       onep82   = onep8 + zep02
       onep84   = onep8 + zep04     
       onep88   = onep8 + zep08
       onep833  = onep8 + zep03 + threeem3
       onep8688 = onep8 + zep06 + eightem3 + eightem4
C
       onep9   = one + zep9
       onep92  = onep9 + zep02
       onep96  = onep9 + zep06
       onep98  = onep9 + zep08
C     
       twop03 = two + zep03
       twop05 = two + zep05
       twop06 = two + zep06        
       twop07 = two + zep07
       twop08 = two + zep08
       twop1  = two + zep1    
       twop2  = two + zep2
       twop30 = two + zep3
       twop4  = two + zep4
       twop5  = five / two
       twop60 = two + zep6
       twop8  = two + zep8
C
       twop11 = twop1 + zep01
       twop12 = twop1 + zep02
       twop13 = twop1 + zep03
       twop16 = twop1 + zep06    
       twop19 = twop1 + zep09
       twop23 = twop2 + zep03
       twop24 = twop2 + zep04
       twop25 = twop2 + zep05
       twop33 = twop30 + zep03
       twop38 = twop30 + zep08
       twop45 = twop4 + zep05
       twop47 = twop4 + zep07
       twop52 = twop5 + zep02
       twop83 = twop8 + zep03        
       twop84 = twop8 + zep04        
       twop85 = twop8 + zep05
C
       twop185 = twop1 + zep08 + fiveem3  
       twop261 = twop2 + zep06 + em03    
       twop645 = twop60 + zep04 + fiveem3
C
       threep02   = three + zep02
       threep06   = three + zep06
       threep14   = three + zep14
       threep15   = three + zep15
       threep25   = three + zep2 + zep05
       threep27   = three + zep2 + zep07
       threep44   = three + zep4 + zep04
       threep46   = three + zep4 + zep06
       threep52   = three + zep5 + zep02
       threep72   = three + zep7 + zep02
       threep76   = three + zep7 + zep06
       threep8    = three + zep8
       threep94   = three + zep9 + zep04
       threep266  = three + zep2 + zep06 + sixem3
       threep464  = three + zep4 + zep06 + fourem3
       threep4641 = threep464 + em04
C
       fourp186 = four + zep1 + zep08 + sixem3
       fourp40  = four + zep4
       fourp76  = four + zep7 + zep06
C
       fivep5  = five + zep5
       fivep56 = five + zep5 + zep06
       fivep76 = five + zep7 + zep06
C
       sip2832 = six + zep2 + zep08 + threeem3 + twoem4 
C       
       ninep24  = nine + zep2 + zep04
       ninep38  = nine + zep3 + zep08
       ninep793 = nine + zep7 + zep09 + threeem3
 
       elevenp5   = eleven + zep5
       twelvep6   = twelve + zep6
       fourteenp2 = fourteen + zep2
       sixteenp5  = sixteen + zep5 
       o88p9844   = 88. + zep9 + zep08 + fourem3 + fourem4
C
       onep3em5     = onep3 * em05
       onep001em32  = 1001. * em35
       ninep637em5  = 9637. * em08
       eightp314em5 = 8314. * em08        
C
       onep1ep30 = onep1 * ep30
       onep2ep30 = onep2 * ep30 
C
       zep999ep31 = zep999 * ep31
C--------------------------------------
C      Decimaux -> fractions
C--------------------------------------
       zep426667      = hundred28 / three100
       zep222         = ten / fourty5
       zep00666666667 = two / three100
       zep01666666667 = five / three100
       zep444         = hundred / twohundred25
       onep222        = one + zep222
       onep8333       = one + twenty5 / thirty
       twop444        = two + zep444                  
       twop6666666667 = eight / three 
       twop6667       = twop6666666667
       fivep333       = sixteen / three   
C--------------------------------------
C--------------------------------------
       zep54      = one - zep46
       zep1381966 = zep13 + eightem3 + em04 + nineem5 + sixem6 + six*em07  
       zep5854102 = zep5 + zep08 + fiveem3 + fourem4 + em05 + two*em07      
C--------------------------------------
C       Constantes physiques
C--------------------------------------
       r_igc     = 8.314472
       stefboltz = 5.6704 / ep08
C--------------------------------------
C       Constantes mathematiques
C--------------------------------------
       sqr2    = sqrt(two)
       sqr3    = sqrt(three)
       sqr5    = sqrt(five)
       sqr6    = sqrt(six)
       sqr7    = sqrt(seven)
       sqr3_2  = sqrt(three_half)
       sqr4_3  = two * sqrt(third)
       sqr16_3 = two * sqr4_3
       pi      = atan2(zero,-one)
       invpi   = one / pi
C--------------------------------------
C      alias
C--------------------------------------
       ep9     = ep09
       zep80   = zep8
       zep90   = zep9
       unp20   = onep2 
       unp30   = onep3
       unp70   = onep7   
       unp90   = onep9
       unem3   = em03
       unem4   = em04
       unem5   = em05
       unem6   = em06
       em3     = em03
       em4     = em04
       em5     = em05
       em6     = em06
       em7     = em07
       em8     = em08
       em9     = em09
       twoem2  = zep02
       fiveem2 = zep05
       nineem2 = zep09
C--------------------------------------
       RETURN    
       END