ttypesetting.mac - pism - [fork] customized build of PISM, the parallel ice sheet model (tillflux branch)
  HTML git clone git://src.adamsgaard.dk/pism
   DIR Log
   DIR Files
   DIR Refs
   DIR LICENSE
       ---
       ttypesetting.mac (1588B)
       ---
            1 texput(E, "E^{n+1}")$
            2 texput(Eold, "E^{n}")$
            3 texput(rho, "\\rho_{i}")$
            4 texput(w, "w^{n}")$
            5 texput(K, "K^{n}")$
            6 texput(R, "R^{n}")$
            7 texput(Phi, "\\Phi^{n}")$
            8 texput(one_over_dz, "\\frac{1}{\\Delta z}")$
            9 texput(half, "\\frac{1}{2}")$
           10 texput(nounify(Up), "\\operatorname{Up}")$
           11 texput(k_s, "k_{s}");
           12 texput(nounify(Delta), "\\Delta")$
           13 texput(nounify(delta_p), "\\mathbf{\\delta}_{+}")$
           14 texput(nounify(delta_m), "\\mathbf{\\delta}_{-}")$
           15 
           16 load("neumann_bc.mac")$
           17 
           18 file : openw("formulas.tex")$
           19 
           20 print_eq(name, equation) := printf(file, "\\newcommand{\\~a}{~a}~%", name, tex1(equation))$
           21 
           22 print_eq(R, rhs(R_eq[k]))$
           23 print_eq(mufactor, rhs(mu_eq))$
           24 
           25 print_eq(discretization, eq)$
           26 
           27 printf(file, "% Values of E on the grid, for the generic equation.~%")$
           28 print_eq(E, E)$
           29 print_eq(El, E[k-1])$
           30 print_eq(Eu, E[k+1])$
           31 
           32 printf(file, "% Matrix entries, w > 0.~%")$
           33 print_eq(Lp, l[p])$
           34 print_eq(Dp, d[p])$
           35 print_eq(Up, u[p])$
           36 print_eq(Bp, b[p])$
           37 
           38 printf(file, "% Matrix entries, w < 0.~%")$
           39 print_eq(Lm, l[m])$
           40 print_eq(Dm, d[m])$
           41 print_eq(Um, u[m])$
           42 print_eq(Bm, b[m])$
           43 
           44 printf(file, "% Neumann B. C. at the base.~%")$
           45 print_eq(neumannb, neumann[0])$
           46 print_eq(Dpb, d_base[p])$
           47 print_eq(Upb, u_base[p])$
           48 print_eq(Bpb, b_base[p])$
           49 print_eq(Dmb, d_base[m])$
           50 print_eq(Umb, u_base[m])$
           51 print_eq(Bmb, b_base[m])$
           52 
           53 printf(file, "% Neumann B. C. at the surface.~%")$
           54 print_eq(neumanns, neumann[k_s])$
           55 print_eq(Lps, l_surface[p])$
           56 print_eq(Dps, d_surface[p])$
           57 print_eq(Bps, b_surface[p])$
           58 print_eq(Lms, l_surface[m])$
           59 print_eq(Dms, d_surface[m])$
           60 print_eq(Bms, b_surface[m])$
           61 
           62 /* Close the file */
           63