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 |
12 |
use dimens_m, only: iim, jjm, llm |
use dimens_m, only: iim, jjm, llm |
13 |
use dimphy, only: klon |
use dimphy, only: klon |
14 |
use temps, only: itau_phy, day_ref, annee_ref |
use disvert_m, only: presnivs |
15 |
use clesphys, only: ecrit_ins |
use dynetat0_m, only: day_ref, annee_ref |
|
use indicesol, only: nbsrf, clnsurf |
|
|
USE calendar, only: ymds2ju |
|
16 |
USE histbeg_totreg_m, ONLY : histbeg_totreg |
USE histbeg_totreg_m, ONLY : histbeg_totreg |
17 |
USE histdef_m, ONLY : histdef |
USE histdef_m, ONLY : histdef |
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 |
21 |
use phyetat0_m, only: rlon, rlat |
use phyetat0_m, only: rlon, rlat |
22 |
use comvert, only: presnivs |
use temps, only: itau_phy |
23 |
|
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 |
28 |
|
|
29 |
REAL zx_lon(iim, jjm + 1), zx_lat(iim, jjm + 1) |
REAL zx_lon(iim, jjm + 1), zx_lat(iim, jjm + 1) |
30 |
real zjulian, zsto, zout |
real zjulian, zsto, zout |
31 |
integer i, nhori, nvert, idayref, nsrf |
integer i, nhori, nvert, nsrf |
32 |
|
|
33 |
!------------------------------------------------------------------- |
!------------------------------------------------------------------- |
34 |
|
|
37 |
zsto = dtime * ecrit_ins |
zsto = dtime * ecrit_ins |
38 |
zout = dtime * ecrit_ins |
zout = dtime * ecrit_ins |
39 |
|
|
40 |
idayref = day_ref |
CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian) |
|
CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) |
|
41 |
|
|
42 |
CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlon, zx_lon) |
CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlon, zx_lon) |
43 |
DO i = 1, iim |
DO i = 1, iim |
49 |
jjm + 1, itau_phy, zjulian, dtime, nhori, nid_ins) |
jjm + 1, itau_phy, zjulian, dtime, nhori, nid_ins) |
50 |
write(*, *)'Inst ', itau_phy, zjulian |
write(*, *)'Inst ', itau_phy, zjulian |
51 |
CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", & |
CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", & |
52 |
llm, presnivs/100., nvert) |
presnivs/100., nvert) |
53 |
|
|
54 |
CALL histdef(nid_ins, "phis", "Surface geop. height", "-", & |
CALL histdef(nid_ins, "phis", "Surface geop. height", "-", & |
55 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
217 |
CALL histdef(nid_ins, "albs", "Surface albedo", "-", & |
CALL histdef(nid_ins, "albs", "Surface albedo", "-", & |
218 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
219 |
"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) |
|
220 |
|
|
|
!IM cf. AM 081204 BEG |
|
|
! HBTM2 |
|
221 |
CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m", & |
CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m", & |
222 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
223 |
"inst(X)", zsto, zout) |
"inst(X)", zsto, zout) |