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