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	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			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