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

Last change on this file since 102 was 97, checked in by slebonnois, 14 years ago

Serie de modifs SL pour homogeneisation des phytitan et phyvenus
Ca touche aussi aux liens phy/dyn (surtout a propos de clesphy0),
a verifier avec les autres, donc...
File size: 7.9 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
98         CALL histdef(nid_ins, "duvdf", "Boundary-layer dU", "m/s2",
99     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
100     .                "inst(X)", zsto,zout)
101c
102         CALL histdef(nid_ins, "dudyn", "Dynamics dU", "m/s2",
103     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
104     .                "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     .                "ave(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     .                "ave(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      ENDIF !lev_histday.GE.3
157c
158c-------------------------------------------------------
159      IF(lev_histday.GE.4) THEN
160c
161         CALL histdef(nid_ins, "dtdyn", "Dynamics dT", "K/s",
162     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
163     .                "inst(X)", zsto,zout)
164c
165         CALL histdef(nid_ins, "dtphy", "Physics dT", "K/s",
166     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
167     .                "inst(X)", zsto,zout)
168c
169         CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s",
170     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
171     .                "inst(X)", zsto,zout)
172c
173         CALL histdef(nid_ins, "dtajs", "Dry adjust. dT", "K/s",
174     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
175     .                "inst(X)", zsto,zout)
176c
177         CALL histdef(nid_ins, "dtswr", "SW radiation dT", "K/s",
178     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
179     .                "inst(X)", zsto,zout)
180c
181         CALL histdef(nid_ins, "dtlwr", "LW radiation dT", "K/s",
182     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
183     .                "inst(X)", zsto,zout)
184c
185         CALL histdef(nid_ins, "dtec", "Cinetic dissip dT", "K/s",
186     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
187     .                "inst(X)", zsto,zout)
188c
189c        CALL histdef(nid_ins, "dvvdf", "Boundary-layer dV", "m/s2",
190c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
191c    .                "inst(X)", zsto,zout)
192c
193      ENDIF !lev_histday.GE.4
194c
195c-------------------------------------------------------
196      IF(lev_histday.GE.5) THEN
197c
198c
199c        call histdef(nid_ins, "taux", 
200c    $         "Zonal wind stress", "Pa", 
201c    $         iim,jjmp1,nhori, 1,1,1, -99, 32,
202c    $         "inst(X)", zsto,zout)
203c
204c        call histdef(nid_ins, "tauy", 
205c    $         "Meridional xind stress", "Pa", 
206c    $         iim,jjmp1,nhori, 1,1,1, -99, 32,
207c    $         "inst(X)", zsto,zout)
208c
209c        CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-",
210c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
211c    .                "inst(X)", zsto,zout)
212c
213c        CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-",
214c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
215c    .                "inst(X)", zsto,zout)
216c
217      ENDIF !lev_histday.GE.5
218c-------------------------------------------------------
219
220         CALL histend(nid_ins)
221c
222         ndex2d = 0
223         ndex3d = 0
224c
225      ENDIF
Note: See TracBrowser for help on using the repository browser.