--- trunk/libf/phylmd/ini_histins.f90 2012/09/20 13:00:41 66 +++ trunk/Sources/phylmd/ini_histins.f 2015/04/29 15:47:56 134 @@ -8,18 +8,19 @@ ! From phylmd/ini_histins.h, v 1.2 2005/05/25 13:10:09 + use clesphys, only: ecrit_ins use dimens_m, only: iim, jjm, llm use dimphy, only: klon - use temps, only: itau_phy, day_ref, annee_ref - use clesphys, only: ecrit_ins - use indicesol, only: nbsrf, clnsurf - USE calendar, only: ymds2ju + use disvert_m, only: presnivs + use dynetat0_m, only: day_ref, annee_ref USE histbeg_totreg_m, ONLY : histbeg_totreg USE histdef_m, ONLY : histdef USE histend_m, ONLY : histend USE histvert_m, ONLY : histvert + use indicesol, only: nbsrf, clnsurf use phyetat0_m, only: rlon, rlat - use disvert_m, only: presnivs + use temps, only: itau_phy + USE ymds2ju_m, only: ymds2ju REAL, intent(in):: dtime ! pas temporel de la physique (s) logical, intent(in):: ok_instan @@ -27,7 +28,7 @@ REAL zx_lon(iim, jjm + 1), zx_lat(iim, jjm + 1) real zjulian, zsto, zout - integer i, nhori, nvert, idayref, nsrf + integer i, nhori, nvert, nsrf !------------------------------------------------------------------- @@ -36,8 +37,7 @@ zsto = dtime * ecrit_ins zout = dtime * ecrit_ins - idayref = day_ref - CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) + CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian) CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlon, zx_lon) DO i = 1, iim @@ -49,7 +49,7 @@ jjm + 1, itau_phy, zjulian, dtime, nhori, nid_ins) write(*, *)'Inst ', itau_phy, zjulian CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", & - llm, presnivs/100., nvert) + presnivs/100., nvert) CALL histdef(nid_ins, "phis", "Surface geop. height", "-", & iim, (jjm + 1), nhori, 1, 1, 1, -99, &