source: trunk/libf/phytitan/ini_histins.h @ 124

Last change on this file since 124 was 107, checked in by slebonnois, 14 years ago

SL : juste correction d'une erreur dans sorties des opacites...
File size: 9.8 KB
Line 
1!
2! $Header: /home/cvsroot/LMDZ4/libf/phylmd/ini_histins.h,v 1.1.1.1 2004/05/19 12:53:08 lmdzadmin Exp $
3!
4      IF (ok_instan) THEN
5c
6          zsto1= dtime
7          zsto = dtime * FLOAT(ecrit_ins)
8          zout = dtime * FLOAT(ecrit_ins)
9c
10         idayref = day_ref
11         CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
12c
13         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlond,zx_lon)
14         DO i = 1, iim
15            zx_lon(i,1) = rlond(i+1)
16            zx_lon(i,jjmp1) = rlond(i+1)
17         ENDDO
18         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlatd,zx_lat)
19         CALL histbeg("histins", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
20     .                 1,iim,1,jjmp1, itau_phy, zjulian, zsto1,
21     .                 nhori, nid_ins)
22         write(*,*)'Inst ', itau_phy, zjulian
23         CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb",
24     .                 klev, presnivs, nvert)
25
26c-------------------------------------------------------
27
28      IF(lev_histday.GE.1) THEN
29c
30ccccccccccccc 2D fields, invariables
31c
32         CALL histdef(nid_ins, "phis", "Surface geop. height", "-",
33     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
34     .                "once",  zsto,zout)
35c
36         CALL histdef(nid_ins, "aire", "Grid area", "-",
37     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
38     .                "once",  zsto,zout)
39c
40ccccccc axe Ls
41         CALL histdef(nid_ins, "ls", "Solar longitude", "degrees",
42     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
43     .                "inst(X)", zsto,zout)
44c
45ccccccccccccc 2D fields, variables
46c
47         CALL histdef(nid_ins, "tsol", "Surface Temperature", "K",
48     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
49     .                "inst(X)", zsto,zout)
50c
51         CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa",
52     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
53     .                "inst(X)", zsto,zout)
54c
55c        CALL histdef(nid_ins, "ue", "Zonal energy transport", "-",
56c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
57c    .                "inst(X)", zsto,zout)
58c
59c        CALL histdef(nid_ins, "ve", "Merid energy transport", "-",
60c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
61c    .                "inst(X)", zsto,zout)
62c
63      ENDIF !lev_histday.GE.1
64c
65c-------------------------------------------------------
66      IF(lev_histday.GE.2) THEN
67c
68ccccccccccccc 3D fields, basics
69c
70         CALL histdef(nid_ins, "temp", "Air temperature", "K",
71     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
72     .                "inst(X)", zsto,zout)
73c
74         CALL histdef(nid_ins, "pres", "Air pressure", "Pa",
75     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
76     .                "inst(X)", zsto,zout)
77c
78         CALL histdef(nid_ins, "geop", "Geopotential height", "m",
79     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
80     .                "inst(X)", zsto,zout)
81c
82         CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s",
83     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
84     .                "inst(X)", zsto,zout)
85c
86         CALL histdef(nid_ins, "vitv", "Meridional wind", "m/s",
87     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
88     .                "inst(X)", zsto,zout)
89c
90         CALL histdef(nid_ins, "vitw", "Vertical wind", "Pa/s",
91     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
92     .                "inst(X)", zsto,zout)
93c
94         CALL histdef(nid_ins, "tops", "Solar rad. at TOA", "W/m2",
95     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
96     .                "inst(X)", zsto,zout)
97c
98c        CALL histdef(nid_ins, "duvdf", "Boundary-layer dU", "m/s2",
99c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
100c    .                "inst(X)", zsto,zout)
101c
102c        CALL histdef(nid_ins, "dudyn", "Dynamics dU", "m/s2",
103c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
104c    .                "inst(X)", zsto,zout)
105c
106      ENDIF !lev_histday.GE.2
107c
108c-------------------------------------------------------
109      IF(lev_histday.GE.3) THEN
110c
111cccccccccccccccccc  Tracers
112c
113         if (iflag_trac.eq.1) THEN
114          if (microfi.eq.1) then
115           DO iq=1,nmicro
116         CALL histdef(nid_ins, tname(iq), ttext(iq), "n/m2",
117     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
118     .                "ins(X)", zsto,zout)
119           ENDDO
120          endif
121          if (nmicro.lt.nqmax) then
122           DO iq=nmicro+1,nqmax
123         CALL histdef(nid_ins, tname(iq), ttext(iq), "ppm",
124     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
125     .                "ins(X)", zsto,zout)
126           ENDDO
127          endif
128         endif
129c
130cccccccccccccccccc  Radiative transfer
131c
132c 2D
133c
134         CALL histdef(nid_ins, "topl", "IR rad. at TOA", "W/m2",
135     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
136     .                "inst(X)", zsto,zout)
137c
138         CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2",
139     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
140     .                "inst(X)", zsto,zout)
141c
142         CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2",
143     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
144     .                "inst(X)", zsto,zout)
145c
146c 3D
147c
148         CALL histdef(nid_ins, "SWnet", "Net SW flux","W/m2",
149     .                iim,jjmp1,nhori, klev,1,klev,nvert,
150     .                32, "inst(X)", zsto,zout)
151c
152         CALL histdef(nid_ins, "LWnet", "Net LW flux","W/m2",
153     .                iim,jjmp1,nhori, klev,1,klev,nvert,
154     .                32, "inst(X)", zsto,zout)
155c
156         DO k=7,NSPECV,10
157           write(str1,'(i2.2)') k
158         CALL histdef(nid_ins,"thv"//str1,"Haze Opa Vis",
159     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
160     .                "ins(X)",zsto,zout)
161         ENDDO
162c
163         DO k=7,NSPECV,10
164           write(str1,'(i2.2)') k
165         CALL histdef(nid_ins,"khv"//str1,"Haze ext Vis ",
166     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
167     .                "ins(X)",zsto,zout)
168         ENDDO
169c
170         DO k=7,NSPECV,10
171           write(str1,'(i2.2)') k
172         CALL histdef(nid_ins,"tgv"//str1,"Haze Opa Vis",
173     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
174     .                "ins(X)",zsto,zout)
175         ENDDO
176c
177         DO k=7,NSPECV,10
178           write(str1,'(i2.2)') k
179         CALL histdef(nid_ins,"kgv"//str1,"Haze ext Vis ",
180     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
181     .                "ins(X)",zsto,zout)
182         ENDDO
183c
184         DO k=8,NSPECI,10
185           write(str1,'(i2.2)') k
186         CALL histdef(nid_ins,"thi"//str1,"Haze Opa IR",
187     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
188     .                "ins(X)",zsto,zout)
189         ENDDO
190c
191         DO k=8,NSPECI,10
192           write(str1,'(i2.2)') k
193         CALL histdef(nid_ins,"khi"//str1,"Haze ext IR ",
194     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
195     .                "ins(X)",zsto,zout)
196         ENDDO
197c
198         DO k=8,NSPECI,10
199           write(str1,'(i2.2)') k
200         CALL histdef(nid_ins,"tgi"//str1,"Haze Opa IR",
201     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
202     .                "ins(X)",zsto,zout)
203         ENDDO
204c
205         DO k=8,NSPECI,10
206           write(str1,'(i2.2)') k
207         CALL histdef(nid_ins,"kgi"//str1,"Haze ext IR ",
208     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
209     .                "ins(X)",zsto,zout)
210         ENDDO
211c
212      ENDIF !lev_histday.GE.3
213c
214c-------------------------------------------------------
215      IF(lev_histday.GE.4) THEN
216c
217         CALL histdef(nid_ins, "dtdyn", "Dynamics dT", "K/s",
218     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
219     .                "inst(X)", zsto,zout)
220c
221         CALL histdef(nid_ins, "dtphy", "Physics dT", "K/s",
222     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
223     .                "inst(X)", zsto,zout)
224c
225         CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s",
226     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
227     .                "inst(X)", zsto,zout)
228c
229         CALL histdef(nid_ins, "dtajs", "Dry adjust. dT", "K/s",
230     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
231     .                "inst(X)", zsto,zout)
232c
233         CALL histdef(nid_ins, "dtswr", "SW radiation dT", "K/s",
234     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
235     .                "inst(X)", zsto,zout)
236c
237         CALL histdef(nid_ins, "dtlwr", "LW radiation dT", "K/s",
238     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
239     .                "inst(X)", zsto,zout)
240c
241c        CALL histdef(nid_ins, "dtec", "Cinetic dissip dT", "K/s",
242c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
243c    .                "inst(X)", zsto,zout)
244c
245c        CALL histdef(nid_ins, "dvvdf", "Boundary-layer dV", "m/s2",
246c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
247c    .                "inst(X)", zsto,zout)
248c
249      ENDIF !lev_histday.GE.4
250c
251c-------------------------------------------------------
252      IF(lev_histday.GE.5) THEN
253c
254c
255c        call histdef(nid_ins, "taux", 
256c    $         "Zonal wind stress", "Pa", 
257c    $         iim,jjmp1,nhori, 1,1,1, -99, 32,
258c    $         "inst(X)", zsto,zout)
259c
260c        call histdef(nid_ins, "tauy", 
261c    $         "Meridional xind stress", "Pa", 
262c    $         iim,jjmp1,nhori, 1,1,1, -99, 32,
263c    $         "inst(X)", zsto,zout)
264c
265c        CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-",
266c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
267c    .                "inst(X)", zsto,zout)
268c
269c        CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-",
270c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
271c    .                "inst(X)", zsto,zout)
272c
273      ENDIF !lev_histday.GE.5
274c-------------------------------------------------------
275
276         CALL histend(nid_ins)
277c
278         ndex2d = 0
279         ndex3d = 0
280c
281      ENDIF
Note: See TracBrowser for help on using the repository browser.