6 |
|
|
7 |
contains |
contains |
8 |
|
|
9 |
subroutine ini_histins(dtime, ok_newmicro) |
subroutine ini_histins(ok_newmicro) |
10 |
|
|
11 |
! From phylmd/ini_histins.h, version 1.2, 2005/05/25 13:10:09 |
! From phylmd/ini_histins.h, version 1.2, 2005/05/25 13:10:09 |
12 |
|
|
13 |
use clesphys, only: ecrit_ins, ok_instan |
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 comconst, only: dtphys |
16 |
|
use dimensions, only: iim, jjm, llm, nqmx |
17 |
use disvert_m, only: presnivs |
use disvert_m, only: presnivs |
18 |
use dynetat0_m, only: day_ref, annee_ref, rlatu, rlonv |
use dynetat0_m, only: day_ref, annee_ref, rlatu, rlonv |
19 |
USE histbeg_totreg_m, ONLY : histbeg_totreg |
USE histbeg_totreg_m, ONLY : histbeg_totreg |
26 |
use phyetat0_m, only: itau_phy |
use phyetat0_m, only: itau_phy |
27 |
USE ymds2ju_m, only: ymds2ju |
USE ymds2ju_m, only: ymds2ju |
28 |
|
|
|
REAL, intent(in):: dtime ! pas temporel de la physique (s) |
|
29 |
logical, intent(in):: ok_newmicro |
logical, intent(in):: ok_newmicro |
30 |
|
|
31 |
! Local: |
! Local: |
37 |
print *, 'Call sequence information: ini_histins' |
print *, 'Call sequence information: ini_histins' |
38 |
|
|
39 |
test_ok_instan: IF (ok_instan) THEN |
test_ok_instan: IF (ok_instan) THEN |
40 |
zsto = dtime * ecrit_ins |
zsto = dtphys * ecrit_ins |
41 |
zout = dtime * ecrit_ins |
zout = dtphys * ecrit_ins |
42 |
CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian) |
CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian) |
43 |
CALL histbeg_totreg("histins", rlonv(:iim) / pi * 180., & |
CALL histbeg_totreg("histins", rlonv(:iim) / pi * 180., & |
44 |
rlatu / pi * 180., 1, iim, & |
rlatu / pi * 180., 1, iim, & |
45 |
1, jjm + 1, itau_phy, zjulian, dtime, nhori, nid_ins) |
1, jjm + 1, itau_phy, zjulian, dtphys, nhori, nid_ins) |
46 |
print *, 'itau_phy = ', itau_phy |
print *, 'itau_phy = ', itau_phy |
47 |
print *, "zjulian = ", zjulian |
print *, "zjulian = ", zjulian |
48 |
CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", & |
CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", & |
132 |
"inst(X)", zsto, zout) |
"inst(X)", zsto, zout) |
133 |
CALL histdef(nid_ins, "msnow", "surface snow amount", "kg/m2", & |
CALL histdef(nid_ins, "msnow", "surface snow amount", "kg/m2", & |
134 |
iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout) |
iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout) |
135 |
|
CALL histdef(nid_ins, "zxfqcalving", "ice calving", "kg m-2 s-1", & |
136 |
|
iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout) |
137 |
|
|
138 |
DO nsrf = 1, nbsrf |
DO nsrf = 1, nbsrf |
139 |
call histdef(nid_ins, "pourc_"//clnsurf(nsrf), & |
call histdef(nid_ins, "pourc_"//clnsurf(nsrf), & |
198 |
CALL histdef(nid_ins, "s_lcl", "Condensation level", "m", & |
CALL histdef(nid_ins, "s_lcl", "Condensation level", "m", & |
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) |
201 |
CALL histdef(nid_ins, "s_capCL", "Conv avlbl pot ener for ABL", "J/m2", & |
CALL histdef(nid_ins, "s_capCL", & |
202 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
"Convective available potential energy for atmospheric boundary " & |
203 |
|
// "layer", "J/m2", iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
204 |
"inst(X)", zsto, zout) |
"inst(X)", zsto, zout) |
205 |
CALL histdef(nid_ins, "s_oliqCL", "Liq Water in BL", "kg/m2", & |
CALL histdef(nid_ins, "s_oliqCL", "Liq Water in BL", "kg/m2", & |
206 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
211 |
CALL histdef(nid_ins, "s_therm", "Exces du thermique", "K", & |
CALL histdef(nid_ins, "s_therm", "Exces du thermique", "K", & |
212 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
213 |
"inst(X)", zsto, zout) |
"inst(X)", zsto, zout) |
|
CALL histdef(nid_ins, "s_trmb1", "deep_cape(HBTM2)", "J/m2", & |
|
|
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
|
|
"inst(X)", zsto, zout) |
|
|
CALL histdef(nid_ins, "s_trmb2", "inhibition (HBTM2)", "J/m2", & |
|
|
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
|
|
"inst(X)", zsto, zout) |
|
|
CALL histdef(nid_ins, "s_trmb3", "Point Omega (HBTM2)", "m", & |
|
|
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
|
|
"inst(X)", zsto, zout) |
|
214 |
CALL histdef(nid_ins, "qsurf", "Surface Air humidity", "", & |
CALL histdef(nid_ins, "qsurf", "Surface Air humidity", "", & |
215 |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
iim, (jjm + 1), nhori, 1, 1, 1, -99, & |
216 |
"inst(X)", zsto, zout) |
"inst(X)", zsto, zout) |