/[escript]/branches/symbolic_from_3470/ripley/test/python/gravity_all.py
ViewVC logotype

Diff of /branches/symbolic_from_3470/ripley/test/python/gravity_all.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3876 by gross, Fri Mar 16 01:12:26 2012 UTC revision 3877 by gross, Mon Mar 19 00:52:00 2012 UTC
# Line 85  beta=1/1000000. Line 85  beta=1/1000000.
85  #  #
86  #  where do we know the gravity:  #  where do we know the gravity:
87  #  #
 chi=1.  
88  x=Function(domain).getX()  x=Function(domain).getX()
89  dz=H/NE_H  dz=H/NE_H
90  H_earth=int(H_earth/dz)*dz # lock to grid  H_earth=int(H_earth/dz)*dz # lock to grid
91  chi=whereNegative(abs(x[DIM-1]-(H_earth+dz/2))-dz/2)  chi=whereNegative(abs(x[DIM-1]-(H_earth+dz/2))-dz/2)
 chi=1  
92  #  #
93  #   normalize  g_hat (data):  #   normalize  g_hat (data):
94  #  #
# Line 161  print "psi =",Lsup(psi_v-psi_l)/Lsup(psi Line 159  print "psi =",Lsup(psi_v-psi_l)/Lsup(psi
159  print "Differences to Data :"  print "Differences to Data :"
160  print "rho =",Lsup(rho_ref-rho_l*g0/LL**2)/Lsup(rho_ref)  print "rho =",Lsup(rho_ref-rho_l*g0/LL**2)/Lsup(rho_ref)
161  print "psi =",Lsup(psi_ref-psi_l*g0)/Lsup(psi_ref)  print "psi =",Lsup(psi_ref-psi_l*g0)/Lsup(psi_ref)
162  print "g =",Lsup(grad(psi_ref-psi_l*g0)[DIM-1])/Lsup(grad(psi_ref)[DIM-1])  print "g =",Lsup(chi*grad(psi_ref-psi_l*g0)[DIM-1])/Lsup(chi*grad(psi_ref)[DIM-1])
163    
164    
165  #saveVTK("u.vtu", rho_ref=rho_ref, psi_ref=psi_ref, g=grad(psi_ref)[DIM-1], rho=rho_v, psi=psi_v, chi=chi)  #saveVTK("u.vtu", rho_ref=rho_ref, psi_ref=psi_ref, g=grad(psi_ref)[DIM-1], rho=rho_v, psi=psi_v, chi=chi)

Legend:
Removed from v.3876  
changed lines
  Added in v.3877

  ViewVC Help
Powered by ViewVC 1.1.26