| Name | Description |
|---|---|
A | |
| acosweight | x=cosweight(sst, nya, oro); |
| advect_backwards | u,v ve u,v |
| airt_compare | |
| airt_gr_compare | Load GR pattern and time series |
| annave | [x, clim] = annave(y); |
| annular_mode | cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/HTRAE/EARTH/atm/netcdf |
| AO_regression_maps | %%%%%%%%%%%%%%%%%%%%%%%%% Color plot of AO %%%%%%%%%%%%%%% |
| ar | AR Computes AR-models of signals using various approaches. |
| armax | ARMAX Computes the prediction error estimate of an ARMAX model. |
| arx | ARX Computes LS-estimates of ARX-models. |
| arx2th | ARX2TH Constructs the THETA-format for an ARX model. |
| arxstruc | ARXSTRUC Computes loss functions for families of ARX-models. |
| atlev | out1 = atlev(dat, lev, ps, hyam, hybm, P0); |
| auxvar | AUXVAR Auxiliary variables for the parameter estimation schemes. |
| axis_limits | axis_limits(num) |
B | |
| back | |
| backwards_quicklook | Plot data |
| Barsugli_Battisti1 | Reproduce figure 4 from BB98 |
| battisti_spectrum | tfield1 = tfield1_10000; |
| battisti_spectrum2 | tfield1 = tfield1_10000; |
| bj | BJ Computes the prediction error estimate of a Box-Jenkins model. |
| blt | BLT Returns FALSE if Lines and Text objects in Figure should be printed in black. |
| bodeplot | BODEPLOT Plots the Bode diagram of a transfer function or spectrum. |
C | |
| ct_hgt_nmc | ctstar = ctstar(1:396); |
| ceof_regress_on_all_2 | Do complex eof: |
| compare_cts | tim = [1200:12000]; |
| compare_cts3 | ctlim = [-.01 360 -90 90]; |
| calc_eq_tend_terms | Get lower layer |
| cartoon2 | tx1 = regress_eof(txann, lpcs, 0); |
| color_shadem | color_shade(x, lev, color) |
| color_davet | cmap = (0.5+cmap)/1.5; |
| colorbar2 | COLORBAR Display color bar (color scale). |
| clo | CLO Clear object |
| calc_clim | Check the curl(tau) method |
| calc_eq_tend_terms | Get lower layer |
| calc_gr_once_again | Load landmask |
| calc_new_gr | Load and configure NMC.REANAL data |
| calc_pacbasin_gr | Load and configure NMC.REANAL data |
| calc_seasonal_gr_pat | Load landmask |
| cancel | Cancel Method to interupt a print job and finish it right now. |
| canform | CANFORM Constructs canonical form model structures. |
| canstart | CANSTART State-space model with initial parameter estimates. |
| cartoon | h = plot3(par2, lat, -100*ones(size(lat)), 'r'); |
| cartoon2 | tx1 = regress_eof(txann, lpcs, 0); |
| catpreview | CATPREVIEW Creates an EPS/TIFF preview file. |
| cdir | |
| ceof_reg_oce_current | Load CPCs |
| ceof_regress_on_all_2 | Get 20deg data from saved data |
| ceof_regress_velocities | Plot regressions on CPC1 of T20' |
| ceof_regress_velocity | Get PCs |
| check_atl_trop_vwnd | Get ATL CTI |
| check_nh_sh_eq_corr | Notes: |
| check_old_new_stats | c = |
| check_pac_trop_vwnd | Define regions, etc |
| clabel | CLABEL Contour plot elevation labels. |
| clean | keep('XAX', 'YAX', 'FRAME', 'CURDIR'); |
| clf | CLF Clear current figure. |
| clim_nmc | Note that December corresponds to month 1, and a new year. |
| clim_plot | ind1 = [1 2 12]; % Winter |
| climatology | sst = sst(360:1800,:,:); |
| climave | [x, clim] = climave(y); |
| clo | CLO Clear object |
| coads_sst_calc_GR | cd /home/disk/tao/data/coads/coads1 |
| color_4levs | |
| color_davet | |
| color_int | |
| color_jeff | [h, c] = color_shade(x, lev, color) |
| color_publish | |
| color_shade | [h, c] = color_shade(x, lev, color) |
| color_shade_dave | [h, c] = color_shade(x, y, dat, lev, color) |
| color_shadem | color_shade(x, lev) |
| color_shadem2 | color_shade(x, lev, color) |
| color_trans | |
| color_under | h = color_under(lat, lon, data, levels, pm) |
| colorbar2 | COLORBAR Display color bar (color scale). |
| colorbar3 | cb = colorbar3(loc, 'force'); |
| colorbar4 | cb = colorbar4(loc, children, 'force'); |
| cometnew | COMET Comet-like trajectory. |
| comp_eof_sst_regress_on_all | ctlim = [115 285 -30 30]; |
| compare | COMPARE Compares the simulated/predicted output with the measured output. |
| compare_coup_stoc | |
| compare_cts | rho = (corr(pcs3(:,1), pcs3(:,1), 1));% + sqrt(corr(ctann, ctann, 2))) / 2 ;% + ... |
| compare_cts2 | temp = nc{'temp'}(tim, yk, xk); |
| compare_cts3 | |
| compare_data | compare(var1, var2, lev, nlev); |
| compare_eof_regions | |
| compare_fields | Start off with the typical variables. |
| compare_prec | |
| compare_qflux | Now check out the surface heat flux budgets. |
| compare_som_fixed | SLP: |
| compare_various_grs | fid = fopen('reg_sst.bin', 'r', 'b'); |
| complex_eof | Complex (Hilbert) EOF: |
| complex_eof_heat_content | Area average heat so process is quicker |
| composite | dc2(lm); |
| composite_1 | Load various maps (warm and cold); |
| cont1d_stereo | contorm(YAX,XAX2,aray2); |
| contin | CONTIN Converts a model to continuous-time form. USE RATHER THD2THC! |
| contourf2 | CONTOURF Filled contour plot. |
| contourfm2 | CONTOURFM Filled contour map. |
| contth | CONTTH Is an auxiliary function for eta2ss and thss2th |
| convert_sigma_to_pres | out = convert_sigma_to_pres(in, ps, newlevs); |
| convert_sigma_to_pres2 | out = convert_sigma_to_pres(in, ps, newlevs); |
| corr | r=corr(vect1, vect2, lag); |
| corr2 | r=corr(vect1, data, lag); |
| corr_ct_mlcurl | tim = 551:1000; |
| corr_indices | tim = 101:550; |
| corr_LP_RAW | |
| corr_nan | function c = covar_nan ( mat1 , mat2 ) ; |
| corr_nan2 | function c = covar_nan ( mat1 , mat2 ) ; |
| corr_nan3 | function c = covar_nan ( mat1 , mat2 , tol , show ) ; |
| corr_sig | [ccoef, sig_pat] = corr_sig(x, y, sig_lev); |
| cosp | function prod=cosp(x,y,nya,minperc) |
| cospectrum | Note: a positive phase implies var1 leads var2 |
| cosweight | x=cosweight(sst, nya); |
| coup_run_stats | Average by season, to make this a bit more manageable |
| coup_windstress | use c = 60m/s, or rossby radius ~10degrees ==> 2c=12000. |
| cov5 | For the 2x2 matrix, only returns a single value... |
| covar_gr | |
| covar_nan | function c = covar_nan ( mat1 , mat2 ) ; |
| covar_nan2 | function c = covar_nan ( mat1 , mat2 ) ; |
| covar_nan3 | function c = covar_nan ( mat1 , mat2 , tol, show ) ; |
| covf | COVF Computes the covariance function estimate for a data matrix. |
| covf2 | COVF2 Computes covariance function estimates |
| cra | CRA Performs correlation analysis to estimate impulse response. |
| create_tenlevel_temperature | Open netcdf file. First execute the command (on hayes): |
| create_tenlevel_uocean | |
| cs | CASE STUDIES IN SYSTEM IDENTIFICATION |
| cs1 | CASE STUDY NUMBER 1 : A GLASS TUBE DRAWING FACTORY |
| cs2 | CASE STUDY # 2 ANALYSING A SIGNAL FROM A TRANSFORMER |
| csirod | |
| csirodi | |
| ct_coads_hgt_nmc | |
| CT_contour_plots | This file contains commands that will generate the figures in |
| ct_gr_z500 | get CTstar: |
| ct_hgt_nmc | |
| ct_lagged_correlation | |
| ct_lpct_regs | Load Data |
| CT_plots | Start with the CT index |
| ct_seasonal | compress data |
| ct_spectrum | Look at power spectra of successive 1000 year chunks |
| ct_spectrum2 | Plot the following: |
| ct_sst | Get climatology for October through March |
| ct_statistics | Monthly |
| ct_uzon_nmc | Note: In producing plots on 20 OCT, 1998, I accidently |
| ct_warmcold | Just a side note here. Recall that the SSTA pattern |
| curlt_eof | Do complex eof: |
D | |
| dcm | function hh = drawcoasts ( lcst, lgrd, lbox, gridx, gridy ) ; |
| define_globals | set graphics defaults |
| define_globals | Define global variables that are commonly used. This |
| drawmapm | drawmapm |
| data | |
| dc | function hh = drawcoasts ( lcst, lgrd, lbox, gridx, gridy ) ; |
| dc2 | dc2(data, colr, pn); |
| dc3 | function hh = drawcoasts ( elev, colr ) ; |
| dc_color | dc2(data); |
| dc_nogrid | function hh = drawcoasts ( lcst, lgrd, lbox, gridx, gridy ) ; |
| dcm | |
| dcm2 | function hh = drawcoasts ( lcst, lgrd, lbox, gridx, gridy ) ; |
| dcmfill | hh = dcm2 ( lev , colr ) ; |
| dcv | dcv(data); |
| default_global | xk = 1:length(lon); |
| default_global1 | xk = 1:length(lon); |
| default_global2 | xk = 1:length(lon); |
| define_1950-93_clim | |
| define_global | |
| define_globals | Define global variables that are commonly used. This |
| define_sst | |
| define_SST_clim | This program contains code to look at SST differences between |
| define_TPAC_yearly_ct | Load landmask |
| define_trop_CT | Edit time: |
| define_trop_GR | [wsst] = getnc('gr_warm_sst.nc', 'SST'); |
| define_weight_netcdf | Define weight |
| delayed_oscillator | ctlim = [115 285 -30 30]; |
| delnew | function [delta, deltaccm] = delnew(ecc, epsi, per) |
| demos | DEMOS Returns demo information to the MATLAB Demo. |
| detrend | DETREND Remove a linear trend from a vector, usually for FFT processing. |
| detrend_NaN | DETREND Remove a linear trend from a vector, usually for FFT processing, |
| detrend_NaN2 | function Y = detrend_NaN ( Xdat , tol , show ) ; |
| detrend_NaN_slow | out = detrend_NaN(in); |
| dg | dg(lat, lon); |
| dg2 | dg(lat, lon); |
| disabledberror | DISABLEDBERROR turns off "dbstop if error" and returns previous |
| dof_corr | dof = dof_corr(y, x); |
| draw_landmap | drawmapm |
| draw_landmap2 | drawmapm |
| drawbox | h = drawbox(lims, sty); |
| drawboxm | h = drawbox(lims, sty); |
| drawmap2 | dc2(data); |
| drawmapm | drawmapm |
| dtemp_dt | Now get tprime |
| dtrend | DTREND Removes trends from data sets. |
E | |
| eq_tend_terms | Get lower layer |
| eq_tend_terms2 | Load PCS: |
| eq_tend_terms | Plot the data; |
| eq_tend_terms2 | Get windstress if necessary |
| eq_tend_terms | Plot all the terms |
| eq_tend_terms2 | Get windstress if necessary |
| eq_tend_terms | Plot the data; |
| eddy_time_plots | save eddy_interped_grsom.mat vtw vtc kew kec newlev lat lon |
| eddy_time_plots2 | Read Variables from wgr and cgr runs |
| eigenvalues | |
| eigplot | eigplot(s1, s2, s3, er1, er2); |
| enabledberror | ENABLEDBERROR turns on "dbstop if error" if hasDBError is true. |
| enso_composites | Define ENSO years somehow here. |
| eof2 | function [pc,eof_reg,sing]=eof2(input_array,n,nya) |
| eof_20deg_depth | load depth_20deg_ave_yr101-550.mat |
| eof_20deg_isotherm | tim = [551:1000]; |
| eof_5093 | cd /home/disk/hayes2/werner/matlab |
| eof_dan | [lam, lds, pcs, per] = eof_dan(data); |
| eof_dan2 | [lam, lds, per] = eof_dan(c, isrot, nkp); |
| eof_dan_NaN | [lam, lds, pcs, per] = eof_dan(data); |
| eof_dan_noweight | [lam, lds, pcs, per] = eof_dan(data); |
| eof_heat_content | Filter Heat |
| eof_height | nc = netcdf('sal_gauss.cdf', 'nowrite'); |
| eof_nan | function [pc,eof_reg,sing]=eof2(input_array,n,nya) |
| eof_routine | Normalized eof routine: |
| eof_slp | |
| eof_sst | lims = [-0.1 360 -90 90]; |
| eof_sst_slp | First, non-standardized |
| epflux_clim | |
| epflux_t31 | |
| eq_depth_terms | |
| eq_ind_corr | title(['Lagged correlations between 10yr LP filtered,'... |
| eq_tend_terms | Load PCS: |
| eq_tend_terms2 | Load PCS: |
| eq_tend_terms3 | Get windstress if necessary |
| eta2ss | ETA2SS Auxiliary routine to TH2SS. |
| etfe | ETFE Computes the Empirical Transfer Function Estimate and Periodogram. |
F | |
| figure1 | ct = getheat(lims, 1:7, tim, 'temp'); |
| figure_landscape | figure_landscape(fig_num); |
| figure_tall | figure_tall(fig_num); |
| fill_landmap | hh = fill_landmap ( lev , colr ) ; |
| fact | fac = fact(num); |
| fcst_2_CCM | fcst_2_CCM.m |
| ffplot | BODEPLOT Plots the Bode diagram of a transfer function or spectrum. |
| ffsdcal | FFSDCAL Auxiliary function to TH2FF. |
| fiactha | FIACTHA Finds indices of handles HH that are active uicontrols ('Value'==1). |
| fiacthad | FIACTHAD Finds the selected data sets. |
| fiactham | FIACTHA Finds the selected models. |
| fig1 | hc = getheat(lims, 1:7, tim); |
| Fig_spect | ctann = rpcs(1:450,num); |
| figure1 | ct = getheat(lims, 1:7, tim, 'temp'); |
| figure_landscape | |
| figure_orient | |
| figure_tall | |
| fill_landmap | hh = fill_landmap ( lev , colr ) ; |
| fill_landmap2 | hh = fill_landmap ( lev , colr ) ; |
| fill_oceanmap | hh = fill_oceanmap(lev, colr); |
| findmax | maxind = findmax(x); |
| findmin | minind = findmin(x); |
| findX | function X0=findX(F,F0,X); |
| findzeros_of_dot | tau = 0.5; nt = length(tau); |
| finish | Finish Method to end a PrintJob at its current page. |
| fixpar | FIXPAR Fixes parameters in state-space and ARX model structures. |
| fixxx | FIXXX Subroutine for fixpar |
| fl | |
| flux_gr_ice | [swl1, lwl1, lhl1, shl1] = getnc('clim_5093.nc', v1, v2, v3, v4); |
| flux_sst_corr | Look at the point by point correlation maps between SST and different |
| fo | |
| form_ticklabels | function tick_labels = form_ticklabels ( tick_labels, nd ); |
| fpdo | |
| frame | function frame ( arg1, arg2 ) ; |
| framem | FRAMEM Toggle and control the display of the map frame |
| framem2 | FRAMEM Toggle and control the display of the map frame |
| framem_orig | FRAMEM Toggle and control the display of the map frame |
| freqfunc | FREQFUNC Frequency functions are created by SPA, ETFE, and TH2FF. |
| fstab | FSTAB FSTAB(A) stabilizes a MONIC polynomial with respect to the |
G | |
| get_ubar_gradtprime | [ub_dhpdx, vb_dhpdy, lat_out, lon_out, depth_out] = ... |
| get_uprime_gradtbar | [up_dhbardx, vp_dhbardy, lat_out, lon_out, depth_out] = ... |
| get_wbar_dtprimedz | [wdhdz, lat_out, lon_out, depth_out] = ... |
| get_wprime_dtbardz | [wdhdz, lat_out, lon_out, depth_out] = ... |
| getheat | [heat, lat, lon, depth, middepth] = getheat(lev, tim, ctlim); |
| getheat | [heat, lat, lon, depth, middepth] = getheat(lims, lev, tim, varn); |
| getheat | [heat, lat, lon, depth, middepth] = getheat(lev, tim, ctlim); |
| get_hctend_terms | Get ubar_dHCprimedx and vbar_dHCprimedy |
| get_hp_hctend_terms | Get ubar_dHCprimedx and vbar_dHCprimedy |
| get_hctend_terms | lims = [89 310 -15 15]; |
| get_hctend_terms | Get ubar_dHCprimedx and vbar_dHCprimedy |
| gu_philander | save gu_philand.mat vbtp wbtp vptb wptb ubtp uptb lims tim lags |
| getheat | [heat, lat, lon, depth, middepth] = getheat(lims, lev, tim, varn); |
| global_axes | global_axes(hsz, vsz, midh, midv, tmarg): |
| global_latlon | global_latlon(lat, lon, frame); |
| get_varn | theVar = get_varn(nc); |
| getnc2 | [data, lat, lon, level, time] = getnc2(filin, varn, lims, lev, tim); |
| get_global | |
| global_axes | global_figs(hsz, vsz, midh, midv, tmarg): |
| global_xy | global_xy(lat, lon); |
| gcont | function [h, c]=gcont(var, clev) |
| get_20deg_isodepth | tim = [551:1000]; |
| get_20deg_vel | First, load temperature field to get climatology |
| get_annual_timeseries | Get landmask |
| get_ao_regs | Perform EOF |
| get_cdtem | |
| get_contour_int | [clev, cint] = get_contour_int(dat, nlev); |
| get_coord | theCoord = get_coord(nc); |
| get_ct_regs | Remove annual cycle |
| get_cti | |
| get_dhcdt | [dhcdt, lon_out, lat_out, tim] = get_dhcdt(pcs, lims, nfrm, tim, lev); |
| get_dof | dofx = dof(data); |
| get_dutdx | [dubartdx, dtbarudx, lon_out, lat_out] = get_dutdx(nfrm, pcs, lims, lev, tim); |
| get_dvtdy | [dvbartdy, dtbarvdy, lon_out, lat_out] = get_dvtdy(nfrm, pcs, lims, lev, tim); |
| get_flux_terms | evp : evaporation (mm/day) |
| get_global | |
| get_hctend_terms | Get ubar_dHCprimedx and vbar_dHCprimedy |
| get_hp_flux_terms | evp : evaporation (mm/day) |
| get_hp_hctend_terms | Get ubar_dHCprimedx and vbar_dHCprimedy |
| get_MIX_heat_flux | |
| get_MIX_heat_flux_mon | |
| get_month1 | outcell = get_month_name1(monind); |
| get_month2 | outcell = get_month_name1(monind); |
| get_monthly_data | cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/HTRAE/HTRAE/atm/netcdf |
| get_monthly_slp | Load SLP data from disk |
| get_NAO | Save data |
| get_nclatlon | [lat, lon, yk, xk] = get_nclatlon(lims, nc); |
| get_nclatlon2 | [lat, lon, yk, xk] = get_nclatlon(lims, nc); |
| get_nmc | [out1, lat, lon] = getnmc(varn1, ... , xylim, lev, tim); |
| get_recvar | theVar = get_recvar(nc); |
| get_taux_std | Start with FSU |
| get_time | tind = get_time(yr1, yr2, startyr); |
| get_ubar_gradtprime | [ub_dhpdx, vb_dhpdy, lat_out, lon_out, depth_out] = ... |
| get_ubtp | [ubtp, vbtp, lat_out, lon_out, depth_out] = ... |
| get_udtdx | [ubtp, uptb, lat_out, lon_out, depth_out] = ... |
| get_uprime_gradtbar | [up_dhbardx, vp_dhbardy, lat_out, lon_out, depth_out] = ... |
| get_uptb | [uptb, vptb, lat_out, lon_out, depth_out] = ... |
| get_varn | theVar = get_varn(nc); |
| get_vdtdy | [vbtp, vptb, lat_out, lon_out, depth_out] = ... |
| get_w_from_divu | Load u and v data; |
| get_wbar_dtprimedz | [wdhdz, lat_out, lon_out, depth_out] = ... |
| get_wbtp | [wbtp, lat_out, lon_out, depth_out] = ... |
| get_wdtdz | [wbtp, wptb, lat_out, lon_out, depth_out] = ... |
| get_windstress | Currently, this is the model wind stress: |
| get_wprime_dtbardz | [wdhdz, lat_out, lon_out, depth_out] = ... |
| get_wptb | [wptb, lat_out, lon_out, depth_out] = ... |
| get_zonal_wind | |
| getargth | GETARGTH Gets the auxiliary argument for a state-space model |
| getatt | varargout = getatt(filin, varargin); |
| getct | |
| getctm | |
| getff | GETFF Selects the frequency function (for own plotting). |
| getflx | [out1, ...] = getc(var1, ... , xylim, tim); |
| getget | GETGET Do a GET on Handle Graphics object or a GET_PARAM on Simulink object. |
| getheat | [heat, lat, lon, depth, middepth] = getheat(lims, lev, tim, varn); |
| getlev | |
| getll | [lat, lon, depth, lm] = getllc(varname, ctlim); |
| getll_ccm | [lat, lon, lev, hyam, hybm, P0] = getll(filin); |
| getmfth | GETMFTH Gets the name of the .m-file that defines the model structure |
| getml | [out1, ...] = getc(var1, ... , xylim, tim); |
| getnc | [out1, ...] = getc(var1, ... , xylim, lev, tim); |
| getnc2 | [data, lat, lon, level, time] = getnc2(filin, varn, lims, lev, tim); |
| getnc_ccm | [varn_out] = getnc(filin, varn, lims, lev, tim); |
| getncap | GETNCAP Gets the number of data points and the number of parameters for a model. |
| getpp | GETPP Get all PaperProperties. |
| gett | GETT Gets the sampling interval for a model. |
| getzp | GETZP Extracts the zeros and the poles from the ZEPO-format created by TH2ZP.% [Z,P] = GETZP(ZEPO,KU,KY) |
| ghostscript | GHOSTSCRIPT Function to convert a PostScript file to another format. |
| global_axes | global_figs(hsz, vsz, midh, midv, tmarg): |
| global_figs | global_figs(hsz, vsz, midh, midv, tmarg): |
| global_latlon | global_latlon(lat, lon, frame, wrap_grenw); |
| global_xy | global_xy(lat, lon); |
| globll | |
| gn | GN Computes the Gauss-Newton direction for PE-criteria |
| gnns | GNNS Computes the Gauss-Newton direction for PE-criteria. |
| go2d | [out, dims] = go2d(in); |
| gquiv | function [h, hx] = gquiv(var1, var2, clev, nthin, units) |
| gquiv2 | function [h, hx] = gquiv(var1, var2, clev, nthin, units) |
| gr_ct_tpac_compare | mquiv(top1, top2, 1, 'stereo', [90 270]); |
| gr_ice_noice | Plot out patterns regressed on LP filtered GR time series: |
| gr_nmc | Load and configure NMC.REANAL data |
| GR_no_ENSO | |
| greyshade | CONTOURF Filled contour plot. |
| greyshd | greyshd(x, lev) |
| gridm5 | |
| gshade | function h=gshade(var, clev) |
| gshade2 | greyshd(x, lev) |
| gu_philander | Look at heat content regressions |
| gu_philander2 | %%%%%%%%%%%%%%%%%%%% Try a Deser et. al. plot |
| gu_philander3 | Plot negative lags |
| gu_philander4 | save gu_philand.mat vbtp wbtp vptb wptb ubtp uptb lims tim lags |
| gu_philander5 | lag1 = -15:-10; lag2 = -9:-4; lp = length(lag1); |
| gu_philander_20degiso | Note: everything in kp2 is also in kp1. However, kp1 has two |
| gu_philander_HP | save gu_philand.mat vbtp wbtp vptb wptb ubtp uptb lims tim lags |
| gu_philander_thesis | save gu_philand2.mat vbtp wbtp vptb wptb ubtp uptb lims tim lags |
| gvect | function h=gcont(var1, var2, flim, clev) |
| gvect2 | function h=gcont(var1, var2, flim, clev) |
H | |
| heat_ceof | ctlim = [110 265 20 65]; % North Pacific Region |
| heat_ceof | tim = 551:1000; |
| hc_eof | lims = [110 300 -60 60]; |
| heat_ceof | tim = 101:1000; |
| heat_ceof | ctlim = [110 300 -30 30]; |
| heat_ceof | tim = 551:1000; |
| HC_clim | print -dps2 HC0-80_80-270.ps |
| hc_eof | save heat_pcs.mat rlam rlds rpcs rper lims lev lat lon ... |
| hc_eofs | |
| hc_sst_lagged_regression | cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data |
| hc_sst_lagged_regression_color | |
| heat_ceof | tim = 551:1000; |
| heat_content_clim | cwat = 4.218e3; % heat capacity of liquid water, J/(kg K) |
| Heatflux | print -dps2 HFLX.ps |
| hermite | y = hermite(x, n); |
| hgt_gr_compare | Load GR pattern and time series |
| hgt_nmc_ccm | Get climatology |
| hist_PC1_0-270HC | |
| hline | function h = hline ( y, linetype, x, zdat ) ; |
| hoffmoeller_corr_plot | Plot the lagged correlation coefficient for various depths |
| hoffmoeller_plot | cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data |
| hofmoeller_depth_temp | ctlim = [-.1 360 -6 6]; |
| hofmoeller_depth_temp2 | ctlim = [-.1 360 -6 6]; |
| hofmoeller_depth_temp3 | Load level 1 and level 5 monthly temperatures |
| hofmoeller_depth_temp4 | First, look depth sections across EQ |
| HP10_eq_tend_terms | Get windstress if necessary |
| HP10_LP9_wreg | |
| HP10_spatial_tendterms | prec = getflx('rnd', lims, tim); |
| HP_HC_evolution_30s30n | |
| hp_lp_SST_hoffmoeller | if i == 1; |
| hpgl | HPGL Method to add control characters to beginning of HPGL files. |
I | |
| isscalar | function iss = isscalar ( u ) ; |
| isvector | function isv = isvector ( u ) ; |
| Ice | print -dps2 ICO.ps |
| idarxstr | IDARXSTR Creates a matrix of arxorders based on a string definition |
| idbuildw | IDBUILDW This function handles the creation of all the ident VIEW windows. |
| idbwtext | IDBWTEXT This function builds the INFO TEXT dialog. |
| idcmdld | IDCMDLD Set up System Identification demos for MATLAB Demo. |
| idconfcp | IDCONFCP Asks if confidence intervals should be computed, and does so if desired. |
| iddatfig | IDDATFIG Generates all view curves for ident data objects. |
| iddemo | The SYSTEM IDENTIFICATION TOOLBOX is an analysis module |
| iddemo1 | This case study concerns data collected from a laboratory scale |
| iddemo2 | In this example we compare several different methods of |
| iddemo3 | This case study is based on the same "hairdyer" data as study # 1. |
| iddemo4 | iddemo4 |
| iddemo5 | IDDEMO5 |
| iddemo6 | The command SEGMENT segments data that are generated from |
| iddemo7 | In this demo we shall demonstrate how to use several the SITB to |
| iddmtab | IDDMTAB Creates additional data/model Summary Board tables. |
| ident | IDENT Starts the Graphical User Interface for the System Identification Toolbox. |
| idfilt | IDFILT Filters data. |
| idgenfig | IDGENFIG Generates all view curves for ident. |
| idgtws | SCRIPT FILE IDGTWS |
| idinput | IDINPUT Generates input signals for identification. |
| idinseva | IDINSEVA Inserts the working or validation data set. |
| idlayout | IDLAYOUT Sets certain layout-defaults for the Graphical User Interface |
| idlaytab | IDLAYTAB Handles closing of ident, and the layout values. |
| idmhit | IDMHIT Determines whether figure is over an object of a certain type. |
| idmodred | IDMODRED Reduces the order of models. |
| idmsize | IDMSIZE Sets default value for the variable maxsize (see help AUXVAR) |
| idmspop | IDMSPOP Returns the model structure popup string. |
| idmwwb | IDMWWB Handles the window button callback in the main ident window. |
| idnextw | IDNEXTW Finds the next available data or model axes. |
| idnonzer | IDNONZER Returns subset of input vector that are valid handles. |
| idopttog | IDOPTTOG Toggles checked options. |
| idparest | IDPAREST The basic callback function for parametric estimation. |
| idplot | IDPLOT Plots input - output data. |
| idresamp | IDRESAMP Resamples data by decimation and interpolation. |
| idsim | IDSIM Simulates a given dynamic system. |
| idsimsd | IDSIMSD Illustrates the uncertainty in simulated model responses. |
| idsimss | IDSIMSS simulates a given system |
| idstrip | IDSTRIP Converts a string of numbers, separated by spaces, to numbers. |
| iduiarx | IDUIARX Handles everything about ARX estimation in the ident GUI. |
| iduiaxes | IDUIAXES Handles the axes options. |
| iduiaxis | IDUIAXIS Sets up the axes limits dialog (axlimdlg) for ident plots. |
| iduibn | IDUIBN Handles model estimation for model that are defined By Name. |
| iduicalc | IDUICALC Unpacks the (vectorized) theta model and provides help to idgenfig. |
| iduiclpw | IDUICLPW Clears lines in the plot windows |
| iduiconf | IDUICONF Manages the confidence interval lines. |
| iduicra | IDUICRA Handles the correlation analysis dialog. |
| iduidemo | IDUIDEMO Runs the ident demo. |
| iduidrop | IDUIDROP Manages dropping one axes on another in ident window. |
| iduiedit | IDUIEDIT Handles various edit functions. |
| iduifile | IDUIFILE Handles file operations in ident. |
| iduifilt | IDUIFILT Handles prefiltering of data for ident. |
| iduifscr | SCRIPT FILE IDUIFSCR |
| iduigco | IDUIGCO Handle of current object. |
| iduigetd | IDUIGETD Gets the desired active data set |
| iduigetp | IDUIGETP Help function for managing XIDplow handles. |
| iduihelp | IDUIHELP Wrapper function to hthelp. |
| iduiinsd | IDUIINSD Inserts data into the Data Summary Board. |
| iduiinsm | IDUIINSM Handles the insertion of models into the Model Summary Board |
| iduiio | IDUIIO Handles estimation of parametric models in input-output form. |
| iduiiono | IDUIIONO Handles everything related to input-output numbers |
| iduiiter | IDUIITER Handles all iteration control for parameter estimation methods. |
| iduikeyp | IDUIKEYP The ident keypress function callback. |
| iduilay | IDUILAY Compute positions for nicely laid out uicontrols. |
| iduilay1 | IDUILAY1 Help function to compute positions for nicely laid out uicontrols. |
| iduilay2 | IDUILAY2 Computes window width based on the number of buttons. |
| iduimbcb | IDUIMBCB Handles the Mouse Button Callbacks for ident plot windows. |
| iduimod | IDUIMOD Main manager of the plot windows. |
| iduims | IDUIMS Enables the right popups in the Orders Editor dialog. |
| iduinpar | IDUINPAR Performs non-parametric estimation for ident. |
| iduiopt | IDUIOPT Handles the setting of all options. |
| iduipoin | IDUIPOIN Sets and resets the window pointer to watch and arrow |
| iduipop | IDUIPOP Manages the callbacks from ident's pop-up menus. |
| iduipw | IDUIPW Callback for all plot windows in the ident GUI. |
| iduiqs | IDUIQS Quickstart operation for ident. |
| iduisel | IDUISEL Manages the selection of input/output data variables and data portions. |
| iduisess | IDUISESS Handles load and save session as well as renaming of session. |
| iduispa | IDUISPA Handles the spectral analysis dialog. |
| iduiss | IDUISS Handles estimation of models in state-space form. |
| iduistat | IDUISTAT Manages the status line in main ident window. |
| iduital | IDUITAL Callback for 'Title and Labels' menu item. |
| iduivis | IDUIVIS Sets the visibility ONOFF for the valid handles in the list HANDLES. |
| iduiwast | IDUIWAST Shows the contents of the wastebasket. |
| iduiwok | IDUIWOK Checks if window number k in the ident GUI exists. |
| idunlink | IDUNLINK Performs the unlinking of figure with no wino. |
| idvmenus | IDVMENUS Sets the submenus for the different view windows. |
| india_sst | |
| ini_atrest | |
| inival | INIVAL Computes initial values for the prediction error method. |
| inputcheck | INPUTCHECK Method to validate input arguments to PRINT. |
| instructions | Instructions for making contour plots in MATLAB |
| interp_hgt_domsom | a = 'hi' |
| interp_oi2CCM | interp2CCM.m |
| intraseasonal_NAO | Average into 5x5 squares |
| iselement | function ise = iselement ( value, set ) ; |
| isfigure | ISFIGURE True for Figure handles. |
| ishghandle | ISHGHANDLE True for Handle Graphics object handles. |
| isinrange | function aa = isinrange ( x, xlim ) |
| isnum | function out1 = isnum( in1 ) |
| isscalar | function iss = isscalar ( u ) ; |
| isslhandle | ISSLHANDLE True for Simulink object handles for models or subsystem. |
| isthss | ISTHSS Tests if the model structure is of state-space type |
| isvector | function isv = isvector ( u ) ; |
| iv | IV Computes instrumental variable estimates for single output ARX-models. |
| iv4 | IV4 Computes approximately optimal IV-estimates for ARX-models. |
| iv4mv | IV4MV Computes approximately optimal IV-estimates for multivariate ARX-models |
| ivar | IVAR Computes IV-estimates for the AR-part of a scalar time series. |
| ivstruc | IVSTRUC Computes the output error fit for families of single-output ARX-models. |
| ivx | IVX Computes instrumental variable estimates for ARX-models. |
| ivxmv | IVXMV Estimates (multivariate) ARX-models using instrumental variables |
J | |
| Jozef_help | For large data sets - we'll assume that the data is in a matrix called |
K | |
| keep_var | [xkeep, ykeep] = keep_var(lim, x, y); |
| keep_var | [xkeep, ykeep] = keep_var(lim, x, y); |
| keep_var | [xkeep, ykeep] = keep_var(lim, x, y); |
| keep | out = keep(varargin); |
| keep_latlon | [lat, lon] = keep_latlon(lim, x, y); |
| keep_latlon2 | [lat, lon] = keep_latlon(lim, y, x); |
| keep_var | [xkeep, ykeep] = keep_var(lim, x, y); |
| keep_var2 | [xkeep, ykeep] = keep_var(lim, x, y); |
| keep_var3 | [xkeep, ykeep] = keep_var(lim, x, y); |
| keep_varm | [xkeep, ykeep] = keep_var(lim, x, y); |
L | |
| lagged_correlation | Load PCS: |
| lagged_regressions | |
| link_toolboxes | This file should be run only at the beginning |
| label_eastwest | function lab = label_eastwest ( xe, nd ); |
| lag_regress_on_ct | cint = [0.1 0.1 0.1]; |
| lagged_1000yr_regression | |
| lagged_corr_pcs | |
| lagged_corr_slpx_ct | load MIX_slpx_DYN_ATM_pslpx_LP9.mat; dofilt = 1; |
| lagged_correlation | Load PCS: |
| lagged_filt_ct_regressions | ctlim = [110 300 -75 65]; |
| lagged_regression | Load the pcs: |
| lagged_regressions | |
| landfill | h = landfill(elev); |
| landshade | h = landshade(elev); |
| laplacian | lap = laplacian(y, x, order, wrap); |
| later_years | clear |
| ldeo_sst_calc_GR | Load and configure NMC.REANAL data |
| lin_remove | lin_remove: linearly remove a time series from data |
| lin_remove_NaN | function Y = lin_remove_NaN ( Xdat , Xtim , tol , show ) ; |
| lin_remove_NaN_defunct | lin_remove_NaN: linearly remove a time series from data |
| link_alexis_scripts | test |
| link_alexis_scripts2 | Start by defining global variables TRUE and FALSE: |
| link_toolboxes | This file should be run only at the beginning |
| lm_over | function hout=landmask(lc,oc,kc,le,oe,ke,lwdth,lgrd,grdx,grdy) |
| lm_under | function hout=landmask(lc,oc,kc,le,oe,ke,lwdth,lgrd,grdx,grdy) |
| loadct | |
| loadpcs | |
| look_3dceof_temp_flux | Get regressions of fluxes: |
| look_3dceof_vort_div | Get regressions of fluxes: |
| look_at_coup_run | cd /home/disk/hayes2/dvimont/ocean/coup/stoc |
| look_at_stoc_run | cd /home/disk/hayes2/dvimont/ocean/coup/coup |
| look_dutdx | Load CPCs |
| look_hc_terms | lims2 = [180 270 -3 3]; tit = 'EQ'; % Equatorial index |
| look_prec_ITCZ_run | |
| LP10_T_zonal_150to180 | !PS-Adobe-3.0 |
| LP9_eq_tend_terms | Get windstress if necessary |
| LP9_spatial_tendterms | prec = getflx('rnd', lims, tim); |
| LP_HC_evolution_30s30n | Look at movie: |
| LPPC1_regression | Plot data |
| ls | |
| lynn_commands | height=load('geop500mb98111000.txt'); |
M | |
| make_time_plots | Look at CT runs: |
| ml_coup_slp_eof | load ML_ANN_sst.mat; slp = sst2; |
| ml_slp_pcs_ppcs_2 | Get ppcs |
| mshade | mshade(data, lev, colr); |
| mean2 | function [ mx, nsum, sx] = mean2 ( x ) ; |
| maptrimp2 | MAPTRIMP Trims a patch map to a specified region |
| mapfill | hh = dcm2 ( lev , colr ) ; |
| maptrimp2 | MAPTRIMP Trims a patch map to a specified region |
| make_5x5_average | Start with SLP: |
| make_annual_heatflux | MIX data |
| make_clim_maps | Define axis properties |
| make_daily_hgts | This program makes one big variable with all the daily |
| make_level5_temp | |
| make_ML_annave | varn = 'psl' |
| make_regressions | The following files are in the directory below: |
| make_som_dom_plots | CONVERGENCE |
| make_time_plots | Load variables directly from netcdf files |
| make_time_plots2 | Read Variables from wgr and cgr runs |
| make_Z500_plots | |
| makebars | MAKEBARS Make data for bar charts. |
| manual_ticks | function manual_ticks ( xyz, tick_positions, tick_labels, ca ); |
| map_axis | hh = map_axis(mstyle, origin, frme); |
| map_contour | [c, h] = map_contour(data, cint, line_style); |
| map_contour_pn | [c, h] = map_contour(data, cint, usezero); |
| map_quiver | h = map_quiver(xdat, ydat, scale, nskip); |
| map_surface | h = map_surface(data, alt, 'PropertyName', ProptertyValue, ...); |
| map_surface_interp | h = map_surface(data, alt, 'PropertyName', ProptertyValue, ...); |
| mapfill | hh = dcm2 ( lev , colr ) ; |
| maptrimp | MAPTRIMP Trims a patch map to a specified region |
| maptrimp2 | MAPTRIMP Trims a patch map to a specified region |
| maptrimp3 | MAPTRIMP Trims a patch map to a specified region |
| maxes | maxes(mstyle, origin); |
| maxmin | [maxind, minind] = maxmin(x); |
| mcolor | mcolor(data, cint, mstyle); |
| mcolor3 | mcont3(data, cint, mstyle, origin); |
| mcolor4 | [h1, h2] = mcolor(data, mstyle, center); |
| mcont | mcont(data, levs, mstyle, origin) |
| mcont2 | [h, c] = mcont2(data, cint, mstyle, zeroline); |
| mcont3 | mcont3(data, cint, mstyle, origin); |
| mean2 | function [ mx, nsum, sx] = mean2 ( x ) ; |
| mf2th | MF2TH Packages user defined model structures into the THETA model format. |
| midprefs | MIDPREFS Choose directory for idprefs.mat, the ident start-up info file. |
| mincrit | MINCRIT Minimizes prediction error criteria |
| MIX_LP_regressions | MIX run |
| MIX_regressions | load SH_MIX_slpx_DYN_ATM_pslpx_RAW.mat; dofilt = 0; |
| mketaarx | MKETAARX constructs an ETA-model structure for (multivariate) ARX-models |
| mktheta | MKTHETA synonymous to POLY2TH |
| ml_coup_slp_eof | [reg3, c3] = regress_eof(sstc, lpseudopcs(:,num), 0); |
| ml_eof_regressions | [reg1, c1] = regress_eof(sst, tim, 0); |
| ml_ppcs_thesis | load ML_SLP_eof_shem_-60to-20.mat |
| ML_projections | First, project rlds and lds onto slp and slpc, respectively |
| ml_run_stats | varn = 'psl' |
| ml_run_stats2 | Average by season, to make this a bit more manageable |
| ml_run_stats3 | Get CT index |
| ml_run_stats4 | for biff = 1:3; |
| ml_slp_pcs_ppcs | Get ppcs |
| ml_slp_pcs_ppcs_2 | Get ppcs |
| model_ekman_upwelling | First, find where the thermocline is |
| modstruc | MODSTRUC Constructs model structures to be used in MS2TH. |
| motor | MOTOR Describes the dc-motor with time-constant t (= par) and |
| movie_ceofs | Try a movie of this |
| movie_hc_25yrlp | |
| movies | Total movie |
| mquiv | mquiv(xdat, ydat, len, mstyle, origin) |
| ms2th | MS2TH Packages standard state-space parameterizations into the THETA format. |
| mshade | mshade(data, levs, mstyle, origin) |
| mshade3 | [c, h] = mshade3(data, lev, colr); |
| mvarx | MVARX Estimates (multivariate) ARX-models |
| myrunning_ave | function y = running_ave ( x, n, usenan ) ; |
| myrunning_sum | function [y, nsum] = running_sum ( x, n, usenan ) ; |
N | |
| n4sid | N4SID Estimates a state-space model using a sub-space method. |
| name | NAME Method to check or create valid filename. |
| nao_ccm | |
| nc_new_time | Repeat for cold.nc |
| ncep | |
| ncread | |
| newpage | NEWPAGE Method to finish the current page of a PrintJob and start next. |
| newtime | keep = find((newdate ~= 19851001) & (newdate ~= 19851101)); |
| newtime2 | |
| nhem_eofs | file_base = [1:73]; |
| nino3 | x = nn3(sst, nxa, nya, frame); |
| north | er=north(s,dat); |
| noselection | NOSELECTION Select/Deselect all objects in Figure. |
| np_vs_eq_ceof_look | load LPregmap_3deof_totpac.mat |
| npac_sst_eofs | Save the data |
| nuderst | NUDERST Selects the step size for numerical differentiation |
| nyqplot | NYQPLOT plots the Nyquist diagram of a frequency function. |
O | |
| ocean_current_clim | [u, v] = getnc('u', 'v', ctlim, lev, tim); |
| ocean_vel | print -dps2 Oce_currents.ps |
| ocean_velocity_vert_ave | xlabel(['Max = ' num2str(round(c(2)*100)/100) ' cm s^-^1']) |
| oe | OE Computes the prediction error estimate of an output-error model. |
P | |
| plot_nmc | |
| plot_60s60n_terms | Load PCS: |
| plot_lp_pac_regressions | Plot the data |
| plot_pac_terms | order = ['ubtp'; 'vbtp'; 'wbtp'; 'uptb'; 'vptb'; 'wptb']; |
| plot_regressions | |
| plot_30s30n_terms | Get ubar_dHCprimedx and vbar_dHCprimedy |
| plot_30s30n_EOF_regressions | load gr_pcs.mat lam lds pcs per lims cthp gr |
| plot_30s30n_regressions | Plot the data |
| plot_30s30n_terms | Get ubar_dHCprimedx and vbar_dHCprimedy |
| plot_ceof_heat_sft | [temp, lat, lon, depth, middepth] = getheat(lims, lev, tim); |
| plot_pac_terms | prec = getflx('rnd', lims, tim); |
| plot_wbtp_tcline_reg | Plot wbT'_z |
| plot_ceof_heat_sft | [temp, lat, lon, depth, middepth] = getheat(lims, lev, tim); |
| pdo_statistics | subplot(2,2,2*i-1); |
| printjob | PRINTJOB Constructor for PrintJob objects. |
| pncontm | [c, h] = pncont(xa, ya, x, clev, lab, s); |
| pac_atl_3dCEOF_look | Look at lagged regressions |
| pac_setup1 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
| pac_setup_old | Average sst to make computations a little more reasonable |
| pac_vs_global | mcont(bot, clevb, 'giso', [0 180]); |
| pc_spectrum | ctlim = [110 300 -30 30]; |
| pc_statistics | |
| pd | |
| pdo_calc | look at SST regressions, to determine what modes 1 and 2 are |
| PDO_looks | Look at partial autocorrelation function |
| pdo_statistics | |
| pdo_statistics2 | |
| pe | PE Computes prediction errors. |
| pefilt | PEFILT Auxiliary routine to PE |
| pem | PEM Computes the prediction error estimate of a general linear model. |
| pemss | PEMSS Minimizes prediction error criteria |
| phase | PHASE Computes the phase of a complex vector |
| plot_15n65n_regressions | Plot the data |
| plot_15n65n_terms | Get ubar_dHCprimedx and vbar_dHCprimedy |
| plot_30s30n_EOF_regressions | load LP9_detrend_L1-7_EOF_yr101-1000.mat |
| plot_30s30n_EOF_terms | Get ubar_dHCprimedx and vbar_dHCprimedy |
| plot_30s30n_regressions | Plot the data |
| plot_30s30n_terms | Load PCS: |
| plot_3d | lim(i) = max(max(max(max(tcor(i,:,:,:))))); |
| plot_60s60n_terms | Load PCS: |
| plot_adv | Load PCS: |
| plot_ceof_heat_sft | [temp, lat, lon, depth, middepth] = getheat(lev, tim, lims); |
| plot_climatology | Plot the climatology for the backward rotating earth run. |
| plot_ct_index_quick | |
| plot_daily_1daycorr | Start with earth |
| plot_daily_hgts | Round 2 |
| plot_eddy | top = squeeze(mean(vtw500-vta500)); |
| plot_eke | lims = [30 150 0 60]; |
| plot_EOF1_EOF2_0-270mHC | |
| plot_eq_currents | tser = lpcs(:,1); lags = -8:2:6; biff = 2; |
| plot_hctend_terms | |
| plot_hgt | |
| plot_lag_regress2 | Load PCS: |
| plot_lag_regressions | color_shade(squeeze(tcoef(i+ind,:,:).^2), 0.1, [.9 .9 .9]); |
| plot_lag_regressions_color | Heat content and ocean currents |
| plot_lp_pac_regressions | Plot the data |
| plot_lp_pac_terms | %%%%%%%%%%%%%%%% Plot dynamic heat content terms |
| plot_mmc | Start with EARTH |
| plot_monthly_ct_time | |
| plot_nmc | title('PRECIPITATION: NMC NDJFMA regressed on GR index'); |
| plot_off_cdrom | |
| plot_olr | Get oi_clim data |
| plot_pac | EOF's |
| plot_pac_regressions | 11111111111111111111111111111111111111 Plot the data |
| plot_pac_terms | prec = getflx('rnd', lims, tim); |
| plot_pcs | tim(timind), pcs(timind, 2), '--k'); |
| plot_PREC_for_SUN | print -dps2 CT_OMEGA_PRECIP.ps |
| plot_precip | |
| plot_reg_HEAT_SFC | load LP10_detrend_L1-7_yr101-1000.mat; pcs3 = pcs; |
| plot_reg_TAUX_z250 | num2str(round(per(num))) ... |
| plot_regressions | Old stuff |
| plot_regressions_HP_LP | load RAW_detrend_L1-7_EOF_yr101-1000.mat; rpc = -1*pcs; |
| plot_regressions_HP_LP_thesis | varn = 'temp'; tit = 'SST'; units = 'K'; cint = 0.05; pn = 1; cint2 = 0.25; |
| plot_som_hard | |
| plot_sst_anom | |
| plot_sst_z | |
| plot_stereo2 | plot_stereo(x, lat, lon, proj,orig); |
| plot_tendencies | Get pcs |
| plot_TPac_vs_NMC | clhgt = atlev(hgt3, lev1, ps2, hyam, hybm, P0); |
| plot_tpgr | top = squeeze(mean(tem1(:,ind,:,:)-tem3(:,ind,:,:))); |
| plot_variance_explained | temp = getnc('temp', lims, lev, tim); |
| plot_velocity_vertical | Look at zonal cross sections |
| plot_vort | Try some plots removing the zonal average. |
| plot_wbtp_tcline_reg | temp2 = squeeze(shiftdim(mean2(shiftdim(temp, 2)), 2)); |
| plot_wdtdz_terms | Load PCS: |
| plot_wpgr | |
| plot_z | |
| plot_zerolag_terms | load LP9_detrend_L1-7_EOF_yr101-1000.mat |
| pn_reg_ct | [ppat, npat, pat, clim] = reg_ct(filin, var_nam, lev); |
| PNA_index | |
| PNA_index_regs | Get data |
| pncont | [h, c] = pncont(xa, ya, x, clev, lab, s); |
| pncont_old | pncont(xa, ya, x, clev, lab); |
| pncontf | pncont(xa, ya, x, clev, lab, s); |
| pncontm | pncont(xa, ya, x, clev, lab); |
| pncontm2 | pncont(ya, xa, x, clev, lab, s); |
| pnpcolorm | PCOLORM Projected matrix map in the z = 0 plane |
| pnsemilogy | [h, c] = pncont(xa, ya, x, clev, lab, s); |
| poly2th | POLY2TH Constructs a "theta-matrix" from given polynomials. |
| polyform | POLYFORM computes the polynomials associated with a given model |
| pop_anal_tropics | Get data |
| pop_anal_tropics2 | Next, SST |
| pop_hold | function pop_hold ( arg1 ) ; |
| pop_stack | function [ stack, var ] = pop_stack ( stack ) ; |
| pop_temp | function var = pop_temp ( arg1 ) ; |
| popd | function popd ( silent ) ; |
| pos_vs_neg_LPCT | The following files are in the directory below: |
| positions | POSITIONS Determine the union of the PaperPositions of all objects passed in. |
| prec_compare | |
| prec_GR_CT | print -dps2 PREC_CT_GR_t31_DOM.ps |
| precip_ct_nmc | filin = 'prate.mon.mean.nc'; |
| predict | PREDICT Computes the m-step ahead prediction. |
| predicts | PREDICTS Auxiliary function to PREDICT. |
| prepare | PREPARE Method to modify a Figure or Simulink model for printing. |
| preparehg | PREPAREHG Method to ready a Figure for printing. |
| preparepointers | PREPAREPOINTERS Set Pointers of all open Figures to Watch. |
| prepareui | PREPAREUI Method to draw Uicontrol objects in output. |
| present | PRESENT presents a parametric model on the screen. |
| presenta | PRESENTA Subroutine to present. |
| presents | PRESENTS Presents a parametric model on the screen. |
| print2 | PRINT Print Figure or model. Save to disk as image or M-file. |
| print_clim | This should print a list of ccm2nc commands for the oi_82_98 |
| printccm2nc_dan | for i=1:64 |
| printccm2ncbig | for i=1:64 |
| printjob | PRINTJOB Constructor for PrintJob objects. |
| printnewlink | a = [' 1981-12'] |
| printtau | printtau(tau,fname); |
| privdos | DOS Execute DOS command in save current directory and return result. |
| project_windstress | load heat_pcs_eq.mat |
| ptpreparehg | PREPAREHG Method of PrintTemplate object that formats a Figure for output. |
| ptrestorehg | FORMAT Method that restores a Figure after formating it for output. |
| push_hold | function push_hold |
| push_stack | function stack = push_stack ( val, stack ) ; |
| push_temp | function push_temp ( var ) |
| pushd | function pushd ( directory, silent ) ; |
Q | |
| quiverm2 | QUIVERM Two dimensional quiver plot projected on a map axes |
| quiverm2 | QUIVER Quiver plot. |
| quiver2 | QUIVER Quiver plot. |
| quiver3 | QUIVER3 3-D quiver plot. |
| quiver_barb | QUIVER Quiver plot. |
| quiver_lab | QUIVER_LAB Quiver plot. |
| quiver_label | USAGE: quiver_label(U,scale,xtext,ytext,xarrow,type) |
| quiver_label2 | USAGE: quiver_label(U,scale,xtext,ytext,xarrow,type) |
| quiver_label_giso | USAGE: quiver_label_giso(U,scale,xtext,ytext,type) |
| quiver_labelm | USAGE: quiver_label(U,scale,xtext,ytext,xarrow,type) |
| quiver_lau | QUIVER Quiver plot. |
| quiver_orig | QUIVER Quiver plot. |
| quiver_scale | QUIVER_SCALE Quiver plot. |
| quiver_vert | QUIVER Quiver plot. |
| quiverm | QUIVERM Two dimensional quiver plot projected on a map axes |
| quiverm2 | QUIVER Quiver plot. |
| quiverscale | QUIVERSCALE creates a scale at the bottom of the quiver plot, |
| quivscale2 | quiver(XAX, YAX, a, a, 1); |
R | |
| reemerge1 | Plot the data |
| regress_flux | print -dps2 HP10_Prec_Icon_lag0.ps |
| regress_velocity | |
| regress_flux | tim = [451:900]; |
| rossby_kelvin_structure | Note: equatorial rossby radius of deformation is given by: |
| range | function [ minval, maxval ] = range( array ); |
| random | |
| range | function [ minval, maxval ] = range( array ); |
| rarmax | RARMAX Computes estimates recursively for an ARMAX model. |
| rarx | RARX Computes estimates recursively for an ARX model. |
| rave | function y = rave ( x, n, usenan ) ; |
| rbj | RBJ Computes estimates recursively for a BOX-JENKINS model. |
| rcont | rcont(x, n) |
| rd | help readnc |
| rdhis | reads data from a netcdf file with multiple time steps (monthly data) and |
| rdnc | [p0,hyai,hybi,hyam,hybm,t,u,v,z3,ps,psl,omega,precl,precc] = rdnc(filin); |
| read_clim | [dat, lat, lon, lev] = ... |
| read_hgt | [slpnmc, nxanmc, nyanmc, slph, nxah, nyah, slpc, nxac, nyac]=read_nmc; |
| read_np | [sst, oro, nxa, nya]=read_sst2; |
| read_sst | [sst, oro, nxa, nya]=read_sst2; |
| read_time | out1 = read_time(fdir, lev, var1); |
| read_write | [slpnmc, nxanmc, nyanmc, slph, nxah, nyah, slpc, nxac, nyac]=read_nmc; |
| read_z3 | [z, nxa, nya]=read_z3; |
| Readme | Release Notes System Identification Toolbox 4.0.3 |
| reemerge1 | Try the COADS data |
| reg_ccmsst_on_nmc | Now, get NMC data |
| reg_ct | [pat, clim] = reg_ct(filin, var_nam, lev); |
| reg_gr | [pat, clim] = reg_ct(filin, var_nam, lev); |
| reg_obs_onto_BPppcs | Load variables for regression |
| reg_SST_onto_SLP_variosBPs | Weight by cosine lat: |
| reg_stuff_onto_BPppcs | Load variables for regression |
| regress | clear |
| regress_20dISO_taux_obs | Get 20dISO depth |
| regress_ceof | [out, temtim] = regress_ceof(dat, pcs, nfrm); |
| regress_CPC_windstress | Load CPCs |
| regress_dan | fld=regress(tm,dat); |
| regress_ddt_co2 | Load PC1 |
| regress_eof | [out, ccoef] = regress_eof(dat, pcs, lags); |
| regress_EQ_SST_indices | For fun, let's see how much mass the atmosphere loses. |
| regress_flux | tim = [451:900]; |
| regress_gr_eddy | Look at 850mb VT: |
| regress_hp_lp_ct | cut down on the size of the data |
| regress_justin | function Y = lin_remove_NaN ( Xdat , Xtim , tol , show ) ; |
| regress_nan | out1 = regress_nan ( vect, data ) ; |
| regress_NP_eof | lims = [110 265 20 65]; % North Pacific Region |
| regress_out | DETREND Linearly remove time series from data. |
| regress_sv | Load PCS: |
| regress_taux_hanom | eta/h = 6e-4; |
| regress_taux_hanom2 | Set up coefficients for centered difference |
| regress_thermocline | Load PCS: |
| regress_velocity | Start with CT regressions |
| regress_vertical_temp | [ntim, nlev, nlat, nlon] = size(temp1); |
| regress_wbc | [temp, lat, lon, depth, middepth] = getheat(lims, lev, tim, varn); |
| remove_anncyc | [x, clim] = annave(y); |
| remove_mean | [x, clim] = annave(y); |
| render | RENDER Method to draw a model or Figure on current page of a print job. |
| reshape_ann_to_mon | [outdat, ntim, nlat, nlon] = reshape_ann_to_mon(indat); |
| resid | RESID Computes and tests the residuals associated with a model. |
| residual_term | The vertical diffusivity between the first and second |
| restore | RESTORE Reset a Figure or Simulink model after printing. |
| restorehg | RESTOREHG Reset a Figure after printing. |
| restorepointers | RESTOREPOINTERS Restore Pointers of all open Figures. |
| restoreui | RESTOREUI Remove Images used to mimic user interface controls in output. |
| rm_mean2 | function y = rm_mean2 ( x ) ; |
| roe | ROE Computes estimates recursively for an output error model. |
| rossby_kelvin_structure | Note: equatorial rossby radius of deformation is given by: |
| rossby_kelvin_structure_2 | Note: equatorial rossby radius of deformation is given by: |
| rotate_gr | Load and configure NMC.REANAL data |
| rotvar | ROTVAR computes the standard deviations of roots to polynomials |
| rpem | RPEM Computes estimates recursively for a general model. |
| rplr | RPLR Computes PLR estimates recursively for a general model. |
| rsum | function [y, nsum] = rsum ( x, n, usenan ) ; |
| running_average | function y = running_ave ( x, n, usenan ) ; |
| rvect | rvect(x, y, kp, len) |
| rw_nmc | rw_nmc; |
S | |
| sd | function sdriver(c,dt,mx,maxtime,lat,iplot,distribution) |
| sd2 | function sdriver(c,dt,mx,maxtime,lat,iplot,distribution) |
| subplot2 | this_ax = subplot2(nCol, thisPlot); |
| std2 | function stdx=std2(x) |
| sum2 | function [sx, nsum] = sum2 ( x ) ; |
| sp | h = sp(nCol, thisPlot); |
| sp_old | h = sp(pos); |
| startup | |
| save_20deg_ceof | Do complex eof: |
| screenpos | SCREENPOS Position Figure using Points units but keep it on screen. |
| sd | SUBPLOT Create axes in tiled positions. |
| sd2 | SUBPLOT Create axes in tiled positions. |
| search | SEARCH Searches for a lower value of the criterion function |
| searchax | SEARCHAX searches for lower values of the prediction error criterion. |
| searchbj | SEARCHBJ searches for lower values of the prediction error criterion |
| searchoe | SEARCHOE searches for lower values of the prediction error criterion. |
| seasave | av = seasave(decx, x); |
| seasonal_eof | |
| seasonal_gr | Look at Surface Winds: |
| seasonal_spect | save seas_ct.mat ct1 ct5; |
| segment | SEGMENT Segments data and tracks abruptly changing systems. |
| selstruc | SELSTRUC Selects model structures according to various criteria. |
| send | SEND Send output file to hardcopy device. |
| sess | SEARCH Searches for a lower value of the criterion function |
| setpp | SETPP Set all paper properties. |
| setset | SETSET Do a SET on Handle Graphics object or a SET_PARAM on Simulink object. |
| sett | SETT Sets the sampling interval in TH-structures and frequency functions. |
| setup | SETUP Open the printer setup dialog. |
| sfunid | SFUNID an S-function which performs system identification. |
| shade_ccoef | shade_ccoef(data); |
| shade_solid | shade_solid(data, lev, color) |
| shade_solid2 | shade_solid(data, lev, alt, color) |
| shade_tscore | greyshd(x, lev) |
| shortrun | % Movie |
| sigifun | This function is used in solving for sigma_i in the |
| sitb | The SYSTEM IDENTIFICATION TOOLBOX contains the following M-files: |
| slblocks | SLBLOCKS Defines the block library for a specific Toolbox or Blockset. |
| smooth | dout = smooth(vect, data); |
| som_prec | filin = '36yearly.cdf'; |
| som_setup | tsnew = tsnew(:,keep); |
| som_ts_eof | filin = 'histape.cdf.36'; |
| sp | h = sp(nCol, thisPlot); |
| sp2 | h = sp(pos); |
| sp_old | h = sp(pos); |
| spa | SPA Performs spectral analysis. |
| spback | h = sp(pos); |
| sph_curl1 | function [ vor, lat2, lon2 ] = sph_curl ( u, v, lat, lon, wrap ) ; |
| sph_div1 | [div, lat2, lon2] = sph_div1(u, v, lat, lon, wrap); |
| sph_grad | function [ x, y ] = sph_grad ( z, ax, ay ) ; |
| sph_gradx | function x = sph_gradx ( z, rx, ry ) ; |
| sph_gradx1 | dfdx = sph_gradx1(data, ry, rx, wrap); |
| sph_grady | function y = sph_grady ( z, rx, ry, icos ) ; |
| sph_grady1 | [dfdy, ry2] = sph_grady1(data, ry, rx, icos); |
| sptalk | subplot_talk(nrow, ncol, thisPlot); |
| spthes | h = sp(pos); |
| squeeze_axes | axis_limits(num) |
| ss2th | SS2TH Produces a parameterized state-space model. |
| ssmodel | SSMODEL In case of a userwritten model structure definition, this |
| ssmodx8 | SSMODX8 The standard state-space model |
| ssmodx9 | SSMODX9 The standard state-space model |
| ssssaux | SSSSAUX Auxiliary file to n4sid. |
| sst_ceof | data; save ML_NH_ceof.mat lam lds pcs per lims; back |
| sst_clim_compare | cd /home/disk/tao/data/coads/coads1a |
| sst_clim_diffs | nino_oi.m |
| sst_complex_eof | The following files are in the directory below: |
| sst_eof | The following files are in the directory below: |
| sst_eof2 | |
| sst_eof_heat_regress | save sst_eof_bp_4.5-50yr.mat ld10 pc10 lam per kp lat lon ctlim |
| sst_fcst | Load data from gr_nmc |
| sst_tau_hc_color_fig | ct = getnc('temp', [180 270 -6 6], 1, tim); |
| sst_tau_hc_fig | ct = getnc('temp', [180 270 -6 6], 1, tim); |
| standardize | function y = standardize ( x ) |
| standardize_NaN | function y = standardize ( x ) |
| start | START Method to start a print job. |
| startup | |
| startup2 | |
| startup_orig | startup file |
| startup_wallace | |
| status | STATUS Return structure of status on PrintJob. |
| std2 | function stdx=std2(x) |
| std_ct | |
| std_NaN | stdx = std_NaN(x); |
| stereo2 | STEREO Stereographic Azimuthal Projection |
| storm_1pt_corr_map | Look at correlation maps |
| struc | STRUC generates typical structure matrices for ARXSTRUC and IVSTRUC. |
| subplot | SUBPLOT Create axes in tiled positions. |
| subplot2 | this_ax = subplot2(nCol, thisPlot); |
| subplot_dan | SUBPLOT Create axes in tiled positions. |
| subset | x = subset2(y, lat, lon, lims); |
| subset_old | x = subset(y, nxa, nya, frame); |
| sum2 | function [sx, nsum] = sum2 ( x ) ; |
| sum4vms | SUM Sum of the elements. |
| sumboard | SUMBOARD Creates the ident figure. |
| sun_to_alpha_precip | |
| surf_wind_compare | Load GR pattern and time series |
| surfacem2 | SURFACEM Display a matrix map warped to a projected graticule |
| svd_intraseasonal | Get landmask |
| SVD_routine | Weight data accordingly |
| svd_slp_taux | Load data |
| svdplot | set(gca,'YTick',[-.5:.25:.5]) |
T | |
| tendency_terms | Load PCS: |
| tendency_terms | Load PCS: |
| tcline_regression | |
| thin | out1 = thin(in1, inc); |
| T42_t31_compare | [clevt, cintt] = get_contour_int(top, 12); |
| table_correlation | ct = getheat([180 270 -6 6], 1:3, tim); |
| tao_windstress_std | 1.6199 : -6 to 6; |
| TAU_curl | print -dps2 TAU_prec.ps |
| tau_to_strf | [strf, lat_out, lon_out] = tau_to_strf(taux, tauy, lat, lon, bc, H); |
| tau_to_strf2 | [strf, lat_out, lon_out] = tau_to_strf(taux, tauy, lat, lon, bc, H); |
| tcline | print -dps2 Tcline.ps |
| tcline_regression | subplot(7,2,2*i-1); |
| temp_3d_ceof_look | for ind = 1:10; |
| temp_3d_eof | ctlim = [-0.1 17.5 47.5 80]; |
| temp_epflux | |
| tendency_terms | Load PCS: |
| tendency_terms2 | Get pcs |
| tendency_terms_lreg | Load PCS: |
| test_AR2_spect | cd ~/model/enso |
| test_eof_routine_AO | |
| th2arx | TH2ARX converts a THETA-format model to an ARX-model. |
| th2ff | TH2FF Computes a model's frequency function,along with its standard deviation |
| th2par | TH2PAR converts the theta-format to parameters and covariance matrix. |
| th2poly | TH2POLY computes the polynomials associated with a given model. |
| th2ss | TH2SS Transforms a model in THETA-format to state-space. |
| th2ssaux | TH2SSAUX An auxiliary routine to TH2SS |
| th2tf | TH2TF Transforms from the THETA-format to transfer functions. |
| th2zp | TH2ZP Computes zeros, poles, static gains and their standard deviations. |
| thc2thd | THC2THD Converts a continuous time model to discrete time. |
| thd2thc | THD2THC converts a model to continuous-time form. |
| theor_coup_amp | Define stuff |
| thermocline_structure | ctlim = [-.1 360 -35 -20]; |
| thes | |
| theta | THETA is a matrix containing information about model structure, estimated |
| thin | out1 = thin(in1, inc); |
| thinit | THINIT Initial parameter values for iterative estimation. |
| thss | THETA(SS) A structure for defining general linar state-space models |
| thss2th | THSS2TH Converts an internal TH(SS)-format to the standard multi-input- |
| tightmap | TIGHTMAP removes whitespace around a map |
| tightmap2 | TIGHTMAP removes whitespace around a map |
| tightmap_orig | TIGHTMAP removes whitespace around a map |
| time_freq_anal | Plot hoffmoeller diagram of heat content and SST |
| time_series | |
| timeseries | r = rand(1, 100000)-0.5; |
| title2 | title2 Graph Title on upper right corner |
| topography | tem2 = mean(tem2); |
| tpac_global_compare | cnum = 10; |
| trf | TRF Computes a model's frequency function. |
| trfcont | TRFCONT Computes the continuous time transfer function of a model. |
| trfsaux | TRFSAUX Auxiliary function to TRFSS |
| trfsd | TRFSD Computes a model's frequency function,along with its standard deviation |
| trfss | TRFSS Auxiliary routine to TRF |
| trfsssd | TRFSSSD Auxiliary routine to TH2FF |
| ts_eof2 | Look at regressions |
| tscore | score = tscore(dof, tail_per); |
| ttest | tstat = ttest(x1, x2); |
| ttest_plot | win = sort([2:12:ntim 3:12:ntim 4:12:ntim]); |
U | |
| ubar_grad_dHCprime | Loop through HP LP and RAW |
| ubar_gradtprime | Start with ubardt'dx and vbardt'dy |
| uncosweight | x=cosweight(sst, nya); |
| unfixpar | UNFIXPAR Unfixes parameters in state-space and ARX model structures. |
| unfixxx | UNFIXXX Subroutine for UNFIXPAR |
| uprime_gradHCbar | Start with u'dtbardx, v'dtbardy |
| uprime_gradtbar | Start with u'dtbardx and v'dtbardy |
V | |
| validate | VALIDATE Method to check state of PrintJob object. |
| var_nan | function v = var_nan ( mat1 ) ; |
| varmaxt | USE: [Frot,AT,Cscor,Vrot,h] = varmaxt(Fm,L,norm,A) |
| vectsig | [fval, dof1, dof2] = vectsig(c1, c2); |
| vel_to_strfnc | [strf, lat_out, lon_out] = vel_to_strfnc(vg, lat, lon); |
| velocity_plots | %%%%%%%%%%% |
| vert_ubtp_uptb | [ubtp, uptb] = vert_ubtp_uptb(pcs, lims, lags); |
| vert_vdtdy_wdtdz | [vdtdy, wdtdz] = ... |
| vert_vptb_wptb | [vdtdy, wdtdz] = ... |
| vertical_ocean_structure | ctlim = [-.1 360 -6 6]; |
| vi | function vi ( arg1, arg2, arg3, arg4, arg5, arg6, arg7, arg8, arg9, arg10 ) |
| vline | function h = vline ( x, linetype, y, zdat ) ; |
| vphi2 | this function makes a vector plot, and plots the |
| vphi4 | this function makes a vector plot, and plots the |
W | |
| wbar_dtpdz_hc0-80 | Load PCS: |
| w_intbar_dHCprimedz | First, get ubar and vbar |
| w_intbar_dtprimedz | First, get ubar and vbar |
| w_intprime_dHCbardz | ctlim = [108 302 -62 62]; |
| w_intprime_dtbardz | First get tbar |
| walker_circ | o = owd; v = vwd; u = uwd; |
| walker_eof | |
| warm_WPAC.define | |
| wave_speed | Kelvin wave |
| wbar_dtpdz_hc0-80 | Load PCS: |
| wbar_dtprimedz | First, get wbar |
| we_bar_dtprimedz | First, get wbar |
| where_sig | [locmax, locmin] = where_sig(x); |
| wind_CT_GR | |
| wind_stress | load LP10_L1-7_CEOF.mat; tit = 'LP10'; |
| windstress_regcoef | First, project rlds and lds onto slp and slpc, respectively |
| work_dan | Step 1: Calculate an SST anomaly data set for the 4 by 6 |
| worldmap2 | WORLDMAP maps a country or region using the WORLDLO atlas data |
| write_covar | cd /home/disk/hayes2/dvimont/nmc/ugrd.vgrd |
| write_list | function list = write_list ( y, ndigit ) ; |
| write_netcdf | Define dimensions |
| wstp_upwelling | First look at mean: |
| wtrsfcairtemppna | This program draws regression and one-pt. correlation maps corresponding to |
Z | |
| zero_fillnan | function x = zero_fillnan ( x ) ; |
| z500_eof | z5 = getnc('NH_500mb_DJF.nc', 'Z3'); |
| z500_setup | pc = pc * diag(1./sqrt(lam)); |
| zero_fillnan | function x = zero_fillnan ( x ) ; |
| zonal_depth_eq | print -dps2 eq_tcline_neglag.ps |
| zonal_eddy | |
| zonal_u | Just a side note here. Recall that the SSTA pattern |
| zp | ZP Computes zeros, poles and static gains associated with a model |
| zpform | ZPFORM Formats zero-pole matrices from different models for use in ZPPLOT |
| zpplot | ZPPLOT Plots zeros and poles. |
| zpsd | ZPSD See TH2ZP |
| zpsdpl | ZPSDPL Plots standard deviations in zero-pole plots. |
| zpss | ZPSS An auxiliary function to ZP. |
| zpssnosp | ZPSS An auxiliary function to ZP |
| zpsssd | ZPSSSD An auxiliary function to TH2ZP |