1 | CCC $Header$ |
---|
2 | CDIR$ LIST |
---|
3 | SUBROUTINE h3cint(kt) |
---|
4 | C |
---|
5 | #if defined key_trc_hamocc3 |
---|
6 | C |
---|
7 | CDIR$ NOLIST |
---|
8 | USE oce_trc |
---|
9 | USE trp_trc |
---|
10 | USE sms |
---|
11 | IMPLICIT NONE |
---|
12 | CDIR$ LIST |
---|
13 | INTEGER kt |
---|
14 | INTEGER ji, jj |
---|
15 | INTEGER iman |
---|
16 | INTEGER ipdtant, ipdtmo |
---|
17 | INTEGER nvit1t,nvit2t |
---|
18 | REAL zpdtan, zman, zpdtmo, zdemi |
---|
19 | REAL zt |
---|
20 | C |
---|
21 | ipdtant = (12*730*60*60)/int(rdttra(1)) |
---|
22 | iman = ilecvit |
---|
23 | ipdtmo = ipdtant / iman |
---|
24 | C |
---|
25 | zpdtan = float ( ipdtant) |
---|
26 | zman = float ( iman) |
---|
27 | zpdtmo = zpdtan / zman |
---|
28 | zdemi = zpdtmo / 2. |
---|
29 | zt = ( float ( kt ) + zdemi) / zpdtmo |
---|
30 | |
---|
31 | C recherche de l'indice des enregistrements |
---|
32 | C du modele dynamique encadrant le pas de temps kt. |
---|
33 | C -------------------------------------------------- |
---|
34 | C |
---|
35 | xtvit = zt - int ( zt) |
---|
36 | nvit1t = int (( float ( kt) + zdemi)/ zpdtmo) |
---|
37 | nvit2t = nvit1t+1 |
---|
38 | nvit1t = MOD ( nvit1t, iman) |
---|
39 | IF ( nvit1t .EQ. 0 ) nvit1t = iman |
---|
40 | nvit2t = MOD ( nvit2t, iman) |
---|
41 | IF ( nvit2t .EQ. 0 ) nvit2t = iman |
---|
42 | C |
---|
43 | C vertical slab |
---|
44 | C =============== |
---|
45 | C |
---|
46 | IF (icice.NE.3.and.icice.ne.1) THEN |
---|
47 | DO jj = 1,jpj |
---|
48 | C |
---|
49 | c --------------------------------------------------------- |
---|
50 | c interpolation de la couverture de glace et du vent |
---|
51 | c --------------------------------------------------------- |
---|
52 | c |
---|
53 | DO ji=1,jpi |
---|
54 | freeze(ji,jj) = (1.-xtvit)*cicemo(ji,jj,nvit1t) |
---|
55 | $ +xtvit*cicemo(ji,jj,nvit2t) |
---|
56 | END DO |
---|
57 | END DO |
---|
58 | ENDIF |
---|
59 | IF (igaswind.NE.3.and.igaswind.ne.1) THEN |
---|
60 | DO jj = 1,jpj |
---|
61 | DO ji=1,jpi |
---|
62 | vatm(ji,jj) = (1.-xtvit)*wsmo(ji,jj,nvit1t) |
---|
63 | $ +xtvit*wsmo(ji,jj,nvit2t) |
---|
64 | ENDDO |
---|
65 | END DO |
---|
66 | ENDIF |
---|
67 | IF (igaswind.EQ.2) THEN |
---|
68 | DO jj = 1,jpj |
---|
69 | DO ji=1,jpi |
---|
70 | kgwanin(ji,jj) = (1.-xtvit)*kgwanmo(ji,jj,nvit1t) |
---|
71 | $ +xtvit*kgwanmo(ji,jj,nvit2t) |
---|
72 | ENDDO |
---|
73 | END DO |
---|
74 | ENDIF |
---|
75 | C |
---|
76 | C |
---|
77 | #endif |
---|
78 | C |
---|
79 | RETURN |
---|
80 | END |
---|
81 | |
---|
82 | |
---|