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	guez	34	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	guez	61	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	guez	34	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	guez	61	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	guez	34	end if
49			CALL histend(nid_day)
50			ENDIF
51
52			end subroutine ini_histday
53
54			end module ini_histday_m