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