HW1-Sil

esults from 1 to 11. didn't start 12 yet.

matlab code

--- Get this [|COARE soundings .nc data file]It is 1726 soundings on a 5mb grid 1000-100 mb. (181 values) Get the IDL or [|Matlab] codes: main ones are esat, mixrat (or saturation mixing ratio), entropy, invert_entropy Read the COARE soundings into memory. Notice array sizes.

% % get data : array [181,1726] % p = nc_varget (fname, 'p'); % P is constant 5mb interval, 1000-100.[1000,995,990,...110,105,100] T = nc_varget (fname, 'T'); RH = nc_varget (fname, 'RH'); % [181,1726] irtemp = nc_varget (fname, 'IRTEMP'); % [1726] %-- Compute q (mixing ratio) = {[0.622 *esat(T) / (p-esat(T)) ] *RH} /100 = mixingratio(esat(T),p) *RH / 100

% #4. compute q : traditionally w and ws are used for q (q is the same as qv) and qsat % ws = 0.622 .* es ./ (press-es); % w=(ws.*rh)/100.; % what this value means? kg/kg

*Figures : mean profile of T, q(and q_sat) and RH Compute the entropy, total water, and pressure of a near-surface parcel. Call these s0, qt0, p0. %===================================================================== % #6. Compute s0, qt0, and p0 %- T0=T_mean(1); p0=p_mean(1); qt0=q_mean(1); s0=entropy(T0, p0, qt0);

% entropy of the actual sounding (qtot=q=qv, b/c no condensate) ent = entropy(T,p,q);

% SATURATION ent of the sounding (qtot=qsat) entsat = entropy(T,p,mixingratio(esat(T),p));

% dry version (qtot=> q=0) entdry = entropy(T,p,1e-7*q);

Figures: mean profiles of Entropy



Same as previous figure, but y axis is height in this figure.







working on #12.