10 |
I pcoefh,yu1,yv1,ftsol,pctsrf, |
I pcoefh,yu1,yv1,ftsol,pctsrf, |
11 |
I frac_impa,frac_nucl, |
I frac_impa,frac_nucl, |
12 |
I pphis,paire,dtime,itap) |
I pphis,paire,dtime,itap) |
13 |
USE ioipsl |
USE histwrite_m |
14 |
USE histcom |
use histcom |
|
|
|
15 |
use dimens_m |
use dimens_m |
16 |
use indicesol |
use indicesol |
17 |
use dimphy |
use dimphy |
35 |
c divers: |
c divers: |
36 |
c ------- |
c ------- |
37 |
c |
c |
38 |
real pdtphys ! pas d'integration pour la physique (seconde) |
real, intent(in):: pdtphys ! pas d'integration pour la physique (seconde) |
39 |
c |
c |
40 |
integer physid, itap |
integer physid |
41 |
|
integer, intent(in):: itap |
42 |
save physid |
save physid |
43 |
integer ndex2d(iim*(jjm+1)),ndex3d(iim*(jjm+1)*klev) |
integer ndex2d(iim*(jjm+1)),ndex3d(iim*(jjm+1)*klev) |
44 |
|
|
54 |
real pt(klon,klev),t(klon,klev) |
real pt(klon,klev),t(klon,klev) |
55 |
c |
c |
56 |
REAL, intent(in):: rlon(klon), rlat(klon) |
REAL, intent(in):: rlon(klon), rlat(klon) |
57 |
real dtime |
real, intent(in):: dtime |
58 |
REAL zx_tmp_3d(iim,jjm+1,klev),zx_tmp_2d(iim,jjm+1) |
REAL zx_tmp_3d(iim,jjm+1,klev),zx_tmp_2d(iim,jjm+1) |
59 |
|
|
60 |
c Couche limite: |
c Couche limite: |
135 |
ndex3d = 0 |
ndex3d = 0 |
136 |
i=itap |
i=itap |
137 |
CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d) |
CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d) |
138 |
CALL histwrite(physid,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex2d) |
CALL histwrite(physid,"phis",i,zx_tmp_2d) |
139 |
c |
c |
140 |
i=itap |
i=itap |
141 |
CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d) |
CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d) |
142 |
CALL histwrite(physid,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex2d) |
CALL histwrite(physid,"aire",i,zx_tmp_2d) |
143 |
|
|
144 |
iadvtr=iadvtr+1 |
iadvtr=iadvtr+1 |
145 |
c |
c |
246 |
|
|
247 |
ccccc |
ccccc |
248 |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, t, zx_tmp_3d) |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, t, zx_tmp_3d) |
249 |
CALL histwrite(physid,"t",itap,zx_tmp_3d, |
CALL histwrite(physid,"t",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
250 |
|
|
251 |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, mfu, zx_tmp_3d) |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, mfu, zx_tmp_3d) |
252 |
CALL histwrite(physid,"mfu",itap,zx_tmp_3d, |
CALL histwrite(physid,"mfu",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
253 |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, mfd, zx_tmp_3d) |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, mfd, zx_tmp_3d) |
254 |
CALL histwrite(physid,"mfd",itap,zx_tmp_3d, |
CALL histwrite(physid,"mfd",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
255 |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, en_u, zx_tmp_3d) |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, en_u, zx_tmp_3d) |
256 |
CALL histwrite(physid,"en_u",itap,zx_tmp_3d, |
CALL histwrite(physid,"en_u",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
257 |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, de_u, zx_tmp_3d) |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, de_u, zx_tmp_3d) |
258 |
CALL histwrite(physid,"de_u",itap,zx_tmp_3d, |
CALL histwrite(physid,"de_u",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
259 |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, en_d, zx_tmp_3d) |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, en_d, zx_tmp_3d) |
260 |
CALL histwrite(physid,"en_d",itap,zx_tmp_3d, |
CALL histwrite(physid,"en_d",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
261 |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, de_d, zx_tmp_3d) |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, de_d, zx_tmp_3d) |
262 |
CALL histwrite(physid,"de_d",itap,zx_tmp_3d, |
CALL histwrite(physid,"de_d",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
263 |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, coefh, zx_tmp_3d) |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, coefh, zx_tmp_3d) |
264 |
CALL histwrite(physid,"coefh",itap,zx_tmp_3d, |
CALL histwrite(physid,"coefh",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
265 |
|
|
266 |
c ajou... |
c ajou... |
267 |
do k=1,klev |
do k=1,klev |
271 |
enddo |
enddo |
272 |
|
|
273 |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, fm_therm1, zx_tmp_3d) |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, fm_therm1, zx_tmp_3d) |
274 |
CALL histwrite(physid,"fm_th",itap,zx_tmp_3d, |
CALL histwrite(physid,"fm_th",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
275 |
c |
c |
276 |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, entr_therm, zx_tmp_3d) |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1, entr_therm, zx_tmp_3d) |
277 |
CALL histwrite(physid,"en_th",itap,zx_tmp_3d, |
CALL histwrite(physid,"en_th",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
278 |
cccc |
cccc |
279 |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1,frac_impa,zx_tmp_3d) |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1,frac_impa,zx_tmp_3d) |
280 |
CALL histwrite(physid,"frac_impa",itap,zx_tmp_3d, |
CALL histwrite(physid,"frac_impa",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
281 |
|
|
282 |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1,frac_nucl,zx_tmp_3d) |
CALL gr_fi_ecrit(klev,klon,iim,jjm+1,frac_nucl,zx_tmp_3d) |
283 |
CALL histwrite(physid,"frac_nucl",itap,zx_tmp_3d, |
CALL histwrite(physid,"frac_nucl",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
284 |
|
|
285 |
CALL gr_fi_ecrit(1, klon,iim,jjm+1, pyu1,zx_tmp_2d) |
CALL gr_fi_ecrit(1, klon,iim,jjm+1, pyu1,zx_tmp_2d) |
286 |
CALL histwrite(physid,"pyu1",itap,zx_tmp_2d,iim*(jjm+1), |
CALL histwrite(physid,"pyu1",itap,zx_tmp_2d) |
|
. ndex2d) |
|
287 |
|
|
288 |
CALL gr_fi_ecrit(1, klon,iim,jjm+1, pyv1,zx_tmp_2d) |
CALL gr_fi_ecrit(1, klon,iim,jjm+1, pyv1,zx_tmp_2d) |
289 |
CALL histwrite(physid,"pyv1",itap,zx_tmp_2d,iim*(jjm+1) |
CALL histwrite(physid,"pyv1",itap,zx_tmp_2d) |
|
. ,ndex2d) |
|
290 |
|
|
291 |
CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol1, zx_tmp_2d) |
CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol1, zx_tmp_2d) |
292 |
CALL histwrite(physid,"ftsol1",itap,zx_tmp_2d, |
CALL histwrite(physid,"ftsol1",itap,zx_tmp_2d) |
|
. iim*(jjm+1),ndex2d) |
|
293 |
CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol2, zx_tmp_2d) |
CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol2, zx_tmp_2d) |
294 |
CALL histwrite(physid,"ftsol2",itap,zx_tmp_2d, |
CALL histwrite(physid,"ftsol2",itap,zx_tmp_2d) |
|
. iim*(jjm+1),ndex2d) |
|
295 |
CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol3, zx_tmp_2d) |
CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol3, zx_tmp_2d) |
296 |
CALL histwrite(physid,"ftsol3",itap,zx_tmp_2d, |
CALL histwrite(physid,"ftsol3",itap,zx_tmp_2d) |
|
. iim*(jjm+1),ndex2d) |
|
297 |
CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol4, zx_tmp_2d) |
CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol4, zx_tmp_2d) |
298 |
CALL histwrite(physid,"ftsol4",itap,zx_tmp_2d, |
CALL histwrite(physid,"ftsol4",itap,zx_tmp_2d) |
|
. iim*(jjm+1),ndex2d) |
|
299 |
|
|
300 |
CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf1, zx_tmp_2d) |
CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf1, zx_tmp_2d) |
301 |
CALL histwrite(physid,"psrf1",itap,zx_tmp_2d, |
CALL histwrite(physid,"psrf1",itap,zx_tmp_2d) |
|
. iim*(jjm+1),ndex2d) |
|
302 |
CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf2, zx_tmp_2d) |
CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf2, zx_tmp_2d) |
303 |
CALL histwrite(physid,"psrf2",itap,zx_tmp_2d, |
CALL histwrite(physid,"psrf2",itap,zx_tmp_2d) |
|
. iim*(jjm+1),ndex2d) |
|
304 |
CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf3, zx_tmp_2d) |
CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf3, zx_tmp_2d) |
305 |
CALL histwrite(physid,"psrf3",itap,zx_tmp_2d, |
CALL histwrite(physid,"psrf3",itap,zx_tmp_2d) |
|
. iim*(jjm+1),ndex2d) |
|
306 |
CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf4, zx_tmp_2d) |
CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf4, zx_tmp_2d) |
307 |
CALL histwrite(physid,"psrf4",itap,zx_tmp_2d, |
CALL histwrite(physid,"psrf4",itap,zx_tmp_2d) |
|
. iim*(jjm+1),ndex2d) |
|
308 |
|
|
309 |
if (ok_sync) call histsync(physid) |
if (ok_sync) call histsync(physid) |
310 |
c if (ok_sync) call histsync |
c if (ok_sync) call histsync |