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

Last change on this file since 5005 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 
1c
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
8c
9       pi = 4.0 * ATAN(1.0)
10       degres = 180./pi
11       rlong = rlonu * degres
12       rlatg = rlatu * degres
13c
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)
19c
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)
24c
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)
29c
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)
34c
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)
39c
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)
44c
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)
49c
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)
54c
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)
59c
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)
64c
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)
69c
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)
74c
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)
79c
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)
84c
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)
89c
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)
94c
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)
99c
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)
104c
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)
109c
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)
114c
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)
119c
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)
124c
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)
129c
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)
134c
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)
139c
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)
144c
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)
149c
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)
154c
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)
159c
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)
164c
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)
169c
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)
174c
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)
179c
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)
184c
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)
189c
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)
194c
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)
199c
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)
204c
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)
209c
210c=================================================================
211c
212         CALL histend(nid_ctesGCM)
213c
214c=================================================================
Note: See TracBrowser for help on using the repository browser.