Changeset 1626 for LMDZ5/trunk


Ignore:
Timestamp:
May 15, 2012, 10:37:02 AM (13 years ago)
Author:
Laurent Fairhead
Message:

Modifications nécessaires sur les dates pour avoir 1+1=2 avec l'ancienne physique

Necessary modifications to ensure 1+1=2 with the old physics package

Location:
LMDZ5/trunk/libf
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/dyn3d/leapfrog.F

    r1625 r1626  
    225225
    226226      jD_cur = jD_ref + day_ini - day_ref +                             &
    227      &          int (itau * dtvr / daysec)
     227     &          itau/day_step
    228228      jH_cur = jH_ref + start_time +                                    &
    229      &          (itau * dtvr / daysec - int(itau * dtvr / daysec))
     229     &          mod(itau,day_step)/float(day_step)
    230230      jD_cur = jD_cur + int(jH_cur)
    231231      jH_cur = jH_cur - int(jH_cur)
    … …  
    382382!     &              (itau * dtvr / daysec - int(itau * dtvr / daysec))
    383383           jD_cur = jD_ref + day_ini - day_ref +                        &
    384      &          int (itau * dtvr / daysec)
     384     &          itau/day_step
    385385           jH_cur = jH_ref + start_time +                               &
    386      &          (itau * dtvr / daysec - int(itau * dtvr / daysec))
     386     &              mod(itau,day_step)/float(day_step)
    387387           jD_cur = jD_cur + int(jH_cur)
    388388           jH_cur = jH_cur - int(jH_cur)

  • LMDZ5/trunk/libf/dyn3dpar/leapfrog_p.F

    r1625 r1626  
    247247
    248248      jD_cur = jD_ref + day_ini - day_ref +                             &
    249      &          int (itau * dtvr / daysec)
     249     &          itau/day_step
    250250      jH_cur = jH_ref + start_time +                                    &
    251      &          (itau * dtvr / daysec - int(itau * dtvr / daysec))
     251     &         mod(itau,day_step)/float(day_step)
    252252      if (jH_cur > 1.0 ) then
    253253        jD_cur = jD_cur +1.
    … …  
    706706c$OMP BARRIER
    707707           jD_cur = jD_ref + day_ini - day_ref
    708      $        + int (itau * dtvr / daysec)
     708     $        + itau/day_step
    709709           jH_cur = jH_ref + start_time +                                &
    710      &              (itau * dtvr / daysec - int(itau * dtvr / daysec))
     710     &              mod(itau,day_step)/float(day_step)
    711711!         call ju2ymds(jD_cur+jH_cur, an, mois, jour, secondes)
    712712           if (jH_cur > 1.0 ) then
Note: See TracChangeset for help on using the changeset viewer.