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 histcom, only: histbeg_totreg, histvert, histend, histdef
    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", &
               pxsize=iim, pysize=jjm+1, phoriid=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 histcom, only: histbeg_totreg, histvert, histend, histdef
15	use phyetat0_m, only: rlon, rlat
16	use clesphys, only: ecrit_day
17	use grid_change, only: gr_phy_write_2d
18	use comvert, only: presnivs
19
20	REAL, intent(in):: dtime ! pas temporel de la physique (s)
21	logical, intent(in):: ok_journe
22	integer, intent(out):: nid_day
23	INTEGER, intent(in):: nq ! nombre de traceurs (y compris vapeur d'eau)
24
25	! Variables local to the procedure:
26	REAL zx_lat(iim, jjm + 1)
27	integer nhori, nvert
28	real zjulian
29
30	!--------------------------------
31
32	IF (ok_journe) THEN
33	CALL ymds2ju(annee_ref, 1, day_ref, 0., zjulian)
34	zx_lat = gr_phy_write_2d(rlat)
35	CALL histbeg_totreg("histday", rlon(2: iim+1), zx_lat(1, :), 1, iim, &
36	1, jjm + 1, itau_phy, zjulian, dtime, nhori, nid_day)
37	CALL histvert(nid_day, "presnivs", "Vertical levels", "mb", &
38	llm, presnivs/100., nvert)
39	if (nq <= 4) then
40	call histdef(nid_day, "Sigma_O3_Royer", &
41	"column-density of ozone, in a cell, from Royer", "DU", &
42	pxsize=iim, pysize=jjm+1, phoriid=nhori, pzsize=llm, &
43	par_oriz=1, par_szz=llm, pzid=nvert, popp="ave(X)", &
44	pfreq_opp=dtime, pfreq_wrt=real(ecrit_day))
45	end if
46	CALL histend(nid_day)
47	ENDIF
48
49	end subroutine ini_histday
50
51	end module ini_histday_m