Skip to content
Snippets Groups Projects
param.f90 1.62 KiB
Newer Older
Victoria Suiazova's avatar
Victoria Suiazova committed
module PARAM
      IMPLICIT REAL (A-H,O-Z)
      real, parameter ::  RO = 1.2
      real, parameter ::  CP =8.4
      real, parameter ::  AKA=.40E0
      real, parameter ::  AP0=.40E0
      real, parameter ::  G=9.81E0
      real, parameter ::  A0=1.15E0
      real, parameter ::  A6=3.5E0
      real, parameter ::  G4=16.0E0
      real, parameter ::  G10=16.0E0
      real, parameter ::  B4=4.7E0
      real, parameter ::  ALFAM=.0144E0
      real, parameter ::  BETAM=.111E0
      real, parameter ::  AN=.000015E0
      real, parameter ::  P4=.71E0
      real, parameter ::  H1=10.0E0
      real, parameter ::  X8=16.3E0
      real, parameter ::  AN1=5.0E0/6.0E0
      real, parameter ::  AN2=.45E0
      real, parameter ::  AL1=AKA*P4
      real, parameter ::  G0=1.2
      real, parameter ::  AL2=(.14E0*(30.0E0**AN2))*(P4**.8E0)
      real, parameter ::  A2=ALOG(H1*(G/ALFAM))
      real, parameter ::  A3=BETAM*AN*A2
      real, parameter ::  R0=.9E0/B4
!C*      real, parameter ::  AN5=(A6/A0)**4
!C*      real, parameter ::  D1=(2.0E0*G10-AN5*G4-SQRT((AN5*G4)**2+4.0E0*AN5*G10*(G10-G4)))/(2.0E0*G10**2)
!C*      real, parameter ::  Y10=(1.0E0-G4*D1)**.25E0
!C*      real, parameter ::  X10=(1.0E0-G10*D1)**.5E0
!C*      real, parameter ::  P1=2.0E0*ATAN(Y10)+ALOG((Y10-1.0E0)/(Y10+1.0E0))
!C*      real, parameter ::  P0=ALOG((X10-1.0E0)/(X10+1.0E0))
      
!      AN5=(A6/A0)**4
!      D1=(2.0E0*G10-AN5*G4-SQRT((AN5*G4)**2+4.0E0*AN5*G10*(G10-G4)))/(2.0E0*G10**2)
!      Y10=(1.0E0-G4*D1)**.25E0
!      X10=(1.0E0-G10*D1)**.5E0
!      P1=2.0E0*ATAN(Y10)+ALOG((Y10-1.0E0)/(Y10+1.0E0))
!      P0=ALOG((X10-1.0E0)/(X10+1.0E0))

end module PARAM