4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
subroutine ini_histins(dtime, ok_instan, nid_ins) |
subroutine ini_histins(dtime, ok_instan, nid_ins, itau_phy) |
8 |
|
|
9 |
! From phylmd/ini_histins.h, v 1.2 2005/05/25 13:10:09 |
! From phylmd/ini_histins.h, v 1.2 2005/05/25 13:10:09 |
10 |
|
|
11 |
use clesphys, only: ecrit_ins |
use clesphys, only: ecrit_ins |
12 |
use dimens_m, only: iim, jjm, llm |
use dimens_m, only: iim, jjm, llm, nqmx |
13 |
use dimphy, only: klon |
use dimphy, only: klon |
14 |
use disvert_m, only: presnivs |
use disvert_m, only: presnivs |
15 |
use dynetat0_m, only: day_ref, annee_ref |
use dynetat0_m, only: day_ref, annee_ref |
18 |
USE histend_m, ONLY : histend |
USE histend_m, ONLY : histend |
19 |
USE histvert_m, ONLY : histvert |
USE histvert_m, ONLY : histvert |
20 |
use indicesol, only: nbsrf, clnsurf |
use indicesol, only: nbsrf, clnsurf |
21 |
|
use iniadvtrac_m, only: tname, ttext |
22 |
use phyetat0_m, only: rlon, rlat |
use phyetat0_m, only: rlon, rlat |
|
use temps, only: itau_phy |
|
23 |
USE ymds2ju_m, only: ymds2ju |
USE ymds2ju_m, only: ymds2ju |
24 |
|
|
25 |
REAL, intent(in):: dtime ! pas temporel de la physique (s) |
REAL, intent(in):: dtime ! pas temporel de la physique (s) |
26 |
logical, intent(in):: ok_instan |
logical, intent(in):: ok_instan |
27 |
integer, intent(out):: nid_ins |
integer, intent(out):: nid_ins |
28 |
|
integer, intent(in):: itau_phy |
29 |
|
|
30 |
|
! Local: |
31 |
REAL zx_lon(iim, jjm + 1), zx_lat(iim, jjm + 1) |
REAL zx_lon(iim, jjm + 1), zx_lat(iim, jjm + 1) |
32 |
real zjulian, zsto, zout |
real zjulian, zsto, zout |
33 |
integer i, nhori, nvert, nsrf |
integer i, nhori, nvert, nsrf, iq, it |
34 |
|
|
35 |
!------------------------------------------------------------------- |
!------------------------------------------------------------------- |
36 |
|
|
37 |
IF (ok_instan) THEN |
IF (ok_instan) THEN |
|
|
|
38 |
zsto = dtime * ecrit_ins |
zsto = dtime * ecrit_ins |
39 |
zout = dtime * ecrit_ins |
zout = dtime * ecrit_ins |
40 |
|
|
218 |
CALL histdef(nid_ins, "albs", "Surface albedo", "-", & |
CALL histdef(nid_ins, "albs", "Surface albedo", "-", & |
219 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
220 |
"inst(X)", zsto, zout) |
"inst(X)", zsto, zout) |
|
CALL histdef(nid_ins, "albslw", "Surface albedo LW", "-", & |
|
|
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
|
|
"inst(X)", zsto, zout) |
|
221 |
|
|
|
!IM cf. AM 081204 BEG |
|
|
! HBTM2 |
|
222 |
CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m", & |
CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m", & |
223 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
224 |
"inst(X)", zsto, zout) |
"inst(X)", zsto, zout) |
292 |
iim, (jjm + 1), nhori, llm, 1, llm, nvert, & |
iim, (jjm + 1), nhori, llm, 1, llm, nvert, & |
293 |
"inst(X)", zsto, zout) |
"inst(X)", zsto, zout) |
294 |
|
|
295 |
|
CALL histdef(nid_ins, "zmasse", "column density of air in cell", & |
296 |
|
"kg m-2", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", & |
297 |
|
zsto, zout) |
298 |
|
|
299 |
|
DO it = 1, nqmx - 2 |
300 |
|
! champ 2D |
301 |
|
iq=it+2 |
302 |
|
CALL histdef(nid_ins, tname(iq), ttext(iq), "U/kga", iim, jjm+1, & |
303 |
|
nhori, llm, 1, llm, nvert, "inst(X)", zsto, zout) |
304 |
|
CALL histdef(nid_ins, "fl"//tname(iq), "Flux "//ttext(iq), & |
305 |
|
"U/m2/s", iim, jjm+1, nhori, llm, 1, llm, nvert, & |
306 |
|
"inst(X)", zsto, zout) |
307 |
|
|
308 |
|
!---Ajout Olivia |
309 |
|
CALL histdef(nid_ins, "d_tr_th_"//tname(iq), & |
310 |
|
"tendance thermique"// ttext(iq), "?", & |
311 |
|
iim, jjm+1, nhori, llm, 1, llm, nvert, & |
312 |
|
"inst(X)", zsto, zout) |
313 |
|
CALL histdef(nid_ins, "d_tr_cv_"//tname(iq), & |
314 |
|
"tendance convection"// ttext(iq), "?", & |
315 |
|
iim, jjm+1, nhori, llm, 1, llm, nvert, & |
316 |
|
"inst(X)", zsto, zout) |
317 |
|
CALL histdef(nid_ins, "d_tr_cl_"//tname(iq), & |
318 |
|
"tendance couche limite"// ttext(iq), "?", & |
319 |
|
iim, jjm+1, nhori, llm, 1, llm, nvert, & |
320 |
|
"inst(X)", zsto, zout) |
321 |
|
!---fin Olivia |
322 |
|
|
323 |
|
ENDDO |
324 |
|
|
325 |
CALL histend(nid_ins) |
CALL histend(nid_ins) |
326 |
ENDIF |
ENDIF |
327 |
|
|