source: LMDZ6/branches/Amaury_dev/libf/phylmd/screenc_mod.F90 @ 5139

Last change on this file since 5139 was 5139, checked in by abarral, 3 months ago

Put nuage.h, flux_arp.h, compbl.h into modules
Move unused phylmd/ini_hist* to obsolete
  • Property svn:executable set to *
File size: 8.6 KB
Line 
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, &
20                         pref, delu, delte, delq, s_pblh, prain, tsol, pat1)
21        USE lmdz_flux_arp, ONLY: fsens, flat, betaevap, ust, tg, ok_flux_surf, ok_prescr_ust, ok_prescr_beta, ok_forc_tsurf
22
23      IMPLICIT NONE
24!-----------------------------------------------------------------------
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)
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
51! okri----input-L- TRUE si on veut tester le nb. Richardson entre la sfce
52!                  et zref par rapport au Ri entre la sfce et la 1ere couche
53! ri1-----input-R- nb. Richardson entre la surface et la 1ere couche
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
61      LOGICAL, INTENT(IN) :: zxli, okri
62      REAL, DIMENSION(klon), INTENT(IN) :: speed, temp, q_zref
63      REAL, INTENT(IN) :: zref
64      REAL, DIMENSION(klon), INTENT(IN) :: ts
65      REAL, DIMENSION(klon), INTENT(IN) :: qsurf, psol
66      REAL, DIMENSION(klon), INTENT(INOUT):: z0m, z0h
67      REAL, DIMENSION(klon), INTENT(IN) :: ustar, testar, qstar, ri1
68
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     
73
74      REAL, DIMENSION(klon), INTENT(OUT) :: pref, delu, delte, delq
75!-----------------------------------------------------------------------
76      include "YOMCST.h"
77
78! Variables locales 
79      INTEGER :: i
80      REAL, DIMENSION(klon) :: cdram, cdrah, cdran, zri1, gref,ycdragm,zri_zero
81
82!-------------------------------------------------------------------------
83      DO i=1, knon
84        gref(i) = zref*RG
85      ENDDO
86
87! Richardson at reference level
88
89!      CALL coefcdrag (klon, knon, nsrf, zxli, &
90!                    speed, temp, q_zref, gref, &
91!                    psol, ts, qsurf, rugos, &
92!                    okri, ri1, &
93!                    cdram, cdrah, cdran, zri1, &
94!                    pref)
95! Fuxing WANG, 04/03/2015, replace the coefcdrag by the merged version: cdrag
96      CALL cdrag (knon, nsrf, &
97                    speed, temp, q_zref, gref, &
98                    psol, s_pblh, ts, qsurf, z0m, z0h, &
99                    zri_zero,0,                &
100                    cdram, cdrah, zri1, pref, prain, tsol, pat1)
101      DO i = 1, knon
102        IF(ok_prescr_ust) THEN
103! La aussi il faut forcer avec ust (FC + MPL 20160210)
104        ycdragm(i) = ust*ust/(1.+speed(i))/speed(i)
105        cdram=ycdragm
106        delu(i) = ust/sqrt(cdram(i))
107        ELSE
108        delu(i) = ustar(i)/sqrt(cdram(i))
109        ENDIF
110        delte(i)= (testar(i)* sqrt(cdram(i)))/ &
111                   cdrah(i)
112        delq(i)= (qstar(i)* sqrt(cdram(i)))/ &
113                  cdrah(i)
114      ENDDO
115
116
117      END SUBROUTINE screenc
118
119      SUBROUTINE screencn(klon, knon, nsrf, zxli, &
120                         speed, temp, q_zref, zref, &
121                         ts, qsurf, z0m, z0h, psol, &
122                         cdrm, cdrh,  okri, &
123                         ri1, iri1, &
124                         pref, delm, delh, zri1, s_pblh, prain, tsol, pat1)
125        USE lmdz_flux_arp, ONLY: fsens, flat, betaevap, ust, tg, ok_flux_surf, ok_prescr_ust, ok_prescr_beta, ok_forc_tsurf
126
127      IMPLICIT NONE
128!-----------------------------------------------------------------------
129
130! Objet : calcul "correcteur" des anomalies du vent, de la temperature
131!         potentielle et de l'humidite relative au niveau de reference zref et
132!         par rapport au 1er niveau (pour u) ou a la surface (pour theta et q)
133!         a partir des equations de Louis.
134
135! Reference : Hess, Colman et McAvaney (1995)
136
137! I. Musat, 01.07.2002
138!-----------------------------------------------------------------------
139
140! klon----input-I- dimension de la grille physique (= nb_pts_latitude X nb_pts_longitude)
141! knon----input-I- nombre de points pour un type de surface
142! nsrf----input-I- indice pour le type de surface; voir indice_sol_mod.F90
143! zxli----input-L- TRUE si calcul des cdrags selon Laurent Li
144! speed---input-R- module du vent au 1er niveau du modele
145! temp----input-R- temperature de l'air au 1er niveau du modele
146! q_zref--input-R- humidite relative au 1er niveau du modele
147! zref----input-R- altitude de reference
148! ts------input-R- temperature de l'air a la surface
149! qsurf---input-R- humidite relative a la surface
150! z0m, z0h---input-R- rugosite
151! psol----input-R- pression au sol
152! ustar---input-R- facteur d'echelle pour le vent
153! testar--input-R- facteur d'echelle pour la temperature potentielle
154! qstar---input-R- facteur d'echelle pour l'humidite relative
155! okri----input-L- TRUE si on veut tester le nb. Richardson entre la sfce
156!                  et zref par rapport au Ri entre la sfce et la 1ere couche
157! ri1-----input-R- nb. Richardson entre la surface et la 1ere couche
158
159! pref----input-R- pression au niveau de reference
160! delu----input-R- anomalie du vent par rapport au 1er niveau
161! delte---input-R- anomalie de la temperature potentielle par rapport a la surface
162! delq----input-R- anomalie de l'humidite relative par rapport a la surface
163
164      INTEGER, INTENT(IN) :: klon, knon, nsrf
165      LOGICAL, INTENT(IN) :: zxli, okri
166      REAL, DIMENSION(klon), INTENT(IN) :: speed, temp, q_zref
167      REAL, INTENT(IN) :: zref
168      REAL, DIMENSION(klon), INTENT(IN) :: ts, qsurf, psol
169      REAL, DIMENSION(klon), INTENT(INOUT) :: z0m, z0h
170      REAL, DIMENSION(klon), INTENT(IN) :: cdrm, cdrh, ri1
171      REAL, DIMENSION(klon), INTENT(INOUT) :: s_pblh
172      REAL, DIMENSION(klon), INTENT(IN) :: prain
173      REAL, DIMENSION(klon), INTENT(IN) :: tsol
174      REAL, DIMENSION(klon), INTENT(IN) :: pat1 !pression premier lev     
175      INTEGER, INTENT(IN)  :: iri1 ! Richardson de la 1ere couche
176
177      REAL, DIMENSION(klon), INTENT(OUT) :: pref, delm, delh, zri1
178!171220     REAL, DIMENSION(klon) :: cdram, cdrah, zri1
179      REAL, DIMENSION(klon) :: cdram, cdrah
180!-----------------------------------------------------------------------
181      include "YOMCST.h"
182
183! Variables locales 
184      INTEGER :: i
185      REAL, DIMENSION(klon) :: cdran, gref,ycdragm
186
187!-------------------------------------------------------------------------
188      DO i=1, knon
189        gref(i) = zref*RG
190      ENDDO
191
192! Richardson at reference level
193
194      CALL cdrag(knon, nsrf, &
195                    speed, temp, q_zref, gref, &
196                    psol, s_pblh, ts, qsurf, z0m, z0h, &
197                    ri1, iri1, &
198                    cdram, cdrah, zri1, pref, prain, tsol, pat1)
199      DO i = 1, knon
200        delm(i) = sqrt(cdrm(i))/sqrt(cdram(i))
201!verifier que temp est un temperat potentielle.
202        delh(i)= (cdrh(i)* sqrt(cdram(i)))/ &
203                   (cdrah(i)*sqrt(cdrm(i)))
204      ENDDO
205
206
207      END SUBROUTINE screencn
208END MODULE screenc_mod
Note: See TracBrowser for help on using the repository browser.