Changeset 762 for trunk/LMDZ.MARS/libf/aeronomars

Timestamp:

Aug 27, 2012, 10:19:09 AM (12 years ago)

Author:

emillour

Message:

Mars GCM:
Bug correction in the varying EUV case; now correctly account for things
as the year cycles.
FGG

File:

• trunk/LMDZ.MARS/libf/aeronomars/jthermcalc_e107.F (modified) (2 diffs)

trunk/LMDZ.MARS/libf/aeronomars/jthermcalc_e107.F

@@ -139 +139 @@
       if(realday.lt.date_e107(1)) then
          e107=e107_tab(1)
+      else if(realday.ge.date_e107(669)) then
+         e107=e107_tab(669)   
       else if(realday.ge.date_e107(1).and.
-     $        realday.le.date_e107(669)) then
-         do i=1,669
+     $        realday.lt.date_e107(669)) then
+         do i=1,668
             if(realday.ge.date_e107(i).and.
      $           realday.lt.date_e107(i+1)) then
                tinf=i
                tsup=i+1
-               e107=e107_tab(tinf)+(zday-date_e107(tinf))*
+               e107=e107_tab(tinf)+(realday-date_e107(tinf))*
      $              (e107_tab(tsup)-e107_tab(tinf))
             endif
          enddo
       endif
+
       !Photoabsorption coefficients at TOA as a function of E10.7
       do j=1,nabs
@@ -159 +162 @@
          enddo
       enddo
-
 ! Interpolation to the tabulated photoabsorption coefficients for each species
 ! in each spectral interval