Last change
on this file since 67 was
4,
checked in by vancop, 8 years ago
|
initial import /Users/ioulianikolskaia/Boulot/CODES/LIM1D/ARCHIVE/TMP/LIM1D_v3.20/
|
File size:
1.2 KB
|
Rev | Line | |
---|
[4] | 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.