1 | ! |
---|
2 | ! $Id$ |
---|
3 | ! |
---|
4 | SUBROUTINE freinage(knon, uu, vv, & |
---|
5 | tt,veget,lai, height,ypaprs,ypplay,drag_pro,d_u,d_v) |
---|
6 | |
---|
7 | !ONLINE: |
---|
8 | use dimphy, only: klon, klev |
---|
9 | ! USE control, ONLY: nvm |
---|
10 | ! USE indice_sol_mod, only : nvm_orch |
---|
11 | |
---|
12 | USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO & |
---|
13 | , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA & |
---|
14 | , R_ecc, R_peri, R_incl & |
---|
15 | , RA, RG, R1SA & |
---|
16 | , RSIGMA & |
---|
17 | , R, RMD, RMV, RD, RV, RCPD & |
---|
18 | , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12 & |
---|
19 | , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w & |
---|
20 | , RCW, RCS & |
---|
21 | , RLVTT, RLSTT, RLMLT, RTT, RATM & |
---|
22 | , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS & |
---|
23 | , RALPD, RBETD, RGAMD |
---|
24 | IMPLICIT NONE |
---|
25 | |
---|
26 | |
---|
27 | |
---|
28 | include "clesphys.h" |
---|
29 | include "YOEGWD.h" |
---|
30 | !FC |
---|
31 | include "dimpft.h" |
---|
32 | include "compbl.h" |
---|
33 | |
---|
34 | ! 0. DECLARATIONS: |
---|
35 | |
---|
36 | ! 0.1 INPUTS |
---|
37 | |
---|
38 | REAL, DIMENSION(klon,klev), INTENT(IN) :: ypplay |
---|
39 | REAL, DIMENSION(klon,klev+1), INTENT(IN) :: ypaprs |
---|
40 | |
---|
41 | |
---|
42 | REAL, DIMENSION(klon, klev), INTENT(IN) :: uu |
---|
43 | REAL, DIMENSION(klon, klev), INTENT(IN) :: vv |
---|
44 | REAL, DIMENSION(klon, klev), INTENT(IN) :: tt |
---|
45 | REAL, DIMENSION(klon,nvm_lmdz), INTENT(IN) :: veget,lai |
---|
46 | REAL, DIMENSION(klon,nvm_lmdz), INTENT(IN) :: height |
---|
47 | |
---|
48 | REAL, DIMENSION(klon,klev) :: wind |
---|
49 | REAL, DIMENSION(klon, klev) :: yzlay |
---|
50 | INTEGER knon |
---|
51 | |
---|
52 | ! 0.2 OUTPUTS |
---|
53 | |
---|
54 | REAL, DIMENSION(klon, klev), INTENT(OUT) :: d_v ! change in v |
---|
55 | REAL, DIMENSION(klon, klev), INTENT(OUT) :: d_u ! change in v |
---|
56 | !knon nombre de points concernes |
---|
57 | REAL, DIMENSION(klon,klev) :: sumveg ! change in v |
---|
58 | |
---|
59 | REAL, DIMENSION(klon,klev), INTENT(OUT) :: drag_pro |
---|
60 | ! (KLON, KLEV) tendencies on winds |
---|
61 | |
---|
62 | |
---|
63 | INTEGER k,jv,i |
---|
64 | |
---|
65 | |
---|
66 | !FCCCC REAL Cd_frein |
---|
67 | |
---|
68 | ! 0.3.1 LOCAL VARIABLE |
---|
69 | |
---|
70 | |
---|
71 | !----------------------------------------------------------------- |
---|
72 | |
---|
73 | ! 1. INITIALISATIONS |
---|
74 | |
---|
75 | |
---|
76 | ! Cd_frein = 7.5E-2 ! (0.075) ! Drag from MASSON 2009 |
---|
77 | !FC ESSAI |
---|
78 | ! Cd_frein = 1.5E-2 ! (0.075) ! Drag from MASSON 2009 |
---|
79 | ! Cd_frein = 0.005 ! (0.075) ! Drag from MASSON 2009 |
---|
80 | |
---|
81 | ! initialisation |
---|
82 | d_u(:,:) =0. |
---|
83 | d_v(:,:) =0. |
---|
84 | drag_pro(:,:) =0. |
---|
85 | sumveg(:,:) =0. |
---|
86 | !! print*, "Cd_frein" , Cd_frein |
---|
87 | |
---|
88 | wind(:,:)= sqrt(uu(:,:)*uu(:,:)+vv(:,:)*vv(:,:)) |
---|
89 | |
---|
90 | yzlay(1:knon,1)= & |
---|
91 | RD*tt(1:knon,1)/(0.5*(ypaprs(1:knon,1)+ypplay(1:knon,1))) & |
---|
92 | *(ypaprs(1:knon,1)-ypplay(1:knon,1))/RG |
---|
93 | DO k=2,klev |
---|
94 | yzlay(1:knon,k)= & |
---|
95 | yzlay(1:knon,k-1)+RD*0.5*(tt(1:knon,k-1)+tt(1:knon,k)) & |
---|
96 | /ypaprs(1:knon,k)*(ypplay(1:knon,k-1)-ypplay(1:knon,k))/RG |
---|
97 | END DO |
---|
98 | |
---|
99 | ! verifier les indexes ..... |
---|
100 | !! print*, " calcul de drag_pro FC " |
---|
101 | |
---|
102 | do k= 1,klev |
---|
103 | |
---|
104 | do jv=2,nvm_lmdz ! (on peut faire 9 ?) |
---|
105 | |
---|
106 | do i=1,knon |
---|
107 | |
---|
108 | sumveg(i,k)= sumveg(i,k)+ veget(i,jv) |
---|
109 | |
---|
110 | ! if ( (height(i,jv) .gt. yzlay(i,k)) .AND. (height(i,jv) .gt. 0.1) .and. LAI(i,jv).gt.0. ) then |
---|
111 | if ( (height(i,jv) .gt. yzlay(i,k)) .AND. (height(i,jv) .gt. 0.1) ) then |
---|
112 | !FC attention veut on le test sur le LAI ? |
---|
113 | if (ifl_pbltree.eq.1) then |
---|
114 | drag_pro(i,k)= drag_pro(i,k)+ & |
---|
115 | veget(i,jv) |
---|
116 | elseif (ifl_pbltree.eq.2) then |
---|
117 | drag_pro(i,k)= drag_pro(i,k)+ & |
---|
118 | 6*LAI(i,jv)*veget(i,jv)*( yzlay(i,k)*(height(i,jv)-yzlay(i,k))/(height(i,jv)*height(i,jv)+ 0.01)) |
---|
119 | elseif (ifl_pbltree.eq.3) then |
---|
120 | drag_pro(i,k)= drag_pro(i,k)+ & |
---|
121 | veget(i,jv)*( yzlay(i,k)*(height(i,jv)-yzlay(i,k))/(height(i,jv)*height(i,jv)+ 0.01)) |
---|
122 | elseif (ifl_pbltree.eq.0) then |
---|
123 | drag_pro(i,k)=0.0 |
---|
124 | endif |
---|
125 | else |
---|
126 | drag_pro(i,k)= drag_pro(i,k) |
---|
127 | endif |
---|
128 | |
---|
129 | |
---|
130 | enddo |
---|
131 | enddo |
---|
132 | enddo |
---|
133 | do k=1,klev |
---|
134 | where (sumveg(1:knon,k) > 0.05 ) |
---|
135 | ! drag_pro(1:knon,k)=Cd_frein*drag_pro(1:knon,k)/sumveg(1:knon,k) |
---|
136 | drag_pro(1:knon,k)=Cd_frein*drag_pro(1:knon,k) |
---|
137 | elsewhere |
---|
138 | drag_pro(1:knon,k)=0.0 |
---|
139 | endwhere |
---|
140 | d_u(1:knon,k) =(-1)*drag_pro(1:knon,k)*uu(1:knon,k)*wind(1:knon,k) |
---|
141 | d_v(1:knon,k) =(-1)*drag_pro(1:knon,k)*vv(1:knon,k)*wind(1:knon,k) |
---|
142 | enddo |
---|
143 | return |
---|
144 | |
---|
145 | END SUBROUTINE freinage |
---|
146 | |
---|