9 subroutine slat_geos5_tamp_v0(im, tau_amp, xlatdeg, tau_gw)
15 real :: tau_amp, xlatdeg(im), tau_gw(im)
16 real :: latdeg, flat_gw, tem
23 latdeg = abs(xlatdeg(i))
24 if (latdeg < 15.3)
then
25 tem = (latdeg-3.0) / 8.0
26 flat_gw = 0.75 * exp(-tem * tem)
27 if (flat_gw < 1.2 .and. latdeg <= 3.0) flat_gw = 0.75
28 elseif (latdeg < 31.0 .and. latdeg >= 15.3)
then
30 elseif (latdeg < 60.0 .and. latdeg >= 31.0)
then
31 tem = (latdeg-60.0) / 23.0
32 flat_gw = 0.50 * exp(- tem * tem)
33 elseif (latdeg >= 60.0)
then
34 tem = (latdeg-60.0) / 70.0
35 flat_gw = 0.50 * exp(- tem * tem)
37 tau_gw(i) = tau_amp*flat_gw
43 subroutine slat_geos5_v0(im, xlatdeg, tau_gw)
52 real,
parameter :: tau_amp = 100.e-3
53 real :: trop_gw, flat_gw
61 if (-15.3 < latdeg .and. latdeg < 15.3)
then
62 flat_gw = trop_gw*exp(-( (abs(latdeg)-3.)/8.0)**2)
63 if (flat_gw < 1.2 .and. abs(latdeg) <= 3.) flat_gw = trop_gw
64 else if (latdeg > -31. .and. latdeg <= -15.3)
then
66 else if (latdeg < 31. .and. latdeg >= 15.3)
then
68 else if (latdeg > -60. .and. latdeg <= -31.)
then
69 flat_gw = 0.50*exp(-((abs(latdeg)-60.)/23.)**2)
70 else if (latdeg < 60. .and. latdeg >= 31.)
then
71 flat_gw = 0.50*exp(-((abs(latdeg)-60.)/23.)**2)
72 else if (latdeg <= -60.)
then
73 flat_gw = 0.50*exp(-((abs(latdeg)-60.)/70.)**2)
74 else if (latdeg >= 60.)
then
75 flat_gw = 0.50*exp(-((abs(latdeg)-60.)/70.)**2)
77 tau_gw(i) = tau_amp*flat_gw