trunk/phylmd/screenc.f

module screenc_m

  IMPLICIT NONE

contains

  SUBROUTINE screenc(klon, knon, nsrf, speed, temp, q_zref, zref, ts, &
       qsurf, rugos, psol, ustar, testar, qstar, pref, delu, delte, delq)
    
    ! From LMDZ4/libf/phylmd/screenc.F90, version 1.1.1.1, 2004/05/19 12:53:09

    ! Objet : calcul "correcteur" des anomalies du vent, de la
    ! temp\'erature potentielle et de l'humidit\'e relative au niveau
    ! de r\'ef\'erence zref et par rapport au 1er niveau (pour u) ou
    ! \`a la surface (pour theta et q) \`a partir des equations de
    ! Louis.

    ! Reference: Hess, Colman et McAvaney (1995)

    ! I. Musat, 01.07.2002

    use coefcdrag_m, only: coefcdrag
    use SUPHEC_M, only: RG

    INTEGER, intent(in):: klon
    ! klon----input-I- dimension de la grille physique (=
    ! nb_pts_latitude X nb_pts_longitude)
    INTEGER, intent(in):: knon
    ! knon----input-I- nombre de points pour un type de surface
    INTEGER, intent(in):: nsrf
    ! nsrf----input-I- indice pour le type de surface; voir indicesol.inc
    REAL, dimension(klon), intent(in):: speed, temp, q_zref
    ! speed---input-R- module du vent au 1er niveau du modele
    ! temp----input-R- temperature de l'air au 1er niveau du modele
    ! q_zref--input-R- humidite relative au 1er niveau du modele
    REAL, intent(in):: zref
    ! zref----input-R- altitude de reference
    REAL, dimension(klon), intent(in):: ts, qsurf, rugos, psol
    ! ts------input-R- temperature de l'air a la surface
    ! qsurf---input-R- humidite relative a la surface
    ! rugos---input-R- rugosite
    ! psol----input-R- pression au sol
    REAL, intent(in):: ustar(:) ! (knon) facteur d'\'echelle pour le vent
    REAL, dimension(klon), intent(in):: testar
    ! testar--input-R- facteur d'echelle pour la temperature potentielle
    REAL, dimension(klon), intent(in):: qstar
    ! qstar---input-R- facteur d'echelle pour l'humidite relative

    REAL, intent(out):: pref(:) ! (knon) pression au niveau de reference
    REAL, dimension(klon), intent(out):: delu
    ! delu----input-R- anomalie du vent par rapport au 1er niveau
    REAL, dimension(klon), intent(out):: delte
    ! delte---input-R- anomalie de la temperature potentielle par
    ! rapport a la surface
    REAL, dimension(klon), intent(out):: delq
    ! delq----input-R- anomalie de l'humidite relative par rapport a la surface

    ! Local:
    INTEGER i
    REAL, dimension(knon):: cdram, cdrah, gref

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

    DO i=1, knon
       gref(i) = zref*RG
    ENDDO

    ! Richardson at reference level

    CALL coefcdrag(nsrf, speed(:knon), temp(:knon), q_zref(:knon), &
         gref(:knon), psol(:knon), ts(:knon), qsurf(:knon), rugos(:knon), &
         cdram, cdrah, pref)

    DO i = 1, knon
       delu(i) = ustar(i) / sqrt(cdram(i))
       delte(i) = testar(i) * sqrt(cdram(i)) / cdrah(i)
       delq(i) = qstar(i) * sqrt(cdram(i)) / cdrah(i)
    ENDDO

  END SUBROUTINE screenc

end module screenc_m
1	module screenc_m
2
3	IMPLICIT NONE
4
5	contains
6
7	SUBROUTINE screenc(klon, knon, nsrf, speed, temp, q_zref, zref, ts, &
8	qsurf, rugos, psol, ustar, testar, qstar, pref, delu, delte, delq)
9
10	! From LMDZ4/libf/phylmd/screenc.F90, version 1.1.1.1, 2004/05/19 12:53:09
11
12	! Objet : calcul "correcteur" des anomalies du vent, de la
13	! temp\'erature potentielle et de l'humidit\'e relative au niveau
14	! de r\'ef\'erence zref et par rapport au 1er niveau (pour u) ou
15	! \`a la surface (pour theta et q) \`a partir des equations de
16	! Louis.
17
18	! Reference: Hess, Colman et McAvaney (1995)
19
20	! I. Musat, 01.07.2002
21
22	use coefcdrag_m, only: coefcdrag
23	use SUPHEC_M, only: RG
24
25	INTEGER, intent(in):: klon
26	! klon----input-I- dimension de la grille physique (=
27	! nb_pts_latitude X nb_pts_longitude)
28	INTEGER, intent(in):: knon
29	! knon----input-I- nombre de points pour un type de surface
30	INTEGER, intent(in):: nsrf
31	! nsrf----input-I- indice pour le type de surface; voir indicesol.inc
32	REAL, dimension(klon), intent(in):: speed, temp, q_zref
33	! speed---input-R- module du vent au 1er niveau du modele
34	! temp----input-R- temperature de l'air au 1er niveau du modele
35	! q_zref--input-R- humidite relative au 1er niveau du modele
36	REAL, intent(in):: zref
37	! zref----input-R- altitude de reference
38	REAL, dimension(klon), intent(in):: ts, qsurf, rugos, psol
39	! ts------input-R- temperature de l'air a la surface
40	! qsurf---input-R- humidite relative a la surface
41	! rugos---input-R- rugosite
42	! psol----input-R- pression au sol
43	REAL, intent(in):: ustar(:) ! (knon) facteur d'\'echelle pour le vent
44	REAL, dimension(klon), intent(in):: testar
45	! testar--input-R- facteur d'echelle pour la temperature potentielle
46	REAL, dimension(klon), intent(in):: qstar
47	! qstar---input-R- facteur d'echelle pour l'humidite relative
48
49	REAL, intent(out):: pref(:) ! (knon) pression au niveau de reference
50	REAL, dimension(klon), intent(out):: delu
51	! delu----input-R- anomalie du vent par rapport au 1er niveau
52	REAL, dimension(klon), intent(out):: delte
53	! delte---input-R- anomalie de la temperature potentielle par
54	! rapport a la surface
55	REAL, dimension(klon), intent(out):: delq
56	! delq----input-R- anomalie de l'humidite relative par rapport a la surface
57
58	! Local:
59	INTEGER i
60	REAL, dimension(knon):: cdram, cdrah, gref
61
62	!-------------------------------------------------------------------------
63
64	DO i=1, knon
65	gref(i) = zref*RG
66	ENDDO
67
68	! Richardson at reference level
69
70	CALL coefcdrag(nsrf, speed(:knon), temp(:knon), q_zref(:knon), &
71	gref(:knon), psol(:knon), ts(:knon), qsurf(:knon), rugos(:knon), &
72	cdram, cdrah, pref)
73
74	DO i = 1, knon
75	delu(i) = ustar(i) / sqrt(cdram(i))
76	delte(i) = testar(i) * sqrt(cdram(i)) / cdrah(i)
77	delq(i) = qstar(i) * sqrt(cdram(i)) / cdrah(i)
78	ENDDO
79
80	END SUBROUTINE screenc
81
82	end module screenc_m