1 | c |
---|
2 | c $Id: ini_bilKP_ave.h 3102 2017-12-03 20:27:42Z acozic $ |
---|
3 | c |
---|
4 | IF (ok_journe) THEN |
---|
5 | c |
---|
6 | zsto = dtime |
---|
7 | zout = ecrit_day |
---|
8 | typeval=tave |
---|
9 | c |
---|
10 | idayref = day_ref |
---|
11 | CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) |
---|
12 | c |
---|
13 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) |
---|
14 | cym DO i = 1, iim |
---|
15 | cym zx_lon(i,1) = rlon(i+1) |
---|
16 | cym zx_lon(i,jjmp1) = rlon(i+1) |
---|
17 | cym ENDDO |
---|
18 | DO ll=1,klev |
---|
19 | znivsig(ll)=REAL(ll) |
---|
20 | ENDDO |
---|
21 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) |
---|
22 | cym write(*,*)'zx_lon = ',zx_lon(:,1) |
---|
23 | cym write(*,*)'zx_lat = ',zx_lat(1,:) |
---|
24 | cym CALL histbeg("histbilKP_ave", iim,zx_lon(:,1), jjmp1, |
---|
25 | cym . zx_lat(1,:), |
---|
26 | cym . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, |
---|
27 | cym . nhori, nid_bilKPave) |
---|
28 | CALL histbeg_phy("histbilKP_ave", itau_phy, zjulian, dtime, |
---|
29 | . nhori, nid_bilKPave) |
---|
30 | |
---|
31 | write(*,*)'Journee ', itau_phy, zjulian |
---|
32 | CALL histvert(nid_bilKPave, "presnivs", |
---|
33 | . "Vertical levels","mb", |
---|
34 | . klev, presnivs/100., nvert) |
---|
35 | c |
---|
36 | c |
---|
37 | c Champs 3D: |
---|
38 | c |
---|
39 | CALL histdef(nid_bilKPave,"ue", |
---|
40 | . "Zonal energy transport","-", |
---|
41 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
42 | . typeval, zsto,zout) |
---|
43 | c |
---|
44 | CALL histdef(nid_bilKPave,"ve", |
---|
45 | . "Merid energy transport","-", |
---|
46 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
47 | . typeval, zsto,zout) |
---|
48 | c |
---|
49 | CALL histdef(nid_bilKPave,"uq", |
---|
50 | . "Zonal humidity transport","-", |
---|
51 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
52 | . typeval, zsto,zout) |
---|
53 | c |
---|
54 | CALL histdef(nid_bilKPave,"vq", |
---|
55 | . "Merid humidity transport","-", |
---|
56 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
57 | . typeval, zsto,zout) |
---|
58 | c |
---|
59 | c Champs 3D: |
---|
60 | c |
---|
61 | CALL histdef(nid_bilKPave,"temp", |
---|
62 | . "Air temperature","K", |
---|
63 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
64 | . typeval, zsto,zout) |
---|
65 | c |
---|
66 | CALL histdef(nid_bilKPave,"ovap", |
---|
67 | . "Specific humidity","Kg/Kg", |
---|
68 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
69 | . typeval, zsto,zout) |
---|
70 | c |
---|
71 | CALL histdef(nid_bilKPave,"geop", |
---|
72 | . "Geopotential height","m", |
---|
73 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
74 | . typeval, zsto,zout) |
---|
75 | c |
---|
76 | CALL histdef(nid_bilKPave,"vitu", |
---|
77 | . "Zonal wind","m/s", |
---|
78 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
79 | . typeval, zsto,zout) |
---|
80 | c |
---|
81 | CALL histdef(nid_bilKPave,"vitv", |
---|
82 | . "Meridional wind","m/s", |
---|
83 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
84 | . typeval, zsto,zout) |
---|
85 | c |
---|
86 | CALL histdef(nid_bilKPave, "vitw", |
---|
87 | . "Vertical wind", "m/s", |
---|
88 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
89 | . typeval, zsto,zout) |
---|
90 | c |
---|
91 | CALL histdef(nid_bilKPave, "pres", |
---|
92 | . "Inter-Layer Air pressure", |
---|
93 | . "Pa", |
---|
94 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
95 | . typeval, zsto,zout) |
---|
96 | c |
---|
97 | CALL histdef(nid_bilKPave, "play", |
---|
98 | . "Mean-Layer Air pressure", |
---|
99 | . "Pa", |
---|
100 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
101 | . typeval, zsto,zout) |
---|
102 | c |
---|
103 | CALL histdef(nid_bilKPave, "oliq", |
---|
104 | . "Liquid water content", |
---|
105 | . "kg/kg", |
---|
106 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
107 | . typeval, zsto,zout) |
---|
108 | c |
---|
109 | CALL histdef(nid_bilKPave, "dtdyn", |
---|
110 | . "Dynamics dT", "K/s", |
---|
111 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
112 | . typeval, zsto,zout) |
---|
113 | c |
---|
114 | CALL histdef(nid_bilKPave, "dqdyn", |
---|
115 | . "Dynamics dQ", "Kg/Kg/s", |
---|
116 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
117 | . typeval, zsto,zout) |
---|
118 | c |
---|
119 | CALL histdef(nid_bilKPave, "dtcon", |
---|
120 | . "Convection dT", "K/s", |
---|
121 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
122 | . typeval, zsto,zout) |
---|
123 | c |
---|
124 | CALL histdef(nid_bilKPave, "ducon", |
---|
125 | . "Convection du", "m/s2", |
---|
126 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
127 | . typeval, zsto,zout) |
---|
128 | c |
---|
129 | CALL histdef(nid_bilKPave, "dvcon", |
---|
130 | . "Convection dv", "m/s2", |
---|
131 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
132 | . typeval, zsto,zout) |
---|
133 | c |
---|
134 | CALL histdef(nid_bilKPave,"dqcon", |
---|
135 | . "Convection dQ","Kg/Kg/s", |
---|
136 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
137 | . typeval, zsto,zout) |
---|
138 | c |
---|
139 | CALL histdef(nid_bilKPave, "dtlsc", |
---|
140 | . "Condensation dT", "K/s", |
---|
141 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
142 | . typeval, zsto,zout) |
---|
143 | c |
---|
144 | CALL histdef(nid_bilKPave,"dqlsc", |
---|
145 | . "Condensation dQ","Kg/Kg/s", |
---|
146 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
147 | . typeval, zsto,zout) |
---|
148 | c |
---|
149 | CALL histdef(nid_bilKPave,"dtvdf", |
---|
150 | . "Boundary-layer dT","K/s", |
---|
151 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
152 | . typeval, zsto,zout) |
---|
153 | c |
---|
154 | CALL histdef(nid_bilKPave, "dqvdf", |
---|
155 | . "Boundary-layer dQ", |
---|
156 | . "Kg/Kg/s", |
---|
157 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
158 | . typeval, zsto,zout) |
---|
159 | c |
---|
160 | CALL histdef(nid_bilKPave,"dtajs", |
---|
161 | . "Ajustement sec dT","K/s", |
---|
162 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
163 | . typeval, zsto,zout) |
---|
164 | c |
---|
165 | CALL histdef(nid_bilKPave, "dqajs", |
---|
166 | . "Ajustement sec dQ", |
---|
167 | . "Kg/Kg/s", |
---|
168 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
169 | . typeval, zsto,zout) |
---|
170 | c |
---|
171 | CALL histdef(nid_bilKPave, "dteva", |
---|
172 | . "Reevaporation dT", "K/s", |
---|
173 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
174 | . typeval, zsto,zout) |
---|
175 | c |
---|
176 | CALL histdef(nid_bilKPave,"dqeva", |
---|
177 | . "Reevaporation dQ", |
---|
178 | . "Kg/Kg/s", |
---|
179 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
180 | . typeval, zsto,zout) |
---|
181 | |
---|
182 | c |
---|
183 | CALL histdef(nid_bilKPave, "dtswr", |
---|
184 | . "SW radiation dT", "K/s", |
---|
185 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
186 | . typeval, zsto,zout) |
---|
187 | c |
---|
188 | CALL histdef(nid_bilKPave, "dtsw0", |
---|
189 | . "SW radiation dT", "K/s", |
---|
190 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
191 | . typeval, zsto,zout) |
---|
192 | c |
---|
193 | CALL histdef(nid_bilKPave, "dtlwr", |
---|
194 | . "LW radiation dT", "K/s", |
---|
195 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
196 | . typeval, zsto,zout) |
---|
197 | c |
---|
198 | CALL histdef(nid_bilKPave, "dtlw0", |
---|
199 | . "LW radiation dT", "K/s", |
---|
200 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
201 | . typeval, zsto,zout) |
---|
202 | c |
---|
203 | CALL histdef(nid_bilKPave,"duvdf", |
---|
204 | . "Boundary-layer dU","m/s2", |
---|
205 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
206 | . typeval, zsto,zout) |
---|
207 | c |
---|
208 | CALL histdef(nid_bilKPave,"dvvdf", |
---|
209 | . "Boundary-layer dV","m/s2", |
---|
210 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
211 | . typeval, zsto,zout) |
---|
212 | c |
---|
213 | IF (ok_orodr) THEN |
---|
214 | IF (ok_orolf) THEN |
---|
215 | CALL histdef(nid_bilKPave, "duoli", |
---|
216 | . "Orography dU", "m/s2", |
---|
217 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
218 | . typeval, zsto,zout) |
---|
219 | c |
---|
220 | CALL histdef(nid_bilKPave, "dvoli", |
---|
221 | . "Orography dV", "m/s2", |
---|
222 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
223 | . typeval, zsto,zout) |
---|
224 | c |
---|
225 | ENDIF |
---|
226 | ENDIF |
---|
227 | C |
---|
228 | CALL histdef(nid_bilKPave, "duphy", |
---|
229 | . "Physiq dU","-", |
---|
230 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
231 | . typeval, zsto,zout) |
---|
232 | C |
---|
233 | CALL histdef(nid_bilKPave, "dvphy", |
---|
234 | . "Physiq dV","-", |
---|
235 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
236 | . typeval, zsto,zout) |
---|
237 | C |
---|
238 | CALL histdef(nid_bilKPave, "dtphy", |
---|
239 | . "Physiq dT","-", |
---|
240 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
241 | . typeval, zsto,zout) |
---|
242 | C |
---|
243 | CALL histdef(nid_bilKPave, "dqphy", |
---|
244 | . "Physiq dQ","-", |
---|
245 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
246 | . typeval, zsto,zout) |
---|
247 | C |
---|
248 | CALL histdef(nid_bilKPave, "dqlphy", |
---|
249 | . "Physiq dQl","-", |
---|
250 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
251 | . typeval, zsto,zout) |
---|
252 | C |
---|
253 | C |
---|
254 | CALL histend(nid_bilKPave) |
---|
255 | c |
---|
256 | ndex2d = 0 |
---|
257 | ndex3d = 0 |
---|
258 | c |
---|
259 | ENDIF ! fin de test sur ok_journe |
---|