Global Index (short | long) | Local contents | Local Index (short | long)
This script calls | |
---|---|
clear cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/ct_t31_ense/data cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data [zw, pw] = getnc('wgr_ice.nc', 'Z3', 'PS'); [zc, pc] = getnc('cgr_ice.nc', 'Z3', 'PS'); [zw, pw] = getnc('wgr_SOM_hard.nc', 'Z3', 'PS'); [zc, pc] = getnc('cgr_SOM_hard.nc', 'Z3', 'PS'); [zw, pw] = getnc('wct_ense.nc', 'Z3', 'PS'); [zc, pc] = getnc('cct_ense.nc', 'Z3', 'PS'); [lat, lon, lev, hyam, hybm, P0] = getll('cct_ense.nc'); get_global default_global if size(zw,1) == 12; zw = zw([1 2 3 4 11 12], :, :, :); pw = pw([1 2 3 4 11 12], :, :); zc = zc([1 2 3 4 11 12], :, :, :); pc = pc([1 2 3 4 11 12], :, :); end lev2 = 500; zw2c = atlev(zw(:,5:18,:,:), lev2, pw, hyam(5:18), hybm(5:18), P0); zc2c = atlev(zc(:,5:18,:,:), lev2, pc, hyam(5:18), hybm(5:18), P0); top = squeeze(mean(zw2c-zc2c))/2; figure(1); sp(1) if ismap(gca); clma; end; mcont(top, [-100:20:100]); title('CT\_ENSEMBLE Run: DJF 500mb Height (WARM - COLD)/2'); xlabel('Contour Interval 20m'); cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data zw = getnc('wgr_SOM_hard.nc', 11, 'Z3'); zc = getnc('cgr_SOM_hard.nc', 11, 'Z3'); zw = getnc('wgr_ice.nc', 11, 'Z3'); zc = getnc('cgr_ice.nc', 11, 'Z3'); zw = squeeze(zw); zc = squeeze(zc); top = squeeze(mean(zw2-zc2))/2; sp(2) if ismap(gca); clma; end; mcont(top, [-100:5:100]); title('GR\_SOM\_Hard Run: 500mb Height (WARM - COLD)/2'); xlabel('Contour Interval 5m'); cd /home/disk/tao/dvimont/pub/david FRAME = [-.1 360 -90 90]; [xk, yk] = keep_var(FRAME, lon, lat); lat2 = lat(yk); lon2 = lon(xk); XAX = lon2; YAX = lat2; var2 = top(yk, xk);