Documentation of plot_nmc

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

Help text

     title('PRECIPITATION:  NMC NDJFMA regressed on GR index');
	title('SLP:  NMC NDJFMA regressed on GR index');

Cross-Reference Information

Listing of script plot_nmc

clear
cd /home/disk/tao/data/nmc.reanalysis/monthly
filin = 'hgt.mon.mean.nc';
[lat, lon, level] = getnc(filin, 'lat', 'lon', 'level');
get_global
default_global
cd /home/disk/tao/dvimont/matlab/CCM/GR
filin = 'hgt.mon.mean.nc';
var_nam = 'hgt';  
cint = 5; clev = [-50:cint:50];
gstyle = 'stereo';  orig = [90 270];
[pat250, clim] = reg_gr(filin, var_nam, 250);
[pat500, clim] = reg_gr(filin, var_nam, 500);
filin = 'prate.mon.mean.nc';
var_nam = 'prate';
gstyle = 'giso';  orig = [0 180]; FRAME = [0 360 -60 60];
[pat250, clim] = reg_gr(filin, var_nam, 1);
pat250 = pat250 * 3600 * 24;
cint = 0.5; clev = [-5.25:cint:5.25];
filin = 'sst.mon.mean.nc';
var_nam = 'air';
[pat500, clim] = reg_gr(filin, var_nam, 1);
filin = 'uwnd.mon.mean.nc';
[u, clim] = reg_gr(filin, 'uwnd', 1000);
filin = 'vwnd.mon.mean.nc';
[v, clim] = reg_gr(filin, 'vwnd', 1000);
filin = 'uwnd.mon.mean.nc';
[pat500, clim] = reg_gr(filin, 'uwnd', 500);
cint = 0.5; clev = [-10:cint:10];
filin = 'slp.mon.mean.nc';
var_nam = 'slp';
[pat250, clim] = reg_gr(filin, var_nam, 1);
gstyle = 'giso';  orig = [0 180]; FRAME = [0 360 -70 70];
cint = 0.5; clev = [-10:cint:10];
cd /home/disk/tao/data/nmc.reanalysis/monthly
[lat, lon] = getll(filin);
default_global
figure(1)
sp(1)
     if ismap; clma; end;
     mcont(pat250, clev, gstyle, orig)
     title('250 mb HEIGHT:  NMC NDJFMA regressed on GR index');
     xlabel(['Contour Interval:  ' num2str(cint) ' m std^-^1']);
sp(2)
     clma
     mcont(pat500, clev, gstyle, orig)
     title('500 mb HEIGHT:  NMC NDJFMA regressed on GR index');
%     title('500 mb ZONAL WIND:  NMC NDJFMA regressed on GR index');
%	title('SURFACE TEMP:  NMC Winter regressed on GR index');
     xlabel(['Contour Interval:  ' num2str(cint) ' m std^-^1']);

cd /home/disk/tao/dvimont/matlab/NMC/Plots2
%print -dps2 250mb_NMC_nwin.ps

gstyle = 'giso';  orig = [0 180]; FRAME = [0 360 -40 40];
XAX = thin(lon, 2);
YAX = flipud(thin(lat, 2));
udat = flipud(thin(u, 2)); vdat = flipud(thin(v, 2));

sp(2)
	clma
	[h, a, c] = mquiv(udat, vdat, 1, 'giso', [0 180]);
   	title('1000 mb WIND:  NMC NDJFMA regressed on GR index');
        xlabel(['Max Wind:  ' num2str(c(2)) ' m s^-^1'])

cd /home/disk/tao/dvimont/matlab/NMC/Plots