help-octave
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Getting Matlab geophysics code running in Octave


From: Lester Anderson
Subject: Re: Getting Matlab geophysics code running in Octave
Date: Mon, 5 Aug 2013 07:33:34 +0100

Hi Ben,

Yes it was the wrong defval - changed that it that sorts some issues! Thanks for spotting that.

octave:1> cd c:\matlab\functions\simmons_matlab
octave:2> coheref2
warning: isstr is obsolete and will be removed from a future version of Octave, please use ischar instead
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 62
62:   keyboard
debug> who
Variables in the current scope:

E     G     T     Te    ans   argn  drho  f2    g     lbd   r     rc    v

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 63
63:   defval('xver',0)
debug> who
Variables in the current scope:

E     G     T     Te    ans   argn  drho  f2    g     lbd   r     rc    v

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 64
64:   if xver==1
debug> who
Variables in the current scope:

E     G     T     Te    ans   argn  drho  f2    g     lbd   r     rc    v     xver

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 64
64:   if xver==1
debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 69
69:   k=2*pi./lbd;
debug> who
Variables in the current scope:

E     G     T     Te    ans   argn  drho  f2    g     lbd   r     rc    v     xver

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 73
73:   D=(E*Te.^3)/(12*(1-v^2));
debug> who
Variables in the current scope:

E     G     T     Te    ans   argn  drho  f2    g     k     lbd   r     rc    v     xver

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 76
76:   if length(Te)>=1 && length(f2)==1 && length(r)==1
debug> who
Variables in the current scope:

D     E     G     T     Te    ans   argn  drho  f2    g     k     lbd   r     rc    v     xver

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 76
76:   if length(Te)>=1 && length(f2)==1 && length(r)==1
debug> who
Variables in the current scope:

D     E     G     T     Te    ans   argn  drho  f2    g     k     lbd   r     rc    v     xver

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 76
76:   if length(Te)>=1 && length(f2)==1 && length(r)==1
debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 80
80:     [LL,FF2]=meshgrid(lbd,f2);
debug> who
Variables in the current scope:

D     E     G     T     Te    ans   argn  drho  f2    g     k     lbd   r     rc    v     xver

debug> whos
Variables in the current scope:

   Attr Name        Size                     Bytes  Class
   ==== ====        ====                     =====  =====
        D           1x1                          8  double
        E           1x1                          8  double
        G           1x1                          8  double
    f   T           1x1                          8  double
    f   Te          1x1                          8  double
        ans         1x1                          8  double
  a     argn        4x6                         24  char
    f   drho        1x1                          8  double
    f   f2          1x5                         40  double
    f   g           1x1                          8  double
        k           1x1000                    8000  double
    f   lbd         1x1000                    8000  double
    f   r           1x1                          8  double
    f   rc          1x1                          8  double
        v           1x1                          8  double
    f   xver        1x1                          8  double

Total is 2041 elements using 16160 bytes

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 81
81:     DD=D; RR=r;
debug> who
Variables in the current scope:

D     E     FF2   G     LL    T     Te    ans   argn  drho  f2    g     k     lbd   r     rc    v     xver

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 81
81:     DD=D; RR=r;
debug> who
Variables in the current scope:

D     DD    E     FF2   G     LL    T     Te    ans   argn  drho  f2    g     k     lbd   r     rc    v     xver

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 88
88:   KK4=(2*pi./LL).^4;
debug> who
Variables in the current scope:

D     E     G     RR    Te    argn  f2    k     r     v
DD    FF2   LL    T     ans   drho  g     lbd   rc    xver

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 89
89:   KK=(2*pi./LL);
debug> who
Variables in the current scope:

D     E     G     LL    T     ans   drho  g     lbd   rc    xver
DD    FF2   KK4   RR    Te    argn  f2    k     r     v

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 92
92:   xai=1+DD.*KK4/drho/g;
debug> who
Variables in the current scope:

D     E     G     KK4   RR    Te    argn  f2    k     r     v
DD    FF2   KK    LL    T     ans   drho  g     lbd   rc    xver

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 93
93:   phi=1+DD.*KK4/rc/g;
debug> who
Variables in the current scope:

D     E     G     KK4   RR    Te    argn  f2    k     r     v     xver
DD    FF2   KK    LL    T     ans   drho  g     lbd   rc    xai

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 94
94:   beta=rc./xai/drho;
debug> who
Variables in the current scope:

D     E     G     KK4   RR    Te    argn  f2    k     phi   rc    xai
DD    FF2   KK    LL    T     ans   drho  g     lbd   r     v     xver

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 99
99:   Ctop=(xai+FF2.*rc^2./drho^2.*phi-RR.*sqrt(FF2)*rc./drho.*[phi.*xai+1]).^2;
debug> who
Variables in the current scope:

D     E     G     KK4   RR    Te    argn  drho  g     lbd   r     v     xver
DD    FF2   KK    LL    T     ans   beta  f2    k     phi   rc    xai

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 100
100:   Cbot1=xai.^2+FF2.*rc^2./drho^2-2*RR.*sqrt(FF2)*rc/drho.*xai;
debug> who
Variables in the current scope:

Ctop  DD    FF2   KK    LL    T     ans   beta  f2    k     phi   rc    xai
D     E     G     KK4   RR    Te    argn  drho  g     lbd   r     v     xver

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 101
101:   Cbot2=1+FF2.*rc^2./drho^2.*phi.^2-2*RR.*sqrt(FF2)*rc/drho.*phi;
debug> who
Variables in the current scope:

Cbot1  D      E      G      KK4    RR     Te     argn   drho   g      lbd    r      v      xver
Ctop   DD     FF2    KK     LL     T      ans    beta   f2     k      phi    rc     xai

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 103
103:   G2b=Ctop./Cbot1./Cbot2;
debug> who
Variables in the current scope:

Cbot1  Ctop   DD     FF2    KK     LL     T      ans    beta   f2     k      phi    rc     xai
Cbot2  D      E      G      KK4    RR     Te     argn   drho   g      lbd    r      v      xver

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 105
105:   if xver==1
debug> who
Variables in the current scope:

Cbot1  Ctop   DD     FF2    G2b    KK4    RR     Te     argn   drho   g      lbd    r      v      xver
Cbot2  D      E      G      KK     LL     T      ans    beta   f2     k      phi    rc     xai

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 105
105:   if xver==1
debug> who
Variables in the current scope:

Cbot1  Ctop   DD     FF2    G2b    KK4    RR     Te     argn   drho   g      lbd    r      v      xver
Cbot2  D      E      G      KK     LL     T      ans    beta   f2     k      phi    rc     xai

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 116
116:   Zb=-2*pi*G*rc*exp(-KK*T).*...
debug> who
Variables in the current scope:

Cbot1  Ctop   DD     FF2    G2b    KK4    RR     Te     argn   drho   g      lbd    r      v      xver
Cbot2  D      E      G      KK     LL     T      ans    beta   f2     k      phi    rc     xai

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 121
121:   Zb=Zb/1e-5;
debug> who
Variables in the current scope:

Cbot1  Ctop   DD     FF2    G2b    KK4    RR     Te     ans    beta   f2     k      phi    rc     xai
Cbot2  D      E      G      KK     LL     T      Zb     argn   drho   g      lbd    r      v      xver

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 124
124:   Zf=Zb+2*pi*G*rc/1e-5;
debug> who
Variables in the current scope:

Cbot1  Ctop   DD     FF2    G2b    KK4    RR     Te     ans    beta   f2     k      phi    rc     xai
Cbot2  D      E      G      KK     LL     T      Zb     argn   drho   g      lbd    r      v      xver

debug> whos
Variables in the current scope:

   Attr Name        Size                     Bytes  Class
   ==== ====        ====                     =====  =====
        Cbot1       5x1000                   40000  double
        Cbot2       5x1000                   40000  double
        Ctop        5x1000                   40000  double
        D           1x1                          8  double
        DD          1x1                          8  double
        E           1x1                          8  double
        FF2         5x1000                   40000  double
        G           1x1                          8  double
        G2b         5x1000                   40000  double
        KK          5x1000                   40000  double
        KK4         5x1000                   40000  double
        LL          5x1000                   40000  double
        RR          1x1                          8  double
    f   T           1x1                          8  double
    f   Te          1x1                          8  double
        Zb          5x1000                   40000  double
        ans         1x1                          8  double
  a     argn        4x6                         24  char
        beta        5x1000                   40000  double
    f   drho        1x1                          8  double
    f   f2          1x5                         40  double
    f   g           1x1                          8  double
        k           1x1000                    8000  double
    f   lbd         1x1000                    8000  double
        phi         5x1000                   40000  double
    f   r           1x1                          8  double
    f   rc          1x1                          8  double
        v           1x1                          8  double
        xai         5x1000                   40000  double
    f   xver        1x1                          8  double

Total is 62043 elements using 496176 bytes

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 126
126:   l=lbd;
debug> who
Variables in the current scope:

Cbot1  Ctop   DD     FF2    G2b    KK4    RR     Te     Zf     argn   drho   g      lbd    r      v      xver
Cbot2  D      E      G      KK     LL     T      Zb     ans    beta   f2     k      phi    rc     xai

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 129
129:   varns={G2b,k,l,Zb,Zf};
debug> who
Variables in the current scope:

Cbot1  Ctop   DD     FF2    G2b    KK4    RR     Te     Zf     argn   drho   g      l      phi    rc     xai
Cbot2  D      E      G      KK     LL     T      Zb     ans    beta   f2     k      lbd    r      v      xver

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 130
130:   varargout=varns(1:nargout);
debug> who
Variables in the current scope:

Cbot1  Ctop   DD     FF2    G2b    KK4    RR     Te     Zf     argn   drho   g      l      phi    rc     varns  xver
Cbot2  D      E      G      KK     LL     T      Zb     ans    beta   f2     k      lbd    r      v      xai

debug> dbstep
stopped in c:\matlab\functions\simmons_matlab\forsyth.m at line 197
197:
debug> who
Variables in the current scope:

Cbot1      DD         G2b        RR         Zf         drho       l          rc         xai
Cbot2      E          KK         T          ans        f2         lbd        v          xver
Ctop       FF2        KK4        Te         argn       g          phi        varargout
D          G          LL         Zb         beta       k          r          varns

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\semilogx.m at line 38
38:   [h, varargin, nargs] = __plt_get_axis_arg__ ("semilogx", varargin{:});
debug> who
Variables in the current scope:

argn      varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\semilogx.m at line 40
40:   if (nargs < 1)  --------> draws a blank figure window here scaled 0-1 in x,y
debug>
debug> who
Variables in the current scope:

argn      h         nargs     varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\semilogx.m at line 44
44:   oldh = gca ();
debug> who
Variables in the current scope:

argn      h         nargs     varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\semilogx.m at line 46
46:     axes (h);
debug> who
Variables in the current scope:

argn      h         nargs     oldh      varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\semilogx.m at line 47
47:     newplot ();
debug> who
Variables in the current scope:

argn      h         nargs     oldh      varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\semilogx.m at line 49
49:     set (h, "xscale", "log");
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\semilogx.m at line 50
50:     if (any( strcmp (get (gca, "nextplot"), {"new", "replace"}))) ---------------> changes xaxis scaling 1e^-1 to1e^+0
debug>
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\semilogx.m at line 51
51:       set (h, "xminortick", "on");
debug> who
Variables in the current scope:

argn      h         nargs     oldh      varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\semilogx.m at line 54
54:     tmp = __plt__ ("semilogx", h, varargin{:});
debug> who
Variables in the current scope:

argn      h         nargs     oldh      varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\semilogx.m at line 56
56:     if (nargout > 0) ---------------> draws plot with semilog X (1e+1 to 1e+4) and linear y and 5 curves
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\semilogx.m at line 57
57:       retval = tmp;
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\semilogx.m at line 60
60:     axes (oldh);
debug> who
Variables in the current scope:

argn      h         nargs     oldh      retval    tmp       varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\semilogx.m at line 63
63: endfunction
debug> who
Variables in the current scope:

argn      h         nargs     oldh      retval    tmp       varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\gca.m at line 42
42:   if (nargin == 0)
debug> who
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\gca.m at line 42
42:   if (nargin == 0)
debug> who
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\gca.m at line 43
43:     h = get (gcf (), "currentaxes");
debug> who
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\gca.m at line 44
44:     if (isempty (h))
debug> who
Variables in the current scope:

h

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\gca.m at line 44
44:     if (isempty (h))
debug> who
Variables in the current scope:

h

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\gca.m at line 51
51: endfunction
debug> who
Variables in the current scope:

h

debug> dbstep
warning: set: allowing rev to match xdir value reverse
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\ylim.m at line 41
41:   ret = __axes_limits__ ("ylim", varargin{:});    xaxis limits switched left to right to give max 1e4 on the left
debug>
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\ylim.m at line 43
43:   if (! isempty (ret))
debug> who
Variables in the current scope:

argn      ret       varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\ylim.m at line 43
43:   if (! isempty (ret))
debug> who
Variables in the current scope:

argn      ret       varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\ylim.m at line 46
46: endfunction
debug> who
Variables in the current scope:

argn      ret       varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\xlim.m at line 45
45:   ret = __axes_limits__ ("xlim", varargin{:});
debug> who
Variables in the current scope:

argn      varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\xlim.m at line 47
47:   if (! isempty (ret))
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\xlim.m at line 47
47:   if (! isempty (ret))
debug> who
Variables in the current scope:

argn      ret       varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\xlim.m at line 50
50: endfunction
debug> who
Variables in the current scope:

argn      ret       varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\xlabel.m at line 36
36:   [h, varargin, nargin] = __plt_get_axis_arg__ ("xlabel", varargin{:});
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\xlabel.m at line 38
38:   if (rem (nargin, 2) != 1)
debug> who
Variables in the current scope:

argn      h         nargin    varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\xlabel.m at line 38
38:   if (rem (nargin, 2) != 1)
debug> who
Variables in the current scope:

argn      h         nargin    varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\xlabel.m at line 42
42:   tmp = __axis_label__ (h, "xlabel", varargin{:},
debug> who
Variables in the current scope:

argn      h         nargin    varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\xlabel.m at line 45
45:   if (nargout > 0)  ------------------> labels the xaxis
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\xlabel.m at line 45
45:   if (nargout > 0)
debug> who
Variables in the current scope:

argn      h         nargin    tmp       varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\xlabel.m at line 45
45:   if (nargout > 0)
debug> who
Variables in the current scope:

argn      h         nargin    tmp       varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\xlabel.m at line 46
46:     retval = tmp;
debug> who
Variables in the current scope:

argn      h         nargin    tmp       varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\xlabel.m at line 49
49: endfunction
debug> who
Variables in the current scope:

argn      h         nargin    retval    tmp       varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\ylabel.m at line 30
30:   [h, varargin, nargin] = __plt_get_axis_arg__ ("ylabel", varargin{:});
debug> who
Variables in the current scope:

argn      varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\ylabel.m at line 32
32:   if (rem (nargin, 2) != 1)
debug> who
Variables in the current scope:

argn      h         nargin    varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\ylabel.m at line 32
32:   if (rem (nargin, 2) != 1)
debug> who
Variables in the current scope:

argn      h         nargin    varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\ylabel.m at line 36
36:   tmp = __axis_label__ (h, "ylabel", varargin{:},
debug> who
Variables in the current scope:

argn      h         nargin    varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\ylabel.m at line 39
39:   if (nargout > 0)  ------------------> labels the yaxis
debug> who
Variables in the current scope:

argn      h         nargin    tmp       varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\ylabel.m at line 39
39:   if (nargout > 0)
debug> who
Variables in the current scope:

argn      h         nargin    tmp       varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\ylabel.m at line 40
40:     retval = tmp;
debug> who
Variables in the current scope:

argn      h         nargin    tmp       varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\ylabel.m at line 43
43: endfunction
debug> who
Variables in the current scope:

argn      h         nargin    retval    tmp       varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 41
41:   [ax, varargin, nargs] = __plt_get_axis_arg__ ("grid", varargin{:});
debug> who
Variables in the current scope:

argn      varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 43
43:   grid_on = (strcmp (get (ax, "xgrid"), "on")
debug> who
Variables in the current scope:

argn      ax        nargs     varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 47
47:   minor_on = (strcmp (get (ax, "xminorgrid"), "on")
debug> who
Variables in the current scope:

argn      ax        grid_on   nargs     varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 51
51:   if (nargs > 2)
debug> who
Variables in the current scope:

argn      ax        grid_on   minor_on  nargs     varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 51
51:   if (nargs > 2)
debug> who
Variables in the current scope:

argn      ax        grid_on   minor_on  nargs     varargin

debug> dbtep
error: 'dbtep' undefined near line 2 column 1
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 51
51:   if (nargs > 2)
debug> who
Variables in the current scope:

argn      ax        grid_on   minor_on  nargs     varargin

debug> who
Variables in the current scope:

argn      ax        grid_on   minor_on  nargs     varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 56
56:     x = varargin{1};
debug> who
Variables in the current scope:

argn      ax        grid_on   minor_on  nargs     varargin

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 57
57:     if (ischar (x))
debug> who
Variables in the current scope:

argn      ax        grid_on   minor_on  nargs     varargin  x

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 57
57:     if (ischar (x))
debug> who
Variables in the current scope:

argn      ax        grid_on   minor_on  nargs     varargin  x

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 58
58:       if (strcmpi (x, "off"))
debug> who
Variables in the current scope:

argn      ax        grid_on   minor_on  nargs     varargin  x

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 58
58:       if (strcmpi (x, "off"))
debug> who
Variables in the current scope:

argn      ax        grid_on   minor_on  nargs     varargin  x

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 58
58:       if (strcmpi (x, "off"))
debug> who
Variables in the current scope:

argn      ax        grid_on   minor_on  nargs     varargin  x

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 61
61:         grid_on = true;
debug> who
Variables in the current scope:

argn      ax        grid_on   minor_on  nargs     varargin  x

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 87
87:   if (grid_on)
debug> who
Variables in the current scope:

argn      ax        grid_on   minor_on  nargs     varargin  x

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 87
87:   if (grid_on)
debug> who
Variables in the current scope:

argn      ax        grid_on   minor_on  nargs     varargin  x

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 88
88:     set (ax, "xgrid", "on", "ygrid", "on", "zgrid", "on");
debug> who
Variables in the current scope:

argn      ax        grid_on   minor_on  nargs     varargin  x

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 89
89:     if (minor_on)
debug> who
Variables in the current scope:

argn      ax        grid_on   minor_on  nargs     varargin  x

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 89
89:     if (minor_on)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 92
92:       set (ax, "xminorgrid", "off", "yminorgrid", "off", "zminorgrid", "off");
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\grid.m at line 99
99: endfunction ----------------> grid drawn on figure x and y
debug> dbstep
warning: get: allowing linew to match line property linewidth --------------> bold curve lines
warning: get: allowing linew to match line property linewidth
warning: get: allowing linew to match line property linewidth
warning: get: allowing linew to match line property linewidth
warning: get: allowing linew to match line property linewidth
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 54
54:   if (nargin > 0 && numel (varargin{1}) == 1 && ishandle (varargin{1})
debug>
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 54
54:   if (nargin > 0 && numel (varargin{1}) == 1 && ishandle (varargin{1})
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 61
61:     ax = gca ();
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 62
62:     fig = gcf ();
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 63
63:     nargs = numel (varargin);
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 66
66:   hold_all = false;
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 67
67:   if (nargs == 0)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 67
67:   if (nargs == 0)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 67
67:   if (nargs == 0)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 70
70:     state = varargin{1};
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 71
71:     if (ischar (state))
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 71
71:     if (ischar (state))
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 72
72:       if (strcmpi (state, "off"))
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 72
72:       if (strcmpi (state, "off"))
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 72
72:       if (strcmpi (state, "off"))
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 72
72:       if (strcmpi (state, "off"))
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 78
78:         turn_hold_off = false;
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 87
87:   if (turn_hold_off)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 87
87:   if (turn_hold_off)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 90
90:     set (ax, "nextplot", "add");
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 91
91:     set (fig, "nextplot", "add");
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 93
93:   set (ax, "__hold_all__", hold_all);
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\hold.m at line 95
95: endfunction
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\gca.m at line 42
42:   if (nargin == 0)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\gca.m at line 42
42:   if (nargin == 0)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\gca.m at line 43
43:     h = get (gcf (), "currentaxes");
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\gca.m at line 44
44:     if (isempty (h))
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\gca.m at line 44
44:     if (isempty (h))
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\gca.m at line 51
51: endfunction
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\axes.m at line 30
30:   if (nargin == 0 || nargin > 1)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\axes.m at line 30
30:   if (nargin == 0 || nargin > 1)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\axes.m at line 32
32:     idx = find (strcmpi (varargin(1:2:end), "parent"), 1, "first");
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\axes.m at line 33
33:     if (! isempty (idx) && length (varargin) >= 2*idx)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\axes.m at line 33
33:     if (! isempty (idx) && length (varargin) >= 2*idx)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\axes.m at line 37
37:       cf = gcf ();
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\axes.m at line 39
39:     tmp = __go_axes__ (cf, varargin{:});
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\axes.m at line 40
40:     if (__is_handle_visible__ (tmp)) --------------> clears the figure ALL BLANK here
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\axes.m at line 41
41:       set (ancestor (cf, "figure"), "currentaxes", tmp);
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\axes.m at line 58
58:   if (nargout > 0)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\axes.m at line 58
58:   if (nargout > 0)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\axes.m at line 59
59:     h = tmp;
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\plot\axes.m at line 62
62: endfunction
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\help\unimplemented.m at line 27
27:   is_matlab_function = true;
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\help\unimplemented.m at line 30
30:   switch (fcn)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\help\unimplemented.m at line 30
30:   switch (fcn)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\help\unimplemented.m at line 30
30:   switch (fcn)
debug> who
Variables in the current scope:

fcn                 is_matlab_function

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\help\unimplemented.m at line 30
30:   switch (fcn)
debug> who
Variables in the current scope:

fcn                 is_matlab_function

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\help\unimplemented.m at line 30
30:   switch (fcn)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\help\unimplemented.m at line 30
30:   switch (fcn)
debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\help\unimplemented.m at line 57
57:     if (ismember (fcn, missing_functions ()))
debug> who
Variables in the current scope:

fcn                 is_matlab_function

debug> dbstep
stopped in C:\Octave\Octave3.6.4_gcc4.6.2\share\octave\3.6.4\m\help\unimplemented.m at line 57
57:     if (ismember (fcn, missing_functions ()))
debug>

Once the figure was cleared it never came back - had to stop at this point.

Looks like a call to axes.m wipes it.

Sorry for the very long mail but took a while to find where it blanks out!

Let me know if this is meaningful.
Cheers
Lester


On 5 August 2013 00:41, Ben Abbott <address@hidden> wrote:
On Aug 4, 2013, at 12:29 PM, Lester Anderson wrote:

> Hi,
>
> I found the defval.m code. From the functions I uploaded, it looks like all it was doing was defining the values of the constant variables, now wehether that is the issue why it was not running correct I do not know. But maybe the solution is more simple.
>
> In the xample I uploaded, all the key things are available, so there are no other toolbox things required. Defval is clearly a Matlab intrinsic function. As I mentioned before, values for 'g' and 'G' can be declared normally without the need for the defval function.
>
> If you can figure out how the plot was created, from those two functions, that would be a big help. If I ignored the defval in the code and replaced it as R=... or G=... etc, it sort of runs but not to generate the output as pected from the image result. Not having worked with Matlab before, I am still learning the procedures with Octave; I have come from a Fortran 95 background. Since a lot of geophysics code is in Matlab, I thought try Octave and see how it copes.
>
> Hopefully I have not gone into the deep end !
>
> Cheers
> Lester

I think you got the wrong defval.m.  The one your looking for is available from the same page you gave a link to earlier.

        http://geoweb.princeton.edu/people/simons/DOTM/defval.m

Ben



reply via email to

[Prev in Thread] Current Thread [Next in Thread]