Skip to content
Snippets Groups Projects
Commit cb880976 authored by SERPE User's avatar SERPE User
Browse files

Adding coeff files for exoplanetary magnetic field models

parent d328b579
No related branches found
No related tags found
No related merge requests found
Pipeline #28825 failed
pro COEFF_ADLeonis_300G,g,h,n,dip,crt,tdip
;#coefficient modele ADLeonis_300G
;#g et h:coeff. of Legendre polynomials
;#n:ordre du developpement
;# n: development order
;# dip:position of the du dipole wrt the rotation axis
;# crt: current sheet model
n=3
G=fltarr(4,4) & H=fltarr(4,4)
G[1,0]=-206.506
G[1,1]=77.9852 & H(1,1)=203.159
G[2,0]=-0.00000
G[2,1]=-0.00000 & H[2,1]=0.00000
G[2,2]=0.0 & H[2,2]=0.0
G[3,0]=-0.00000
G[3,1]=0.00000 & H[3,1]=0.00000
G[3,2]=0.0 & H[3,2]=0.0
G[3,3]=0.0 & H[3,3]=0.0
;#dip = [polar_angle, 360-azimuth_angle]
dip=[46.5000,111.000]
;# tdip :
;# t3d,/reset,rotate=[0, 0, 360-azimuth_angle]
;# -> rotates 360-azimuth_angle degrees counterclockwise about the Z axis.
;# T3D,rotate=[0,-polar_angle,0]
;# -> rotate polar_angle degs clockwise about the Y axis
;# tdip=transpose(!p.t)
;# Thix matrix allows to work in jovimagnetic coordinates
;# Matrix is:
;# [cosa.cosb -sina cosa.sinb]
;# [sina.cosb cosa sina.sinb]
;# [ -sinb 0 cosb]
tdip=[[-0.24668422,-0.64263436,-0.72537437],[0.93358043,-0.35836795,0.0000000],[-0.25995093,-0.67719531,0.68835458]]
;# Current sheet model
crt=[0, 0, 0, 0, 0]
return
end
pro COEFF_ADLeonis_460G,g,h,n,dip,crt,tdip
;#coefficient modele ADLeonis_460G
;#g et h:coeff. of Legendre polynomials
;#n:ordre du developpement
;# n: development order
;# dip:position of the du dipole wrt the rotation axis
;# crt: current sheet model
n=3
G=fltarr(4,4) & H=fltarr(4,4)
G[1,0]=-236.918
G[1,1]=166.637 & H(1,1)=-357.354
G[2,0]=-0.00000
G[2,1]=-0.00000 & H[2,1]=-0.00000
G[2,2]=0.0 & H[2,2]=0.0
G[3,0]=-0.00000
G[3,1]=0.00000 & H[3,1]=-0.00000
G[3,2]=0.0 & H[3,2]=0.0
G[3,3]=0.0 & H[3,3]=0.0
;#dip = [polar_angle, 360-azimuth_angle]
dip=[59.0000,245.000]
;# tdip :
;# t3d,/reset,rotate=[0, 0, 360-azimuth_angle]
;# -> rotates 360-azimuth_angle degrees counterclockwise about the Z axis.
;# T3D,rotate=[0,-polar_angle,0]
;# -> rotate polar_angle degs clockwise about the Y axis
;# tdip=transpose(!p.t)
;# Thix matrix allows to work in jovimagnetic coordinates
;# Matrix is:
;# [cosa.cosb -sina cosa.sinb]
;# [sina.cosb cosa sina.sinb]
;# [ -sinb 0 cosb]
tdip=[[-0.21766450,0.46678302,-0.85716730],[-0.90630779,-0.42261826,0.0000000],[-0.36225455,0.77685740,0.51503807]]
;# Current sheet model
crt=[0, 0, 0, 0, 0]
return
end
pro COEFF_tau_boo_A_axi,g,h,n,dip,crt,tdip
;#coefficient modele tau_boo_A_axi
;#g et h:coeff. of Legendre polynomials
;#n:ordre du developpement
;# n: development order
;# dip:position of the du dipole wrt the rotation axis
;# crt: current sheet model
;# Parameters of the model"
;# dipole_moment: 3.8929639
;# polar_angle: 0
;# azimuth_angle: 0
;# offset: 0.0
;# body_radius: 9.2597664e+08
;# body_flattening: 0
n=3
G=fltarr(4,4) & H=fltarr(4,4)
G[1,0]=3.89296
G[1,1]=0.00000 & H(1,1)=0.00000
G[2,0]=0.00000
G[2,1]=0.00000 & H[2,1]=0.00000
G[2,2]=0.0 & H[2,2]=0.0
G[3,0]=0.00000
G[3,1]=0.00000 & H[3,1]=0.00000
G[3,2]=0.0 & H[3,2]=0.0
G[3,3]=0.0 & H[3,3]=0.0
;#dip = [polar_angle, 360-azimuth_angle]
dip=[0,360.000]
;# tdip :
;# t3d,/reset,rotate=[0, 0, 360-azimuth_angle]
;# -> rotates 360-azimuth_angle degrees counterclockwise about the Z axis.
;# T3D,rotate=[0,-polar_angle,0]
;# -> rotate polar_angle degs clockwise about the Y axis
;# tdip=transpose(!p.t)
;# Thix matrix allows to work in jovimagnetic coordinates
;# Matrix is:
;# [cosa.cosb -sina cosa.sinb]
;# [sina.cosb cosa sina.sinb]
;# [ -sinb 0 cosb]
tdip=[[1.0000000,2.4492936e-16,0.0000000],[-2.4492936e-16,1.0000000,0.0000000],[0.0000000,0.0000000,1.0000000]]
;# Current sheet model
crt=[0, 0, 0, 0, 0]
return
end
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment