Changeset 1200 for trunk

Timestamp:

Mar 4, 2014, 3:34:09 PM (11 years ago)

Author:

msylvestre

Message:

rings.F90 : update of the rings opacities with values from Cassini/UVIS observations

File:

• trunk/LMDZ.GENERIC/libf/phystd/rings.F90 (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/phystd/rings.F90

@@ -1 +1 @@
 SUBROUTINE rings(ngrid, declin, ptime, rad, flat, eclipse)
 ! Calculates Saturn's rings shadowing
+! Includes rings opacities measured by Cassini/UVIS
 ! Authors: M. Sylvestre, M. Capderou, S. Guerlet, A. Spiga
 
@@ -17 +18 @@
     REAL :: rpol   ! polar radius of the planet
     REAL :: e      ! shape excentricity of the planet : (1-e*e) = (1-f)*(1-f)    
-    INTEGER, PARAMETER :: nb_a=1 ! number of subdivisions of the A ring
+    INTEGER, PARAMETER :: nb_a=4 ! number of subdivisions of the A ring
     INTEGER, PARAMETER :: nb_b=3 ! number of subdivisions of the B ring
-    INTEGER, PARAMETER :: nb_c=1 ! number of subdivisions of the C ring
+    INTEGER, PARAMETER :: nb_c=3 ! number of subdivisions of the C ring
+    INTEGER, PARAMETER :: nb_ca = 2 ! number of subdivisions in the Cassini division
     INTEGER :: i
 
@@ -29 +31 @@
     REAL, DIMENSION(nb_c) :: C_Rint ! internal radii of the subdivisions of the C ring
     REAL, DIMENSION(nb_c) :: C_Rext ! external radii of the subdivisions of the C ring 
-                
+    REAL, DIMENSION(nb_ca) :: Ca_Rint ! internal radii of the subdivisions of the Cassini Division
+    REAL, DIMENSION(nb_ca) :: Ca_Rext ! external radii of the subdivisions of the Cassini Division
+
     ! Opacities of the rings : for each one we can give different opacities for each part
     REAL, DIMENSION(nb_a) :: tau_A ! opacity of the A ring
     REAL, DIMENSION(nb_b) :: tau_B ! opacity of the B ring
     REAL, DIMENSION(nb_c) :: tau_C ! opacity of the C ring
-                
+    REAL, DIMENSION(nb_ca) :: tau_Ca ! opacity of the Cassini Division 
+
     ! Parameters used to calculate if a point is under a ring subdivision's shadow
     REAL :: phi_S                             ! subsolar point longitude
@@ -65 +70 @@
     eclipse(:) = 2000.
 
+! Model of the rings with Cassini/UVIS opacities
+
     ! Size of the rings
-    A_Rint(1) = 2.01*rad
-    A_Rext(1) = 2.26*rad
-    B_Rint(1) = 1.55*rad
-    B_Rext(1) = 1.67*rad
-    B_Rint(2) = 1.67*rad
+    A_Rint(1) = 2.03*rad
+    A_Rext(1) = 2.06*rad
+    A_Rint(2) = 2.06*rad
+    A_Rext(2) = 2.09*rad
+    A_Rint(3) = 2.09*rad
+    A_Rext(3) = 2.12*rad
+    A_Rint(4) = 2.12*rad
+    A_Rext(4) = 2.27*rad
+
+    B_Rint(1) = 1.53*rad
+    B_Rext(1) = 1.64*rad
+    B_Rint(2) = 1.64*rad
     B_Rext(2) = 1.83*rad
     B_Rint(3) = 1.83*rad
-    B_Rext(3) = 1.92*rad
-    C_Rint(1) = 1.21*rad
-    C_Rext(1) = 1.53*rad
+    B_Rext(3) = 1.95*rad
     
-    ! Opacities of the rings (TBD: update with most recent values?)
-    tau_A(1) = 0.4
-    tau_B(1) = 0.8
-    tau_B(2) = 2.
-    tau_B(3) = 1.4
-    tau_C(1) = 0.08
+    C_Rint(1) = 1.24*rad
+    C_Rext(1) = 1.29*rad
+    C_Rint(2) = 1.29*rad
+    C_Rext(2) = 1.43*rad
+    C_Rint(3) = 1.43*rad
+    C_Rext(3) = 1.53*rad
+
+    Ca_Rint(1) = 1.95*rad
+    Ca_Rext(1) = 1.99*rad
+    Ca_Rint(2) = 1.99*rad
+    Ca_Rext(2) = 2.03*rad
+
+
+    ! Opacities of the rings
+    tau_A(1) = 1.24
+    tau_A(2) = 0.81
+    tau_A(3) = 0.67
+    tau_A(4) = 0.58
+                
+    tau_B(1) = 1.29
+    tau_B(2) = 5.13  
+    tau_B(3) = 2.84
+    
+    tau_C(1) = 0.06
+    tau_C(2) = 0.10
+    tau_C(3) = 0.14
+
+    tau_Ca(1) = 0.06
+    tau_Ca(2) = 0.24
 
     ! Convert to cartesian coordinates