Changeset 3010


Ignore:
Timestamp:
Jul 24, 2023, 8:59:01 AM (17 months ago)
Author:
llange
Message:

MARS PCM
Fixing a bug introduced in r2963: timeperi was forced to be >0. Consequently, Lsp could only vary between 180° - 360°. Simulations made with Lsp < 180° should not be correct.
Remove the abs().
LL

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.MARS/libf/phymars/iniorbit.F

    r2963 r3010  
    6868
    6969      timeperi=2.*atan(sqrt((1.+e_elips)/(1.-e_elips))*tan(zx0/2.))
    70       if(timeperi.lt.0) timeperi=-timeperi
    7170      PRINT*,'iniorbit: Perihelion solar long. Ls (deg)=',
    7271     &       360.-timeperi*180./pi
Note: See TracChangeset for help on using the changeset viewer.