summaryrefslogtreecommitdiffstats
path: root/plot_rho_rz.py
diff options
context:
space:
mode:
authorVolker Hoffmann <volker@cheleb.net>2013-08-21 15:07:26 +0200
committerVolker Hoffmann <volker@cheleb.net>2013-08-21 15:07:26 +0200
commite0f29d6c30cb1b39b3a3f6276b109c570921d732 (patch)
tree88d7dfe9b596e3a20908351456a6cb8ed4bda9d5 /plot_rho_rz.py
parentb381060940d4dfd7387516f8da756d2d765876a8 (diff)
read,save,load and use rho_rz>0 for clim (so the log10 of a potential negative density does not break the plots) -- @todo: some sort of notification if rho<0 [clearly unphysical!]
Diffstat (limited to 'plot_rho_rz.py')
-rw-r--r--plot_rho_rz.py6
1 files changed, 3 insertions, 3 deletions
diff --git a/plot_rho_rz.py b/plot_rho_rz.py
index 046c544..38b9629 100644
--- a/plot_rho_rz.py
+++ b/plot_rho_rz.py
@@ -27,11 +27,11 @@ for iout in iouts:
disk = Disk(iout)
disk.load_npz_minmax()
if first:
- rho_rz_lo = disk.rho_rz_min; rho_rz_hi = disk.rho_rz_max
+ rho_rz_lo = disk.rho_rz_min_gt0; rho_rz_hi = disk.rho_rz_max_gt0
first = False
else:
- if disk.rho_rz_min < rho_rz_lo: rho_rz_lo = disk.rho_rz_min
- if disk.rho_rz_max > rho_rz_hi: rho_rz_hi = disk.rho_rz_max
+ if disk.rho_rz_min_gt0 < rho_rz_lo: rho_rz_lo = disk.rho_rz_min_gt0
+ if disk.rho_rz_max_gt0 > rho_rz_hi: rho_rz_hi = disk.rho_rz_max_gt0
rho_rz_lim = [np.log10(rho_rz_lo), np.log10(rho_rz_hi)]