source: LMDZ6/branches/Amaury_dev/libf/phylmd/ini_bilKP_ins.h @ 5136

Last change on this file since 5136 was 5116, checked in by abarral, 4 months ago

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