phylmd/Interface_surf/screenc.f90

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 cdrag_m, only: cdrag
    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
    ! speed---input-R- module du vent au 1er niveau du modele
    REAL, dimension(klon), intent(in):: temp
    ! temp----input-R- temperature de l'air au 1er niveau du modele
    REAL, dimension(klon), intent(in):: q_zref
    ! 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
    ! ts------input-R- temperature de l'air a la surface
    REAL, intent(in):: qsurf(:) ! (knon) humidite relative a la surface
    REAL, dimension(klon), intent(in):: rugos
    ! rugos---input-R- rugosite
    REAL, dimension(klon), intent(in):: psol
    ! 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 cdrag(nsrf, speed(:knon), temp(:knon), q_zref(:knon), gref(:knon), &
         psol(:knon), ts(:knon), qsurf, 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	guez	225	module screenc_m
2
3			IMPLICIT NONE
4
5			contains
6
7	guez	227	SUBROUTINE screenc(klon, knon, nsrf, speed, temp, q_zref, zref, ts, &
8			qsurf, rugos, psol, ustar, testar, qstar, pref, delu, delte, delq)
9	guez	225
10	guez	227	! From LMDZ4/libf/phylmd/screenc.F90, version 1.1.1.1, 2004/05/19 12:53:09
11	guez	225
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	guez	275	use cdrag_m, only: cdrag
23	guez	227	use SUPHEC_M, only: RG
24
25	guez	225	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	guez	304	REAL, dimension(klon), intent(in):: speed
33	guez	225	! speed---input-R- module du vent au 1er niveau du modele
34	guez	304	REAL, dimension(klon), intent(in):: temp
35	guez	225	! temp----input-R- temperature de l'air au 1er niveau du modele
36	guez	304	REAL, dimension(klon), intent(in):: q_zref
37	guez	225	! q_zref--input-R- humidite relative au 1er niveau du modele
38			REAL, intent(in):: zref
39			! zref----input-R- altitude de reference
40	guez	304	REAL, dimension(klon), intent(in):: ts
41	guez	225	! ts------input-R- temperature de l'air a la surface
42	guez	304	REAL, intent(in):: qsurf(:) ! (knon) humidite relative a la surface
43			REAL, dimension(klon), intent(in):: rugos
44	guez	225	! rugos---input-R- rugosite
45	guez	304	REAL, dimension(klon), intent(in):: psol
46	guez	225	! psol----input-R- pression au sol
47	guez	246	REAL, intent(in):: ustar(:) ! (knon) facteur d'\'echelle pour le vent
48			REAL, dimension(klon), intent(in):: testar
49	guez	225	! testar--input-R- facteur d'echelle pour la temperature potentielle
50	guez	246	REAL, dimension(klon), intent(in):: qstar
51	guez	225	! qstar---input-R- facteur d'echelle pour l'humidite relative
52
53	guez	272	REAL, intent(out):: pref(:) ! (knon) pression au niveau de reference
54	guez	225	REAL, dimension(klon), intent(out):: delu
55			! delu----input-R- anomalie du vent par rapport au 1er niveau
56			REAL, dimension(klon), intent(out):: delte
57			! delte---input-R- anomalie de la temperature potentielle par
58			! rapport a la surface
59			REAL, dimension(klon), intent(out):: delq
60			! delq----input-R- anomalie de l'humidite relative par rapport a la surface
61
62			! Local:
63			INTEGER i
64	guez	274	REAL, dimension(knon):: cdram, cdrah, gref
65	guez	225
66			!-------------------------------------------------------------------------
67
68			DO i=1, knon
69			gref(i) = zref*RG
70			ENDDO
71
72			! Richardson at reference level
73
74	guez	275	CALL cdrag(nsrf, speed(:knon), temp(:knon), q_zref(:knon), gref(:knon), &
75	guez	304	psol(:knon), ts(:knon), qsurf, rugos(:knon), cdram, cdrah, pref)
76	guez	225
77			DO i = 1, knon
78	guez	227	delu(i) = ustar(i) / sqrt(cdram(i))
79			delte(i) = testar(i) * sqrt(cdram(i)) / cdrah(i)
80			delq(i) = qstar(i) * sqrt(cdram(i)) / cdrah(i)
81	guez	225	ENDDO
82
83			END SUBROUTINE screenc
84
85			end module screenc_m