source: LMDZ5/branches/testing/libf/phylmd/ini_histrac.h @ 5423

Last change on this file since 5423 was 2298, checked in by Laurent Fairhead, 10 years ago

Merged trunk changes -r2237:2291 into testing branch
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 11.2 KB
Line 
1!
2! $Id $
3!
4  IF (ecrit_tra>0.) THEN
5!$OMP MASTER
6!!!     CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
7! correction pour l heure initiale                               !jyg
8!                                                               !jyg
9      CALL ymds2ju(annee_ref, 1, day_ref, hour, zjulian)         !jyg
10
11     CALL histbeg_phy("histrac", itau_phy, zjulian, pdtphys,nhori, nid_tra)
12     CALL histvert(nid_tra, "presnivs", "Vertical levels", "Pa",klev, presnivs, nvert,"down")
13
14     zsto = pdtphys
15     zout = ecrit_tra
16     CALL histdef(nid_tra, "phis", "Surface geop. height", "-",   &
17          iim,jj_nb,nhori, 1,1,1, -99, 32,"once",  zsto,zout)
18     CALL histdef(nid_tra, "aire", "Grid area", "-",              &
19          iim,jj_nb,nhori, 1,1,1, -99, 32,"once",  zsto,zout)
20     CALL histdef(nid_tra, "zmasse", "column density of air in cell", &
21          "kg m-2", iim, jj_nb, nhori, klev, 1, klev, nvert, 32, "ave(X)", &
22          zsto,zout)
23! RomP >>>
24     CALL histdef(nid_tra, "sourceBE", "source 7Be", &
25          "at/kgA/s", iim, jj_nb, nhori, klev, 1, klev, nvert, 32, "ave(X)", &
26          zsto,zout)
27! RomP <<<
28
29!TRACEURS
30!----------------
31     DO it = 1,nbtr
32!!        iiq = niadv(it+2)                                                         ! jyg
33        iiq = niadv(it+nqo)                                                         ! jyg
34
35! CONCENTRATIONS
36        CALL histdef(nid_tra, tname(iiq), ttext(iiq), "U/kga",    &
37             iim,jj_nb,nhori, klev,1,klev,nvert, 32,"ave(X)", zsto,zout)
38
39! TD LESSIVAGE
40        IF (lessivage .AND. aerosol(it)) THEN
41           CALL histdef(nid_tra, "fl"//tname(iiq),"Flux "//ttext(iiq), &
42                "at/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,       &
43                "ave(X)", zsto,zout)
44           CALL histdef(nid_tra, "d_tr_ls_"//tname(iiq),      &
45                "tendance lessivage large scale"// ttext(iiq), "?",&
46                iim,jj_nb,nhori, klev,1,klev,nvert, 32,       &
47                "ave(X)", zsto,zout)
48           CALL histdef(nid_tra, "d_tr_insc_"//tname(iiq),      &
49                "tendance lessivage large scale"// ttext(iiq), "?",&
50                iim,jj_nb,nhori, klev,1,klev,nvert, 32,       &
51                "ave(X)", zsto,zout)
52           CALL histdef(nid_tra, "d_tr_bcscav_"//tname(iiq),      &
53                "tendance lessivage large scale"// ttext(iiq), "?",&
54                iim,jj_nb,nhori, klev,1,klev,nvert, 32,       &
55                "ave(X)", zsto,zout)
56           CALL histdef(nid_tra, "d_tr_evls_"//tname(iiq),      &
57                "tendance lessivage large scale"// ttext(iiq), "?",&
58                iim,jj_nb,nhori, klev,1,klev,nvert, 32,       &
59                "ave(X)", zsto,zout)
60!  Tracer concentration in LS precipitation at surface
61           CALL histdef(nid_tra, "qpr_ls_"//tname(iiq),       &
62                "concentration in LS precip"// ttext(iiq), "at/kgw", &
63                iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
64                "ave(X)", zsto,zout)     
65                 END IF
66
67! TD THERMIQUES
68        IF (iflag_thermals.gt.0) THEN
69           CALL histdef(nid_tra, "d_tr_th_"//tname(iiq),      &
70                "tendance thermique"// ttext(iiq), "?",       &
71                iim,jj_nb,nhori, klev,1,klev,nvert, 32,       &
72                "ave(X)", zsto,zout)
73        ENDIF
74
75! TD CONVECTION
76        IF (iflag_con.GE.2) THEN
77           CALL histdef(nid_tra, "d_tr_cv_"//tname(iiq),   &
78                "tendance convection"// ttext(iiq), "?",   &
79                iim,jj_nb,nhori, klev,1,klev,nvert, 32,    &
80                "ave(X)", zsto,zout)
81        ENDIF
82
83! RomP >>>
84        IF (iflag_con.EQ.30) THEN
85           CALL histdef(nid_tra, "d_tr_cvMA_"//tname(iiq),   &
86                "tendance convection"// ttext(iiq), "?",&
87                iim,jj_nb,nhori, klev,1,klev,nvert, 32,    &
88                "ave(X)", zsto,zout)
89           CALL histdef(nid_tra, "d_tr_trsp_"//tname(iiq),   &
90                "tendance transport "// ttext(iiq), "at/kga",   &
91                iim,jj_nb,nhori, klev,1,klev,nvert, 32,    &
92                "ave(X)", zsto,zout)
93           CALL histdef(nid_tra, "d_tr_sscav_"//tname(iiq),   &
94                "tendance lessivage flux satures "// ttext(iiq), "at/kga",   &
95                iim,jj_nb,nhori, klev,1,klev,nvert, 32,    &
96                "ave(X)", zsto,zout)
97           CALL histdef(nid_tra, "d_tr_sat_"//tname(iiq),   &
98                "tendance flux satures "// ttext(iiq), "at/kga",  &
99                iim,jj_nb,nhori, klev,1,klev,nvert, 32,    &
100                "ave(X)", zsto,zout)
101           CALL histdef(nid_tra, "d_tr_uscav_"//tname(iiq),  &
102                "tendance flux insatures "// ttext(iiq), "at/kga", &
103                iim,jj_nb,nhori, klev,1,klev,nvert, 32,    &
104                "ave(X)", zsto,zout)
105           CALL histdef(nid_tra, "tr_pr_"//tname(iiq),  &
106                "concentration dans precip"// ttext(iiq), "at/kga", &
107                iim,jj_nb,nhori, klev,1,klev,nvert, 32,    &
108                "ave(X)", zsto,zout)
109           CALL histdef(nid_tra, "tr_aa_"//tname(iiq),  &
110                "concentration precip issu AA"// ttext(iiq), "at/kga", &
111                iim,jj_nb,nhori, klev,1,klev,nvert, 32,    &
112                "ave(X)", zsto,zout)
113           CALL histdef(nid_tra, "tr_mel_"//tname(iiq),  &
114                "concentration precip issu melange"// ttext(iiq), "at/kga", &
115                iim,jj_nb,nhori, klev,1,klev,nvert, 32,    &
116                "ave(X)", zsto,zout)
117           CALL histdef(nid_tra, "tr_di_"//tname(iiq),  &
118                "concentration dans descente insaturee"// ttext(iiq), "at/kga", &
119                iim,jj_nb,nhori, klev,1,klev,nvert, 32,    &
120                "ave(X)", zsto,zout)
121           CALL histdef(nid_tra, "tr_trspdi_"//tname(iiq),  &
122                "conc descente insaturee MA"// ttext(iiq), "at/kga", &
123                iim,jj_nb,nhori, klev,1,klev,nvert, 32,    &
124                "ave(X)", zsto,zout)
125           CALL histdef(nid_tra, "zmfd1a_"//tname(iiq),  &
126                "zmfd1a"// ttext(iiq), "_", &
127                iim,jj_nb,nhori, klev,1,klev,nvert, 32,    &
128                "ave(X)", zsto,zout)
129           CALL histdef(nid_tra, "zmfphi2_"//tname(iiq),  &
130                "zmfphi2"// ttext(iiq), "_", &
131                iim,jj_nb,nhori, klev,1,klev,nvert, 32,    &
132                "ave(X)", zsto,zout)
133           CALL histdef(nid_tra, "zmfdam_"//tname(iiq),  &
134                "zmfdam"// ttext(iiq), "_", &
135                iim,jj_nb,nhori, klev,1,klev,nvert, 32,    &
136                "ave(X)", zsto,zout)
137          ENDIF
138! RomP <<<
139           CALL histdef(nid_tra, "dtrdyn_"//tname(iiq),  &
140                "td dyn tra"// ttext(iiq), "at/kga", &
141                iim,jj_nb,nhori, klev,1,klev,nvert, 32,  &
142                "ave(X)", zsto,zout)
143! TD decroissance radioactive
144           CALL histdef(nid_tra, "d_tr_dec_"//tname(iiq),   &
145                "tendance decroi radio "// ttext(iiq), "",  &
146                iim,jj_nb,nhori, klev,1,klev,nvert, 32,  &
147                "ave(X)", zsto,zout)
148
149! TD COUCHE-LIMITE
150      IF (iflag_vdf_trac>=0) THEN
151        CALL histdef(nid_tra, "d_tr_cl_"//tname(iiq),      &
152             "tendance couche limite"// ttext(iiq), "?",   &
153             iim,jj_nb,nhori, klev,1,klev,nvert, 32,       &
154             "ave(X)", zsto,zout)
155!  Dry deposit (1st layer and surface)
156        CALL histdef(nid_tra, "d_tr_dry_"//tname(iiq),       &
157             "tendancy dry deposit"// ttext(iiq), "at/kga/step", &
158             iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
159             "ave(X)", zsto,zout)     
160        CALL histdef(nid_tra, "flux_tr_dry_"//tname(iiq),       &
161             "dry deposit at surf (downward)"// ttext(iiq), "at/m2/step", &
162             iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
163             "ave(X)", zsto,zout)     
164      ENDIF
165     ENDDO
166
167     CALL histdef(nid_tra, "Mint", "Mint","",         &
168          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
169          "inst(X)", zout,zout)
170     CALL histdef(nid_tra, "frac_impa", "frac_impa","",         &
171          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
172          "inst(X)", zout,zout)
173     CALL histdef(nid_tra, "frac_nucl", "frac_nucl","",         &
174          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
175          "inst(X)", zout,zout)
176!---------------   
177!
178! VENT (niveau 1)
179     CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-",      &
180          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
181          "inst(X)",  zout,zout)     
182     CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-",      &
183          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
184          "inst(X)",  zout,zout)
185
186! TEMPERATURE DU SOL
187     CALL histdef(nid_tra, "ftsol1", "temper sol", "-",    &
188          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
189          "inst(X)",  zout,zout)
190     CALL histdef(nid_tra, "ftsol2", "temper sol", "-",    &
191          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
192          "inst(X)",  zout,zout)
193     CALL histdef(nid_tra, "ftsol3", "temper sol", "-",    &
194          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
195          "inst",  zout,zout)
196     CALL histdef(nid_tra, "ftsol4", "temper sol", "-",    &
197          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
198          "inst(X)",  zout,zout)
199
200! NATURE DU SOL
201     CALL histdef(nid_tra, "psrf1", "nature sol", "-",     &
202          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
203          "inst(X)",  zout,zout)
204     CALL histdef(nid_tra, "psrf2", "nature sol", "-",     &
205          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
206          "inst(X)",  zout,zout)
207     CALL histdef(nid_tra, "psrf3", "nature sol", "-",     &
208          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
209          "inst(X)",  zout,zout)
210     CALL histdef(nid_tra, "psrf4", "nature sol", "-",     &
211          iim,jj_nb,nhori, 1,1,1, -99, 32,                 & 
212          "inst(X)",  zout,zout)
213! DIVERS
214     CALL histdef(nid_tra, "pplay", "pressure","-",        &
215          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
216          "inst(X)", zout,zout)
217     CALL histdef(nid_tra, "T", "temperature","K",         &
218          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
219          "inst(X)", zout,zout)
220     CALL histdef(nid_tra, "mfu", "flux u mont","-",       &
221          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
222          "ave(X)", zsto,zout)
223     CALL histdef(nid_tra, "mfd", "flux u decen","-",      &
224          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
225          "ave(X)", zsto,zout)
226     CALL histdef(nid_tra, "en_u", "flux u mont","-",      &
227          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
228          "ave(X)", zsto,zout)
229     CALL histdef(nid_tra, "en_d", "flux u mont","-",      &
230          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
231          "ave(X)", zsto,zout)
232     CALL histdef(nid_tra, "de_d", "flux u mont","-",      &
233          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
234          "ave(X)", zsto,zout)
235     CALL histdef(nid_tra, "de_u", "flux u decen","-",     &
236          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
237          "ave(X)", zsto,zout)
238     CALL histdef(nid_tra, "coefh", "turbulent coef","-",  &
239          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
240          "ave(X)", zsto,zout)   
241     
242     CALL histend(nid_tra)
243!$OMP END MASTER
244  END IF ! ecrit_tra>0.
245 
Note: See TracBrowser for help on using the repository browser.