libf/phylmd/ini_histday.f90

module ini_histday_m

  implicit none

contains

  subroutine ini_histday(dtime, ok_journe, nid_day, nq)

    ! From phylmd/ini_histday.h, v 1.3 2005/05/25 13:10:09

    use dimens_m, only: iim, jjm, llm
    use temps, only: itau_phy, day_ref, annee_ref
    USE calendar, only: ymds2ju
    USE histbeg_totreg_m, ONLY : histbeg_totreg
    USE histdef_m, ONLY : histdef
    USE histend_m, ONLY : histend
    USE histvert_m, ONLY : histvert
    use phyetat0_m, only: rlon, rlat
    use clesphys, only: ecrit_day
    use grid_change, only: gr_phy_write_2d
    use comvert, only: presnivs

    REAL, intent(in):: dtime ! pas temporel de la physique (s)
    logical, intent(in):: ok_journe
    integer, intent(out):: nid_day
    INTEGER, intent(in):: nq ! nombre de traceurs (y compris vapeur d'eau)

    ! Variables local to the procedure:
    REAL zx_lat(iim, jjm + 1)
    integer nhori, nvert
    real zjulian

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

    IF (ok_journe) THEN
       CALL ymds2ju(annee_ref, 1, day_ref, 0., zjulian)
       zx_lat = gr_phy_write_2d(rlat)
       CALL histbeg_totreg("histday", rlon(2: iim+1), zx_lat(1, :), 1, iim, &
            1, jjm + 1, itau_phy, zjulian, dtime, nhori, nid_day)
       CALL histvert(nid_day, "presnivs", "Vertical levels", "mb", &
            llm, presnivs/100., nvert)
       if (nq <= 4) then
          call histdef(nid_day, "Sigma_O3_Royer", &
               "column-density of ozone, in a cell, from Royer", "DU", &
               xsize=iim, ysize=jjm+1, horiid=nhori, pzsize=llm, par_oriz=1, &
               par_szz=llm, pzid=nvert, popp="ave(X)", pfreq_opp=dtime, &
               pfreq_wrt=real(ecrit_day))
       end if
       CALL histend(nid_day)
    ENDIF

  end subroutine ini_histday

end module ini_histday_m
1	module ini_histday_m
2
3	implicit none
4
5	contains
6
7	subroutine ini_histday(dtime, ok_journe, nid_day, nq)
8
9	! From phylmd/ini_histday.h, v 1.3 2005/05/25 13:10:09
10
11	use dimens_m, only: iim, jjm, llm
12	use temps, only: itau_phy, day_ref, annee_ref
13	USE calendar, only: ymds2ju
14	USE histbeg_totreg_m, ONLY : histbeg_totreg
15	USE histdef_m, ONLY : histdef
16	USE histend_m, ONLY : histend
17	USE histvert_m, ONLY : histvert
18	use phyetat0_m, only: rlon, rlat
19	use clesphys, only: ecrit_day
20	use grid_change, only: gr_phy_write_2d
21	use comvert, only: presnivs
22
23	REAL, intent(in):: dtime ! pas temporel de la physique (s)
24	logical, intent(in):: ok_journe
25	integer, intent(out):: nid_day
26	INTEGER, intent(in):: nq ! nombre de traceurs (y compris vapeur d'eau)
27
28	! Variables local to the procedure:
29	REAL zx_lat(iim, jjm + 1)
30	integer nhori, nvert
31	real zjulian
32
33	!--------------------------------
34
35	IF (ok_journe) THEN
36	CALL ymds2ju(annee_ref, 1, day_ref, 0., zjulian)
37	zx_lat = gr_phy_write_2d(rlat)
38	CALL histbeg_totreg("histday", rlon(2: iim+1), zx_lat(1, :), 1, iim, &
39	1, jjm + 1, itau_phy, zjulian, dtime, nhori, nid_day)
40	CALL histvert(nid_day, "presnivs", "Vertical levels", "mb", &
41	llm, presnivs/100., nvert)
42	if (nq <= 4) then
43	call histdef(nid_day, "Sigma_O3_Royer", &
44	"column-density of ozone, in a cell, from Royer", "DU", &
45	xsize=iim, ysize=jjm+1, horiid=nhori, pzsize=llm, par_oriz=1, &
46	par_szz=llm, pzid=nvert, popp="ave(X)", pfreq_opp=dtime, &
47	pfreq_wrt=real(ecrit_day))
48	end if
49	CALL histend(nid_day)
50	ENDIF
51
52	end subroutine ini_histday
53
54	end module ini_histday_m