source: trunk/libf/phytitan/ini_histday.h @ 6

Last change on this file since 6 was 3, checked in by slebonnois, 15 years ago

Creation de repertoires:

  • chantiers : pour communiquer sur nos projets de modifs
  • documentation : pour stocker les docs

Ajout de:

  • libf/phytitan : physique de Titan
  • libf/chimtitan: chimie de Titan
  • libf/phyvenus : physique de Venus
File size: 7.8 KB
RevLine 
[3]1!
2! $Header: /home/cvsroot/LMDZ4/libf/phylmd/ini_histday.h,v 1.2 2005/01/27 10:06:12 fairhead Exp $
3!
4      IF (ok_journe) THEN
5c
6         zsto = dtime
7         zout = dtime * FLOAT(ecrit_day)
8c zsto1: pour des flux radiatifs calcules tous les radpas appels physiq
9         zsto1= dtime * FLOAT(radpas)
10c
11         idayref = day_ref
12         z4julian = zjulian
13         CALL ymds2ju(annee_ref, 1, idayref, zero, z4julian)
14c
15         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlond,zx_lon)
16         DO i = 1, iim
17            zx_lon(i,1) = rlond(i+1)
18            zx_lon(i,jjmp1) = rlond(i+1)
19         ENDDO
20         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlatd,zx_lat)
21         z4_lon = zx_lon
22         z4_lat = zx_lat
23         zpresnivs=presnivs
24         CALL histbeg("histday", iim,z4_lon(:,1), jjmp1,z4_lat(1,:),
25     .                 1,iim,1,jjmp1, itau_phy, z4julian, zsto, 
26     .                 nhori, nid_day)
27         write(*,*)'Journee ', itau_phy, zjulian
28         CALL histvert(nid_day, "presnivs", "Vertical levels", "mb",
29     .                 klev, zpresnivs, nvert)
30
31c-------------------------------------------------------
32      IF(lev_histday.GE.1) THEN
33c
34ccccccccccccc 2D fields, invariables
35c
36         CALL histdef(nid_day, "phis", "Surface geop. height", "-",
37     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
38     .                "once",  zsto,zout)
39c
40         CALL histdef(nid_day, "aire", "Grid area", "-",
41     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
42     .                "once",  zsto,zout)
43c
44ccccccc axe Ls
45         CALL histdef(nid_day, "ls", "Solar longitude", "degrees",
46     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
47     .                "ave(X)", zsto,zout)
48c
49ccccccccccccc 2D fields, variables
50c
51         CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
52     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
53     .                "ave(X)", zsto,zout)
54c
55         CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
56     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
57     .                "ave(X)", zsto,zout)
58c
59c        CALL histdef(nid_day, "ue", "Zonal energy transport", "-",
60c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
61c    .                "ave(X)", zsto,zout)
62c
63c        CALL histdef(nid_day, "ve", "Merid energy transport", "-",
64c     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
65c     .                "ave(X)", zsto,zout)
66c
67      ENDIF !lev_histday.GE.1
68c
69c-------------------------------------------------------
70      IF(lev_histday.GE.2) THEN
71c
72ccccccccccccc 3D fields, basics
73c
74         CALL histdef(nid_day, "temp", "Air temperature", "K",
75     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
76     .                "ave(X)", zsto,zout)
77c
78         CALL histdef(nid_day, "pres", "Air pressure", "Pa",
79     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
80     .                "ave(X)", zsto,zout)
81c
82         CALL histdef(nid_day, "geop", "Geopotential height", "m",
83     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
84     .                "ave(X)", zsto,zout)
85c
86         CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
87     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
88     .                "ave(X)", zsto,zout)
89c
90         CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
91     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
92     .                "ave(X)", zsto,zout)
93c
94         CALL histdef(nid_day, "vitw", "Vertical wind", "Pa/s",
95     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
96     .                "ave(X)", zsto,zout)
97c
98         CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
99     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
100     .                "ave(X)", zsto1,zout)
101c
102         CALL histdef(nid_day, "dudyn", "Dynamics dU", "m/s2",
103     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
104     .                "ave(X)", zsto,zout)
105c
106         CALL histdef(nid_day, "duvdf", "Boundary-layer dU", "m/s2",
107     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
108     .                "ave(X)", zsto,zout)
109c
110         if (iflag_trac.eq.1) THEN
111          if (microfi.eq.1) then
112           DO iq=1,nmicro
113         CALL histdef(nid_day, tname(iq), ttext(iq), "n/m2",
114     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
115     .                "ave(X)", zsto,zout)
116           ENDDO
117          endif
118          if (nmicro.lt.nqmax) then
119           DO iq=nmicro+1,nqmax
120         CALL histdef(nid_day, tname(iq), ttext(iq), "ppm",
121     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
122     .                "ave(X)", zsto,zout)
123           ENDDO
124          endif
125         endif
126c
127      ENDIF !lev_histday.GE.2
128c
129c-------------------------------------------------------
130      IF(lev_histday.GE.3) THEN
131c
132cccccccccccccccccc  Radiative transfer
133c
134c 2D
135c
136         CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
137     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
138     .                "ave(X)", zsto1,zout)
139c
140         CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2",
141     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
142     .                "ave(X)", zsto1,zout)
143c
144         CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2",
145     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
146     .                "ave(X)", zsto1,zout)
147c
148c 3D
149c
150         CALL histdef(nid_day, "SWnet", "Net SW flux","W/m2",
151     .                iim,jjmp1,nhori, klev,1,klev,nvert,
152     .                32, "ave(X)", zsto1,zout)
153c
154         CALL histdef(nid_day, "LWnet", "Net LW flux","W/m2",
155     .                iim,jjmp1,nhori, klev,1,klev,nvert,
156     .                32, "ave(X)", zsto1,zout)
157c
158      ENDIF !lev_histday.GE.3
159c
160c-------------------------------------------------------
161      IF(lev_histday.GE.4) THEN
162c
163         CALL histdef(nid_day, "dtdyn", "Dynamics dT", "K/s",
164     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
165     .                "ave(X)", zsto,zout)
166c
167         CALL histdef(nid_day, "dtphy", "Physics dT", "K/s",
168     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
169     .                "ave(X)", zsto,zout)
170c
171         CALL histdef(nid_day, "dtvdf", "Boundary-layer dT", "K/s",
172     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
173     .                "ave(X)", zsto,zout)
174c
175         CALL histdef(nid_day, "dtajs", "Dry adjust. dT", "K/s",
176     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
177     .                "ave(X)", zsto,zout)
178c
179         CALL histdef(nid_day, "dtswr", "SW radiation dT", "K/s",
180     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
181     .                "ave(X)", zsto,zout)
182c
183         CALL histdef(nid_day, "dtlwr", "LW radiation dT", "K/s",
184     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
185     .                "ave(X)", zsto,zout)
186c
187         CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s",
188     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
189     .                "ave(X)", zsto,zout)
190c
191      ENDIF !lev_histday.GE.4
192c
193c-------------------------------------------------------
194      IF(lev_histday.GE.5) THEN
195c
196c
197c        call histdef(nid_day, "taux", 
198c    $         "Zonal wind stress", "Pa", 
199c    $         iim,jjmp1,nhori, 1,1,1, -99, 32,
200c    $         "ave(X)", zsto,zout)
201c
202c        call histdef(nid_day, "tauy", 
203c    $         "Meridional xind stress", "Pa", 
204c    $         iim,jjmp1,nhori, 1,1,1, -99, 32,
205c    $         "ave(X)", zsto,zout)
206c
207c        CALL histdef(nid_day, "cdrm", "Momentum drag coef.", "-",
208c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
209c    .                "ave(X)", zsto,zout)
210c
211c        CALL histdef(nid_day, "cdrh", "Heat drag coef.", "-",
212c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
213c    .                "ave(X)", zsto,zout)
214c
215      ENDIF !lev_histday.GE.5
216c-------------------------------------------------------
217c
218         CALL histend(nid_day)
219c
220         ndex2d = 0
221         ndex3d = 0
222c
223      ENDIF ! fin de test sur ok_journe
Note: See TracBrowser for help on using the repository browser.