Please find the complete problem description in the attached file.
a=-12 degrees celcious φ is latitude in radians (which should vary from the North to the South pole), b1=40oC, z is the height above sea level (which should vary from the surface to 22km), zT is the approximate tropopause height (you may assume that this is constant at 12km).
Using the thermal wind shear balance equation derive the corresponding zonal wind field in the Northern hemisphere. This will require you to calculate the meridional gradient of the temperature to determine the zonal velocities vertical derivative. Using the assumption that the surface wind is zero everywhere you can then derive the zonal wind field using Euler's method. Remember that thermal wind shear balance is not valid near equator so only derive the zonalwind up to 20oN.© BrainMass Inc. brainmass.com October 10, 2019, 2:58 am ad1c9bdddf
Some explanations are in comments in the attached Matlab script.
fi = [-90:1:90]*pi/180; % latitude in radians
z = [0:0.1:22]; % height in km
[fi2, z2] = meshgrid(fi,z);
%z2 = z' * ones(1,length(fi)); % altitude ...
This solution helps with a thermodynamics that involves plotting temperature, longitude and latitude.