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 |
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 |
integer physid |
41 |
integer, intent(in):: itap |
integer, intent(in):: itap |
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: |
116 |
c====================================================================== |
c====================================================================== |
117 |
|
|
118 |
ok_sync = .true. |
ok_sync = .true. |
|
print*,'Dans phystokenc.F' |
|
|
print*,'iadvtr= ',iadvtr |
|
|
print*,'istphy= ',istphy |
|
|
print*,'istdyn= ',istdyn |
|
119 |
|
|
120 |
IF (iadvtr.eq.0) THEN |
IF (iadvtr.eq.0) THEN |
|
|
|
121 |
CALL initphysto('phystoke', |
CALL initphysto('phystoke', |
122 |
. rlon,rlat,dtime, dtime*istphy,dtime*istphy,nqmx,physid) |
. rlon,rlat,dtime, dtime*istphy,dtime*istphy,nqmx,physid) |
|
|
|
|
write(*,*) 'apres initphysto ds phystokenc' |
|
|
|
|
|
|
|
123 |
ENDIF |
ENDIF |
124 |
c |
c |
125 |
ndex2d = 0 |
ndex2d = 0 |
126 |
ndex3d = 0 |
ndex3d = 0 |
127 |
i=itap |
i=itap |
128 |
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) |
129 |
CALL histwrite(physid,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex2d) |
CALL histwrite(physid,"phis",i,zx_tmp_2d) |
130 |
c |
c |
131 |
i=itap |
i=itap |
132 |
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) |
133 |
CALL histwrite(physid,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex2d) |
CALL histwrite(physid,"aire",i,zx_tmp_2d) |
134 |
|
|
135 |
iadvtr=iadvtr+1 |
iadvtr=iadvtr+1 |
136 |
c |
c |
237 |
|
|
238 |
ccccc |
ccccc |
239 |
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) |
240 |
CALL histwrite(physid,"t",itap,zx_tmp_3d, |
CALL histwrite(physid,"t",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
241 |
|
|
242 |
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) |
243 |
CALL histwrite(physid,"mfu",itap,zx_tmp_3d, |
CALL histwrite(physid,"mfu",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
244 |
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) |
245 |
CALL histwrite(physid,"mfd",itap,zx_tmp_3d, |
CALL histwrite(physid,"mfd",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
246 |
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) |
247 |
CALL histwrite(physid,"en_u",itap,zx_tmp_3d, |
CALL histwrite(physid,"en_u",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
248 |
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) |
249 |
CALL histwrite(physid,"de_u",itap,zx_tmp_3d, |
CALL histwrite(physid,"de_u",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
250 |
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) |
251 |
CALL histwrite(physid,"en_d",itap,zx_tmp_3d, |
CALL histwrite(physid,"en_d",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
252 |
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) |
253 |
CALL histwrite(physid,"de_d",itap,zx_tmp_3d, |
CALL histwrite(physid,"de_d",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
254 |
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) |
255 |
CALL histwrite(physid,"coefh",itap,zx_tmp_3d, |
CALL histwrite(physid,"coefh",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
256 |
|
|
257 |
c ajou... |
c ajou... |
258 |
do k=1,klev |
do k=1,klev |
262 |
enddo |
enddo |
263 |
|
|
264 |
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) |
265 |
CALL histwrite(physid,"fm_th",itap,zx_tmp_3d, |
CALL histwrite(physid,"fm_th",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
266 |
c |
c |
267 |
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) |
268 |
CALL histwrite(physid,"en_th",itap,zx_tmp_3d, |
CALL histwrite(physid,"en_th",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
269 |
cccc |
cccc |
270 |
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) |
271 |
CALL histwrite(physid,"frac_impa",itap,zx_tmp_3d, |
CALL histwrite(physid,"frac_impa",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
272 |
|
|
273 |
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) |
274 |
CALL histwrite(physid,"frac_nucl",itap,zx_tmp_3d, |
CALL histwrite(physid,"frac_nucl",itap,zx_tmp_3d) |
|
. iim*(jjm+1)*klev,ndex3d) |
|
275 |
|
|
276 |
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) |
277 |
CALL histwrite(physid,"pyu1",itap,zx_tmp_2d,iim*(jjm+1), |
CALL histwrite(physid,"pyu1",itap,zx_tmp_2d) |
|
. ndex2d) |
|
278 |
|
|
279 |
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) |
280 |
CALL histwrite(physid,"pyv1",itap,zx_tmp_2d,iim*(jjm+1) |
CALL histwrite(physid,"pyv1",itap,zx_tmp_2d) |
|
. ,ndex2d) |
|
281 |
|
|
282 |
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) |
283 |
CALL histwrite(physid,"ftsol1",itap,zx_tmp_2d, |
CALL histwrite(physid,"ftsol1",itap,zx_tmp_2d) |
|
. iim*(jjm+1),ndex2d) |
|
284 |
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) |
285 |
CALL histwrite(physid,"ftsol2",itap,zx_tmp_2d, |
CALL histwrite(physid,"ftsol2",itap,zx_tmp_2d) |
|
. iim*(jjm+1),ndex2d) |
|
286 |
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) |
287 |
CALL histwrite(physid,"ftsol3",itap,zx_tmp_2d, |
CALL histwrite(physid,"ftsol3",itap,zx_tmp_2d) |
|
. iim*(jjm+1),ndex2d) |
|
288 |
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) |
289 |
CALL histwrite(physid,"ftsol4",itap,zx_tmp_2d, |
CALL histwrite(physid,"ftsol4",itap,zx_tmp_2d) |
|
. iim*(jjm+1),ndex2d) |
|
290 |
|
|
291 |
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) |
292 |
CALL histwrite(physid,"psrf1",itap,zx_tmp_2d, |
CALL histwrite(physid,"psrf1",itap,zx_tmp_2d) |
|
. iim*(jjm+1),ndex2d) |
|
293 |
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) |
294 |
CALL histwrite(physid,"psrf2",itap,zx_tmp_2d, |
CALL histwrite(physid,"psrf2",itap,zx_tmp_2d) |
|
. iim*(jjm+1),ndex2d) |
|
295 |
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) |
296 |
CALL histwrite(physid,"psrf3",itap,zx_tmp_2d, |
CALL histwrite(physid,"psrf3",itap,zx_tmp_2d) |
|
. iim*(jjm+1),ndex2d) |
|
297 |
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) |
298 |
CALL histwrite(physid,"psrf4",itap,zx_tmp_2d, |
CALL histwrite(physid,"psrf4",itap,zx_tmp_2d) |
|
. iim*(jjm+1),ndex2d) |
|
299 |
|
|
300 |
if (ok_sync) call histsync(physid) |
if (ok_sync) call histsync(physid) |
301 |
c if (ok_sync) call histsync |
c if (ok_sync) call histsync |