Changeset 1177 for trunk

Timestamp:

Feb 14, 2014, 9:04:37 AM (11 years ago)

Author:

sglmd

Message:

Corrected a mistake in the geometry of the ring's transmission

File:

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

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

@@ -125 +125 @@
     do i=1, nb_A 
         where(r(:) .ge. A_Rint(i) .and. r(:) .le. A_Rext(i) .and. eclipse(:) .ne. 1000.)
-            eclipse(:) = 1. - exp(-tau_A(i)/cos(declin))
+            eclipse(:) = 1. - exp(-tau_A(i)/abs(sin(declin)))
         end where
     end do 
@@ -131 +131 @@
     do i=1, nb_B 
         where(r(:) .ge. B_Rint(i) .and. r(:) .le. B_Rext(i) .and. eclipse(:) .ne. 1000.)
-            eclipse(:) = 1. - exp(-tau_B(i)/cos(declin))
+            eclipse(:) = 1. - exp(-tau_B(i)/abs(sin(declin)))
         end where
     enddo
@@ -137 +137 @@
     do i=1, nb_C 
         where(r(:) .ge. C_Rint(i) .and. r(:) .le. C_Rext(i) .and. eclipse(:) .ne. 1000.)
-            eclipse(:) = 1. - exp(-tau_C(i)/cos(declin))
+            eclipse(:) = 1. - exp(-tau_C(i)/abs(sin(declin)))
         end where
     enddo