1 | |
---|
2 | subroutine initial_heino() |
---|
3 | |
---|
4 | ! ************************************************** |
---|
5 | ! ** initialisation * |
---|
6 | ! * * |
---|
7 | ! ************************************************** |
---|
8 | |
---|
9 | USE module3D_phy |
---|
10 | USE DEFORM_DECLAR |
---|
11 | ! use sliding_dragging_heino |
---|
12 | |
---|
13 | implicit none |
---|
14 | |
---|
15 | !cdc if (itracebug.eq.1) call tracebug(' Heino: entree dans routine intial_heino') |
---|
16 | |
---|
17 | !====================================== La reponse est 42 =========== |
---|
18 | write(num_rep_42,*) |
---|
19 | write(num_rep_42,*)'Heino : Modification de certains parametres' |
---|
20 | write(num_rep_42,*)'____________________________________________________' |
---|
21 | write(num_rep_42,*) |
---|
22 | write(num_rep_42,*) ' la valeur des parametres suivants est modifiee' |
---|
23 | write(num_rep_42,*) |
---|
24 | !===================================================================== |
---|
25 | |
---|
26 | |
---|
27 | ! définition des pas de temps en fonction du pas d'espace |
---|
28 | ! loigliss=4 |
---|
29 | ! write(42,*)'loi de glissement Heino',loigliss |
---|
30 | |
---|
31 | |
---|
32 | ! write(num_rep_42,*) |
---|
33 | ! write(num_rep_42,*)' ice streams' |
---|
34 | ! write(num_rep_42,*) '-----------' |
---|
35 | ! toblim=5.e5 |
---|
36 | ! write(num_rep_42,*) 'seuil de contrainte basale, toblim=',toblim |
---|
37 | |
---|
38 | ! hwatstream=5000. !meters |
---|
39 | |
---|
40 | ! hwatermax=.2 ! hauteur max d'eau sous glaciaire |
---|
41 | ! write(num_rep_42,*)'hauteur d eau maximum', hwatermax |
---|
42 | |
---|
43 | ! write(num_rep_42,*)' passage en stream si hwater > hwatstream=', & |
---|
44 | ! hwatstream |
---|
45 | ! print*,'test' |
---|
46 | ! Loi de déformation prescrite par Heino |
---|
47 | SF(1)=3. |
---|
48 | TTRANS(1)=-10. |
---|
49 | ! BAT1(1)=SF(1)*7.65e-17 ! Pa-3 a-1 Heino idem Huybrechts ancien |
---|
50 | BAT1(1)=7.645058e-17 ! meilleur ajustement si Tmelt=0 |
---|
51 | Q1(1)=60.e+3 ! energie d'activation T < TTRANS(1) |
---|
52 | |
---|
53 | ! Rappel : valeurs habituelles Ritz BAT1(1)=SF(1)*0.166e-15 |
---|
54 | ! Glen T<Trans Q1(1)=78.2e+3 |
---|
55 | |
---|
56 | ! BAT2(1)=SF(1)*0.28642e-15 ! Pa -3 a-1 ancien Heino |
---|
57 | BAT2(1)=2.8586262e-16 ! meilleur ajustement si Tmelt=0 |
---|
58 | Q2(1)=139.e+3 |
---|
59 | |
---|
60 | ! Rappel : valeurs habituelles Ritz BAT2(1)=SF(1)*0.2e-15 |
---|
61 | ! Glen T>Ttrans Q2(1)=95.45e+3 T > TTRANS(1) |
---|
62 | |
---|
63 | TTRANS(2)=-10. |
---|
64 | ! SF(2)=0. |
---|
65 | SF(2)=0. ! pas de fluidite newtonnienne dans Heino |
---|
66 | ! SF(2)=3. |
---|
67 | ! SF(2)=6. ! coefficient multiplicateur |
---|
68 | |
---|
69 | ! Pour les temperatures inferieures a TTRPHI |
---|
70 | |
---|
71 | BAT1(2)=8.313e-8 ! fluidite Newtonienne |
---|
72 | Q1(2)=40.e+3 ! energie d'activation |
---|
73 | |
---|
74 | ! Pour les temperatures superieures a TTRPHI |
---|
75 | BAT2(2)=0. ! fluidite Newtonienne |
---|
76 | BAT2(2)=BAT1(2) ! attention rajouté pour l'ice shelf |
---|
77 | |
---|
78 | ! Q2(2)=SPHI*8.313e-8 |
---|
79 | Q2(2)=60.e+3 ! energie d'activation |
---|
80 | |
---|
81 | write(num_rep_42,*) |
---|
82 | write(num_rep_42,*) 'Loi de deformation Heino n, sf, ttrans, bat1, Q1, bat2, Q2' |
---|
83 | do iglen=1,2 |
---|
84 | write(num_rep_42,fmt=123) int(glen(iglen)), sf(iglen), ttrans(iglen),bat1(iglen), & |
---|
85 | Q1(iglen), bat2(iglen), Q2(iglen) |
---|
86 | 123 format(i2,1x,f0.2,1x,f0.3,1x,4(es10.3,1x)) |
---|
87 | ! application des sf |
---|
88 | bat1(iglen)=bat1(iglen)*sf(iglen) |
---|
89 | bat2(iglen)=bat2(iglen)*sf(iglen) |
---|
90 | end do |
---|
91 | |
---|
92 | |
---|
93 | |
---|
94 | ! ------------------------ fin loi de déformation -------------- |
---|
95 | |
---|
96 | |
---|
97 | |
---|
98 | |
---|
99 | ndisp=1 |
---|
100 | write(num_rep_42,*)' sorties temporelles tous les',ndisp |
---|
101 | |
---|
102 | ! dtt=0.2, dtmax=0.2 pour Heinrich events |
---|
103 | ! dtt=.25 |
---|
104 | ! dtmax=0.25 |
---|
105 | |
---|
106 | write(num_rep_42,*)'dtmax est maintenant impose ici (avant il y avait mise a 5 ans dans icethick)' |
---|
107 | write(num_rep_42,*)' pas de temps : dtmax=',dtmax,' dtt=',dtt |
---|
108 | |
---|
109 | ! testdiag=0.02 |
---|
110 | ! testdiag=0.005 |
---|
111 | write(num_rep_42,*)'testdiag=',testdiag |
---|
112 | end subroutine initial_heino |
---|