source: trunk/LMDZ.TITAN/libf/phytitan/ini_histmth.h @ 149

Last change on this file since 149 was 106, checked in by slebonnois, 14 years ago

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