Changeset 1204 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

Mar 7, 2014, 4:53:30 PM (11 years ago)

Author:

msylvestre

Message:

rings.F90 : calculation of the shadow under the Cassini Division was omitted in the previous version

File:

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

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

@@ -18 +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=4 ! number of subdivisions of the A ring
-    INTEGER, PARAMETER :: nb_b=3 ! number of subdivisions of the B ring
-    INTEGER, PARAMETER :: nb_c=3 ! number of subdivisions of the C 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 = 3 ! number of subdivisions of the C ring
     INTEGER, PARAMETER :: nb_ca = 2 ! number of subdivisions in the Cassini division
     INTEGER :: i
@@ -109 +109 @@
                 
     tau_B(1) = 1.29
-    tau_B(2) = 5.13  
-    tau_B(3) = 2.84
+    tau_B(2) = 5.13 
+    tau_B(3) = 2.84 
     
     tau_C(1) = 0.06
@@ -176 +176 @@
     enddo
 
+    do i=1, nb_ca
+        where(r(:) .ge. Ca_Rint(i) .and. r(:) .le. Ca_Rext(i) .and. eclipse(:) .ne. 1000.)
+            eclipse(:) = 1. - exp(-tau_Ca(i)/abs(sin(declin)))
+        end where
+    enddo
+
     ! At the other places and the excluded ones, eclipse is 0. 
     where(eclipse(:) .eq. 2000. .or. eclipse(:) .eq. 1000.)