Global Index (short | long) | Local contents | Local Index (short | long)
Try some plots removing the zonal average.
This script calls | |
---|---|
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/ct/data filin = 'wct.nc'; nc = netcdf(filin, 'nowrite'); vtw = nc{'VT'}(1:3,:,:,:); kew = nc{'VVPUU'}(1:3,:,:,:); uw = nc{'U'}(1:3,:,:,:); vw = nc{'V'}(1:3,:,:,:); airw = nc{'AIRT'}(1:3,:,:,:); psw = nc{'PS'}(1:3,:,:); omw = nc{'OMEGA'}(1:3,:,:); zw = nc{'Z3'}(1:3,:,:,:); nc = close(nc); filin = 'cct.nc'; nc = netcdf(filin, 'nowrite'); vtc = nc{'VT'}(1:3,:,:,:); kec = nc{'VVPUU'}(1:3,:,:,:); uc = nc{'U'}(1:3,:,:,:); vc = nc{'V'}(1:3,:,:,:); airc = nc{'AIRT'}(1:3,:,:,:); psc = nc{'PS'}(1:3,:,:); omc = nc{'OMEGA'}(1:3,:,:); zc = nc{'Z3'}(1:3,:,:,:); lat = nc{'lat'}(:); lon = nc{'lon'}(:); hybm = nc{'hybm'}(:); hyam = nc{'hyam'}(:); levs = nc{'lev'}(:); P0 = nc{'P0'}(:); nc = close(nc); lev = 500; uw500 = atlev(uw, lev, psw, hyam, hybm, P0); vw500 = atlev(vw, lev, psw, hyam, hybm, P0); vorw500 = sph_curl(squeeze(mean(uw500)), squeeze(mean(vw500)), lon, lat); uc500 = atlev(uc, lev, psc, hyam, hybm, P0); vc500 = atlev(vc, lev, psc, hyam, hybm, P0); vorc500 = sph_curl(squeeze(mean(uc500)), squeeze(mean(vc500)), lon, lat); zw500 = atlev(zw, lev, psw, hyam, hybm, P0); zc500 = atlev(zc, lev, psc, hyam, hybm, P0); cd /home/disk/tao/dvimont/matlab/CCM/CT/Vars save 500mb_vars.mat uw500 vw500 zw500 uc500 vc500 zc500 lev = 250; uw250 = atlev(uw, lev, psc, hyam, hybm, P0); vw250 = atlev(vw, lev, psc, hyam, hybm, P0); vorw250 = sph_curl(squeeze(mean(uw250)), squeeze(mean(vw250)), lon, lat); uc250 = atlev(uc, lev, psc, hyam, hybm, P0); vc250 = atlev(vc, lev, psc, hyam, hybm, P0); vorc250 = sph_curl(squeeze(mean(uc250)), squeeze(mean(vc250)), lon, lat); zw250 = atlev(zw, lev, psw, hyam, hybm, P0); zc250 = atlev(zc, lev, psc, hyam, hybm, P0); save 250mb_vars.mat uw250 vw250 zw250 uc250 vc250 zc250 global XAX YAX FRAME XAX = lon; YAX = lat; FRAME = [0 360 -90 90]; top = vorw250; bot = vorc250; top = vorw500; bot = vorc500; top = divw250; bot = divc250; top = divw500; bot = divc500; figure(1) range1 = [-30:1:30]*1e-5; sp(1) mcont(top, range1); sp(2) mcont(bot, range1); figure(2) range2 = [-10:.5:10]*1e-5; mcont(top-bot, range2); divw500 = sph_div(squeeze(mean(uw500)),squeeze(mean(vw500)),lon,lat); divc500 = sph_div(squeeze(mean(uc500)),squeeze(mean(vc500)),lon,lat); divw250 = sph_div(squeeze(mean(uw250)),squeeze(mean(vw250)),lon,lat); divc250 = sph_div(squeeze(mean(uc250)),squeeze(mean(vc250)),lon,lat); [ntim, nlat, nlon] = size(zw500); zon = ones(nlon, 1) * reshape(squeeze(mean2(shiftdim(zw500, 2))), 1, ntim*nlat); zon = reshape(shiftdim(zon, 1), ntim, nlat, nlon); zw5star = zw500 - zon; zon = ones(nlon, 1) * reshape(squeeze(mean2(shiftdim(zc500, 2))), 1, ntim*nlat); zon = reshape(shiftdim(zon, 1), ntim, nlat, nlon); zc5star = zc500 - zon; zon = ones(nlon, 1) * reshape(squeeze(mean2(shiftdim(zw250, 2))), 1, ntim*nlat); zon = reshape(shiftdim(zon, 1), ntim, nlat, nlon); zw25star = zw250 - zon; zon = ones(nlon, 1) * reshape(squeeze(mean2(shiftdim(zc250, 2))), 1, ntim*nlat); zon = reshape(shiftdim(zon, 1), ntim, nlat, nlon); zc25star = zc250 - zon; top = squeeze(mean(zw5star)); bot = squeeze(mean(zc5star)); top = squeeze(mean(zw25star)); bot = squeeze(mean(zc25star)); figure(1) range = [-500:50:500]; sp(1) mcont(top, range) sp(2) mcont(bot, range) figure(2) range = [-400:20:400]; mcont(top-bot, range)