FilesExpand file tree

 main

/

fargo2python.py

Copy path

Blame

More file actions

Blame

More file actions

Latest commit

 

History

History

History

126 lines (112 loc) · 4.11 KB

 main

/

fargo2python.py

Copy path

Top

File metadata and controls

• Code
• Blame

126 lines (112 loc) · 4.11 KB

Raw

Copy raw file

Download raw file

Open symbols panel

Edit and raw actions

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

42

43

44

45

46

47

48

49

50

51

52

53

54

55

56

57

58

59

60

61

62

63

64

65

66

67

68

69

70

71

72

73

74

75

76

77

78

79

80

81

82

83

84

85

86

87

88

89

90

91

92

93

94

95

96

97

98

99

100

101

102

103

104

105

106

107

108

109

110

111

112

113

114

115

116

117

118

119

120

121

122

123

124

125

126

# ===================================================================

# FARGO to PYTHON

# code written by Clement Baruteau (CB), with contributions from

# Pablo Benitez-Llambay and Gaylor Wafflard-Fernandez

# ===================================================================

#

# program can run with either Python 2.X or Python 3.X.

#

# NB: logarithmic color scale requires matplotblib version > 3.0 to

# display color bar properly

#

# ===================================================================

# =========================================

# TO DO LIST

# =========================================

# - think of a better handling of 3D case? Right now I can only plot the primitive fields in a

# vertical plane, not the fields built upon primitive ones (like the vorticity etc.).

# - show RH around planet

# - axifield (minor): add 1D visc evolving profiles, add planet position via dashed line

# =========================================

# Plotting parameters

import matplotlib

matplotlib.rcParams.update({'font.size': 24})

matplotlib.rc('font', family='Arial')

# =====================

# 1. IMPORT GLOBAL VARIABLES

# =====================

import par

if par.allfluids == 'Yes':

matplotlib.rcParams.update({'font.size': 20})

# =====================

# 2. DISPLAY 2D FIELDS

# =====================

if par.onedprofile == 'No' and par.plot_field == True:

from plot_twodfield import *

plottwodfield()

# =====================

# 3. DISPLAY AZIMUTHALLY-AVERAGED RADIAL 1D PROFILES

# =====================

if par.onedprofile != 'No' and par.plot_field == True:

if par.onedspacetime == 'Yes':

from plot_spacetimediagram import *

plotspacetimediagram()

else:

from plot_onedprofile import *

plotonedprofile()

# =====================

# 4. DISPLAY TIME EVOLUTION OF DISC TORQUE OR POWER ON PLANET(s)

# =====================

if par.plot_tqwk != 'No':

from plot_tqwk import *

plottqwk()

# =====================

# 5. DISPLAY TIME EVOLUTION OF PLANET ORBITAL ELEMENTS

# =====================

if par.plot_planet != 'No':

from plot_planet import *

plotplanet()

# =====================

# 6. DISPLAY STUFF FOR DUST PARTICLES (only for Dusty-FARGO ADSG)

# =====================

if par.plot_dust != 'No' and par.fargo3d == 'No':

from plot_dust import *

plotdust()

# =====================

# 7. DISC'S CENTER OF MASS

# =====================

if ( ('plot_disccom' in open('paramsf2p.dat').read()) and (par.plot_disccom != 'No') ):

from plot_disccom import *

plotdisccom()

# =====================

# 8. DISC MASS (time)

# =====================

if ( ('plot_discmass' in open('paramsf2p.dat').read()) and (par.plot_discmass != 'No') ):

from plot_discmass import *

plotdiscmass()

# =====================

# 9. ratio of orbit-crossing and librating inverse vortensities (time)

# =====================

if ( ('plot_libcross' in open('paramsf2p.dat').read()) and (par.plot_libcross != 'No') ):

from plot_libcross import *

plotlibcross()

# =====================

# 10. disc turbulent properties

# =====================

if ( ('plot_turb' in open('paramsf2p.dat').read()) and (par.plot_turb != 'No') ):

from plot_turb import *

if par.plot_turb == 'power_spectrum':

plotpowerspectrum()

if par.plot_turb == 'auto_correlation':

plotautocorrelationtimescale()

if par.plot_turb == 'alphas':

plot_alphas()

if par.plot_turb == 'histo_field':

plot_histofield()

if par.plot_turb == 'histo_torque':

plot_histotorque()

if par.plot_turb == 'time_alphareynolds':

plot_time_alphaRey()

if par.plot_turb == 'time_alphagrav':

plot_time_alphaGrav()

if par.plot_turb == 'time_alphagrav_ext':

plot_time_alphaGrav_external()

# =====================

# 11. DISC'S ECCENTRICITY OR PERICENTRE ARGUMENT

# =====================

if ( ('plot_discecc' in open('paramsf2p.dat').read()) and (par.plot_discecc != 'No') ):

from plot_discbinary import *

plotdiscecc()

if ( ('plot_discperarg' in open('paramsf2p.dat').read()) and (par.plot_discperarg != 'No') ):

from plot_discbinary import *

plotdiscperarg()