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, intent(in):: ts(:) ! (knon) 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, intent(in):: psol(:) ! (knon) 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, psol, ts, &
         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	341	REAL, intent(in):: ts(:) ! (knon) temperature de l'air a la surface
41	guez	304	REAL, intent(in):: qsurf(:) ! (knon) humidite relative a la surface
42			REAL, dimension(klon), intent(in):: rugos
43	guez	225	! rugos---input-R- rugosite
44	guez	341	REAL, intent(in):: psol(:) ! (knon) pression au sol
45	guez	246	REAL, intent(in):: ustar(:) ! (knon) facteur d'\'echelle pour le vent
46			REAL, dimension(klon), intent(in):: testar
47	guez	225	! testar--input-R- facteur d'echelle pour la temperature potentielle
48	guez	246	REAL, dimension(klon), intent(in):: qstar
49	guez	225	! qstar---input-R- facteur d'echelle pour l'humidite relative
50
51	guez	272	REAL, intent(out):: pref(:) ! (knon) pression au niveau de reference
52	guez	225	REAL, dimension(klon), intent(out):: delu
53			! delu----input-R- anomalie du vent par rapport au 1er niveau
54			REAL, dimension(klon), intent(out):: delte
55			! delte---input-R- anomalie de la temperature potentielle par
56			! rapport a la surface
57			REAL, dimension(klon), intent(out):: delq
58			! delq----input-R- anomalie de l'humidite relative par rapport a la surface
59
60			! Local:
61			INTEGER i
62	guez	274	REAL, dimension(knon):: cdram, cdrah, gref
63	guez	225
64			!-------------------------------------------------------------------------
65
66			DO i=1, knon
67			gref(i) = zref*RG
68			ENDDO
69
70			! Richardson at reference level
71
72	guez	341	CALL cdrag(nsrf, speed(:knon), temp(:knon), q_zref(:knon), gref, psol, ts, &
73			qsurf, rugos(:knon), cdram, cdrah, pref)
74	guez	225
75			DO i = 1, knon
76	guez	227	delu(i) = ustar(i) / sqrt(cdram(i))
77			delte(i) = testar(i) * sqrt(cdram(i)) / cdrah(i)
78			delq(i) = qstar(i) * sqrt(cdram(i)) / cdrah(i)
79	guez	225	ENDDO
80
81			END SUBROUTINE screenc
82
83			end module screenc_m