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

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)