libf/phylmd/ini_histhf.f90

module ini_histhf_m

  implicit none

contains

  subroutine ini_histhf(dtime, nid_hf, nid_hf3d)

    ! From phylmd/ini_histhf.h, version 1.3 2005/05/25 13:10:09

    use dimens_m, only: iim, jjm, llm
    use temps, only: day_ref, annee_ref, itau_phy
    use dimphy, only: klon
    USE calendar, only: ymds2ju
    USE histbeg_totreg_m, ONLY : histbeg_totreg
    USE histend_m, ONLY : histend
    USE histvert_m, ONLY : histvert
    use phyetat0_m, only: rlon, rlat
    use comvert, only: presnivs
    use ini_histhf3d_m, only: ini_histhf3d

    REAL, intent(in):: dtime ! pas temporel de la physique (s)
    integer, intent(out):: nid_hf, nid_hf3d

    REAL zx_lon(iim, jjm + 1), zx_lat(iim, jjm + 1)
    integer idayref
    real zjulian
    integer i, nhori, nvert

    !-----------------------------------------------

    idayref = day_ref
    CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)

    CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlon, zx_lon)
    DO i = 1, iim
       zx_lon(i, 1) = rlon(i+1)
       zx_lon(i, (jjm + 1)) = rlon(i+1)
    ENDDO

    CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlat, zx_lat)
    CALL histbeg_totreg("histhf", zx_lon(:, 1), zx_lat(1, :), 1, iim, 1, &
         (jjm + 1), itau_phy, zjulian, dtime, nhori, nid_hf)

    CALL histvert(nid_hf, "presnivs", "Vertical levels", "mb", &
         llm, presnivs/100., nvert)

    call ini_histhf3d(dtime, nid_hf3d)
    CALL histend(nid_hf)

  end subroutine ini_histhf

end module ini_histhf_m
1	guez	34	module ini_histhf_m
2	guez	3
3	guez	34	implicit none
4	guez	3
5			contains
6
7	guez	20	subroutine ini_histhf(dtime, nid_hf, nid_hf3d)
8	guez	3
9			! From phylmd/ini_histhf.h, version 1.3 2005/05/25 13:10:09
10
11			use dimens_m, only: iim, jjm, llm
12			use temps, only: day_ref, annee_ref, itau_phy
13			use dimphy, only: klon
14	guez	30	USE calendar, only: ymds2ju
15	guez	61	USE histbeg_totreg_m, ONLY : histbeg_totreg
16			USE histend_m, ONLY : histend
17			USE histvert_m, ONLY : histvert
18	guez	3	use phyetat0_m, only: rlon, rlat
19	guez	20	use comvert, only: presnivs
20	guez	34	use ini_histhf3d_m, only: ini_histhf3d
21	guez	3
22			REAL, intent(in):: dtime ! pas temporel de la physique (s)
23			integer, intent(out):: nid_hf, nid_hf3d
24
25			REAL zx_lon(iim, jjm + 1), zx_lat(iim, jjm + 1)
26			integer idayref
27			real zjulian
28			integer i, nhori, nvert
29
30			!-----------------------------------------------
31
32			idayref = day_ref
33			CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
34
35			CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlon, zx_lon)
36			DO i = 1, iim
37			zx_lon(i, 1) = rlon(i+1)
38			zx_lon(i, (jjm + 1)) = rlon(i+1)
39			ENDDO
40
41			CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlat, zx_lat)
42	guez	15	CALL histbeg_totreg("histhf", zx_lon(:, 1), zx_lat(1, :), 1, iim, 1, &
43			(jjm + 1), itau_phy, zjulian, dtime, nhori, nid_hf)
44	guez	3
45			CALL histvert(nid_hf, "presnivs", "Vertical levels", "mb", &
46			llm, presnivs/100., nvert)
47
48	guez	20	call ini_histhf3d(dtime, nid_hf3d)
49	guez	3	CALL histend(nid_hf)
50
51			end subroutine ini_histhf
52
53	guez	34	end module ini_histhf_m