source: trunk/LMDZ.TITAN.old/libf/phytitan/ini_histday.h @ 3709

Last change on this file since 3709 was 1530, checked in by emillour, 10 years ago

Venus and Titan GCMs:
Updates in the physics to keep up with updates in LMDZ5 (up to
LMDZ5 trunk, rev 2350) concerning dynamics/physics separation:

  • Adapted makelmdz and makelmdz_fcm script to stop if trying to compile 1d model or newstart or start2archive in parallel.
  • got rid of references to "dimensions.h" in physics. Within physics packages, use nbp_lon (=iim), nbp_lat (=jjmp1) and nbp_lev (=llm) from module mod_grid_phy_lmdz (in phy_common) instead. Only partially done for Titan, because of many hard-coded commons; a necessary first step will be to clean these up (using modules).

EM
File size: 13.6 KB
Line 
1      IF (ok_journe) THEN
2
3         zsto = dtime
4         zout = dtime * REAL(ecrit_day)
5c zsto1: pour des flux radiatifs calcules tous les radpas appels physiq
6         zsto1= dtime * REAL(radpas)
7
8         idayref = day_ref
9         CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
10
11         CALL histbeg_phy("histday.nc", itau_phy, zjulian, dtime,
12     .                 nhori, nid_day)
13
14!$OMP MASTER
15         CALL histvert(nid_day, "presnivs", "Vertical levels", "Pa",
16     .                 klev, presnivs, nvert)
17
18c-------------------------------------------------------
19      IF(lev_histday.GE.1) THEN
20
21ccccccccccccc 2D fields, invariables
22
23         CALL histdef(nid_day, "phis", "Surface geop. height", "-",
24     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 
25     .                "once",  zsto,zout)
26
27         CALL histdef(nid_day, "aire", "Grid area", "-",
28     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 
29     .                "once",  zsto,zout)
30
31ccccccc axe Ls
32         CALL histdef(nid_day, "ls", "Solar longitude", "degrees",
33     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 
34     .                "ave(X)", zsto,zout)
35
36ccccccccccccc 2D fields, variables
37
38         CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
39     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 
40     .                "ave(X)", zsto,zout)
41
42         CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
43     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 
44     .                "ave(X)", zsto,zout)
45
46c        CALL histdef(nid_day, "ue", "Zonal energy transport", "-",
47c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 
48c    .                "ave(X)", zsto,zout)
49
50c        CALL histdef(nid_day, "ve", "Merid energy transport", "-",
51c     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 
52c     .                "ave(X)", zsto,zout)
53
54      ENDIF !lev_histday.GE.1
55
56c-------------------------------------------------------
57      IF(lev_histday.GE.2) THEN
58
59ccccccccccccc 3D fields, basics
60
61         CALL histdef(nid_day, "temp", "Air temperature", "K",
62     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
63     .                "ave(X)", zsto,zout)
64
65         CALL histdef(nid_day, "pres", "Air pressure", "Pa",
66     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
67     .                "ave(X)", zsto,zout)
68
69         CALL histdef(nid_day, "geop", "Geopotential height", "m",
70     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
71     .                "ave(X)", zsto,zout)
72
73         CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
74     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
75     .                "ave(X)", zsto,zout)
76
77         CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
78     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
79     .                "ave(X)", zsto,zout)
80
81         CALL histdef(nid_day, "vitw", "Vertical wind", "Pa/s",
82     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
83     .                "ave(X)", zsto,zout)
84
85         CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
86     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 
87     .                "ave(X)", zsto1,zout)
88
89         CALL histdef(nid_day, "duvdf", "Boundary-layer dU", "m/s2",
90     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
91     .                "ave(X)", zsto,zout)
92
93         CALL histdef(nid_day, "dudyn", "Dynamics dU", "m/s2",
94     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
95     .                "ave(X)", zsto,zout)
96
97cccccccccccccccccc  Tracers
98
99         if (iflag_trac.eq.1) THEN
100          if (microfi.ge.1) then
101c           DO iq=1,nmicro
102c             CALL histdef(nid_day, tname(iq), ttext(iq), "n/m2",
103c     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
104c     .                "ave(X)", zsto,zout)
105c           ENDDO
106             CALL histdef(nid_day, "qaer","nb tot aer" , "n/m2",
107     .                    nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
108     .                    "ave(X)", zsto,zout)
109
110            if (clouds.eq.1) then
111             CALL histdef(nid_day, "qnoy","nb tot noy" , "n/m2",
112     .                    nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
113     .                    "ave(X)", zsto,zout)
114             CALL histdef(nid_day, "qgl1","V tot gl1" , "m3/m2",
115     .                    nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
116     .                    "ave(X)", zsto,zout)
117             CALL histdef(nid_day, "qgl2","V tot gl2" , "m3/m2",
118     .                    nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
119     .                    "ave(X)", zsto,zout)
120             CALL histdef(nid_day, "qgl3","V tot gl3" , "m3/m2",
121     .                    nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
122     .                    "ave(X)", zsto,zout)
123c--------------
124c ----- SATURATION ESP NUAGES
125               CALL histdef(nid_day,"ch4sat", "saturation CH4", "--",
126     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
127     .                "ave(X)", zsto,zout)
128               CALL histdef(nid_day,"c2h6sat", "saturation C2H6", "--",
129     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
130     .                "ave(X)", zsto,zout)
131               CALL histdef(nid_day,"c2h2sat", "saturation C2H2", "--",
132     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
133     .                "ave(X)", zsto,zout)
134c --------------
135c ----- RESERVOIR DE SURFACE
136               CALL histdef(nid_day, "reserv", "Reservoir surface","m",
137     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
138     .                "ave(X)", zsto,zout)
139c --------------
140c ----- ECHANGE GAZ SURF/ATM (evaporation)
141               CALL histdef(nid_day, "evapch4", "Evaporation CH4","m",
142     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
143     .                "ave(X)", zsto,zout)
144c --------------
145c ----- PRECIPITATIONS (precipitations cumulatives)
146               CALL histdef(nid_day,"prech4","Precip CH4","m",
147     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
148     .                "ave(X)", zsto,zout)
149               CALL histdef(nid_day,"prec2h6","Precip C2H6",
150     .                "m",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
151     .                "ave(X)", zsto,zout)
152               CALL histdef(nid_day,"prec2h2","Precip C2H2",
153     .                "m",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
154     .                "ave(X)", zsto,zout)
155               CALL histdef(nid_day,"prenoy","Precip NOY",
156     .                "um/s",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
157     .                "ave(X)", zsto,zout)
158               CALL histdef(nid_day,"preaer","Precip AER",
159     .                "um/s",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
160     .                "ave(X)", zsto,zout)
161c --------------
162c ----- FLUX GLACE
163               CALL histdef(nid_day,"flxgl1", "flux gl CH4",
164     .              "kg/m2/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
165     .              "ave(X)", zsto,zout)
166               CALL histdef(nid_day,"flxgl2", "flux gl C2H6",
167     .              "kg/m2/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
168     .              "ave(X)", zsto,zout)
169               CALL histdef(nid_day,"flxgl3", "flux gl C2H2",
170     .              "kg/m2/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
171     .              "ave(X)", zsto,zout)
172c --------------
173c ----- RAYON DES GOUTTES
174               CALL histdef(nid_day,"rcldbar", "rayon moyen goutte",
175     .                "m",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
176     .                "ave(X)", zsto,zout)
177            endif
178          endif
179c --------------
180c ----- TRACEURS CHIMIQUES
181          if (nmicro.lt.nqmax) then
182           DO iq=nmicro+1,nqmax
183         CALL histdef(nid_day, tname(iq), ttext(iq), "ppm",
184     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
185     .                "ave(X)", zsto,zout)
186           ENDDO
187          endif
188         endif
189
190      ENDIF !lev_histday.GE.2
191
192c-------------------------------------------------------
193      IF(lev_histday.GE.3) THEN
194
195cccccccccccccccccc  Radiative transfer
196
197c 2D
198
199         CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
200     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 
201     .                "ave(X)", zsto1,zout)
202
203         CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2",
204     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 
205     .                "ave(X)", zsto1,zout)
206
207         CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2",
208     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 
209     .                "ave(X)", zsto1,zout)
210
211c 3D
212
213         CALL histdef(nid_day, "SWnet", "Net SW flux","W/m2",
214     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
215     .                32, "ave(X)", zsto1,zout)
216
217         CALL histdef(nid_day, "LWnet", "Net LW flux","W/m2",
218     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
219     .                32, "ave(X)", zsto1,zout)
220
221c --------------
222c ----- OPACITE BRUME
223         DO k=7,NSPECV,10
224           write(str2,'(i2.2)') k
225         CALL histdef(nid_day,"thv"//str2,"Haze Opa Vis",
226     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
227     .                "ave(X)",zsto1,zout)
228         ENDDO
229
230         DO k=8,NSPECI,10
231           write(str2,'(i2.2)') k
232         CALL histdef(nid_day,"thi"//str2,"Haze Opa IR",
233     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
234     .                "ave(X)",zsto1,zout)
235         ENDDO
236
237c --------------
238c ----- EXTINCTION BRUME
239         DO k=7,NSPECV,10
240           write(str2,'(i2.2)') k
241         CALL histdef(nid_day,"khv"//str2,"Haze ext Vis ",
242     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
243     .                "ave(X)",zsto1,zout)
244         ENDDO
245
246         DO k=8,NSPECI,10
247           write(str2,'(i2.2)') k
248         CALL histdef(nid_day,"khi"//str2,"Haze ext IR ",
249     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
250     .                "ave(X)",zsto1,zout)
251         ENDDO
252
253c --------------
254c ----- OPACITE GAZ
255         DO k=7,NSPECV,10
256           write(str2,'(i2.2)') k
257         CALL histdef(nid_day,"tgv"//str2,"Gas Opa Vis",
258     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
259     .                "ave(X)",zsto1,zout)
260         ENDDO
261
262         DO k=8,NSPECI,10
263           write(str2,'(i2.2)') k
264         CALL histdef(nid_day,"tgi"//str2,"Gas Opa IR",
265     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
266     .                "ave(X)",zsto1,zout)
267         ENDDO
268
269c --------------
270c ----- EXTINCTION GAZ
271         DO k=7,NSPECV,10
272           write(str2,'(i2.2)') k
273         CALL histdef(nid_day,"kgv"//str2,"Gas ext Vis ",
274     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
275     .                "ave(X)",zsto1,zout)
276         ENDDO
277
278         DO k=8,NSPECI,10
279           write(str2,'(i2.2)') k
280         CALL histdef(nid_day,"kgi"//str2,"Gas ext IR ",
281     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
282     .                "ave(X)",zsto1,zout)
283         ENDDO
284
285c --------------
286c ----- OPACITE NUAGES
287         if (clouds.eq.1) then
288           CALL histdef(nid_day,"tcld","Cld Opa proxy",
289     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
290     .                "ave(X)",zsto,zout)
291
292c --------------
293c ----- EXTINCTION NUAGES
294           CALL histdef(nid_day,"kcld","Cld Ext proxy",
295     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
296     .                "ave(X)",zsto,zout)
297         endif
298
299      ENDIF !lev_histday.GE.3
300
301c-------------------------------------------------------
302      IF(lev_histday.GE.4) THEN
303
304         CALL histdef(nid_day, "dtdyn", "Dynamics dT", "K/s",
305     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
306     .                "ave(X)", zsto,zout)
307
308         CALL histdef(nid_day, "dtphy", "Physics dT", "K/s",
309     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
310     .                "ave(X)", zsto,zout)
311
312         CALL histdef(nid_day, "dtvdf", "Boundary-layer dT", "K/s",
313     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
314     .                "ave(X)", zsto,zout)
315
316         CALL histdef(nid_day, "dtajs", "Dry adjust. dT", "K/s",
317     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
318     .                "ave(X)", zsto,zout)
319
320         CALL histdef(nid_day, "dtswr", "SW radiation dT", "K/s",
321     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
322     .                "ave(X)", zsto,zout)
323
324         CALL histdef(nid_day, "dtlwr", "LW radiation dT", "K/s",
325     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
326     .                "ave(X)", zsto,zout)
327
328c        CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s",
329c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
330c    .                "ave(X)", zsto,zout)
331
332      ENDIF !lev_histday.GE.4
333
334c-------------------------------------------------------
335      IF(lev_histday.GE.5) THEN
336
337
338c        call histdef(nid_day, "taux", 
339c    $         "Zonal wind stress", "Pa", 
340c    $         nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
341c    $         "ave(X)", zsto,zout)
342
343c        call histdef(nid_day, "tauy", 
344c    $         "Meridional xind stress", "Pa", 
345c    $         nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
346c    $         "ave(X)", zsto,zout)
347
348c        CALL histdef(nid_day, "cdrm", "Momentum drag coef.", "-",
349c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 
350c    .                "ave(X)", zsto,zout)
351
352c        CALL histdef(nid_day, "cdrh", "Heat drag coef.", "-",
353c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 
354c    .                "ave(X)", zsto,zout)
355
356      ENDIF !lev_histday.GE.5
357c-------------------------------------------------------
358
359         CALL histend(nid_day)
360
361      ENDIF ! fin de test sur ok_journe
Note: See TracBrowser for help on using the repository browser.