Documentation of seasonal_gr

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

Help text

  Look at Surface Winds:

Cross-Reference Information

Listing of script seasonal_gr

cd /home/disk/tao/dvimont/matlab/CCM/GR
filin = 'prate.mon.mean.nc';
var_nam = 'prate';
[prategr, pratecl] = reg_gr(filin, var_nam, 1);
prategrt = prategr * (3600 * 24 * 1e3 / 1e3);
cd /home/disk/tao/data/nmc.reanalysis/monthly
filin = 'prate.mon.mean.nc';
nc = netcdf(filin, 'nowrite');
  nmclev = nc{'level'}(:);
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
nc = close(nc);
cd /home/disk/tao/dvimont/matlab/CCM/GR
get_global
define_global
FRAME = [0 360 -45 45];
figure(1)
sp(1)
     mcont(prategrt, [-5:.5:5], 'gortho', [0 180]);
     yticks = 0; xticks = 0;
     manual_ticks('x',xticks,['']); manual_ticks('y',yticks,['EQ']);
     title('NMC:  SON Precipitation Regressed on GR')
     xlabel('Contour Interval:  0.5 mm (day)^-^1 (std)^-^1')
cd /home/disk/tao/dvimont/matlab/CCM/GR
filin = 'air.mon.mean.srfc.nc';
var_nam = 'air';
[sfctgr, sfctcl] = reg_gr(filin, var_nam, 1);
cd /home/disk/tao/data/nmc.reanalysis/monthly
nc = netcdf(filin, 'nowrite');
  nmclev = nc{'level'}(:);
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
nc = close(nc);
cd /home/disk/tao/data/nmc.reanalysis
filin = 'landmask.nc';
nc = netcdf(filin, 'nowrite');
  lm = nc{'lsmask'}(:);
nc = close(nc);
cd /home/disk/tao/dvimont/matlab/CCM/GR
define_global
top = zeros(size(sfctgr));
top(find(lm == 1)) = sfctgr(find(lm == 1));
sp(2)
     mcont(top, [-5:.1:5], 'giso', [0 180]);
     yticks = 0; xticks = 0;
     manual_ticks('x',xticks,['']); manual_ticks('y',yticks,['EQ']);
     title('NMC:  SON Temperature Regressed on GR')
     xlabel('Contour Interval:  0.1 K (std)^-^1')
cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots


cd /home/disk/tao/dvimont/matlab/CCM/GR

lev = 1000;
filin = 'uwnd.mon.mean.nc';
var_nam = 'uwnd';
[u1000gr, u1000cl] = reg_gr(filin, var_nam, lev);
filin = 'vwnd.mon.mean.nc';
var_nam = 'vwnd';
[v1000gr, v1000cl] = reg_gr(filin, var_nam, lev);

top1 = thin(u1000gr,2);
top2 = thin(v1000gr,2);
FRAME = [0 360 -50 50];

figure(1)
sp(1)
XAX = thin(lon,2);
YAX = thin(lat,2);
     mquiv(top1, top2, 1.5, 'gortho', [0 180])
     yticks = 0; xticks = 0;
     manual_ticks('x',xticks,['']); manual_ticks('y',yticks,['EQ']);
     title('JJA 1000mb Winds Regressed on GR')