T_HFBULKTC: a program to test hfbulktc.m using COARE test data T_HFBULKTC is a program to test hfbulktc.m using COARE test data (file test2_5b.dat). NOTE: Make sure that you use the COARE parameters located in the files as_consts.m and cool_skin.m when conducting the test. VARIBLES: ur = wind speed [m/s] measured at height zr [m] Ta = air temperature [C] measured at height zt [m] rh = relative humidity [%] measured at height zq [m] Pa = air pressure [mb] Ts = sea surface temperature [C] sal = salinity [psu (PSS-78)] dlw = downwelling (INTO water) longwave radiation [W/m^2] dsw = measured insolation [W/m^2] nsw = net shortwave radiation INTO the water [W/m^2] rain = rain rate [mm/hr]
0001 % T_HFBULKTC: a program to test hfbulktc.m using COARE test data 0002 % T_HFBULKTC is a program to test hfbulktc.m using COARE test data 0003 % (file test2_5b.dat). NOTE: Make sure that you use the COARE 0004 % parameters located in the files as_consts.m and cool_skin.m when 0005 % conducting the test. 0006 % 0007 % VARIBLES: 0008 % 0009 % ur = wind speed [m/s] measured at height zr [m] 0010 % Ta = air temperature [C] measured at height zt [m] 0011 % rh = relative humidity [%] measured at height zq [m] 0012 % Pa = air pressure [mb] 0013 % Ts = sea surface temperature [C] 0014 % sal = salinity [psu (PSS-78)] 0015 % dlw = downwelling (INTO water) longwave radiation [W/m^2] 0016 % dsw = measured insolation [W/m^2] 0017 % nsw = net shortwave radiation INTO the water [W/m^2] 0018 % rain = rain rate [mm/hr] 0019 0020 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0021 % 4/7/99: version 1.2 (contributed by AA) 0022 % 8/5/99: version 2.0 0023 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0024 0025 % change directory to yours 0026 load c:\flux_fsu\test\test2_5b.dat 0027 0028 % parse data 0029 ur = test2_5b(:,2); 0030 zr = 15; 0031 Ta = test2_5b(:,4); 0032 zt = 15; 0033 Pa = 1008*ones(size(ur)); 0034 q = test2_5b(:,5); 0035 rh = q./qsat(Ta,Pa)/10; 0036 zq = 15; 0037 Ts = test2_5b(:,14); 0038 sal = 30*ones(size(ur)); 0039 dsw = test2_5b(:,9); 0040 nsw = 0.945*dsw; 0041 dlw = test2_5b(:,10); 0042 rain= test2_5b(:,11); 0043 0044 i=[1:length(ur)]; 0045 % i=[1:10]; 0046 % NO cool-skin; compare to COARE output in file no_skin.out 0047 A1 = hfbulktc(ur(i),zr,Ta(i),zt,rh(i),zq,Pa(i),Ts(i)); 0048 0049 % YES cool-skin; compare to COARE output file yes_skin.out 0050 A2 = hfbulktc(ur(i),zr,Ta(i),zt,rh(i),zq,Pa(i),Ts(i), ... 0051 sal(i),dlw(i),dsw(i),nsw(i)); 0052