~fluidity-core/fluidity/sea-ice-branch

« back to all changes in this revision

Viewing changes to examples/stokes_square_convection/Plot_RMS.py

• browse files at revision 3647
• compare with another revision
• download diff
• download tarball
• view history from revision 3647

Show diffs side-by-side

added

removed

examples/stokes_square_convection/Plot_RMS.py

 

1

#!/usr/bin/python

 

2

 

 

3

# This sript plots up the RMS velocity for both the 24x24 and 48x48 case

 

4

import pylab

 

5

from fluidity_tools import stat_parser as stat

 

6

 

 

7

# Stafiles:

 

8

statfile24="stokes-sc-Ra1e5-24.stat"

 

9

statfile48="stokes-sc-Ra1e5-48.stat"

 

10

 

 

11

# First plot 24x24 case:

 

12

pylab.plot(stat(statfile24)["CoordinateMesh"]["nodes"][-1],

 

13

      stat(statfile24)["Fluid"]["Velocity%magnitude"]["l2norm"][-1],

 

14

      linestyle='None', marker='o', markerfacecolor='1.0')

 

15

 

 

16

# Next plot 48x48 case:

 

17

pylab.plot(stat(statfile48)["CoordinateMesh"]["nodes"][-1],

 

18

      stat(statfile48)["Fluid"]["Velocity%magnitude"]["l2norm"][-1],

 

19

      linestyle='None', marker='o', markerfacecolor='1.0')

 

20

 

 

21

# Plot benchmark value as line for comparison:

 

22

pylab.plot([100,8e4],[193.214,193.214],'k--',lw=0.4)

 

23

 

 

24

pylab.xlabel(r"Vertices")

 

25

pylab.ylabel(r"RMS Velocity")

 

26

pylab.xlim(100,5e3)

 

27

pylab.ylim(192.0,195.0)

 

28

 

 

29

pylab.savefig("RMS_1e5.pdf")

 

30

 

• Older »

Loggerhead is a web-based interface for Breezy
Version: 2.0.1