Context Navigation

source: LMDZ4/branches/LMDZ4-dev-20091210/libf/dyn3d/ini_paramLMDZ_dyn.h @ 3406

Last change on this file since 3406 was 1201, checked in by Laurent Fairhead, 15 years ago
Modifications nécessaires a l'inclusion d'un calendrier réaliste. La date courante est calculée dans leapfrog.F et exprimée en Jour Julien (modifié). On en a profité pour faire un peu de ménage dans la gestion des dates du modèle. Dans la physique, on utilise les routines de passages entre calendrier Julien et Gregorien incluses dans IOIPSL pour calculer le nombre de jours écoulés depuis le 1er janvier (pour les conditions aux limites) ou l'equinoxe (pour le calcul de la longitude solaire). Le calcul de l'orbite reprend celui du gcm planétaire (codé par FH) On décide du calendrier à utiliser à l'aide du paramètre calend du run.def. Par défaut celui-ci est à earth_360d LF
Property svn:eol-style set to `native` Property svn:executable set to ``* Property svn:keywords set to `Author Date Id Revision`
File size: 8.7 KB

Line
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, "idissip",
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=================================================================

Note: See TracBrowser for help on using the repository browser.

Download in other formats: