source: LMDZ5/branches/IPSLCM5A2.1_ISO/libf/phyiso/write_bilKP_ave.h

Last change on this file was 3331, checked in by acozic, 7 years ago

Add modification for isotopes

  • Property svn:executable set to *
File size: 5.6 KB
Line 
1c
2c $Header$
3c
4      IF (ok_journe) THEN
5c
6      ndex2d = 0
7      ndex3d = 0
8c
9c Champs 2D:
10c
11      itau_w = itau_phy + itap + start_time * day_step / iphysiq
12c
13cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, ue_lay,zx_tmp_3d)
14      CALL histwrite_phy(nid_bilKPave,"ue",itau_w,ue_lay)
15c
16cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, ve_lay,zx_tmp_3d)
17      CALL histwrite_phy(nid_bilKPave,"ve",itau_w,ve_lay)
18c
19cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, uq_lay,zx_tmp_3d)
20      CALL histwrite_phy(nid_bilKPave,"uq",itau_w,uq_lay)
21c
22cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, vq_lay,zx_tmp_3d)
23      CALL histwrite_phy(nid_bilKPave,"vq",itau_w,vq_lay)
24c
25c Champs 3D:
26C
27cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
28      CALL histwrite_phy(nid_bilKPave,"temp",itau_w,t_seri)
29c
30cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
31      CALL histwrite_phy(nid_bilKPave,"ovap",itau_w,qx(:,:,ivap))
32c
33cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
34      CALL histwrite_phy(nid_bilKPave,"geop",itau_w,zphi)
35c
36cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
37      CALL histwrite_phy(nid_bilKPave,"vitu",itau_w,u_seri)
38c
39cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
40      CALL histwrite_phy(nid_bilKPave,"vitv",itau_w,v_seri)
41c
42cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
43      CALL histwrite_phy(nid_bilKPave,"vitw",itau_w,omega)
44c
45cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
46      CALL histwrite_phy(nid_bilKPave,"pres",itau_w,pplay)
47c
48cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, paprs, zx_tmp_3d)
49      CALL histwrite_phy(nid_bilKPave,"play",itau_w,paprs)
50c
51cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldliq, zx_tmp_3d)
52      CALL histwrite_phy(nid_bilKPave,"oliq",itau_w,cldliq)
53c
54cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
55      CALL histwrite_phy(nid_bilKPave,"dtdyn",itau_w,d_t_dyn)
56c
57cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
58      CALL histwrite_phy(nid_bilKPave,"dqdyn",itau_w,d_q_dyn)
59c
60cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_con, zx_tmp_3d)
61      CALL histwrite_phy(nid_bilKPave,"dtcon",itau_w,d_t_con)
62c
63cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_con, zx_tmp_3d)
64      CALL histwrite_phy(nid_bilKPave,"ducon",itau_w,d_u_con)
65c
66cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_con, zx_tmp_3d)
67      CALL histwrite_phy(nid_bilKPave,"dvcon",itau_w,d_v_con)
68c
69cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_con, zx_tmp_3d)
70      CALL histwrite_phy(nid_bilKPave,"dqcon",itau_w,d_q_con)
71c
72cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_lsc, zx_tmp_3d)
73      CALL histwrite_phy(nid_bilKPave,"dtlsc",itau_w,d_t_lsc)
74c
75cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_lsc, zx_tmp_3d)
76      CALL histwrite_phy(nid_bilKPave,"dqlsc",itau_w,d_q_lsc)
77c
78cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_vdf, zx_tmp_3d)
79      CALL histwrite_phy(nid_bilKPave,"dtvdf",itau_w,d_t_vdf)
80c
81cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d)
82      CALL histwrite_phy(nid_bilKPave,"dqvdf",itau_w,d_q_vdf)
83c
84cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_ajs, zx_tmp_3d)
85      CALL histwrite_phy(nid_bilKPave,"dtajs",itau_w,d_t_ajs)
86c
87cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_ajs, zx_tmp_3d)
88      CALL histwrite_phy(nid_bilKPave,"dqajs",itau_w,d_q_ajs)
89c
90cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_eva, zx_tmp_3d)
91      CALL histwrite_phy(nid_bilKPave,"dteva",itau_w,d_t_eva)
92c
93cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_eva, zx_tmp_3d)
94      CALL histwrite_phy(nid_bilKPave,"dqeva",itau_w,d_q_eva)
95c
96cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat, zx_tmp_3d)
97      CALL histwrite_phy(nid_bilKPave,"dtswr",itau_w,heat)
98c
99cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat0, zx_tmp_3d)
100      CALL histwrite_phy(nid_bilKPave,"dtsw0",itau_w,heat0)
101c
102cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool, zx_tmp_3d)
103      CALL histwrite_phy(nid_bilKPave,"dtlwr",itau_w,cool)
104c
105cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool0, zx_tmp_3d)
106      CALL histwrite_phy(nid_bilKPave,"dtlw0",itau_w,cool0)
107c
108cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
109      CALL histwrite_phy(nid_bilKPave,"duvdf",itau_w,d_u_vdf)
110c
111cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d)
112      CALL histwrite_phy(nid_bilKPave,"dvvdf",itau_w,d_v_vdf)
113c
114      IF (ok_orodr) THEN
115      IF (ok_orolf) THEN
116c
117      DO k = 1, klev
118      DO i = 1, klon
119        d_u_oli(i,k) = d_u_oro(i,k) + d_u_lif(i,k)
120        d_v_oli(i,k) = d_v_oro(i,k) + d_v_lif(i,k)
121      ENDDO
122      ENDDO
123c
124cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_oli, zx_tmp_3d)
125      CALL histwrite_phy(nid_bilKPave,"duoli",d_u_oli)
126c
127cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_oli, zx_tmp_3d)
128      CALL histwrite_phy(nid_bilKPave,"dvoli",itau_w,d_v_oli)
129c
130      ENDIF
131      ENDIF
132C
133cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u, zx_tmp_3d)
134      CALL histwrite_phy(nid_bilKPave,"duphy",itau_w,d_u)
135c
136cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v, zx_tmp_3d)
137      CALL histwrite_phy(nid_bilKPave,"dvphy",itau_w,d_v)
138c
139cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
140      CALL histwrite_phy(nid_bilKPave,"dtphy",itau_w,d_t)
141c
142cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,1), 
143cymf     .zx_tmp_3d)
144      CALL histwrite_phy(nid_bilKPave,"dqphy",itau_w,d_qx(:,:,1))
145c
146cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,2), 
147cym     .zx_tmp_3d)
148      CALL histwrite_phy(nid_bilKPave,"dqlphy",itau_w,d_qx(:,:,2))
149c
150C
151      if (ok_sync) then
152        call histsync(nid_bilKPave)
153      endif
154       ENDIF
155
Note: See TracBrowser for help on using the repository browser.