source: trunk/LMDZ.TITAN/libf/phytitan/ini_histins.h @ 963

Last change on this file since 963 was 902, checked in by slebonnois, 12 years ago

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