Changeset 762 for trunk/LMDZ.MARS/libf


Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/libf/aeronomars/jthermcalc_e107.F

    r705 r762  
    139139      if(realday.lt.date_e107(1)) then
    140140         e107=e107_tab(1)
     141      else if(realday.ge.date_e107(669)) then
     142         e107=e107_tab(669)   
    141143      else if(realday.ge.date_e107(1).and.
    142      $        realday.le.date_e107(669)) then
    143          do i=1,669
     144     $        realday.lt.date_e107(669)) then
     145         do i=1,668
    144146            if(realday.ge.date_e107(i).and.
    145147     $           realday.lt.date_e107(i+1)) then
    146148               tinf=i
    147149               tsup=i+1
    148                e107=e107_tab(tinf)+(zday-date_e107(tinf))*
     150               e107=e107_tab(tinf)+(realday-date_e107(tinf))*
    149151     $              (e107_tab(tsup)-e107_tab(tinf))
    150152            endif
    151153         enddo
    152154      endif
     155
    153156      !Photoabsorption coefficients at TOA as a function of E10.7
    154157      do j=1,nabs
     
    159162         enddo
    160163      enddo
    161 
    162164! Interpolation to the tabulated photoabsorption coefficients for each species
    163165! in each spectral interval
Note: See TracChangeset for help on using the changeset viewer.