1 | !> \file accum7.f |
---|
2 | !! Calcule de l'accumulation |
---|
3 | !< |
---|
4 | |
---|
5 | !> SUBROUTINE: Accum7 |
---|
6 | !! Calcul l'accumulation de glace anuelle |
---|
7 | !! \author Catherine |
---|
8 | !! \date Avril 99 |
---|
9 | !! @note Accumulation heminord pour la deglaciation |
---|
10 | !! @note Used modules: |
---|
11 | !! @note - use module3D_phy |
---|
12 | !! @note - use climat_forcage_mod |
---|
13 | !! |
---|
14 | !< |
---|
15 | |
---|
16 | subroutine ACCUM7() |
---|
17 | |
---|
18 | ! -------------------------------------------------------------------- |
---|
19 | ! Accumulation heminord pour la deglaciation |
---|
20 | ! |
---|
21 | ! Catherine Avril 99 |
---|
22 | !--------------------------------------------------------------------- |
---|
23 | |
---|
24 | USE module3D_phy |
---|
25 | USE CLIMAT_FORCAGE_MOD |
---|
26 | ! USE CLIMAT_PROFIL_MOD |
---|
27 | |
---|
28 | implicit none |
---|
29 | real SLMIN |
---|
30 | |
---|
31 | |
---|
32 | ! calcul de Tjuly et et Tann |
---|
33 | |
---|
34 | do J=1,NY |
---|
35 | do I=1,NX |
---|
36 | ZS(I,J)=max(SEALEVEL,S(I,J)) |
---|
37 | |
---|
38 | Tann(i,j)= - TEMPGRAD * (Zs(i,j)-S0(i,j)) ! correction d'altitude |
---|
39 | & + delTatime(i,j)+TA0(I,J) |
---|
40 | |
---|
41 | Tjuly(i,j)=- TEMPGRJUL * (Zs(i,j)-S0(i,j)) ! correction d'altitude |
---|
42 | & + delTjtime(i,j)+TJ0(I,J) |
---|
43 | end do |
---|
44 | end do |
---|
45 | |
---|
46 | |
---|
47 | |
---|
48 | c ******* temperature et duree de l'hiver ***** |
---|
49 | PYY=2.*PI/50. |
---|
50 | do J=1,NY |
---|
51 | do I=1,NX |
---|
52 | FT(I,J)=0. |
---|
53 | do K=1,50 |
---|
54 | TEMP=TANN(I,J)-(TJULY(I,J)-TANN(I,J))*cos(PYY*K) |
---|
55 | if (TEMP.le.PSOLID) then |
---|
56 | FT(I,J)=FT(I,J)+1. |
---|
57 | end if |
---|
58 | end do |
---|
59 | FT(I,J)=FT(I,J)/50. |
---|
60 | end do |
---|
61 | end do |
---|
62 | |
---|
63 | c ------------- Precipitation ------------------------------ |
---|
64 | c fonction de exp(0.05 T) |
---|
65 | do J=1,NY |
---|
66 | do I=1,NX |
---|
67 | |
---|
68 | Precip(i,j)=precip0(i,j)*rapactime(i,j)* |
---|
69 | & exp(0.05*(Tann(i,j)-Ta0(i,j))) |
---|
70 | |
---|
71 | c Precip(i,j)=precip0(i,j)*rapactime(i,j) |
---|
72 | |
---|
73 | ACC(I,J)=PRECIP(I,J)*FT(i,j) |
---|
74 | |
---|
75 | end do |
---|
76 | end do |
---|
77 | |
---|
78 | |
---|
79 | |
---|
80 | end |
---|