Global Index (short | long) | Local contents | Local Index (short | long)
This script calls | |
---|---|
clear cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/som_5093/clim_5093/data filin = ['EP_FLUX_SOM_clim.nc']; [uv, u, v, airt, ps] = getnc(filin, ... 'UV', 'U', 'V', 'AIRT', 'PS'); [lat, lon, lev, hyam, hybm, p0] = getll(filin); [ntim, nlev, nlat, nlon] = size(u); win = sort([2:12:120 3:12:120 4:12:120]); ps = reshape(ps, ntim, nlat*nlon); airt = reshape(airt, ntim, nlev, nlat*nlon); vt = reshape(airt, ntim, nlev, nlat*nlon); press = p0 * hyam * ones(120, nlat*nlon) + hybm * ps; ... Unfinished. This is a very memory - intensive job. It should probably be done only one season at a time to reduce memory and time.