Documentation of make_som_dom_plots

Global Index (short | long) | Local contents | Local Index (short | long)

Help text

  CONVERGENCE

Cross-Reference Information

Listing of script make_som_dom_plots

clear
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
filin1 = 'wgr_5093.nc';
filin2 = 'cgr_5093.nc';
filin3 = 'wgr_SOM_hard.nc';
filin4 = 'cgr_SOM_hard.nc';
var1 = 'TAUX'; var2 = 'TAUY';
[dw1, dw2] = getnc(filin1, var1, var2);
[dc1, dc2] = getnc(filin2, var1, var2);
[sw1, sw2] = getnc(filin3, var1, var2);
[sc1, sc2] = getnc(filin4, var1, var2);
[lat, lon] = getll(filin1);
get_global; default_global;
tim = [1 2 3 4 11 12];

u = squeeze(mean(dw1(tim,:,:)));
v = squeeze(mean(dw2(tim,:,:)));
dw = sph_curl(u, v, lon, lat);
u = squeeze(mean(dc1(tim,:,:)));
v = squeeze(mean(dc2(tim,:,:)));
dc = sph_curl(u, v, lon, lat);
u = squeeze(mean(sw1(tim,:,:)));
v = squeeze(mean(sw2(tim,:,:)));
sw = sph_curl(u, v, lon, lat);
u = squeeze(mean(sc1(tim,:,:)));
v = squeeze(mean(sc2(tim,:,:)));
sc = sph_curl(u, v, lon, lat);

lab = 'SFC WIND STRESS CONVERGENCE';
cint = 1; clevs = [-20:cint:20];
FRAME = [0 360 -60 60];

i = 1;

if i == 1;
  top = (dw - dc)*1e8;
  bot = (sw - sc)*1e8;
elseif i == 2;
  top = (sw - dw)*1e8;
  bot = (sc - dc)*1e8;
end

sp(1)
     mcont(top, clevs, 'giso', [0 180]);
     if i == 1;
       title(['NDJFMA ' lab ': (WGR\_DOM - CGR\_DOM)']);
     elseif i == 2;
       title(['NDJFMA ' lab ': (WGR\_SOM - WGR\_DOM)']);
     end
     xlabel(['Contour Interval:  ' num2str(cint) 'e^-^8 s^-^1'])
sp(2)
     mcont(bot, clevs, 'giso', [0 180]);
     if i == 1;
       title(['NDJFMA ' lab ': (WGR\_SOM - CGR\_SOM)']);
     elseif i == 2;
       title(['NDJFMA ' lab ': (CGR\_SOM - CGR\_DOM)']);
     end
     xlabel(['Contour Interval:  ' num2str(cint) 'e^-^8 s^-^1'])