1 | module declar_map_recul |
---|
2 | |
---|
3 | implicit none |
---|
4 | integer, parameter :: nx = 387 |
---|
5 | integer, parameter :: ny = 374 |
---|
6 | integer :: i,j ! pour les boucles |
---|
7 | |
---|
8 | ! donnees d'entree |
---|
9 | real,dimension(nx,ny) :: H0 !< epaisseur |
---|
10 | real,dimension(nx,ny) :: H !< courant |
---|
11 | real,dimension(nx,ny) :: Bsoc !< socle |
---|
12 | real,dimension(nx,ny) :: B !< base |
---|
13 | real,dimension(nx,ny) :: S !< surface |
---|
14 | real,dimension(nx) :: X_15 !< dimension X |
---|
15 | real,dimension(ny) :: Y_15 !< dimension Y |
---|
16 | |
---|
17 | ! utilise dans les calculs |
---|
18 | real,dimension(nx,ny) :: delHdt_mx !< variation d'epaisseur imposee venant de x |
---|
19 | real,dimension(nx,ny) :: delHdt_my !< variation d'epaisseur imposee venant de y |
---|
20 | real,dimension(nx,ny) :: delHdt !< variation d'epaisseur imposee venant de y |
---|
21 | real,dimension(nx,ny) :: H_float !< epaisseur de flottaison |
---|
22 | real,dimension(nx,ny) :: time_float !< temps auquel le noeud se met a flotter |
---|
23 | real,dimension(nx,ny) :: Deltatot_H0 !< variation d'epaisseur totale au moment ou le noeud |
---|
24 | !< commence a decroitre |
---|
25 | |
---|
26 | real,dimension(4) :: dH_voisin ! quand on balaye les voisins |
---|
27 | |
---|
28 | real,dimension(nx,ny) :: time_dep_mx !< temps a partir duquel le recul est autorise (noeuds >) |
---|
29 | real,dimension(nx,ny) :: time_dep_my !< temps a partir duquel le recul est autorise (noeuds ^) |
---|
30 | |
---|
31 | real,dimension(nx,ny) :: x_gl_mx !< position sous maille de la grounding line selon x |
---|
32 | real,dimension(nx,ny) :: y_gl_my !< position sous maille de la grounding line selon y |
---|
33 | !< sur les mailles mineures |
---|
34 | |
---|
35 | integer,dimension(nx,ny) :: mk_traiter !< masque des noeuds a traiter |
---|
36 | integer,dimension(nx,ny) :: mk_gr !< masque des noeuds flottants |
---|
37 | integer,dimension(nx,ny) :: newfloat !< masque des noeuds qui viennent de flotter |
---|
38 | integer,dimension(nx,ny) :: i_delta !< index qui dit si un noeud est en cours de decroissance |
---|
39 | |
---|
40 | real,dimension(nx,ny) :: retreat0_x !< vitesse maxi autorisee selon x : valeur initiale |
---|
41 | real,dimension(nx,ny) :: retreat0_y !< vitesse maxi autorisee selon y : valeur initiale |
---|
42 | real,dimension(nx,ny) :: retreat_x !< vitesse maxi autorisee selon x |
---|
43 | real,dimension(nx,ny) :: retreat_y !< surface maxi autorisee selon y |
---|
44 | |
---|
45 | |
---|
46 | |
---|
47 | real :: dx !< pas d'espace |
---|
48 | real :: dt !< pas de temps |
---|
49 | real :: time !< temps |
---|
50 | real :: retreat_rate !< taux de retrait |
---|
51 | real :: coef_Bflot !< pour ice shelves B = coef_Bflot * H |
---|
52 | integer :: som_voisins !< pour des tests |
---|
53 | integer :: i_time !< pour boucle temps |
---|
54 | integer :: n_time !< pour boucle temps |
---|
55 | |
---|
56 | |
---|
57 | |
---|
58 | |
---|
59 | ! pour les io netcdf |
---|
60 | !------------------------- |
---|
61 | |
---|
62 | integer :: ncid ! numero d'unite |
---|
63 | integer :: status !< operations io netcdf |
---|
64 | character(len=100) :: file_ncdf !< fichier netcdf |
---|
65 | character(len=100) :: fileout_ncdf !< fichier netcdf |
---|
66 | character(len=20),dimension(2) :: dimname ! pour les noms de dimension |
---|
67 | Character*20,dimension(:),pointer :: tabvarname ! pour les noms de variable |
---|
68 | real*8, dimension(:,:), pointer :: work_tab !< tableau 2d real ecrit dans le fichier |
---|
69 | real*8,dimension(:),pointer :: Tabx ! pointeur vers le tableau des x |
---|
70 | real*8,dimension(:),pointer :: Taby ! pointeur vers le tableau des y |
---|
71 | |
---|
72 | end module declar_map_recul |
---|