Changeset 3504 for trunk/LMDZ.PLUTO/libf/phypluto/iniorbit.F
- Timestamp:
- Nov 8, 2024, 10:57:14 AM (13 days ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.PLUTO/libf/phypluto/iniorbit.F
r3184 r3504 1 1 SUBROUTINE iniorbit 2 2 $ (papoastr,pperiastr,pyear_day,pperi_day,pobliq) 3 3 4 4 USE planete_mod, only: apoastr, periastr, year_day, obliquit, 5 5 & peri_day, e_elips, p_elips, timeperi 6 6 use comcstfi_mod, only: pi 7 use mod_phys_lmdz_para, only : is_master 7 8 IMPLICIT NONE 8 9 … … 32 33 peri_day=pperi_day 33 34 34 PRINT*,'iniorbit: Periastron in AU ',periastr 35 PRINT*,'iniorbit: Apoastron in AU ',apoastr 36 PRINT*,'iniorbit: Obliquity in degrees :',obliquit 37 35 if (is_master) then 36 PRINT*,'iniorbit: Periastron in AU ',periastr 37 PRINT*,'iniorbit: Apoastron in AU ',apoastr 38 PRINT*,'iniorbit: Obliquity in degrees :',obliquit 39 endif 38 40 39 41 e_elips=(apoastr-periastr)/(periastr+apoastr) 40 42 p_elips=0.5*(periastr+apoastr)*(1-e_elips*e_elips) 41 43 42 print*,'iniorbit: e_elips',e_elips 43 print*,'iniorbit: p_elips',p_elips 44 44 if (is_master) then 45 print*,'iniorbit: e_elips',e_elips 46 print*,'iniorbit: p_elips',p_elips 47 endif 45 48 !----------------------------------------------------------------------- 46 49 ! compute polar angle and distance to the Sun:
Note: See TracChangeset
for help on using the changeset viewer.