source: LMDZ4/tags/LMDZ4_V2/libf/phylmd/ini_bilKP_ins.h @ 741

Last change on this file since 741 was 741, checked in by (none), 17 years ago

This commit was manufactured by cvs2svn to create tag 'LMDZ4_V2'.
  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 10.4 KB
Line 
1c
2c $Header$
3c
4      IF (ok_journe) THEN
5c
6         zsto = dtime
7         zout = dtime
8         typeval=tinst
9c
10         idayref = day_ref
11         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
12c
13         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
14         DO i = 1, iim
15            zx_lon(i,1) = rlon(i+1)
16            zx_lon(i,jjmp1) = rlon(i+1)
17         ENDDO
18         DO ll=1,klev
19            znivsig(ll)=float(ll)
20         ENDDO
21         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
22         write(*,*)'zx_lon = ',zx_lon(:,1)
23         write(*,*)'zx_lat = ',zx_lat(1,:)
24c
25cIM 280405 BEG
26c
27cIM cf. AM 081204 BEG region
28          imin_ins=1
29          imax_ins=iim
30          jmin_ins=1
31          jmax_ins=jjmp1
32          do i=1,iim-1
33             if(zx_lon(i,1).lt.lonmin_ins) imin_ins=i
34             if(zx_lon(i,1).le.lonmax_ins) imax_ins=i+1
35          enddo
36          do j=1,jjmp1-1
37             if(zx_lat(1,j).ge.latmin_ins) jmax_ins=j
38             if(zx_lat(1,j).gt.latmax_ins) jmin_ins=j
39          enddo
40c
41          print*,'On stoke le fichier bilKP instantanne sur ',
42     s   imin_ins,imax_ins,jmin_ins,jmax_ins
43          print*,'On stoke le fichier bilKP instantanne sur ',
44     s   zx_lon(imin_ins,1),zx_lon(imax_ins,1),
45     s   zx_lat(1,jmin_ins),zx_lat(1,jmax_ins)
46cIM cf. AM 081204 END region
47c
48cIM 280405 END
49c
50         IF(1.EQ.0) THEN
51         CALL histbeg("histbilKP_ins", iim,zx_lon(:,1), jjmp1,
52     .                zx_lat(1,:),
53     .                1,iim,1,jjmp1, itau_phy, zjulian, dtime,
54     .                nhori, nid_bilKPins)
55         ENDIF
56c
57cIM 280405 BEG
58c
59cIM cf. AM 081204 BEG region
60         CALL histbeg("histbilKP_ins", iim,zx_lon(:,1), 
61     .                 jjmp1,zx_lat(1,:),
62     .                 imin_ins,imax_ins-imin_ins+1,
63     .                 jmin_ins,jmax_ins-jmin_ins+1,
64     .                 itau_phy, zjulian, dtime,
65     .                 nhori, nid_bilKPins)
66cIM 081204 END
67c
68cIM 280405 END
69c
70         write(*,*)'Journee ', itau_phy, zjulian
71         CALL histvert(nid_bilKPins, "presnivs",
72     .                "Vertical levels","mb",
73     .                 klev, presnivs/100., nvert)
74c
75c Champs 3D:
76c
77         CALL histdef(nid_bilKPins,"ue",
78     .   "Zonal energy transport","-",
79     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32, 
80     .                typeval, zsto,zout)
81c
82         CALL histdef(nid_bilKPins,"ve",
83     .   "Merid energy transport","-",
84     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32, 
85     .                typeval, zsto,zout)
86c
87         CALL histdef(nid_bilKPins,"uq",
88     .   "Zonal humidity transport","-",
89     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32, 
90     .                typeval, zsto,zout)
91c
92         CALL histdef(nid_bilKPins,"vq",
93     .   "Merid humidity transport","-",
94     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32, 
95     .                typeval, zsto,zout)
96c
97c Champs 3D:
98c
99         CALL histdef(nid_bilKPins, "temp",
100     .   "Air temperature", "K",
101     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
102     .                typeval, zsto,zout)
103c
104         CALL histdef(nid_bilKPins,"ovap",
105     .   "Specific humidity","Kg/Kg",
106     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
107     .                typeval, zsto,zout)
108c
109         CALL histdef(nid_bilKPins,"geop",
110     .   "Geopotential height", "m",
111     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
112     .                typeval, zsto,zout)
113c
114         CALL histdef(nid_bilKPins,"vitu", 
115     .   "Zonal wind", "m/s",
116     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
117     .                typeval, zsto,zout)
118c
119         CALL histdef(nid_bilKPins,"vitv", 
120     .   "Meridional wind", "m/s",
121     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
122     .                typeval, zsto,zout)
123c
124         CALL histdef(nid_bilKPins, "vitw",
125     .   "Vertical wind", "m/s",
126     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
127     .                typeval, zsto,zout)
128c
129         CALL histdef(nid_bilKPins, "pres",
130     .   "Inter-Layer Air pressure",
131     .                "Pa",
132     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
133     .                typeval, zsto,zout)
134c
135         CALL histdef(nid_bilKPins, "play",
136     .   "Mean-Layer Air pressure",
137     .                "Pa",
138     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
139     .                typeval, zsto,zout)
140c
141         CALL histdef(nid_bilKPins, "oliq",
142     .   "Liquid water content", 
143     .                "kg/kg",
144     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
145     .                typeval, zsto,zout)
146c
147         CALL histdef(nid_bilKPins, "dtdyn",
148     .   "Dynamics dT", "K/s",
149     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
150     .                typeval, zsto,zout)
151c
152         CALL histdef(nid_bilKPins, "dqdyn",
153     .   "Dynamics dQ", "Kg/Kg/s",
154     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
155     .                typeval, zsto,zout)
156c
157         CALL histdef(nid_bilKPins, "dtcon",
158     .   "Convection dT", "K/s",
159     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
160     .                typeval, zsto,zout)
161c
162         CALL histdef(nid_bilKPins, "ducon",
163     .   "Convection du", "m/s2",
164     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
165     .                typeval, zsto,zout)
166c
167         CALL histdef(nid_bilKPins, "dvcon",
168     .   "Convection dv", "m/s2",
169     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
170     .                typeval, zsto,zout)
171c
172         CALL histdef(nid_bilKPins,"dqcon",
173     .   "Convection dQ","Kg/Kg/s",
174     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
175     .                typeval, zsto,zout)
176c
177         CALL histdef(nid_bilKPins, "dtlsc",
178     .   "Condensation dT", "K/s",
179     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
180     .                typeval, zsto,zout)
181c
182         CALL histdef(nid_bilKPins,"dqlsc",
183     .   "Condensation dQ","Kg/Kg/s",
184     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
185     .                typeval, zsto,zout)
186c
187         CALL histdef(nid_bilKPins,"dtvdf",
188     .   "Boundary-layer dT","K/s",
189     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
190     .                typeval, zsto,zout)
191c
192         CALL histdef(nid_bilKPins, "dqvdf", 
193     .   "Boundary-layer dQ", 
194     .               "Kg/Kg/s",
195     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
196     .                typeval, zsto,zout)
197c
198         CALL histdef(nid_bilKPins,"dtajs",
199     .   "Ajustement sec dT","K/s",
200     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
201     .                typeval, zsto,zout)
202c
203         CALL histdef(nid_bilKPins,"dqajs",
204     .   "Ajustement sec dQ", 
205     .               "Kg/Kg/s",
206     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
207     .                typeval, zsto,zout)
208c
209         CALL histdef(nid_bilKPins,"dteva",
210     .   "Reevaporation dT","K/s",
211     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
212     .                typeval, zsto,zout)
213c
214         CALL histdef(nid_bilKPins,"dqeva",
215     .   "Reevaporation dQ",
216     .                "Kg/Kg/s",
217     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
218     .                typeval, zsto,zout)
219
220c
221         CALL histdef(nid_bilKPins, "dtswr", 
222     .   "SW radiation dT", "K/s",
223     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
224     .                typeval, zsto,zout)
225c
226         CALL histdef(nid_bilKPins, "dtsw0", 
227     .   "SW radiation dT", "K/s",
228     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
229     .                typeval, zsto,zout)
230c
231         CALL histdef(nid_bilKPins, "dtlwr", 
232     .   "LW radiation dT", "K/s",
233     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
234     .                typeval, zsto,zout)
235c
236         CALL histdef(nid_bilKPins, "dtlw0", 
237     .   "LW radiation dT", "K/s",
238     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
239     .                typeval, zsto,zout)
240c
241         CALL histdef(nid_bilKPins,"duvdf",
242     .   "Boundary-layer dU","m/s2",
243     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
244     .                typeval, zsto,zout)
245c
246         CALL histdef(nid_bilKPins,"dvvdf",
247     .   "Boundary-layer dV","m/s2",
248     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
249     .                typeval, zsto,zout)
250c
251         IF (ok_orodr) THEN
252         IF (ok_orolf) THEN
253         CALL histdef(nid_bilKPins, "duoli", 
254     .   "Orography dU", "m/s2",
255     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
256     .                typeval, zsto,zout)
257c
258         CALL histdef(nid_bilKPins, "dvoli", 
259     .   "Orography dV", "m/s2",
260     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
261     .                typeval, zsto,zout)
262c
263         ENDIF
264         ENDIF
265C
266         CALL histdef(nid_bilKPins, "duphy",
267     .   "Physiq dU","-",
268     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
269     .                typeval, zsto,zout)
270C
271         CALL histdef(nid_bilKPins, "dvphy",
272     .   "Physiq dV","-",
273     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
274     .                typeval, zsto,zout)
275C
276         CALL histdef(nid_bilKPins, "dtphy",
277     .   "Physiq dT","-",
278     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
279     .                typeval, zsto,zout)
280C
281         CALL histdef(nid_bilKPins, "dqphy",
282     .   "Physiq dQ","-",
283     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
284     .                typeval, zsto,zout)
285C
286         CALL histdef(nid_bilKPins, "dqlphy",
287     .   "Physiq dQl","-",
288     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
289     .                typeval, zsto,zout)
290cIM 280405 BEG
291c
292c Champs 2D:
293c
294c u850, v850
295c        DO k=1, nlevSTD
296         DO k=1, 12
297c
298         IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
299         IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
300c
301         IF(bb2.EQ."850") THEN
302c
303          CALL histdef(nid_bilKPins, "u"//bb2,
304     .                 "Zonal wind "//bb2//"mb","m/s",
305     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
306     .                typeval, zsto,zout)
307c
308          CALL histdef(nid_bilKPins, "v"//bb2,
309     .                 "Meridional wind "//bb2//"mb","m/s",
310     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
311     .                typeval, zsto,zout)
312c
313         ENDIF !(bb2.EQ."850") 
314c
315         ENDDO !k=1, 12
316c
317cIM 280405 END
318c
319         CALL histend(nid_bilKPins)
320c
321         ndex2d = 0
322         ndex3d = 0
323c
324      ENDIF ! fin de test sur ok_journe
Note: See TracBrowser for help on using the repository browser.