source:
tags/nemo_dev_x9/NEMO/TOP_SRC/SMS/p4zday.F
@
5089
Last change on this file since 5089 was 186, checked in by opalod, 20 years ago | |
---|---|
|
|
File size: 1.9 KB |
Line | |
---|---|
1 | CDIR$ LIST |
2 | SUBROUTINE p4zday |
3 | #if defined key_passivetrc && defined key_trc_pisces |
4 | CCC--------------------------------------------------------------------- |
5 | CCC |
6 | CCC ROUTINE p4zday : PISCES MODEL |
7 | CCC ***************************** |
8 | CCC |
9 | CCC PURPOSE : |
10 | CCC --------- |
11 | CCC PISCES : compute the day length depending on latitude |
12 | CCC and the day |
13 | CCC |
14 | CC INPUT : |
15 | CC ----- |
16 | CC argument |
17 | CC ktask : task identificator |
18 | CC common |
19 | CC all the common defined in opa |
20 | CC |
21 | CC |
22 | CC OUTPUT : : no |
23 | CC ------ |
24 | CC |
25 | CC EXTERNAL : |
26 | CC -------- |
27 | CC None |
28 | CC |
29 | CC MODIFICATIONS: |
30 | CC -------------- |
31 | CC original : E. Maier-Reimer (GBC 1993) |
32 | CC additions : C. Le Quere (1999) |
33 | CC modifications : O. Aumont (2004) |
34 | CC---------------------------------------------------------------------- |
35 | CC parameters and commons |
36 | CC ====================== |
37 | USE oce_trc |
38 | USE trp_trc |
39 | USE sms |
40 | IMPLICIT NONE |
41 | CC---------------------------------------------------------------------- |
42 | CC local declarations |
43 | CC ================== |
44 | INTEGER ji, jj, iyy |
45 | REAL rum, delta, codel, phi, argu |
46 | C |
47 | C Get year |
48 | C -------- |
49 | C |
50 | iyy = ndastp/10000 |
51 | |
52 | write(numout,*) |
53 | write(numout,*) 'p4zday - Julian day ', nday_year |
54 | write(numout,*) |
55 | |
56 | |
57 | IF (nleapy.EQ.1 .AND. MOD(iyy,4).EQ.0) THEN |
58 | rum = FLOAT(nday_year-80)/366. |
59 | ELSE |
60 | rum = FLOAT(nday_year-80)/365. |
61 | ENDIF |
62 | |
63 | C |
64 | delta = SIN(rum*rpi*2.)*sin(rpi*23.5/180.) |
65 | codel = asin(delta) |
66 | |
67 | DO jj = 1,jpj |
68 | DO ji = 1,jpi |
69 | phi = gphit(ji,jj)*rpi/180. |
70 | argu = tan(codel)*tan(phi) |
71 | strn(ji,jj) = 0. |
72 | argu=min(1.,argu) |
73 | argu=max(-1.,argu) |
74 | strn(ji,jj)=24.-2.*acos(argu)*180./rpi/15. |
75 | strn(ji,jj)=max(strn(ji,jj),0.) |
76 | END DO |
77 | END DO |
78 | C |
79 | #endif |
80 | RETURN |
81 | END |
82 |
Note: See TracBrowser
for help on using the repository browser.