Ignore:
Timestamp:
Feb 24, 2020, 3:36:36 PM (5 years ago)
Author:
jvatant
Message:

Get rid of the old 'sparadrap' in iniorbit.F assuming Mkm instead of AU if periastre or apostre gt 25, (otherwise Neptune is a pb)
Everybody should be using AU, that's all.
--JVO

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.GENERIC/libf/phystd/iniorbit.F

    r1384 r2243  
    3232      peri_day=pperi_day
    3333
    34 
    35       !!!! SPARADRAP TEMPORAIRE !!!!
    36       !!!! SPARADRAP TEMPORAIRE !!!!
    37       !!!! We hope that all cases are above 25 Mkm [OK with Gliese 581d]
    38       IF ( apoastr .gt. 25.) THEN
    39         PRINT*,'!!!!! WARNING !!!!!'
    40         PRINT*,'!!!!! YOU ARE ABOUT TO WITNESS A DIRT HACK !!!!!'
    41         PRINT*,'This must be an old start file.'
    42         PRINT*,'The code changed 19/03/2012: we now use AU.'
    43         PRINT*,'So I am assuming units are in Mkm here'
    44         PRINT*,'and I am performing a conversion towards AU.'
    45         periastr = periastr / 149.598 ! Mkm to AU
    46         apoastr = apoastr / 149.598 ! Mkm to AU
    47       ENDIF
    48       !!!! SPARADRAP TEMPORAIRE !!!!
    49       !!!! SPARADRAP TEMPORAIRE !!!!
    50 
    51  
    5234      PRINT*,'iniorbit: Periastron in AU  ',periastr
    5335      PRINT*,'iniorbit: Apoastron in AU  ',apoastr
Note: See TracChangeset for help on using the changeset viewer.