25 real,
parameter ::
s1=0.30731408,
s2=6.707e-9,
s3=-8.1899e-5,
sa=3.227,&
27 real,
parameter ::
t1=0.29585798,
t2=1.002e-9,
t3=-3.1658e-5,
ta=3.426,&
35 type(ocean_profile_type),
intent(inout) :: station
39 real :: dpth_orig, dpth_new,tdrop
42 if (.not. station%accepted)
return 44 fix_depth = int(station%fix_depth)
45 select case(fix_depth)
53 dpth_orig = station%depth(k)
55 dpth_new = (1.0417*dpth_orig) - (75.096*(1.0-((1.0-(0.0002063*dpth_orig)))**0.5))
56 station%depth(k)=dpth_new
59 station%fix_depth=-1.0
63 dpth_orig = station%depth(k)
65 if (dpth_orig .le. 250.0)
then 66 dpth_new = dpth_orig*0.9572
67 else if (dpth_orig .le. 500.)
then 68 dpth_new = dpth_orig*0.9565
69 else if (dpth_orig .le. 750.0)
then 70 dpth_new = dpth_orig*0.9558
71 else if (dpth_orig .le. 1000.)
then 72 dpth_new = dpth_orig*0.9550
73 else if (dpth_orig .le. 1250.0)
then 74 dpth_new = dpth_orig*0.9542
75 else if (dpth_orig .le. 1500.0)
then 76 dpth_new = dpth_orig*0.9533
78 dpth_new = dpth_orig*0.9524
80 station%depth(k)=dpth_new
83 station%fix_depth=-1.0
86 dpth_orig = station%depth(k)
88 tdrop=(
s1*dpth_orig +
s2) -
s3 89 dpth_new =
sa*tdrop +
sb*tdrop*tdrop
90 station%depth(k)=dpth_new
93 station%fix_depth=-1.0
96 dpth_orig = station%depth(k)
98 tdrop=(
t1*dpth_orig +
t2) -
t3 99 dpth_new =
ta*tdrop +
tb*tdrop*tdrop
100 station%depth(k)=dpth_new
103 station%fix_depth=-1.0
subroutine xbt_drop_rate_adjust(station)
real, parameter, public missing_value