1 | !> \file dragging_hwatermax_0.2_mod.f90 |
---|
2 | !! Defini les zones de stream avec un critere |
---|
3 | !< |
---|
4 | |
---|
5 | !> \namespace dragging_hwatstream |
---|
6 | !! Defini les zones de stream avec 1 critere |
---|
7 | !! \author ... |
---|
8 | !! \date ... |
---|
9 | !! @note Le critere est: |
---|
10 | !! @note * un critere sur la hauteur d'eau |
---|
11 | !! @note Used module |
---|
12 | !! @note - use module3D_phy |
---|
13 | !< |
---|
14 | module dragging_hwatstream |
---|
15 | ! Defini les zones de stream avec un critere sur la hauteur d'eau |
---|
16 | |
---|
17 | USE module3D_phy |
---|
18 | |
---|
19 | |
---|
20 | |
---|
21 | implicit none |
---|
22 | real :: tostick |
---|
23 | |
---|
24 | contains |
---|
25 | !------------------------------------------------------------------------------- |
---|
26 | !> SUBROUTINE: init_dragging |
---|
27 | !!Cette routine fait l'initialisation du dragging. |
---|
28 | !> |
---|
29 | subroutine init_dragging ! Cette routine fait l'initialisation du dragging. |
---|
30 | |
---|
31 | implicit none |
---|
32 | |
---|
33 | namelist/drag_hwat/hwatstream,cf,betamax,toblim,tostick |
---|
34 | |
---|
35 | if (itracebug.eq.1) call tracebug(' Dragging avec hwatermax') |
---|
36 | |
---|
37 | |
---|
38 | ! formats pour les ecritures dans 42 |
---|
39 | 428 format(A) |
---|
40 | |
---|
41 | ! lecture des parametres du run block draghwat |
---|
42 | !-------------------------------------------------------------------- |
---|
43 | |
---|
44 | rewind(num_param) ! pour revenir au debut du fichier param_list.dat |
---|
45 | read(num_param,drag_hwat) |
---|
46 | |
---|
47 | write(num_rep_42,428)'!___________________________________________________________' |
---|
48 | write(num_rep_42,428) '&drag_hwat ! nom du bloc hwatermax ' |
---|
49 | write(num_rep_42,*) |
---|
50 | write(num_rep_42,*) 'hwatsream = ',hwatstream |
---|
51 | write(num_rep_42,*) 'cf = ',cf |
---|
52 | write(num_rep_42,*) 'betamax = ', betamax |
---|
53 | write(num_rep_42,*) 'toblim = ', toblim |
---|
54 | write(num_rep_42,*)'/' |
---|
55 | write(num_rep_42,428) '! hwatstream (m) : critere de passage en stream si hwater > hwatstream ' |
---|
56 | write(num_rep_42,428) '! cf coefficient de la loi de frottement fonction Neff' |
---|
57 | write(num_rep_42,428) '! betamax : (Pa) frottement maxi ' |
---|
58 | write(num_rep_42,428) '! toblim : (Pa) tes pour les iles ' |
---|
59 | write(num_rep_42,*) |
---|
60 | |
---|
61 | tostick=1.e5 ! valeurs pour les points non flgzmx |
---|
62 | |
---|
63 | !------------------------------------------------------------------- |
---|
64 | ! masque stream |
---|
65 | |
---|
66 | mstream_mx(:,:)=1 |
---|
67 | mstream_my(:,:)=1 |
---|
68 | |
---|
69 | |
---|
70 | ! coefficient permettant de modifier le basal drag. |
---|
71 | drag_mx(:,:)=1. |
---|
72 | drag_my(:,:)=1. |
---|
73 | |
---|
74 | |
---|
75 | return |
---|
76 | end subroutine init_dragging |
---|
77 | !________________________________________________________________________________ |
---|
78 | |
---|
79 | !> SUBROUTINE: dragging |
---|
80 | !!Defini la localisation des streams et le frottement basal |
---|
81 | !> |
---|
82 | !------------------------------------------------------------------------- |
---|
83 | subroutine dragging ! defini la localisation des streams et le frottement basal |
---|
84 | |
---|
85 | |
---|
86 | ! les iles n'ont pas de condition neff mais ont la condition toblim |
---|
87 | ! (ce bloc etait avant dans flottab) |
---|
88 | |
---|
89 | do j=2,ny |
---|
90 | do i=2,nx |
---|
91 | ilemx(i,j)=ilemx(i,j).and.(abs(tobmx(i,j)).lt.toblim) |
---|
92 | ilemy(i,j)=ilemy(i,j).and.(abs(tobmy(i,j)).lt.toblim) |
---|
93 | end do |
---|
94 | end do |
---|
95 | |
---|
96 | |
---|
97 | gzm_x: do j=1,ny ! 2,ny |
---|
98 | do i=2,nx ! -1 |
---|
99 | |
---|
100 | ! gzmx pour les hauteurs d'eau superieures a un seuil |
---|
101 | ! attention, il semble qu'il y avait une erreur dans la precedente |
---|
102 | ! formulation car le test etait seulement sur hwater(i,j). |
---|
103 | ! le gzmx initial a ete calculé dans flottab pour les points a cheval sur la grounding line |
---|
104 | |
---|
105 | gzmx(i,j)=gzmx(i,j).or. & |
---|
106 | (((hwater(i,j)+hwater(i-1,j))*0.5.gt.hwatstream) & |
---|
107 | .and..not.flotmx(i,j)) |
---|
108 | |
---|
109 | |
---|
110 | ! calcul du frottement basal (ce bloc etait avant dans neffect) |
---|
111 | if (gzmx(i,j)) then ! stream |
---|
112 | betamx(i,j)=cf*neffmx(i,j) |
---|
113 | else if (flotmx(i,j).or.ilemx(i,j)) then ! flottant ou ile |
---|
114 | betamx(i,j)=0. |
---|
115 | else ! grounded, SIA |
---|
116 | betamx(i,j)=tostick ! frottement glace posee (1 bar) |
---|
117 | endif |
---|
118 | |
---|
119 | end do |
---|
120 | end do gzm_x |
---|
121 | |
---|
122 | gzm_y: do j=2,ny !-1 |
---|
123 | do i=1,nx ! 2,nx |
---|
124 | |
---|
125 | ! gzmy pour les hauteurs d'eau superieures a un seuil |
---|
126 | ! attention, il semble qu'il y avait une erreur dans la precedente |
---|
127 | ! formulation car le test etait seulement sur hwater(i,j). |
---|
128 | |
---|
129 | gzmy(i,j)=gzmy(i,j).or. & |
---|
130 | (((hwater(i,j)+hwater(i,j-1))*0.5.gt.hwatstream) & |
---|
131 | .and..not.flotmy(i,j)) |
---|
132 | |
---|
133 | ! calcul du frottement basal (ce bloc etait avant dans neffect) |
---|
134 | if (gzmy(i,j)) then ! stream |
---|
135 | betamy(i,j)=cf*neffmy(i,j) |
---|
136 | else if (flotmy(i,j).or.ilemy(i,j)) then ! flottant ou ile |
---|
137 | betamy(i,j)=0. |
---|
138 | else ! grounded, SIA |
---|
139 | betamy(i,j)=tostick ! frottement glace posee |
---|
140 | endif |
---|
141 | |
---|
142 | end do |
---|
143 | end do gzm_y |
---|
144 | |
---|
145 | return |
---|
146 | end subroutine dragging |
---|
147 | |
---|
148 | end module dragging_hwatstream |
---|
149 | |
---|