NCEPLIBS-w3emc  2.9.3
w3fa04.f
Go to the documentation of this file.
1 C> @file
2 C> @brief Compute standard pressure, temp, pot temp.
3 C> @author James McDonell @date 1974-06-01
4 
5 C> Computes the standard pressure, temperature, and poten-
6 C> tial temperature given the height in meters (<32 km). For
7 C> the pressure and temperature the results duplicate the values in
8 C> the U.S. standard atmosphere (1962), which is the icao standard
9 C> atmosphere to 54.7487 mb (20 km) and the proposed extension to
10 C> 8.68 mb (32 km). For potential temperature a value of 2/7 is
11 C> used for rd/cp.
12 C>
13 C> Program history log:
14 C> - James McDonell 1974-06-01
15 C> - Ralph Jones 1984-07-05 Change to ibm vs fortran.
16 C> - Ralph Jones 1990-04-27 Change to cray cft77 fortran.
17 C>
18 C> @param[in] HEIGHT Height in meters.
19 C> @param[out] PRESS Standard pressure in millibars.
20 C> @param[out] TEMP Temperature in degrees kelvin.
21 C> @param[out] THETA Potential temperature in degrees kelvin.
22 C>
23 C> @note Not valid for heights greater than 32 km. declare all parameters
24 C> as type real*4.
25 C>
26 C> @author James McDonell @date 1974-06-01
27 
28  SUBROUTINE w3fa04(HEIGHT,PRESS,TEMP,THETA)
29 C
30  REAL M0
31 C
32  DATA
33  *g /9.80665/,
34  *rstar /8314.32/,
35  *m0 /28.9644/,
36  *piso /54.7487/,
37  *ziso /20000./,
38  *salp /-.0010/,
39  *tstr /216.65/,
40  *ptrop /226.321/,
41  *alp /.0065/,
42  *t0 /288.15/,
43  *pzero /1013.25/
44 C
45  rovcp = 2.0 / 7.0
46  r = rstar/m0
47  IF (height.GT.ziso) GO TO 100
48  IF (height.GT.11000.) GO TO 200
49 C
50 C COMPUTE IN TROPOSPHERE
51 C
52  temp = t0 - height * alp
53  press = pzero * ((1.0 - ((alp/t0) * height)) ** (g/(alp * r)))
54  GO TO 300
55 C
56 C COMPUTE LAPSE RATE = -.0010 CASES
57 C
58  100 CONTINUE
59  d = height - ziso
60  press = piso * ((1.-(( salp /tstr) * d )) ** (g/( salp * r)))
61  temp = tstr - d * salp
62  GO TO 300
63 C
64 C COMPUTE ISOTHERMAL CASES
65 C
66  200 CONTINUE
67  d = exp((height - 11000.0) / ((r / g) * tstr))
68  press = ptrop / d
69  temp = tstr
70 C
71  300 CONTINUE
72  theta = temp * ((1000.0 / press) ** rovcp)
73  RETURN
74  END
w3fa04
subroutine w3fa04(HEIGHT, PRESS, TEMP, THETA)
Computes the standard pressure, temperature, and poten- tial temperature given the height in meters (...
Definition: w3fa04.f:29