9 |
|
|
10 |
! From phylmd/ini_histrac.h, version 1.10 2006/02/21 08:08:30 |
! From phylmd/ini_histrac.h, version 1.10 2006/02/21 08:08:30 |
11 |
|
|
12 |
|
use clesphys, only: ecrit_tra |
13 |
use dimens_m, only: iim, jjm, llm |
use dimens_m, only: iim, jjm, llm |
14 |
USE calendar, only: ymds2ju |
use disvert_m, only: presnivs |
15 |
USE histbeg_totreg_m, ONLY : histbeg_totreg |
use dimphy, only: klon |
16 |
|
use grid_change, only: gr_phy_write_2d |
17 |
|
USE histbeg_totreg_m, ONLY: histbeg_totreg |
18 |
USE histdef_m, ONLY : histdef |
USE histdef_m, ONLY : histdef |
19 |
USE histend_m, ONLY : histend |
USE histend_m, ONLY : histend |
20 |
USE histvert_m, ONLY : histvert |
USE histvert_m, ONLY : histvert |
21 |
use temps, only: annee_ref, day_ref, itau_phy |
use iniadvtrac_m, only: tnom, ttext |
|
use iniadvtrac_m, only: niadv, tnom, ttext |
|
|
use dimphy, only: klon |
|
|
use clesphys, only: ecrit_tra |
|
|
use grid_change, only: gr_phy_write_2d |
|
22 |
use phyetat0_m, only: rlon, rlat |
use phyetat0_m, only: rlon, rlat |
23 |
use disvert_m, only: presnivs |
use temps, only: annee_ref, day_ref, itau_phy |
24 |
|
USE ymds2ju_m, only: ymds2ju |
25 |
|
|
26 |
INTEGER, intent(out):: nid_tra |
INTEGER, intent(out):: nid_tra |
27 |
real, intent(in):: pdtphys ! pas d'integration pour la physique (s) |
real, intent(in):: pdtphys ! pas d'integration pour la physique (s) |
37 |
REAL zx_lat(iim, jjm+1) |
REAL zx_lat(iim, jjm+1) |
38 |
INTEGER nhori, nvert |
INTEGER nhori, nvert |
39 |
REAL zsto, zout |
REAL zsto, zout |
40 |
integer it, iq, iiq |
integer it, iq |
41 |
|
|
42 |
!--------------------------------------------------------- |
!--------------------------------------------------------- |
43 |
|
|
63 |
DO it = 1, nq_phys |
DO it = 1, nq_phys |
64 |
! champ 2D |
! champ 2D |
65 |
iq=it+2 |
iq=it+2 |
66 |
iiq=niadv(iq) |
CALL histdef(nid_tra, tnom(iq), ttext(iq), "U/kga", iim, jjm+1, & |
|
CALL histdef(nid_tra, tnom(iq), ttext(iiq), "U/kga", iim, jjm+1, & |
|
67 |
nhori, llm, 1, llm, nvert, "ave(X)", zsto, zout) |
nhori, llm, 1, llm, nvert, "ave(X)", zsto, zout) |
68 |
if (lessivage) THEN |
if (lessivage) THEN |
69 |
CALL histdef(nid_tra, "fl"//tnom(iq), "Flux "//ttext(iiq), & |
CALL histdef(nid_tra, "fl"//tnom(iq), "Flux "//ttext(iq), & |
70 |
"U/m2/s", iim, jjm+1, nhori, llm, 1, llm, nvert, & |
"U/m2/s", iim, jjm+1, nhori, llm, 1, llm, nvert, & |
71 |
"ave(X)", zsto, zout) |
"ave(X)", zsto, zout) |
72 |
endif |
endif |
73 |
|
|
74 |
!---Ajout Olivia |
!---Ajout Olivia |
75 |
CALL histdef(nid_tra, "d_tr_th_"//tnom(iq), & |
CALL histdef(nid_tra, "d_tr_th_"//tnom(iq), & |
76 |
"tendance thermique"// ttext(iiq), "?", & |
"tendance thermique"// ttext(iq), "?", & |
77 |
iim, jjm+1, nhori, llm, 1, llm, nvert, & |
iim, jjm+1, nhori, llm, 1, llm, nvert, & |
78 |
"ave(X)", zsto, zout) |
"ave(X)", zsto, zout) |
79 |
CALL histdef(nid_tra, "d_tr_cv_"//tnom(iq), & |
CALL histdef(nid_tra, "d_tr_cv_"//tnom(iq), & |
80 |
"tendance convection"// ttext(iiq), "?", & |
"tendance convection"// ttext(iq), "?", & |
81 |
iim, jjm+1, nhori, llm, 1, llm, nvert, & |
iim, jjm+1, nhori, llm, 1, llm, nvert, & |
82 |
"ave(X)", zsto, zout) |
"ave(X)", zsto, zout) |
83 |
CALL histdef(nid_tra, "d_tr_cl_"//tnom(iq), & |
CALL histdef(nid_tra, "d_tr_cl_"//tnom(iq), & |
84 |
"tendance couche limite"// ttext(iiq), "?", & |
"tendance couche limite"// ttext(iq), "?", & |
85 |
iim, jjm+1, nhori, llm, 1, llm, nvert, & |
iim, jjm+1, nhori, llm, 1, llm, nvert, & |
86 |
"ave(X)", zsto, zout) |
"ave(X)", zsto, zout) |
87 |
!---fin Olivia |
!---fin Olivia |