9 |
SUBROUTINE orbite(xjour, longi, dist) |
SUBROUTINE orbite(xjour, longi, dist) |
10 |
|
|
11 |
USE yomcst, ONLY : r_ecc, r_peri |
USE yomcst, ONLY : r_ecc, r_peri |
12 |
use numer_rec, only: pi |
use nr_util, only: pi |
13 |
|
|
14 |
! Auteur(s): Z.X. Li (LMD/CNRS) |
! Auteur(s): Z.X. Li (LMD/CNRS) |
15 |
! Date: 1993/08/18 |
! Date: 1993/08/18 |
63 |
USE dimphy, ONLY : klon |
USE dimphy, ONLY : klon |
64 |
USE yomcst, ONLY : r_incl |
USE yomcst, ONLY : r_incl |
65 |
USE phyetat0_m, ONLY : rlat, rlon |
USE phyetat0_m, ONLY : rlat, rlon |
66 |
USE comconst, ONLY : pi |
use nr_util, only: assert, pi |
|
use numer_rec, only: assert |
|
67 |
|
|
68 |
! Author : O. Boucher (LMD/CNRS), d'après les routines "zenith" et |
! Author : O. Boucher (LMD/CNRS), d'après les routines "zenith" et |
69 |
! "angle" de Z.X. Li |
! "angle" de Z.X. Li |
109 |
|
|
110 |
!---------------------------------------------------------------------- |
!---------------------------------------------------------------------- |
111 |
|
|
112 |
call assert((/size(pmu0), size(frac)/) == klon, "zenang") |
if (present(frac)) call assert((/size(pmu0), size(frac)/) == klon, & |
113 |
|
"zenang") |
114 |
|
|
115 |
deux_pi = 2*pi |
deux_pi = 2*pi |
116 |
|
|
117 |
lat_sun = asin(sin(longi * pi / 180.) * sin(r_incl * pi / 180.)) |
lat_sun = asin(sin(longi * pi / 180.) * sin(r_incl * pi / 180.)) |