Global Index (short | long) | Local contents | Local Index (short | long)
This script calls | |
---|---|
clear cd /cdrom var1 = 'UU'; var2 = 'VV'; var3 = 'AIRT'; var4 = 'V'; [uuw1, vvw1, uw1, vw1] = getnc('warm/wgr_5093.nc', ... var1, var2, var3, var4); [uuc1, vvc1, uc1, vc1] = getnc('cold/cgr_5093.nc', ... var1, var2, var3, var4); [uuw2, vvw2, uw2, vw2] = getnc('warm/wgr_SOM_hard.nc', ... var1, var2, var3, var4); [uuc2, vvc2, uc2, vc2] = getnc('cold/cgr_SOM_hard.nc', ... var1, var2, var3, var4); [lat, lon, lev] = getll('cold/cgr_SOM_hard.nc'); get_global; default_global; varw1 = uuw1 - uw1.*vw1; varc1 = uuc1 - uc1.*vc1; varw2 = uuw2 - uw2.*vw2; varc2 = uuc2 - uc2.*vc2; tim = [1 2 3 4 11 12]; figure(1) top = varw1 - varc1; bot = varw2 - varc2; top = vvw1-vvc1; bot = vvw2-vvc2; top = squeeze(mean(top(tim,:,:)))*(1000*3600*24); bot = squeeze(mean(bot(tim,:,:)))*(1000*3600*24); sp(1) mcont(top,[-1:.1:1], 'giso', [0 180]); sp(2) mcont(bot,[-1:.1:.1], 'giso', [0 180]);