4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
|
|
8 |
subroutine ini_histrac(nid_tra, pdtphys, nq_phys, lessivage) |
subroutine ini_histrac(nid_tra, pdtphys, nq_phys, lessivage) |
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 dimens_m, only: iim, jjm, llm |
use dimens_m, only: iim, jjm, llm |
13 |
USE calendar, only: ymds2ju |
USE calendar, only: ymds2ju |
14 |
use histcom, only: histbeg_totreg, histvert, histend, histdef |
USE histbeg_totreg_m, ONLY : histbeg_totreg |
15 |
|
USE histdef_m, ONLY : histdef |
16 |
|
USE histend_m, ONLY : histend |
17 |
|
USE histvert_m, ONLY : histvert |
18 |
use temps, only: annee_ref, day_ref, itau_phy |
use temps, only: annee_ref, day_ref, itau_phy |
19 |
use iniadvtrac_m, only: niadv, tnom, ttext |
use iniadvtrac_m, only: niadv, tnom, ttext |
20 |
use dimphy, only: klon |
use dimphy, only: klon |
21 |
use clesphys, only: ecrit_tra |
use clesphys, only: ecrit_tra |
22 |
use grid_change, only: gr_phy_write_2d |
use grid_change, only: gr_phy_write_2d |
23 |
use phyetat0_m, only: rlon, rlat |
use phyetat0_m, only: rlon, rlat |
24 |
use comvert, only: presnivs |
use disvert_m, only: presnivs |
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) |
45 |
zx_lat(:, :) = gr_phy_write_2d(rlat) |
zx_lat(:, :) = gr_phy_write_2d(rlat) |
46 |
CALL histbeg_totreg("histrac", rlon(2:iim+1), zx_lat(1, :), & |
CALL histbeg_totreg("histrac", rlon(2:iim+1), zx_lat(1, :), & |
47 |
1, iim, 1, jjm+1, itau_phy, zjulian, pdtphys, nhori, nid_tra) |
1, iim, 1, jjm+1, itau_phy, zjulian, pdtphys, nhori, nid_tra) |
48 |
CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb", llm, & |
CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb", presnivs, nvert) |
|
presnivs, nvert) |
|
49 |
|
|
50 |
zsto = pdtphys |
zsto = pdtphys |
51 |
zout = pdtphys * REAL(ecrit_tra) |
zout = pdtphys * REAL(ecrit_tra) |