source: trunk/LMDZ.VENUS/libf/phyvenus/ini_histins.h @ 1442

Last change on this file since 1442 was 1442, checked in by slebonnois, 10 years ago

SL: update of the Venus GCM, + corrections on routines used for newstart/start2archive for Titan and Venus, + some modifications on tools
File size: 12.8 KB
RevLine 
[3]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
[892]5
[902]6         zsto = dtime * ecrit_ins
7         zout = dtime * ecrit_ins
[892]8
[3]9         idayref = day_ref
[97]10         CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
[892]11
12         call histbeg_phy("histins.nc",itau_phy,
[902]13     .                    zjulian,dtime,nhori,nid_ins)
[892]14
15!$OMP MASTER
16         CALL histvert(nid_ins, "presnivs", "Vertical levels", "Pa",
[97]17     .                 klev, presnivs, nvert)
[3]18
19c-------------------------------------------------------
[1305]20      IF(lev_histins.GE.1) THEN
[888]21c
22ccccccccccccc 2D fields, basics
23c
24         CALL histdef(nid_ins, "phis", "Surface geop. height", "-",
[892]25     .                iim,jj_nb,nhori, 1,1,1, nvert, 32, 
[888]26     .                "once",  zsto,zout)
27c
28         CALL histdef(nid_ins, "aire", "Grid area", "-",
[892]29     .                iim,jj_nb,nhori, 1,1,1, nvert, 32, 
[888]30     .                "once",  zsto,zout)
31c
32         CALL histdef(nid_ins, "tsol", "Surface Temperature", "K",
[892]33     .                iim,jj_nb,nhori, 1,1,1, nvert, 32, 
[888]34     .                "ins(X)", zsto,zout)
35c
36         CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa",
[892]37     .                iim,jj_nb,nhori, 1,1,1, nvert, 32, 
[888]38     .                "ins(X)", zsto,zout)
39c
40c        CALL histdef(nid_ins, "ue", "Zonal energy transport", "-",
[892]41c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32, 
[888]42c    .                "ins(X)", zsto,zout)
43c
44c        CALL histdef(nid_ins, "ve", "Merid energy transport", "-",
[892]45c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32, 
[888]46c    .                "ins(X)", zsto,zout)
47c
[902]48c        CALL histdef(nid_ins, "cdragh", "Drag coef on T", "-",
49c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32, 
50c    .                "ins(X)", zsto,zout)
[888]51c
[902]52c        CALL histdef(nid_ins, "cdragm", "Drag coef on U", "-",
53c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32, 
54c    .                "ins(X)", zsto,zout)
[888]55c
[1305]56      ENDIF !lev_histins.GE.1
[888]57c
[3]58c-------------------------------------------------------
[1305]59      IF(lev_histins.GE.2) THEN
[888]60c
61ccccccccccccc 3D fields, basics
62c
63         CALL histdef(nid_ins, "temp", "Air temperature", "K",
[892]64     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]65     .                "ins(X)", zsto,zout)
66c
67         CALL histdef(nid_ins, "pres", "Air pressure", "Pa",
[892]68     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]69     .                "ins(X)", zsto,zout)
70c
71         CALL histdef(nid_ins, "geop", "Geopotential height", "m",
[892]72     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]73     .                "ins(X)", zsto,zout)
74c
75         CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s",
[892]76     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]77     .                "ins(X)", zsto,zout)
78c
79         CALL histdef(nid_ins, "vitv", "Meridional wind", "m/s",
[892]80     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]81     .                "ins(X)", zsto,zout)
82c
83         CALL histdef(nid_ins, "vitw", "Vertical wind", "Pa/s",
[892]84     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]85     .                "ins(X)", zsto,zout)
86c
87         CALL histdef(nid_ins, "dudyn", "Dynamics dU", "m/s2",
[892]88     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]89     .                "ins(X)", zsto,zout)
90c
91         CALL histdef(nid_ins, "duvdf", "Boundary-layer dU", "m/s2",
[892]92     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]93     .                "ins(X)", zsto,zout)
94c
95c        CALL histdef(nid_ins, "mang", "Angular momentum", "kg m2/s",
[892]96c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]97c    .                "ins(X)", zsto,zout)
98c
[902]99c        CALL histdef(nid_ins, "Kz", "vertical diffusion coef", "m2/s",
100c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
101c    .                "ins(X)", zsto,zout)
[888]102c
103c plusieurs traceurs
[1442]104          if (ok_chem) THEN
[888]105            DO iq=1,nqmax
[1442]106c            DO iq=1,10
[888]107             IF (iq.LE.99) THEN
108          WRITE(str2,'(i2.2)') iq
[1442]109          CALL histdef(nid_ins, tname(iq), ttext(iq), "mol/mol",
[892]110     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]111     .                "ins(X)", zsto,zout)
112             ELSE
113          PRINT*, "Trop de traceurs"
114          CALL abort
115             ENDIF
116            ENDDO
[1442]117          CALL histdef(nid_mth, "d_qmoldif CO2", "Dif molec" , "kg/kg",
118     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
119     .                "ave(X)", zsto,zout)
120          CALL histdef(nid_mth, "d_qmoldif O3p", "Dif molec" , "kg/kg",
121     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
122     .                "ave(X)", zsto,zout)
123          CALL histdef(nid_mth, "d_qmoldif N2", "Dif molec" , "kg/kg",
124     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
125     .                "ave(X)", zsto,zout)
126
127          endif
[888]128c
129         CALL histdef(nid_ins, "tops", "Solar rad. at TOA", "W/m2",
[892]130     .                iim,jj_nb,nhori, 1,1,1, nvert, 32, 
[888]131     .                "ins(X)", zsto,zout)
132c
[1442]133          if (ok_cloud) THEN
134
135          if (nb_mode.GE.1) THEN
136           
137c
138         CALL histdef(nid_ins, "NBRTOTm1", "Nbr total droplet",
139     .                "#/cm3", iim,jj_nb,nhori, klev,1,klev, nvert, 32, 
[1305]140     .                "ins(X)", zsto,zout)
[1442]141c
142
143c
144c         CALL histdef(nid_ins, "R_MEDIANm1", "Median radius
145c     .    for log normal distribution" ,
146c     .                "fraction",
147c     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32, 
148c     .                "ins(X)", zsto,zout)
149c
150
151c
152c         CALL histdef(nid_ins, "STDDEVm1", "Std Deviation
153c     .    for log normal distribution",
154c     .                "fraction",
155c     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32, 
156c     .                "ins(X)", zsto,zout)
157c
158
159          if (nb_mode.GE.2) THEN
160
161c
162         CALL histdef(nid_ins, "NBRTOTm2", "Nbr total droplet",
163     .                "#/cm3", iim,jj_nb,nhori, klev,1,klev, nvert, 32, 
[1305]164     .                "ins(X)", zsto,zout)
[1442]165c
166
167c
168c         CALL histdef(nid_ins, "R_MEDIANm2", "Median radius
169c     .    for log normal distribution" ,
170c     .                "fraction",
171c     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32, 
172c     .                "ins(X)", zsto,zout)
173c
174
175c
176c         CALL histdef(nid_ins, "STDDEVm2", "Std Deviation
177c     .    for log normal distribution",
178c     .                "fraction",
179c     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32, 
180c     .                "ins(X)", zsto,zout)
181c
182
183          if (nb_mode.GE.3) THEN
184         
185c
186         CALL histdef(nid_ins, "NBRTOTm3", "Nbr total droplet", "#/cm3",
187     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32, 
[1305]188     .                "ins(X)", zsto,zout)
[1442]189c
190
191c
192c         CALL histdef(nid_ins, "R_MEDIANm3", "Median radius
193c     .   for log normal distribution" ,
194c     .                "fraction",
195c     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32, 
196c     .                "ins(X)", zsto,zout)
197c
198
199c
200c         CALL histdef(nid_ins, "STDDEVm3", "Std Deviation
201c     .    for log normal distribution",
202c     .                "fraction",
203c     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32, 
204c     .                "ins(X)", zsto,zout)
205c
206
207         ENDIF
208         ENDIF
209         ENDIF
210         
211c
212         CALL histdef(nid_ins, "WH2SO4", "Weight fraction H2SO4",
213     .                "fraction",
214     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32, 
[1305]215     .                "ins(X)", zsto,zout)
[1442]216c
217
218c
219         CALL histdef(nid_ins, "rho_droplet", "density cloud droplet",
220     .                "kg.m-3",
221     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32, 
[1305]222     .                "ins(X)", zsto,zout)
[1442]223c
[1305]224
[1442]225                ENDIF
226               
227          if (ok_sedim) THEN
228c
229         CALL histdef(nid_ins, "d_tr_sed_H2SO4", "var mmr from sedim",
230     .                "kg/kg",
231     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32, 
[1305]232     .                "ins(X)", zsto,zout)
[1442]233c
234
235c
236         CALL histdef(nid_ins, "d_tr_sed_H2O", "var mmr from sedim",
237     .                "kg/kg",
238     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32, 
[1305]239     .                "ins(X)", zsto,zout)
[1442]240c
241
242c
243         CALL histdef(nid_ins, "F_sedim", "tendency from sedim",
244     .                "kg.m-2.s-1",
245     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32, 
[1305]246     .                "ins(X)", zsto,zout)
[1442]247c
248                ENDIF
[1305]249      ENDIF !lev_histins.GE.2
[888]250c
251c-------------------------------------------------------
[1305]252      IF(lev_histins.GE.3) THEN
[888]253c
254cccccccccccccccccc  Radiative transfer
255c
256c 2D
257c
258         CALL histdef(nid_ins, "topl", "IR rad. at TOA", "W/m2",
[892]259     .                iim,jj_nb,nhori, 1,1,1, nvert, 32, 
[888]260     .                "ins(X)", zsto,zout)
261c
262         CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2",
[892]263     .                iim,jj_nb,nhori, 1,1,1, nvert, 32, 
[888]264     .                "ins(X)", zsto,zout)
265c
266         CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2",
[892]267     .                iim,jj_nb,nhori, 1,1,1, nvert, 32, 
[888]268     .                "ins(X)", zsto,zout)
269c
270c 3D
271c
272         CALL histdef(nid_ins, "SWnet", "Net SW flux","W/m2",
[892]273     .                iim,jj_nb,nhori, klev,1,klev,nvert,
[888]274     .                32, "ins(X)", zsto,zout)
275c
276         CALL histdef(nid_ins, "LWnet", "Net LW flux","W/m2",
[892]277     .                iim,jj_nb,nhori, klev,1,klev,nvert,
[888]278     .                32, "ins(X)", zsto,zout)
279c
280         CALL histdef(nid_ins, "fluxvdf", "PBL net flux","W/m2",
[892]281     .                iim,jj_nb,nhori, klev,1,klev,nvert,
[888]282     .                32, "ins(X)", zsto,zout)
283c
284         CALL histdef(nid_ins, "fluxdyn", "Dyn. net flux","W/m2",
[892]285     .                iim,jj_nb,nhori, klev,1,klev,nvert,
[888]286     .                32, "ins(X)", zsto,zout)
287c
288         CALL histdef(nid_ins, "fluxajs", "Dry adj. net flux","W/m2",
[892]289     .                iim,jj_nb,nhori, klev,1,klev,nvert,
[888]290     .                32, "ins(X)", zsto,zout)
291c
292c        CALL histdef(nid_ins, "fluxec", "Cin. net flux","W/m2",
[892]293c    .                iim,jj_nb,nhori, klev,1,klev,nvert,
[888]294c    .                32, "ins(X)", zsto,zout)
295c
[1305]296      ENDIF !lev_histins.GE.3
[888]297c
298c-------------------------------------------------------
[1305]299      IF(lev_histins.GE.4) THEN
[888]300c
301         CALL histdef(nid_ins, "dtdyn", "Dynamics dT", "K/s",
[892]302     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]303     .                "ins(X)", zsto,zout)
304c
305c        CALL histdef(nid_ins, "dtphy", "Physics dT", "K/s",
[892]306c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]307c    .                "ins(X)", zsto,zout)
308c
309         CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s",
[892]310     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]311     .                "ins(X)", zsto,zout)
312c
313         CALL histdef(nid_ins, "dtajs", "Dry adjust. dT", "K/s",
[892]314     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]315     .                "ins(X)", zsto,zout)
316c
317         CALL histdef(nid_ins, "dtswr", "SW radiation dT", "K/s",
[892]318     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]319     .                "ins(X)", zsto,zout)
320c
321         CALL histdef(nid_ins, "dtlwr", "LW radiation dT", "K/s",
[892]322     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]323     .                "ins(X)", zsto,zout)
324c
325c        CALL histdef(nid_ins, "dtec", "Cinetic dissip dT", "K/s",
[892]326c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]327c    .                "ins(X)", zsto,zout)
328c
329         CALL histdef(nid_ins, "duajs", "Dry convection dU", "m/s2",
[892]330     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]331     .                "ins(X)", zsto,zout)
332c
333         CALL histdef(nid_ins, "dugwo", "GW oro dU", "m/s2",
[892]334     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]335     .                "ins(X)", zsto,zout)
336c
337         CALL histdef(nid_ins, "dugwno", "GW non-oro dU", "m/s2",
[892]338     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]339     .                "ins(X)", zsto,zout)
340c
341c        CALL histdef(nid_ins, "dvvdf", "Boundary-layer dV", "m/s2",
[892]342c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]343c    .                "ins(X)", zsto,zout)
344c
[1305]345      ENDIF !lev_histins.GE.4
[888]346c
347c-------------------------------------------------------
[1305]348      IF(lev_histins.GE.5) THEN
[888]349c
350c        call histdef(nid_ins, "taux", 
351c    $         "Zonal wind stress", "Pa", 
[892]352c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
[888]353c    $         "ins(X)", zsto,zout)
354c
355c        call histdef(nid_ins, "tauy", 
356c    $         "Meridional xind stress", "Pa", 
[892]357c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
[888]358c    $         "ins(X)", zsto,zout)
359c
360c        CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-",
[892]361c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32, 
[888]362c    .                "ins(X)", zsto,zout)
363c
364c        CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-",
[892]365c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32, 
[888]366c    .                "ins(X)", zsto,zout)
367c
[1305]368      ENDIF !lev_histins.GE.5
[888]369c-------------------------------------------------------
370c
[3]371         CALL histend(nid_ins)
[892]372!$OMP END MASTER
373
[3]374      ENDIF
Note: See TracBrowser for help on using the repository browser.