source: LMDZ6/trunk/libf/phylmd/screenc_mod.f90 @ 5309

Last change on this file since 5309 was 5301, checked in by abarral, 9 days ago

Turn tsoilnudge.h fcg_gcssold.h flux_arp.h into module
  • Property svn:executable set to *
File size: 8.4 KB
RevLine 
[3817]1!
2MODULE screenc_mod
3!
4! This module contains some procedures for calculation of the correction
5! of temperature, specific humidity and wind at a reference level
6!
7  USE cdrag_mod
8  IMPLICIT NONE
9
10CONTAINS
11!
12!****************************************************************************************
13!
14!r original routine svn3623
15!
16      SUBROUTINE screenc(klon, knon, nsrf, zxli, &
17                         speed, temp, q_zref, zref, &
18                         ts, qsurf, z0m, z0h, psol, &
19                         ustar, testar, qstar, okri, ri1, &
[4722]20                         pref, delu, delte, delq, s_pblh, prain, tsol, pat1)
[5285]21      USE yomcst_mod_h
[5301]22      USE flux_arp_mod_h
[5274]23IMPLICIT NONE
[3817]24!-----------------------------------------------------------------------
[5274]25!
26! Objet : calcul "correcteur" des anomalies du vent, de la temperature
27!         potentielle et de l'humidite relative au niveau de reference zref et
28!         par rapport au 1er niveau (pour u) ou a la surface (pour theta et q)
[3817]29!         a partir des equations de Louis.
30!
31! Reference : Hess, Colman et McAvaney (1995)
32!
33! I. Musat, 01.07.2002
34!-----------------------------------------------------------------------
35!
36! klon----input-I- dimension de la grille physique (= nb_pts_latitude X nb_pts_longitude)
37! knon----input-I- nombre de points pour un type de surface
38! nsrf----input-I- indice pour le type de surface; voir indice_sol_mod.F90
39! zxli----input-L- TRUE si calcul des cdrags selon Laurent Li
40! speed---input-R- module du vent au 1er niveau du modele
41! temp----input-R- temperature de l'air au 1er niveau du modele
42! q_zref--input-R- humidite relative au 1er niveau du modele
43! zref----input-R- altitude de reference
44! ts------input-R- temperature de l'air a la surface
45! qsurf---input-R- humidite relative a la surface
46! z0m, z0h---input-R- rugosite
47! psol----input-R- pression au sol
48! ustar---input-R- facteur d'echelle pour le vent
49! testar--input-R- facteur d'echelle pour la temperature potentielle
50! qstar---input-R- facteur d'echelle pour l'humidite relative
[5274]51! okri----input-L- TRUE si on veut tester le nb. Richardson entre la sfce
[3817]52!                  et zref par rapport au Ri entre la sfce et la 1ere couche
[5274]53! ri1-----input-R- nb. Richardson entre la surface et la 1ere couche
[3817]54!
55! pref----input-R- pression au niveau de reference
56! delu----input-R- anomalie du vent par rapport au 1er niveau
57! delte---input-R- anomalie de la temperature potentielle par rapport a la surface
58! delq----input-R- anomalie de l'humidite relative par rapport a la surface
59!
60      INTEGER, intent(in) :: klon, knon, nsrf
[5274]61      LOGICAL, intent(in) :: zxli, okri
[3817]62      REAL, dimension(klon), intent(in) :: speed, temp, q_zref
63      REAL, intent(in) :: zref
[4722]64      REAL, dimension(klon), intent(IN) :: ts
65      REAL, dimension(klon), intent(in) :: qsurf, psol
66      REAL, dimension(klon), intent(inout):: z0m, z0h
[5274]67      REAL, dimension(klon), intent(in) :: ustar, testar, qstar, ri1
[4722]68
[5274]69      REAL, dimension(klon), intent(inout) :: s_pblh
70      REAL, dimension(klon), intent(in) :: prain
71      REAL, dimension(klon), intent(in) :: tsol
72      REAL, DIMENSION(klon), INTENT(IN)    :: pat1 !pression premier lev
[3817]73!
[5274]74      REAL, dimension(klon), intent(out) :: pref, delu, delte, delq
[3817]75!-----------------------------------------------------------------------
76!
77! Variables locales 
78      INTEGER :: i
[4478]79      REAL, dimension(klon) :: cdram, cdrah, cdran, zri1, gref,ycdragm,zri_zero
[3817]80!
81!-------------------------------------------------------------------------
82      DO i=1, knon
83        gref(i) = zref*RG
84      ENDDO
85!
86! Richardson at reference level
87!
88!      CALL coefcdrag (klon, knon, nsrf, zxli, &
89!                    speed, temp, q_zref, gref, &
90!                    psol, ts, qsurf, rugos, &
91!                    okri, ri1, &
92!                    cdram, cdrah, cdran, zri1, &
93!                    pref)
94! Fuxing WANG, 04/03/2015, replace the coefcdrag by the merged version: cdrag
95      CALL cdrag (knon, nsrf, &
96                    speed, temp, q_zref, gref, &
[4722]97                    psol, s_pblh, ts, qsurf, z0m, z0h, &
[4478]98                    zri_zero,0,                &
[4722]99                    cdram, cdrah, zri1, pref, prain, tsol, pat1)
[3817]100      DO i = 1, knon
101        IF(ok_prescr_ust) THEN
102! La aussi il faut forcer avec ust (FC + MPL 20160210)
103        ycdragm(i) = ust*ust/(1.+speed(i))/speed(i)
104        cdram=ycdragm
105        delu(i) = ust/sqrt(cdram(i))
106        ELSE
107        delu(i) = ustar(i)/sqrt(cdram(i))
108        ENDIF
109        delte(i)= (testar(i)* sqrt(cdram(i)))/ &
110                   cdrah(i)
111        delq(i)= (qstar(i)* sqrt(cdram(i)))/ &
112                  cdrah(i)
113      ENDDO
114!
115      RETURN
116      END SUBROUTINE screenc
117!
118      SUBROUTINE screencn(klon, knon, nsrf, zxli, &
119                         speed, temp, q_zref, zref, &
120                         ts, qsurf, z0m, z0h, psol, &
121                         cdrm, cdrh,  okri, &
122                         ri1, iri1, &
[4722]123                         pref, delm, delh, zri1, s_pblh, prain, tsol, pat1)
[5285]124      USE yomcst_mod_h
[5301]125      USE flux_arp_mod_h
[5274]126IMPLICIT NONE
[3817]127!-----------------------------------------------------------------------
[5274]128!
129! Objet : calcul "correcteur" des anomalies du vent, de la temperature
130!         potentielle et de l'humidite relative au niveau de reference zref et
131!         par rapport au 1er niveau (pour u) ou a la surface (pour theta et q)
[3817]132!         a partir des equations de Louis.
133!
134! Reference : Hess, Colman et McAvaney (1995)
135!
136! I. Musat, 01.07.2002
137!-----------------------------------------------------------------------
138!
139! klon----input-I- dimension de la grille physique (= nb_pts_latitude X nb_pts_longitude)
140! knon----input-I- nombre de points pour un type de surface
141! nsrf----input-I- indice pour le type de surface; voir indice_sol_mod.F90
142! zxli----input-L- TRUE si calcul des cdrags selon Laurent Li
143! speed---input-R- module du vent au 1er niveau du modele
144! temp----input-R- temperature de l'air au 1er niveau du modele
145! q_zref--input-R- humidite relative au 1er niveau du modele
146! zref----input-R- altitude de reference
147! ts------input-R- temperature de l'air a la surface
148! qsurf---input-R- humidite relative a la surface
149! z0m, z0h---input-R- rugosite
150! psol----input-R- pression au sol
151! ustar---input-R- facteur d'echelle pour le vent
152! testar--input-R- facteur d'echelle pour la temperature potentielle
153! qstar---input-R- facteur d'echelle pour l'humidite relative
[5274]154! okri----input-L- TRUE si on veut tester le nb. Richardson entre la sfce
[3817]155!                  et zref par rapport au Ri entre la sfce et la 1ere couche
[5274]156! ri1-----input-R- nb. Richardson entre la surface et la 1ere couche
[3817]157!
158! pref----input-R- pression au niveau de reference
159! delu----input-R- anomalie du vent par rapport au 1er niveau
160! delte---input-R- anomalie de la temperature potentielle par rapport a la surface
161! delq----input-R- anomalie de l'humidite relative par rapport a la surface
162!
163      INTEGER, intent(in) :: klon, knon, nsrf
[5274]164      LOGICAL, intent(in) :: zxli, okri
[3817]165      REAL, dimension(klon), intent(in) :: speed, temp, q_zref
166      REAL, intent(in) :: zref
[4722]167      REAL, dimension(klon), intent(in) :: ts, qsurf, psol
168      REAL, dimension(klon), intent(inout) :: z0m, z0h
[5274]169      REAL, dimension(klon), intent(in) :: cdrm, cdrh, ri1
170      REAL, dimension(klon), intent(inout) :: s_pblh
171      REAL, dimension(klon), intent(in) :: prain
172      REAL, dimension(klon), intent(in) :: tsol
173      REAL, DIMENSION(klon), INTENT(IN) :: pat1 !pression premier lev
[3817]174      INTEGER, INTENT(IN)  :: iri1 ! Richardson de la 1ere couche
175!
176      REAL, dimension(klon), intent(out) :: pref, delm, delh, zri1
[5274]177!171220     REAL, dimension(klon) :: cdram, cdrah, zri1
[3817]178      REAL, dimension(klon) :: cdram, cdrah
179!-----------------------------------------------------------------------
[5301]180      !
[3817]181! Variables locales 
182      INTEGER :: i
183      REAL, dimension(klon) :: cdran, gref,ycdragm
184!
185!-------------------------------------------------------------------------
186      DO i=1, knon
187        gref(i) = zref*RG
188      ENDDO
189!
190! Richardson at reference level
191!
[4478]192      CALL cdrag(knon, nsrf, &
[3817]193                    speed, temp, q_zref, gref, &
[4722]194                    psol, s_pblh, ts, qsurf, z0m, z0h, &
[3817]195                    ri1, iri1, &
[4722]196                    cdram, cdrah, zri1, pref, prain, tsol, pat1)
[3817]197      DO i = 1, knon
198        delm(i) = sqrt(cdrm(i))/sqrt(cdram(i))
199!verifier que temp est un temperat potentielle.
200        delh(i)= (cdrh(i)* sqrt(cdram(i)))/ &
201                   (cdrah(i)*sqrt(cdrm(i)))
202      ENDDO
203!
204      RETURN
205      END SUBROUTINE screencn
206END MODULE screenc_mod
Note: See TracBrowser for help on using the repository browser.