source: LMDZ4/branches/pre_V3/libf/phylmd/write_bilKP_ins.h

Last change on this file was 644, checked in by Laurent Fairhead, 20 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.