File: ew_gfs.f

package info (click to toggle)
flextra 5.0-2.1
  • links: PTS
  • area: main
  • in suites: jessie, jessie-kfreebsd
  • size: 860 kB
  • ctags: 402
  • sloc: fortran: 6,987; makefile: 55; sh: 17
file content (23 lines) | stat: -rw-r--r-- 703 bytes parent folder | download | duplicates (7)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
      REAL FUNCTION EW(X)
C     ****************************************************************
C     SAETTIGUNGSDAMPFDRUCK UEBER WASSER IN PA. X IN KELVIN.
C     NACH DER GOFF-GRATCH-FORMEL.
C     ****************************************************************
      EW=0.
      IF(X.LE.0.) then
        WRITE(*,*) 'TEMP: ',X
        STOP 'SORRY: T NOT IN [K]'
      ENDIF
      Y=373.16/X
      A=-7.90298*(Y-1.)
      A=A+(5.02808*0.43429*ALOG(Y))
      C=(1.-(1./Y))*11.344
      C=-1.+(10.**C)
      C=-1.3816*C/(10.**7)
      D=(1.-Y)*3.49149
      D=-1.+(10.**D)
      D=8.1328*D/(10.**3)
      Y=A+C+D
      EW=1013.246*(10.**Y)*100.       ! Saettigungsdampfdruck in Pa
      RETURN
      END