dyn3d/Guide/conf_guide.f

module conf_guide_m

  IMPLICIT NONE

  ! Constantes de rappel, en jours :
  REAL:: tau_min_u = 0.03
  REAL:: tau_max_u = 10.
  REAL:: tau_min_v = 0.03
  REAL:: tau_max_v = 10.
  REAL:: tau_min_t = 0.03
  REAL:: tau_max_t = 10.
  REAL:: tau_min_q = 0.03
  REAL:: tau_max_q = 10.
  REAL:: tau_min_p = 0.03
  REAL:: tau_max_p = 10.

  LOGICAL:: ncep = .false. ! Coordonnee vert NCEP ou ECMWF
  LOGICAL:: ini_anal = .false. ! Initial = analyse
  LOGICAL:: guide_u = .false. ! guidage de u
  LOGICAL:: guide_v = .false. ! gvidage de v
  LOGICAL:: guide_t = .false. ! guidage de T
  LOGICAL:: guide_q = .false. ! guidage de q

  logical:: online = .true. ! controle du guide
  ! hors-ligne: x=x_rea

  ! Dans le cas où on n'a les analyses que sur une bande de latitudes :
  REAL, save:: lat_min_guide ! minimum latitude for nudging, in rad
  real, save:: lat_max_guide ! maximum latitude for nudging, in rad

  logical, save:: ok_guide ! guidage
  REAL, save:: factt ! pas de temps entre deux appels au guidage, en jours

contains

  SUBROUTINE conf_guide

    ! From LMDZ4/libf/dyn3d/conf_guide.F, version 1.1.1.1 2004/05/19 12:53:07
    !  Parametres de controle du run:

    use abort_gcm_m, only: abort_gcm
    use comconst, only: daysec, dtvr
    use conf_gcm_m, only: day_step, iperiod
    use dynetat0_m, only: grossismx, grossismy
    use nr_util, only: assert, pi
    use unit_nml_m, only: unit_nml

    ! Local:

    REAL:: lat_min_guide_deg = -90. ! in degrees
    real:: lat_max_guide_deg = 90. ! in degrees

    namelist /conf_guide_nml/ ncep, ini_anal, guide_u, guide_v, guide_t, &
         guide_q, online, tau_min_u, tau_max_u, tau_min_v, tau_max_v, &
         tau_min_t, tau_max_t, tau_min_q, tau_max_q, tau_min_p, tau_max_p, &
         lat_min_guide_deg, lat_max_guide_deg

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

    print *, "Call sequence information: conf_guide"

    print *, "Enter namelist 'conf_guide_nml'."
    read(unit=*, nml=conf_guide_nml)
    write(unit_nml, nml=conf_guide_nml)

    lat_min_guide = lat_min_guide_deg / 180. * pi
    lat_max_guide = lat_max_guide_deg / 180. * pi

    ok_guide = any((/guide_u, guide_v, guide_t, guide_q/))
    if (ok_guide .and. mod(day_step, 4 * iperiod) /= 0) call &
         abort_gcm("conf_guide", 'ok_guide day_step iperiod')

    if (ok_guide .and. online) then
       factt = dtvr * iperiod / daysec
       print *, "factt = ", factt
       if (abs(grossismx - 1.) >= 0.1 .and. abs(grossismy - 1.) >= 0.1) then
          if (guide_u) call assert(factt / tau_min_u < 1, &
               "conf_guide tau_min_u")
          if (guide_v) call assert(factt / tau_min_v < 1, &
               "conf_guide tau_min_v")
          if (guide_t) call assert(factt / tau_min_t < 1, &
               "conf_guide tau_min_t")
          if (guide_q) call assert(factt / tau_min_q < 1, &
               "conf_guide tau_min_q")
       end if
    end if

  end SUBROUTINE conf_guide

end module conf_guide_m
1	guez	44	module conf_guide_m
2	guez	37
3	guez	44	IMPLICIT NONE
4	guez	3
5	guez	115	! Constantes de rappel, en jours :
6			REAL:: tau_min_u = 0.03
7	guez	83	REAL:: tau_max_u = 10.
8	guez	115	REAL:: tau_min_v = 0.03
9	guez	83	REAL:: tau_max_v = 10.
10	guez	115	REAL:: tau_min_t = 0.03
11	guez	83	REAL:: tau_max_t = 10.
12	guez	115	REAL:: tau_min_q = 0.03
13	guez	83	REAL:: tau_max_q = 10.
14	guez	115	REAL:: tau_min_p = 0.03
15	guez	83	REAL:: tau_max_p = 10.
16	guez	3
17	guez	83	LOGICAL:: ncep = .false. ! Coordonnee vert NCEP ou ECMWF
18			LOGICAL:: ini_anal = .false. ! Initial = analyse
19	guez	115	LOGICAL:: guide_u = .false. ! guidage de u
20			LOGICAL:: guide_v = .false. ! gvidage de v
21	guez	107	LOGICAL:: guide_t = .false. ! guidage de T
22			LOGICAL:: guide_q = .false. ! guidage de q
23	guez	83
24	guez	102	logical:: online = .true. ! controle du guide
25			! hors-ligne: x=x_rea
26
27	guez	140	! Dans le cas où on n'a les analyses que sur une bande de latitudes :
28			REAL, save:: lat_min_guide ! minimum latitude for nudging, in rad
29			real, save:: lat_max_guide ! maximum latitude for nudging, in rad
30	guez	83
31	guez	115	logical, save:: ok_guide ! guidage
32			REAL, save:: factt ! pas de temps entre deux appels au guidage, en jours
33
34	guez	44	contains
35	guez	3
36	guez	44	SUBROUTINE conf_guide
37	guez	3
38	guez	44	! From LMDZ4/libf/dyn3d/conf_guide.F, version 1.1.1.1 2004/05/19 12:53:07
39			! Parametres de controle du run:
40	guez	3
41	guez	115	use abort_gcm_m, only: abort_gcm
42			use comconst, only: daysec, dtvr
43			use conf_gcm_m, only: day_step, iperiod
44	guez	139	use dynetat0_m, only: grossismx, grossismy
45	guez	140	use nr_util, only: assert, pi
46	guez	83	use unit_nml_m, only: unit_nml
47	guez	3
48	guez	140	! Local:
49
50			REAL:: lat_min_guide_deg = -90. ! in degrees
51			real:: lat_max_guide_deg = 90. ! in degrees
52
53	guez	83	namelist /conf_guide_nml/ ncep, ini_anal, guide_u, guide_v, guide_t, &
54	guez	115	guide_q, online, tau_min_u, tau_max_u, tau_min_v, tau_max_v, &
55			tau_min_t, tau_max_t, tau_min_q, tau_max_q, tau_min_p, tau_max_p, &
56	guez	140	lat_min_guide_deg, lat_max_guide_deg
57	guez	83
58	guez	44	!-----------------------------------------------------------------------
59	guez	3
60	guez	44	print *, "Call sequence information: conf_guide"
61	guez	3
62	guez	83	print *, "Enter namelist 'conf_guide_nml'."
63			read(unit=*, nml=conf_guide_nml)
64			write(unit_nml, nml=conf_guide_nml)
65	guez	3
66	guez	140	lat_min_guide = lat_min_guide_deg / 180. * pi
67			lat_max_guide = lat_max_guide_deg / 180. * pi
68
69	guez	115	ok_guide = any((/guide_u, guide_v, guide_t, guide_q/))
70			if (ok_guide .and. mod(day_step, 4 * iperiod) /= 0) call &
71	guez	171	abort_gcm("conf_guide", 'ok_guide day_step iperiod')
72	guez	115
73			if (ok_guide .and. online) then
74			factt = dtvr * iperiod / daysec
75			print *, "factt = ", factt
76			if (abs(grossismx - 1.) >= 0.1 .and. abs(grossismy - 1.) >= 0.1) then
77			if (guide_u) call assert(factt / tau_min_u < 1, &
78			"conf_guide tau_min_u")
79			if (guide_v) call assert(factt / tau_min_v < 1, &
80			"conf_guide tau_min_v")
81			if (guide_t) call assert(factt / tau_min_t < 1, &
82			"conf_guide tau_min_t")
83			if (guide_q) call assert(factt / tau_min_q < 1, &
84			"conf_guide tau_min_q")
85			end if
86			end if
87
88	guez	44	end SUBROUTINE conf_guide
89
90			end module conf_guide_m