source: LMDZ5/branches/LF-private/libf/dyn3d/ini_paramLMDZ_dyn.h @ 2109

Last change on this file since 2109 was 1502, checked in by jghattas, 14 years ago

Modifications for variable idissip :

  • Changed name of variable idissip to dissip_period everywhere to be compatible with old .def files.
  • This variable was before read from physiq.def but the value was overwritten by a calculation in inidissip. Now, if dissip_period=0 calculation is done as before(default). Else the value from physiq.def is used directly.
  • leapfrog : added "AND NOT forward" at line 284 (dyn3d) and line 363(dyn3dpar) necessare if dissip_period not a multiple by iperiod.
  • iniconst : removed calculation of dtdiss (calculation done in inidissip)

FC, JG
  • 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, "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)
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.