Last change
on this file since 26 was
6,
checked in by vancop, 8 years ago
|
initial import of v3.20 /Users/ioulianikolskaia/Boulot/CODES/LIM1D/ARCHIVE/TMP/LIM1D_v3.20/
|
File size:
1.2 KB
|
Line | |
---|
1 | c |
---|
2 | c ******************** |
---|
3 | function pdecli(i,j) |
---|
4 | c ******************** |
---|
5 | c |
---|
6 | c THIS FUNCTION COMPUTES THE SOLAR DECLINATION FOR THE DAY |
---|
7 | c j (IN DECIMAL DEGREES). |
---|
8 | c j = DAY OF THE YEAR (j=1 ON JANUARY 1) |
---|
9 | c i = -1, 0, 1 FOR ODD, NORMAL AND LEAP YEARS, |
---|
10 | c RESPECTIVELY |
---|
11 | c |
---|
12 | c--- |
---|
13 | c Ccpl [Ccp0] => ligne specifique a la version avec [sans] couplage . |
---|
14 | c--- |
---|
15 | include 'type.com' |
---|
16 | include 'const.com' |
---|
17 | xj = j |
---|
18 | a0 = 0.39507671 |
---|
19 | a1 = 22.85684301 |
---|
20 | a2 = -0.38637317 |
---|
21 | a3 = 0.15096535 |
---|
22 | a4 = -0.00961411 |
---|
23 | b1 = -4.29692073 |
---|
24 | b2 = 0.05702074 |
---|
25 | b3 = -0.09028607 |
---|
26 | b4 = 0.00592797 |
---|
27 | pi = 3.1415927 |
---|
28 | if (i) 1,2,3 |
---|
29 | 1 xj = xj-0.5 |
---|
30 | goto 4 |
---|
31 | 3 xj = xj-1. |
---|
32 | goto 4 |
---|
33 | 2 xj = j |
---|
34 | 4 p = pi*(2.0*xj-367.0) |
---|
35 | p = p/yeaday |
---|
36 | p2 = p*2.0 |
---|
37 | p3 = p*3.0 |
---|
38 | p4 = p*4.0 |
---|
39 | dcl1 = a1*cos(p)+a2*cos(p2)+a3*cos(p3)+a4*cos(p4) |
---|
40 | dcl2 = b1*sin(p)+b2*sin(p2)+b3*sin(p3)+b4*sin(p4) |
---|
41 | pdecli = a0+dcl1+dcl2 |
---|
42 | c |
---|
43 | return |
---|
44 | c---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
---|
45 | c- Fin de la routine pdecli - |
---|
46 | end |
---|
Note: See
TracBrowser
for help on using the repository browser.