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 |
48 |
CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", & |
CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", & |
49 |
presnivs/100., nvert) |
presnivs/100., nvert) |
50 |
|
|
51 |
CALL histdef(nid_ins, "phis", "Surface geop. height", "-", & |
CALL histdef(nid_ins, "phis", "surface geopotential", "m2 s-2", & |
52 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
53 |
"once", zsto, zout) |
"once", zsto, zout) |
54 |
CALL histdef(nid_ins, "aire", "Grid area", "-", & |
CALL histdef(nid_ins, "aire", "Grid area", "-", & |
170 |
"rugosite "//clnsurf(nsrf), "-", & |
"rugosite "//clnsurf(nsrf), "-", & |
171 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
172 |
"inst(X)", zsto, zout) |
"inst(X)", zsto, zout) |
173 |
|
call histdef(nid_ins, "u10m_"//clnsurf(nsrf), & |
174 |
|
"zonal wind 10 m "//clnsurf(nsrf), "m s-1", & |
175 |
|
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
176 |
|
"inst(X)", zsto, zout) |
177 |
|
call histdef(nid_ins, "v10m_"//clnsurf(nsrf), & |
178 |
|
"meridional wind 10 m "//clnsurf(nsrf), "m s-1", & |
179 |
|
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
180 |
|
"inst(X)", zsto, zout) |
181 |
END DO |
END DO |
182 |
|
|
183 |
CALL histdef(nid_ins, "rugs", "rugosity", "-", & |
CALL histdef(nid_ins, "rugs", "rugosity", "-", & |
208 |
CALL histdef(nid_ins, "s_therm", "Exces du thermique", "K", & |
CALL histdef(nid_ins, "s_therm", "Exces du thermique", "K", & |
209 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
210 |
"inst(X)", zsto, zout) |
"inst(X)", zsto, zout) |
211 |
CALL histdef(nid_ins, "s_trmb1", "deep_cape(HBTM2)", "J/m2", & |
CALL histdef(nid_ins, "qsurf", "Surface Air humidity", "", & |
|
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
|
|
"inst(X)", zsto, zout) |
|
|
CALL histdef(nid_ins, "s_trmb2", "inhibition (HBTM2)", "J/m2", & |
|
|
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
|
|
"inst(X)", zsto, zout) |
|
|
CALL histdef(nid_ins, "s_trmb3", "Point Omega (HBTM2)", "m", & |
|
212 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
213 |
"inst(X)", zsto, zout) |
"inst(X)", zsto, zout) |
214 |
|
|
261 |
"K/s", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", & |
"K/s", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", & |
262 |
zsto, zout) |
zsto, zout) |
263 |
|
|
264 |
|
if (ok_newmicro) then |
265 |
|
CALL histdef(nid_ins, "re", "cloud droplet effective radius", & |
266 |
|
"micrometer", iim, jjm + 1, nhori, llm, 1, llm, nvert, & |
267 |
|
"inst(X)", zsto, zout) |
268 |
|
CALL histdef(nid_ins, "fl", & |
269 |
|
"denominator of Cloud droplet effective radius", "", iim, & |
270 |
|
jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", zsto, zout) |
271 |
|
end if |
272 |
|
|
273 |
DO it = 1, nqmx - 2 |
DO it = 1, nqmx - 2 |
274 |
! champ 2D |
! champ 2D |
275 |
iq=it+2 |
iq=it+2 |