source: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/write_bilKP_ins.h @ 5023

Last change on this file since 5023 was 644, checked in by Laurent Fairhead, 19 years ago

Synchronisation avec tous les diagnostiques de Ionela IM
Inclusion du slab ocean IM
LF
  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 8.4 KB
Line 
1c
2c $Header$
3c
4      IF (ok_journe) THEN
5c
6      ndex2d = 0
7      ndex3d = 0
8c
9      itau_w = itau_phy + itap
10c
11c Champs 3D:
12c
13      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, ue_lay,zx_tmp_3d)
14      CALL histwrite(nid_bilKPins,"ue",itau_w,zx_tmp_3d,
15     .               iim*jjmp1*klev,ndex3d)
16c
17      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, ve_lay,zx_tmp_3d)
18      CALL histwrite(nid_bilKPins,"ve",itau_w,zx_tmp_3d,
19     .               iim*jjmp1*klev,ndex3d)
20c
21      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, uq_lay,zx_tmp_3d)
22      CALL histwrite(nid_bilKPins,"uq",itau_w,zx_tmp_3d,
23     .               iim*jjmp1*klev,ndex3d)
24c
25      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, vq_lay,zx_tmp_3d)
26      CALL histwrite(nid_bilKPins,"vq",itau_w,zx_tmp_3d,
27     .               iim*jjmp1*klev,ndex3d)
28c
29c Champs 3D:
30C
31      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
32      CALL histwrite(nid_bilKPins,"temp",itau_w,zx_tmp_3d,
33     .                                   iim*jjmp1*klev,ndex3d)
34c
35      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
36      CALL histwrite(nid_bilKPins,"ovap",itau_w,zx_tmp_3d,
37     .                                   iim*jjmp1*klev,ndex3d)
38c
39      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
40      CALL histwrite(nid_bilKPins,"geop",itau_w,zx_tmp_3d,
41     .                                   iim*jjmp1*klev,ndex3d)
42c
43      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
44      CALL histwrite(nid_bilKPins,"vitu",itau_w,zx_tmp_3d,
45     .                                   iim*jjmp1*klev,ndex3d)
46c
47      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
48      CALL histwrite(nid_bilKPins,"vitv",itau_w,zx_tmp_3d,
49     .                                   iim*jjmp1*klev,ndex3d)
50c
51      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
52      CALL histwrite(nid_bilKPins,"vitw",itau_w,zx_tmp_3d,
53     .                                   iim*jjmp1*klev,ndex3d)
54c
55      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
56      CALL histwrite(nid_bilKPins,"pres",itau_w,zx_tmp_3d,
57     .                                   iim*jjmp1*klev,ndex3d)
58c
59      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, paprs, zx_tmp_3d)
60      CALL histwrite(nid_bilKPins,"play",itau_w,zx_tmp_3d,
61     .                                   iim*jjmp1*klev,ndex3d)
62c
63      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldliq, zx_tmp_3d)
64      CALL histwrite(nid_bilKPins,"oliq",itau_w,zx_tmp_3d,
65     .                                   iim*jjmp1*klev,ndex3d)
66c
67      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
68      CALL histwrite(nid_bilKPins,"dtdyn",itau_w,zx_tmp_3d,
69     .                                   iim*jjmp1*klev,ndex3d)
70c
71      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
72      CALL histwrite(nid_bilKPins,"dqdyn",itau_w,zx_tmp_3d,
73     .                                   iim*jjmp1*klev,ndex3d)
74c
75      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_con, zx_tmp_3d)
76      CALL histwrite(nid_bilKPins,"dtcon",itau_w,zx_tmp_3d,
77     .                                   iim*jjmp1*klev,ndex3d)
78c
79      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_con, zx_tmp_3d)
80      CALL histwrite(nid_bilKPins,"ducon",itau_w,zx_tmp_3d,
81     .                                   iim*jjmp1*klev,ndex3d)
82c
83      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_con, zx_tmp_3d)
84      CALL histwrite(nid_bilKPins,"dvcon",itau_w,zx_tmp_3d,
85     .                                   iim*jjmp1*klev,ndex3d)
86c
87      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_con, zx_tmp_3d)
88      CALL histwrite(nid_bilKPins,"dqcon",itau_w,zx_tmp_3d,
89     .                                   iim*jjmp1*klev,ndex3d)
90c
91      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_lsc, zx_tmp_3d)
92      CALL histwrite(nid_bilKPins,"dtlsc",itau_w,zx_tmp_3d,
93     .                                   iim*jjmp1*klev,ndex3d)
94c
95      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_lsc, zx_tmp_3d)
96      CALL histwrite(nid_bilKPins,"dqlsc",itau_w,zx_tmp_3d,
97     .                                   iim*jjmp1*klev,ndex3d)
98c
99      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_vdf, zx_tmp_3d)
100      CALL histwrite(nid_bilKPins,"dtvdf",itau_w,zx_tmp_3d,
101     .                                   iim*jjmp1*klev,ndex3d)
102c
103      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d)
104      CALL histwrite(nid_bilKPins,"dqvdf",itau_w,zx_tmp_3d,
105     .                                   iim*jjmp1*klev,ndex3d)
106c
107      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_ajs, zx_tmp_3d)
108      CALL histwrite(nid_bilKPins,"dtajs",itau_w,zx_tmp_3d,
109     .                                   iim*jjmp1*klev,ndex3d)
110c
111      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_ajs, zx_tmp_3d)
112      CALL histwrite(nid_bilKPins,"dqajs",itau_w,zx_tmp_3d,
113     .                                   iim*jjmp1*klev,ndex3d)
114c
115      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_eva, zx_tmp_3d)
116      CALL histwrite(nid_bilKPins,"dteva",itau_w,zx_tmp_3d,
117     .                                   iim*jjmp1*klev,ndex3d)
118c
119      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_eva, zx_tmp_3d)
120      CALL histwrite(nid_bilKPins,"dqeva",itau_w,zx_tmp_3d,
121     .                                   iim*jjmp1*klev,ndex3d)
122c
123      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat, zx_tmp_3d)
124      CALL histwrite(nid_bilKPins,"dtswr",itau_w,zx_tmp_3d,
125     .                                   iim*jjmp1*klev,ndex3d)
126c
127      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat0, zx_tmp_3d)
128      CALL histwrite(nid_bilKPins,"dtsw0",itau_w,zx_tmp_3d,
129     .                                   iim*jjmp1*klev,ndex3d)
130c
131      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool, zx_tmp_3d)
132      CALL histwrite(nid_bilKPins,"dtlwr",itau_w,zx_tmp_3d,
133     .                                   iim*jjmp1*klev,ndex3d)
134c
135      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool0, zx_tmp_3d)
136      CALL histwrite(nid_bilKPins,"dtlw0",itau_w,zx_tmp_3d,
137     .                                   iim*jjmp1*klev,ndex3d)
138c
139      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
140      CALL histwrite(nid_bilKPins,"duvdf",itau_w,zx_tmp_3d,
141     .                                   iim*jjmp1*klev,ndex3d)
142c
143      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d)
144      CALL histwrite(nid_bilKPins,"dvvdf",itau_w,zx_tmp_3d,
145     .                                   iim*jjmp1*klev,ndex3d)
146c
147      IF (ok_orodr) THEN
148      IF (ok_orolf) THEN
149c
150      DO k = 1, klev
151      DO i = 1, klon
152        d_u_oli(i,k) = d_u_oro(i,k) + d_u_lif(i,k)
153        d_v_oli(i,k) = d_v_oro(i,k) + d_v_lif(i,k)
154      ENDDO
155      ENDDO
156c
157      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_oli, zx_tmp_3d)
158      CALL histwrite(nid_bilKPins,"duoli",itau_w,zx_tmp_3d,
159     .                                   iim*jjmp1*klev,ndex3d)
160c
161      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_oli, zx_tmp_3d)
162      CALL histwrite(nid_bilKPins,"dvoli",itau_w,zx_tmp_3d,
163     .                                   iim*jjmp1*klev,ndex3d)
164c
165      ENDIF
166      ENDIF
167C
168      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u, zx_tmp_3d)
169      CALL histwrite(nid_bilKPins,"duphy",itau_w,zx_tmp_3d,
170     .                                   iim*jjmp1*klev,ndex3d)
171c
172      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v, zx_tmp_3d)
173      CALL histwrite(nid_bilKPins,"dvphy",itau_w,zx_tmp_3d,
174     .                                   iim*jjmp1*klev,ndex3d)
175c
176      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
177      CALL histwrite(nid_bilKPins,"dtphy",itau_w,zx_tmp_3d,
178     .                                   iim*jjmp1*klev,ndex3d)
179c
180      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,1), 
181     .zx_tmp_3d)
182      CALL histwrite(nid_bilKPins,"dqphy",itau_w,zx_tmp_3d,
183     .                                   iim*jjmp1*klev,ndex3d)
184c
185      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,2), 
186     .zx_tmp_3d)
187      CALL histwrite(nid_bilKPins,"dqlphy",itau_w,zx_tmp_3d,
188     .                                   iim*jjmp1*klev,ndex3d)
189c
190cIM 280405 BEG
191c
192c Champs 2D:
193c
194c   Ecriture de champs dynamiques sur des niveaux de pression
195c     DO k=1, nlevSTD
196      DO k=1, 12
197c
198       IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
199       IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
200c
201       IF(bb2.EQ."850") THEN
202c
203        CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,1),zx_tmp_2d)
204        CALL histwrite(nid_bilKPins,"u"//bb2,itau_w,zx_tmp_2d,
205     $       iim*jjmp1,ndex2d)
206c
207        CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,1),zx_tmp_2d)
208        CALL histwrite(nid_bilKPins,"v"//bb2,itau_w,zx_tmp_2d,
209     $       iim*jjmp1,ndex2d)
210c
211       ENDIF !(bb2.EQ."850")
212c
213       ENDDO !k=1, 12
214c
215cIM 280405 END
216C
217      if (ok_sync) then
218        call histsync(nid_bilKPins)
219      endif
220       ENDIF
221
Note: See TracBrowser for help on using the repository browser.