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 clesphys2, only: conv_emanuel |
13 |
use dimens_m, only: iim, jjm, llm, nqmx |
use dimens_m, only: iim, jjm, llm, nqmx |
14 |
use dimphy, only: klon |
use dimphy, only: klon |
15 |
use disvert_m, only: presnivs |
use disvert_m, only: presnivs |
16 |
use dynetat0_m, only: day_ref, annee_ref |
use dynetat0_m, only: day_ref, annee_ref |
17 |
|
use gr_fi_ecrit_m, only: gr_fi_ecrit |
18 |
USE histbeg_totreg_m, ONLY : histbeg_totreg |
USE histbeg_totreg_m, ONLY : histbeg_totreg |
19 |
USE histdef_m, ONLY : histdef |
USE histdef_m, ONLY : histdef |
20 |
USE histend_m, ONLY : histend |
USE histend_m, ONLY : histend |
22 |
use indicesol, only: nbsrf, clnsurf |
use indicesol, only: nbsrf, clnsurf |
23 |
use iniadvtrac_m, only: tname, ttext |
use iniadvtrac_m, only: tname, ttext |
24 |
use phyetat0_m, only: rlon, rlat |
use phyetat0_m, only: rlon, rlat |
|
use temps, only: itau_phy |
|
25 |
USE ymds2ju_m, only: ymds2ju |
USE ymds2ju_m, only: ymds2ju |
26 |
|
|
27 |
REAL, intent(in):: dtime ! pas temporel de la physique (s) |
REAL, intent(in):: dtime ! pas temporel de la physique (s) |
28 |
logical, intent(in):: ok_instan |
logical, intent(in):: ok_instan |
29 |
integer, intent(out):: nid_ins |
integer, intent(out):: nid_ins |
30 |
|
integer, intent(in):: itau_phy |
31 |
|
|
32 |
|
! Local: |
33 |
REAL zx_lon(iim, jjm + 1), zx_lat(iim, jjm + 1) |
REAL zx_lon(iim, jjm + 1), zx_lat(iim, jjm + 1) |
34 |
real zjulian, zsto, zout |
real zjulian, zsto, zout |
35 |
integer i, nhori, nvert, nsrf, iq, it |
integer i, nhori, nvert, nsrf, iq, it |
262 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
263 |
"inst(X)", zsto, zout) |
"inst(X)", zsto, zout) |
264 |
|
|
265 |
!IM cf. AM 081204 END |
if (conv_emanuel) CALL histdef(nid_ins, "ptop", "cloud top pressure", & |
266 |
|
"Pa", iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout) |
267 |
|
|
268 |
! Champs 3D: |
! Champs 3D: |
269 |
|
|
299 |
"kg m-2", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", & |
"kg m-2", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", & |
300 |
zsto, zout) |
zsto, zout) |
301 |
|
|
302 |
|
CALL histdef(nid_ins, "rhum", "Relative humidity", & |
303 |
|
"", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", & |
304 |
|
zsto, zout) |
305 |
|
|
306 |
DO it = 1, nqmx - 2 |
DO it = 1, nqmx - 2 |
307 |
! champ 2D |
! champ 2D |
308 |
iq=it+2 |
iq=it+2 |