source: branches/iLoveclim/SOURCES/Recul_force_grounding_line/declar_proto_recul.f90 @ 30

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

initial import GRISLI trunk
File size: 4.9 KB
Line 
1module declar_proto_recul
2
3use proto_3D_grisli
4use declar_io_ncdf_proto
5
6implicit none
7
8! utilise dans les calculs specifiques au recul force
9
10integer,parameter            :: nb_regions = 14   ! nombre de regions
11
12
13character(len=100) :: number_data_B_file          ! nombre de donnees socle dans la maille
14character(len=100) :: region_file                 ! regions pour les temps de demarrage
15character(len=100) :: file_exp                    ! nom du fichier experiences 
16real               :: retreat_1                   !< taux de retrait 1 si beta < beta_lim en km
17real               :: retreat_2                   !< taux de retrait 2 si beta > beta_lim
18real               :: beta_lim_ret                !< limite valeur beta
19real               :: bed_sigma_1                 !< pour le sigma du socle quand il y a des donnees
20real               :: bed_sigma_2                 !< pour le sigma du socle quand il n'y a pas de donnees
21
22
23integer, dimension(nx,ny)          :: map_region  !< regions
24real, dimension(nx,ny)             :: Bsoc_sigma  !< incertitude sur Bsoc : en entree nombre de points de mesure
25real, dimension(nx,ny)             :: Btest       !< pour le test sur la pente du socle
26 
27real,dimension(nx,ny)              :: delHdx_mx   !< variation d'epaisseur imposee venant de x
28real,dimension(nx,ny)              :: delHdy_my   !< variation d'epaisseur imposee venant de y
29real,dimension(nx,ny)              :: delHdx_SW   !< variation d'epaisseur imposee venant de diag SW-NE
30real,dimension(nx,ny)              :: delHdx_SE   !< variation d'epaisseur imposee venant de diag SE-NW (noeuds \)
31
32real,dimension(nx,ny)              :: delHdt      !< variation d'epaisseur totale
33real,dimension(nx,ny)              :: delHdt_sanity !< variation d'epaisseur maxi autorise par le sanity check
34real,dimension(nx,ny)              :: H_float     !< epaisseur de flottaison
35real,dimension(nx,ny)              :: time_float  !< temps auquel le noeud se met a flotter
36
37
38! modulation des temps de depart en fonction de l'ocean
39
40real,dimension(0:nb_regions)       :: time_region !< temps a partir duquel le recul est autorise par region
41                                                  !< la region 0 est l'ocean
42real,dimension(nx,ny)              :: time_dep    !< temps a partir duquel le recul est autorise (noeuds centres)
43real,dimension(nx,ny)              :: time_dep_mx !< temps a partir duquel le recul est autorise (noeuds >)
44real,dimension(nx,ny)              :: time_dep_my !< temps a partir duquel le recul est autorise (noeuds ^)
45real,dimension(nx,ny)              :: time_dep_SW !< temps a partir duquel le recul est autorise (noeuds /)
46real,dimension(nx,ny)              :: time_dep_SE !< temps a partir duquel le recul est autorise (noeuds \)
47
48
49!< position grounding line sur les mailles mineures
50
51real,dimension(nx,ny)              :: x_gl_mx     !< position sous maille de la grounding line selon x
52real,dimension(nx,ny)              :: y_gl_my     !< position sous maille de la grounding line selon y
53real,dimension(nx,ny)              :: SW_gl_m     !< position sous maille de la grounding line selon SW-NE (noeuds /)
54real,dimension(nx,ny)              :: SE_gl_m     !< position sous maille de la grounding line selon SE-NW (noeuds \)
55                           
56integer,dimension(nx,ny)           :: mk_traiter  !< masque des noeuds a traiter
57integer,dimension(nx,ny)           :: mk_traiter0 !< masque des noeuds a traiter initial
58integer,dimension(nx,ny)           :: mk_gr       !< masque des noeuds flottants
59
60! taux de retrait initial (fonction du beta)
61real,dimension(nx,ny)              :: retreat0_x  !< vitesse maxi autorisee selon x : valeur initiale
62real,dimension(nx,ny)              :: retreat0_y  !< vitesse maxi autorisee selon y : valeur initiale
63real,dimension(nx,ny)              :: retreat0_SW !< vitesse maxi autorisee / SW-NE : valeur initiale  (noeuds /)
64real,dimension(nx,ny)              :: retreat0_SE !< vitesse maxi autorisee \ SE-NW : valeur initiale  (noeuds \)
65
66! taux de retrait  fonction du temps
67real,dimension(nx,ny)              :: retreat_x   !< vitesse maxi autorisee selon x : module par le temps
68real,dimension(nx,ny)              :: retreat_y   !< vitesse maxi autorisee selon x : module par le temps
69real,dimension(nx,ny)              :: retreat_SW  !< vitesse maxi autorisee / SW-NE : valeur initiale  (noeuds /)
70real,dimension(nx,ny)              :: retreat_SE  !< vitesse maxi autorisee \ SE-NW : module par le temps  (noeuds \)
71
72
73
74! tableaux de travail toutes directions
75real,dimension(nx,ny)              :: travail_centre !< valeur sur les noeuds centres
76real,dimension(nx,ny)              :: travail_mx     !< valeur sur noeuds >
77real,dimension(nx,ny)              :: travail_my     !< valeur sur noeuds ^
78real,dimension(nx,ny)              :: travail_SW     !< valeur sur noeuds /
79real,dimension(nx,ny)              :: travail_SE     !< valeur sur noeuds \
80
81
82real,dimension(nx,ny)              :: epsmax         !< taux maxi de deformation
83
84end module declar_proto_recul
Note: See TracBrowser for help on using the repository browser.