6 |
|
|
7 |
contains |
contains |
8 |
|
|
9 |
subroutine ini_histins(dtime) |
subroutine ini_histins(dtime, ok_newmicro) |
10 |
|
|
11 |
! From phylmd/ini_histins.h, version 1.2, 2005/05/25 13:10:09 |
! From phylmd/ini_histins.h, version 1.2, 2005/05/25 13:10:09 |
12 |
|
|
26 |
USE ymds2ju_m, only: ymds2ju |
USE ymds2ju_m, only: ymds2ju |
27 |
|
|
28 |
REAL, intent(in):: dtime ! pas temporel de la physique (s) |
REAL, intent(in):: dtime ! pas temporel de la physique (s) |
29 |
|
logical, intent(in):: ok_newmicro |
30 |
|
|
31 |
! Local: |
! Local: |
32 |
real zjulian, zsto, zout |
real zjulian, zsto, zout |
33 |
integer nhori, nvert, nsrf, iq, it |
integer nhori, nvert, nsrf, iq, it |
36 |
|
|
37 |
print *, 'Call sequence information: ini_histins' |
print *, 'Call sequence information: ini_histins' |
38 |
|
|
39 |
IF (ok_instan) THEN |
test_ok_instan: IF (ok_instan) THEN |
40 |
zsto = dtime * ecrit_ins |
zsto = dtime * ecrit_ins |
41 |
zout = dtime * ecrit_ins |
zout = dtime * ecrit_ins |
42 |
CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian) |
CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian) |
130 |
CALL histdef(nid_ins, "dtsvdfi", "Boundary-layer dTs(g)", "K/s", & |
CALL histdef(nid_ins, "dtsvdfi", "Boundary-layer dTs(g)", "K/s", & |
131 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
132 |
"inst(X)", zsto, zout) |
"inst(X)", zsto, zout) |
133 |
|
CALL histdef(nid_ins, "msnow", "surface snow amount", "kg/m2", & |
134 |
|
iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout) |
135 |
|
|
136 |
DO nsrf = 1, nbsrf |
DO nsrf = 1, nbsrf |
137 |
call histdef(nid_ins, "pourc_"//clnsurf(nsrf), & |
call histdef(nid_ins, "pourc_"//clnsurf(nsrf), & |
258 |
CALL histdef(nid_ins, "dtlw0", "CS LW radiation dT", & |
CALL histdef(nid_ins, "dtlw0", "CS LW radiation dT", & |
259 |
"K/s", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", & |
"K/s", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", & |
260 |
zsto, zout) |
zsto, zout) |
261 |
|
if (ok_newmicro) then |
262 |
|
CALL histdef(nid_ins, "re", "cloud droplet effective radius", & |
263 |
|
"micrometer", iim, jjm + 1, nhori, llm, 1, llm, nvert, & |
264 |
|
"inst(X)", zsto, zout) |
265 |
|
CALL histdef(nid_ins, "fl", & |
266 |
|
"denominator of Cloud droplet effective radius", "", iim, & |
267 |
|
jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", zsto, zout) |
268 |
|
end if |
269 |
|
|
270 |
DO it = 1, nqmx - 2 |
DO it = 1, nqmx - 2 |
271 |
! champ 2D |
! champ 2D |
290 |
ENDDO |
ENDDO |
291 |
|
|
292 |
CALL histend(nid_ins) |
CALL histend(nid_ins) |
293 |
ENDIF |
ENDIF test_ok_instan |
294 |
|
|
295 |
end subroutine ini_histins |
end subroutine ini_histins |
296 |
|
|