2 |
|
|
3 |
implicit none |
implicit none |
4 |
|
|
5 |
|
integer, save:: nid_ins |
6 |
|
|
7 |
contains |
contains |
8 |
|
|
9 |
subroutine ini_histins(dtime, ok_instan, nid_ins, itau_phy) |
subroutine ini_histins(dtime) |
10 |
|
|
11 |
! 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 |
12 |
|
|
13 |
use clesphys, only: ecrit_ins |
use clesphys, only: ecrit_ins, ok_instan |
14 |
use clesphys2, only: conv_emanuel |
use clesphys2, only: conv_emanuel |
15 |
use dimens_m, only: iim, jjm, llm, nqmx |
use dimens_m, only: iim, jjm, llm, nqmx |
|
use dimphy, only: klon |
|
16 |
use disvert_m, only: presnivs |
use disvert_m, only: presnivs |
17 |
use dynetat0_m, only: day_ref, annee_ref, rlatu, rlonv |
use dynetat0_m, only: day_ref, annee_ref, rlatu, rlonv |
18 |
USE histbeg_totreg_m, ONLY : histbeg_totreg |
USE histbeg_totreg_m, ONLY : histbeg_totreg |
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 nr_util, only: pi |
use nr_util, only: pi |
25 |
|
use phyetat0_m, only: itau_phy |
26 |
USE ymds2ju_m, only: ymds2ju |
USE ymds2ju_m, only: ymds2ju |
27 |
|
|
28 |
REAL, intent(in):: dtime ! pas temporel de la physique (s) |
REAL, intent(in):: dtime ! pas temporel de la physique (s) |
|
logical, intent(in):: ok_instan |
|
|
integer, intent(out):: nid_ins |
|
|
integer, intent(in):: itau_phy |
|
29 |
|
|
30 |
! Local: |
! Local: |
31 |
real zjulian, zsto, zout |
real zjulian, zsto, zout |
154 |
"inst(X)", zsto, zout) |
"inst(X)", zsto, zout) |
155 |
|
|
156 |
DO nsrf = 1, nbsrf |
DO nsrf = 1, nbsrf |
|
|
|
157 |
call histdef(nid_ins, "pourc_"//clnsurf(nsrf), & |
call histdef(nid_ins, "pourc_"//clnsurf(nsrf), & |
158 |
"% "//clnsurf(nsrf), "%", & |
"% "//clnsurf(nsrf), "%", & |
159 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
198 |
"rugosite "//clnsurf(nsrf), "-", & |
"rugosite "//clnsurf(nsrf), "-", & |
199 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
200 |
"inst(X)", zsto, zout) |
"inst(X)", zsto, zout) |
|
!XXX |
|
201 |
END DO |
END DO |
202 |
|
|
203 |
CALL histdef(nid_ins, "rugs", "rugosity", "-", & |
CALL histdef(nid_ins, "rugs", "rugosity", "-", & |
204 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
205 |
"inst(X)", zsto, zout) |
"inst(X)", zsto, zout) |
249 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
250 |
"inst(X)", zsto, zout) |
"inst(X)", zsto, zout) |
251 |
|
|
252 |
if (conv_emanuel) CALL histdef(nid_ins, "ptop", "cloud top pressure", & |
if (conv_emanuel) then |
253 |
"Pa", iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout) |
CALL histdef(nid_ins, "ptop", "cloud top pressure", & |
254 |
|
"Pa", iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout) |
255 |
|
CALL histdef(nid_ins, "dnwd0", "unsaturated downdraft", & |
256 |
|
"kg/m2/s", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", & |
257 |
|
zsto, zout) |
258 |
|
end if |
259 |
|
|
260 |
! Champs 3D: |
! Champs 3D: |
261 |
|
|