1 | c |
---|
2 | c $Header$ |
---|
3 | c |
---|
4 | IF (ok_journe) THEN |
---|
5 | c |
---|
6 | ndex2d = 0 |
---|
7 | ndex3d = 0 |
---|
8 | c |
---|
9 | itau_w = itau_phy + itap + start_time * day_step / iphysiq |
---|
10 | c |
---|
11 | c Champs 3D: |
---|
12 | c |
---|
13 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, ue_lay,zx_tmp_3d) |
---|
14 | CALL histwrite_phy(nid_bilKPins,"ue",itau_w,ue_lay) |
---|
15 | c |
---|
16 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, ve_lay,zx_tmp_3d) |
---|
17 | CALL histwrite_phy(nid_bilKPins,"ve",itau_w,ve_lay) |
---|
18 | c |
---|
19 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, uq_lay,zx_tmp_3d) |
---|
20 | CALL histwrite_phy(nid_bilKPins,"uq",itau_w,uq_lay) |
---|
21 | c |
---|
22 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, vq_lay,zx_tmp_3d) |
---|
23 | CALL histwrite_phy(nid_bilKPins,"vq",itau_w,vq_lay) |
---|
24 | c |
---|
25 | c Champs 3D: |
---|
26 | C |
---|
27 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) |
---|
28 | CALL histwrite_phy(nid_bilKPins,"temp",itau_w,t_seri) |
---|
29 | c |
---|
30 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d) |
---|
31 | CALL histwrite_phy(nid_bilKPins,"ovap",itau_w,qx(:,:,ivap)) |
---|
32 | c |
---|
33 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) |
---|
34 | CALL histwrite_phy(nid_bilKPins,"geop",itau_w,zphi) |
---|
35 | c |
---|
36 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) |
---|
37 | CALL histwrite_phy(nid_bilKPins,"vitu",itau_w,u_seri) |
---|
38 | c |
---|
39 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) |
---|
40 | CALL histwrite_phy(nid_bilKPins,"vitv",itau_w,v_seri) |
---|
41 | c |
---|
42 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d) |
---|
43 | CALL histwrite_phy(nid_bilKPins,"vitw",itau_w,omega) |
---|
44 | c |
---|
45 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) |
---|
46 | CALL histwrite_phy(nid_bilKPins,"pres",itau_w,pplay) |
---|
47 | c |
---|
48 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, paprs, zx_tmp_3d) |
---|
49 | CALL histwrite_phy(nid_bilKPins,"play",itau_w,paprs) |
---|
50 | c |
---|
51 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, radocond, zx_tmp_3d) |
---|
52 | CALL histwrite_phy(nid_bilKPins,"oliq",itau_w,cldliq) |
---|
53 | c |
---|
54 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d) |
---|
55 | CALL histwrite_phy(nid_bilKPins,"dtdyn",itau_w,d_t_dyn) |
---|
56 | c |
---|
57 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d) |
---|
58 | CALL histwrite_phy(nid_bilKPins,"dqdyn",itau_w,d_q_dyn) |
---|
59 | c |
---|
60 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_con, zx_tmp_3d) |
---|
61 | CALL histwrite_phy(nid_bilKPins,"dtcon",itau_w,d_t_con) |
---|
62 | c |
---|
63 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_con, zx_tmp_3d) |
---|
64 | CALL histwrite_phy(nid_bilKPins,"ducon",itau_w,d_u_con) |
---|
65 | c |
---|
66 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_con, zx_tmp_3d) |
---|
67 | CALL histwrite_phy(nid_bilKPins,"dvcon",itau_w,d_v_con) |
---|
68 | c |
---|
69 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_con, zx_tmp_3d) |
---|
70 | CALL histwrite_phy(nid_bilKPins,"dqcon",itau_w,d_q_con) |
---|
71 | c |
---|
72 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_lsc, zx_tmp_3d) |
---|
73 | CALL histwrite_phy(nid_bilKPins,"dtlsc",itau_w,d_t_lsc) |
---|
74 | c |
---|
75 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_lsc, zx_tmp_3d) |
---|
76 | CALL histwrite_phy(nid_bilKPins,"dqlsc",itau_w,d_q_lsc) |
---|
77 | c |
---|
78 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_vdf, zx_tmp_3d) |
---|
79 | CALL histwrite_phy(nid_bilKPins,"dtvdf",itau_w,d_t_vdf) |
---|
80 | c |
---|
81 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d) |
---|
82 | CALL histwrite_phy(nid_bilKPins,"dqvdf",itau_w,d_q_vdf) |
---|
83 | c |
---|
84 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_ajs, zx_tmp_3d) |
---|
85 | CALL histwrite_phy(nid_bilKPins,"dtajs",itau_w,d_t_ajs) |
---|
86 | c |
---|
87 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_ajs, zx_tmp_3d) |
---|
88 | CALL histwrite_phy(nid_bilKPins,"dqajs",itau_w,d_q_ajs) |
---|
89 | c |
---|
90 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_eva, zx_tmp_3d) |
---|
91 | CALL histwrite_phy(nid_bilKPins,"dteva",itau_w,d_t_eva) |
---|
92 | c |
---|
93 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_eva, zx_tmp_3d) |
---|
94 | CALL histwrite_phy(nid_bilKPins,"dqeva",itau_w,d_q_eva) |
---|
95 | c |
---|
96 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat, zx_tmp_3d) |
---|
97 | CALL histwrite_phy(nid_bilKPins,"dtswr",itau_w,heat) |
---|
98 | c |
---|
99 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat0, zx_tmp_3d) |
---|
100 | CALL histwrite_phy(nid_bilKPins,"dtsw0",itau_w,heat0) |
---|
101 | c |
---|
102 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool, zx_tmp_3d) |
---|
103 | CALL histwrite_phy(nid_bilKPins,"dtlwr",itau_w,cool) |
---|
104 | c |
---|
105 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool0, zx_tmp_3d) |
---|
106 | CALL histwrite_phy(nid_bilKPins,"dtlw0",itau_w,cool0) |
---|
107 | c |
---|
108 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d) |
---|
109 | CALL histwrite_phy(nid_bilKPins,"duvdf",itau_w,d_u_vdf) |
---|
110 | c |
---|
111 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d) |
---|
112 | CALL histwrite_phy(nid_bilKPins,"dvvdf",itau_w,d_v_vdf) |
---|
113 | c |
---|
114 | IF (ok_orodr) THEN |
---|
115 | IF (ok_orolf) THEN |
---|
116 | c |
---|
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 |
---|
123 | c |
---|
124 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_oli, zx_tmp_3d) |
---|
125 | CALL histwrite_phy(nid_bilKPins,"duoli",itau_w,d_u_oli) |
---|
126 | c |
---|
127 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_oli, zx_tmp_3d) |
---|
128 | CALL histwrite_phy(nid_bilKPins,"dvoli",itau_w,d_v_oli) |
---|
129 | c |
---|
130 | ENDIF |
---|
131 | ENDIF |
---|
132 | C |
---|
133 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u, zx_tmp_3d) |
---|
134 | CALL histwrite_phy(nid_bilKPins,"duphy",itau_w,d_u) |
---|
135 | c |
---|
136 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v, zx_tmp_3d) |
---|
137 | CALL histwrite_phy(nid_bilKPins,"dvphy",itau_w,d_v) |
---|
138 | c |
---|
139 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d) |
---|
140 | CALL histwrite_phy(nid_bilKPins,"dtphy",itau_w,d_t) |
---|
141 | c |
---|
142 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,1), |
---|
143 | cym .zx_tmp_3d) |
---|
144 | CALL histwrite_phy(nid_bilKPins,"dqphy",itau_w,d_qx(:,:,1)) |
---|
145 | c |
---|
146 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,2), |
---|
147 | cym .zx_tmp_3d) |
---|
148 | CALL histwrite_phy(nid_bilKPins,"dqlphy",itau_w,d_qx(:,:,2)) |
---|
149 | c |
---|
150 | cIM 280405 BEG |
---|
151 | c |
---|
152 | c Champs 2D: |
---|
153 | c |
---|
154 | c Ecriture de champs dynamiques sur des niveaux de pression |
---|
155 | c DO k=1, nlevSTD |
---|
156 | DO k=1, 12 |
---|
157 | c |
---|
158 | IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k) |
---|
159 | IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k) |
---|
160 | c |
---|
161 | IF(bb2.EQ."850") THEN |
---|
162 | c |
---|
163 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,1),zx_tmp_2d) |
---|
164 | CALL histwrite_phy(nid_bilKPins,"u"//bb2,itau_w,usumSTD(:,k,1)) |
---|
165 | c |
---|
166 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,1),zx_tmp_2d) |
---|
167 | CALL histwrite_phy(nid_bilKPins,"v"//bb2,itau_w,vsumSTD(:,k,1)) |
---|
168 | c |
---|
169 | ENDIF !(bb2.EQ."850") |
---|
170 | c |
---|
171 | ENDDO !k=1, 12 |
---|
172 | c |
---|
173 | cIM 280405 END |
---|
174 | C |
---|
175 | if (ok_sync) then |
---|
176 | call histsync(nid_bilKPins) |
---|
177 | endif |
---|
178 | ENDIF |
---|
179 | |
---|