WRF-Yumin

T=getnc(file,'T'); theta = T + 300; PB=getnc(file,'PB'); pres = P + PB; temp = theta.*(pres./p1000mb).^kappa; PH=getnc(file,'PH'); PHB=getnc(file,'PHB'); phtotal=PH+PHB;

PSFC=getnc(file,'PSFC'); MU=getnc(file,'MU'); MUB=getnc(file,'MUB'); mu=MU+MUB;

ZNU=getnc(file,'ZNU'); ZNW=getnc(file,'ZNW'); FNM=getnc(file,'FNM'); FNP=getnc(file,'FNP');

% loadcvars.m CF1=getnc(file,'CF1'); CF2=getnc(file,'CF2'); CF3=getnc(file,'CF3'); P_TOP=getnc(file,'P_TOP');

% compute_alt.m

% geopotential phi=zeros(nz,ny,nx);

for kk=1:nz phi(kk,:,:) = 0.5*( ( squeeze(phtotal(kk,:,:)) - ... Rd.*squeeze(temp(kk,:,:)).*log( (mu*ZNU(kk)+P_TOP)./(mu*ZNW(kk) +P_TOP) ) ) + ... ( squeeze(phtotal(kk+1,:,:)) - ... Rd.*squeeze(temp(kk,:,:)).*log( (mu*ZNU(kk)+P_TOP)./(mu*ZNW(kk+1)+P_TOP) ) ) ); z(kk,:,:) = phi(kk,:,:) / g; alt(kk) = mean(mean(z(kk,:,:))); end