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

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