source: trunk/SOURCES/Recul_force_grounding_line/declar_map_recul_mod.f90 @ 159

Last change on this file since 159 was 4, checked in by dumas, 10 years ago

initial import GRISLI trunk

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