My Project
 All Classes Files Functions Variables Macros
ini_paramLMDZ_dyn.h
Go to the documentation of this file.
1 c
2  dt_cum = dtvr*day_step
3 
4 ! zan = annee_ref
5 ! dayref = day_ref
6 ! CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)
7  tau0 = itau_dyn
8 c
9  pi = 4.0 * ATAN(1.0)
10  degres = 180./pi
11  rlong = rlonu * degres
12  rlatg = rlatu * degres
13 c
14  CALL histbeg("paramLMDZ_dyn.nc",
15  . iip1,rlong, jjp1,rlatg,
16  . 1,1,1,1,
17  . tau0, jD_ref+jH_ref , dt_cum,
18  . thoriid, nid_ctesGCM)
19 c
20  CALL histdef(nid_ctesGCM, "prt_level",
21  . "Niveau impression debuggage dynamique",
22  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
23  . "once", dt_cum,dt_cum)
24 c
25  CALL histdef(nid_ctesGCM, "dayref",
26  . "Jour de l etat initial ( = 350 si 20 Decembre par ex.)",
27  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
28  . "once", dt_cum,dt_cum)
29 c
30  CALL histdef(nid_ctesGCM, "anneeref",
31  . "Annee de l etat initial",
32  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
33  . "once", dt_cum,dt_cum)
34 c
35  CALL histdef(nid_ctesGCM, "anneelim",
36  . "Annee du fichier limitxxxx.nc si ok_limitvrai =y",
37  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
38  . "once", dt_cum,dt_cum)
39 c
40  CALL histdef(nid_ctesGCM, "raz_date",
41  . "Remise a zero (raz) date init.: 0 pas de raz;1=date gcm.def",
42  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
43  . "once", dt_cum,dt_cum)
44 c
45  CALL histdef(nid_ctesGCM, "nday",
46  . "Nombre de jours d integration",
47  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
48  . "once", dt_cum,dt_cum)
49 c
50  CALL histdef(nid_ctesGCM, "day_step",
51  . "nombre de pas par jour pour dt = 1 min",
52  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
53  . "once", dt_cum,dt_cum)
54 c
55  CALL histdef(nid_ctesGCM, "iperiod",
56  . "periode pour le pas Matsuno (en pas de temps)",
57  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
58  . "once", dt_cum,dt_cum)
59 c
60  CALL histdef(nid_ctesGCM, "iapp_tracvl",
61  . "frequence du groupement des flux (en pas de temps)",
62  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
63  . "once", dt_cum,dt_cum)
64 c
65  CALL histdef(nid_ctesGCM, "iconser",
66  . "periode de sortie des variables de controle (en pas de temps)",
67  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
68  . "once", dt_cum,dt_cum)
69 c
70  CALL histdef(nid_ctesGCM, "iecri",
71  . "periode d ecriture du fichier histoire (en jour)",
72  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
73  . "once", dt_cum,dt_cum)
74 c
75  CALL histdef(nid_ctesGCM, "periodav",
76  . "periode de stockage fichier histmoy (en jour)",
77  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
78  . "once", dt_cum,dt_cum)
79 c
80  CALL histdef(nid_ctesGCM, "dissip_period",
81  . "periode de la dissipation (en pas) ... a completer",
82  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
83  . "once", dt_cum,dt_cum)
84 c
85  CALL histdef(nid_ctesGCM, "lstardis",
86  . "choix de l operateur de dissipation: 1= star,0=non-star ??",
87  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
88  . "once", dt_cum,dt_cum)
89 c
90  CALL histdef(nid_ctesGCM, "nitergdiv",
91  . "nombre d iterations de l operateur de dissipation gradiv",
92  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
93  . "once", dt_cum,dt_cum)
94 c
95  CALL histdef(nid_ctesGCM, "nitergrot",
96  . "nombre d iterations de l operateur de dissipation nxgradrot",
97  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
98  . "once", dt_cum,dt_cum)
99 c
100  CALL histdef(nid_ctesGCM, "niterh",
101  . "nombre d iterations de l operateur de dissipation divgrad",
102  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
103  . "once", dt_cum,dt_cum)
104 c
105  CALL histdef(nid_ctesGCM, "tetagdiv",
106  ."temps dissipation des + petites long. d ondes pour u,v (gradiv)",
107  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
108  . "once", dt_cum,dt_cum)
109 c
110  CALL histdef(nid_ctesGCM, "tetagrot",
111  ."temps diss. des + petites long. d ondes pour u,v (nxgradrot)",
112  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
113  . "once", dt_cum,dt_cum)
114 c
115  CALL histdef(nid_ctesGCM, "tetatemp",
116  ."temps diss. des + petites long. d ondes pour h (divgrad)",
117  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
118  . "once", dt_cum,dt_cum)
119 c
120  CALL histdef(nid_ctesGCM, "coefdis",
121  ."coefficient pour gamdissip",
122  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
123  . "once", dt_cum,dt_cum)
124 c
125  CALL histdef(nid_ctesGCM, "purmats",
126  ."Choix schema integration temporel: 1=Matsuno,0=Matsuno-leapfrog",
127  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
128  . "once", dt_cum,dt_cum)
129 c
130  CALL histdef(nid_ctesGCM, "ok_guide",
131  ."Guidage: 1=true ,0=false",
132  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
133  . "once", dt_cum,dt_cum)
134 c
135  CALL histdef(nid_ctesGCM, "true_calendar",
136  ."Choix du calendrier",
137  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
138  . "once", dt_cum,dt_cum)
139 c
140  CALL histdef(nid_ctesGCM, "guide_calend",
141  ."Guidage calendrier gregorien: 1=oui ,0=non",
142  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
143  . "once", dt_cum,dt_cum)
144 c
145  CALL histdef(nid_ctesGCM, "iflag_phys",
146  ."Permet de faire tourner le modele sans physique: 1=avec ,0=sans",
147  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
148  . "once", dt_cum,dt_cum)
149 c
150  CALL histdef(nid_ctesGCM, "iphysiq",
151  ."Periode de la physique en pas de temps de la dynamique",
152  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
153  . "once", dt_cum,dt_cum)
154 c
155  CALL histdef(nid_ctesGCM, "clon",
156  ."longitude en degres du centre du zoom",
157  . "deg",iip1,jjp1,thoriid, 1,1,1, -99, 32,
158  . "once", dt_cum,dt_cum)
159 c
160  CALL histdef(nid_ctesGCM, "clat",
161  ."latitude en degres du centre du zoom",
162  . "deg",iip1,jjp1,thoriid, 1,1,1, -99, 32,
163  . "once", dt_cum,dt_cum)
164 c
165  CALL histdef(nid_ctesGCM, "grossismx",
166  ."facteur de grossissement du zoom, selon la longitude",
167  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
168  . "once", dt_cum,dt_cum)
169 c
170  CALL histdef(nid_ctesGCM, "grossismy",
171  ."facteur de grossissement du zoom, selon la latitude",
172  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
173  . "once", dt_cum,dt_cum)
174 c
175  CALL histdef(nid_ctesGCM, "fxyhypb",
176  ."Fonction f(y) hyperbolique si true=1, sinusoidale si false=0",
177  . "-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
178  . "once", dt_cum,dt_cum)
179 c
180  CALL histdef(nid_ctesGCM, "dzoomx",
181  ."extension en longitude de la zone du zoom (fraction zone totale)"
182  . ,"-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
183  . "once", dt_cum,dt_cum)
184 c
185  CALL histdef(nid_ctesGCM, "dzoomy",
186  ."extension en latitude de la zone du zoom (fraction zone totale)"
187  . ,"-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
188  . "once", dt_cum,dt_cum)
189 c
190  CALL histdef(nid_ctesGCM, "taux",
191  ."raideur du zoom en X"
192  . ,"-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
193  . "once", dt_cum,dt_cum)
194 c
195  CALL histdef(nid_ctesGCM, "tauy",
196  ."raideur du zoom en Y"
197  . ,"-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
198  . "once", dt_cum,dt_cum)
199 c
200  CALL histdef(nid_ctesGCM, "ysinus",
201  ."ysinus=1: Ftion f(y) avec y=Sin(latit.)/ ysinus=0: y = latit"
202  . ,"-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
203  . "once", dt_cum,dt_cum)
204 c
205  CALL histdef(nid_ctesGCM, "ip_ebil_dyn",
206  ."PRINTlevel for energy conservation diag.; 0/1= pas de print,
207  . 2= print","-",iip1,jjp1,thoriid, 1,1,1, -99, 32,
208  . "once", dt_cum,dt_cum)
209 c
210 c=================================================================
211 c
212  CALL histend(nid_ctesGCM)
213 c
214 c=================================================================