| Directory: | ./ |
|---|---|
| File: | phys/climb_wind_mod.f90 |
| Date: | 2022-01-11 19:19:34 |
| Exec | Total | Coverage | |
|---|---|---|---|
| Lines: | 79 | 101 | 78.2% |
| Branches: | 136 | 212 | 64.2% |
| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | ! | ||
| 2 | MODULE climb_wind_mod | ||
| 3 | ! | ||
| 4 | ! Module to solve the verctical diffusion of the wind components "u" and "v". | ||
| 5 | ! | ||
| 6 | USE dimphy | ||
| 7 | |||
| 8 | IMPLICIT NONE | ||
| 9 | |||
| 10 | SAVE | ||
| 11 | PRIVATE | ||
| 12 | |||
| 13 | REAL, DIMENSION(:), ALLOCATABLE :: alf1, alf2 | ||
| 14 | !$OMP THREADPRIVATE(alf1,alf2) | ||
| 15 | REAL, DIMENSION(:,:), ALLOCATABLE :: Kcoefm | ||
| 16 | !$OMP THREADPRIVATE(Kcoefm) | ||
| 17 | REAL, DIMENSION(:,:), ALLOCATABLE :: Ccoef_U, Dcoef_U | ||
| 18 | !$OMP THREADPRIVATE(Ccoef_U, Dcoef_U) | ||
| 19 | REAL, DIMENSION(:,:), ALLOCATABLE :: Ccoef_V, Dcoef_V | ||
| 20 | !$OMP THREADPRIVATE(Ccoef_V, Dcoef_V) | ||
| 21 | REAL, DIMENSION(:), ALLOCATABLE :: Acoef_U, Bcoef_U | ||
| 22 | !$OMP THREADPRIVATE(Acoef_U, Bcoef_U) | ||
| 23 | REAL, DIMENSION(:), ALLOCATABLE :: Acoef_V, Bcoef_V | ||
| 24 | !$OMP THREADPRIVATE(Acoef_V, Bcoef_V) | ||
| 25 | LOGICAL :: firstcall=.TRUE. | ||
| 26 | !$OMP THREADPRIVATE(firstcall) | ||
| 27 | |||
| 28 | |||
| 29 | PUBLIC :: climb_wind_down, climb_wind_up | ||
| 30 | |||
| 31 | CONTAINS | ||
| 32 | ! | ||
| 33 | !**************************************************************************************** | ||
| 34 | ! | ||
| 35 | 29963097 | SUBROUTINE climb_wind_init | |
| 36 | |||
| 37 | INTEGER :: ierr | ||
| 38 | CHARACTER(len = 20) :: modname = 'climb_wind_init' | ||
| 39 | |||
| 40 | !**************************************************************************************** | ||
| 41 | ! Allocation of global module variables | ||
| 42 | ! | ||
| 43 | !**************************************************************************************** | ||
| 44 | |||
| 45 |
4/8✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
|
1 | ALLOCATE(alf1(klon), stat=ierr) |
| 46 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (ierr /= 0) CALL abort_physic(modname,'Pb in allocate alf1',1) |
| 47 | |||
| 48 |
4/8✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
|
1 | ALLOCATE(alf2(klon), stat=ierr) |
| 49 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (ierr /= 0) CALL abort_physic(modname,'Pb in allocate alf2',1) |
| 50 | |||
| 51 |
6/12✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
|
2 | ALLOCATE(Kcoefm(klon,klev), stat=ierr) |
| 52 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (ierr /= 0) CALL abort_physic(modname,'Pb in allocate Kcoefm',1) |
| 53 | |||
| 54 |
6/12✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
|
2 | ALLOCATE(Ccoef_U(klon,klev), stat=ierr) |
| 55 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (ierr /= 0) CALL abort_physic(modname,'Pb in allocate Ccoef_U',1) |
| 56 | |||
| 57 |
6/12✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
|
2 | ALLOCATE(Dcoef_U(klon,klev), stat=ierr) |
| 58 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (ierr /= 0) CALL abort_physic(modname,'Pb in allocation Dcoef_U',1) |
| 59 | |||
| 60 |
6/12✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
|
2 | ALLOCATE(Ccoef_V(klon,klev), stat=ierr) |
| 61 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (ierr /= 0) CALL abort_physic(modname,'Pb in allocation Ccoef_V',1) |
| 62 | |||
| 63 |
6/12✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
|
2 | ALLOCATE(Dcoef_V(klon,klev), stat=ierr) |
| 64 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (ierr /= 0) CALL abort_physic(modname,'Pb in allocation Dcoef_V',1) |
| 65 | |||
| 66 |
16/32✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 1 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 1 times.
✓ Branch 14 taken 1 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 1 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 1 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 1 times.
✓ Branch 22 taken 1 times.
✗ Branch 23 not taken.
✓ Branch 24 taken 1 times.
✗ Branch 25 not taken.
✓ Branch 26 taken 1 times.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 1 times.
✓ Branch 30 taken 1 times.
✗ Branch 31 not taken.
|
1 | ALLOCATE(Acoef_U(klon), Bcoef_U(klon), Acoef_V(klon), Bcoef_V(klon), STAT=ierr) |
| 67 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF ( ierr /= 0 ) PRINT*,' pb in allloc Acoef_U and Bcoef_U, ierr=', ierr |
| 68 | |||
| 69 | 1 | firstcall=.FALSE. | |
| 70 | |||
| 71 | 1 | END SUBROUTINE climb_wind_init | |
| 72 | ! | ||
| 73 | !**************************************************************************************** | ||
| 74 | ! | ||
| 75 | 1920 | SUBROUTINE climb_wind_down(knon, dtime, coef_in, pplay, paprs, temp, delp, u_old, v_old, & | |
| 76 | !!! nrlmd le 02/05/2011 | ||
| 77 | Ccoef_U_out, Ccoef_V_out, Dcoef_U_out, Dcoef_V_out, & | ||
| 78 | Kcoef_m_out, alf_1_out, alf_2_out, & | ||
| 79 | !!! | ||
| 80 | Acoef_U_out, Acoef_V_out, Bcoef_U_out, Bcoef_V_out) | ||
| 81 | ! | ||
| 82 | ! This routine calculates for the wind components u and v, | ||
| 83 | ! recursivly the coefficients C and D in equation | ||
| 84 | ! X(k) = C(k) + D(k)*X(k-1), X=[u,v], k=[1,klev] is the vertical layer. | ||
| 85 | ! | ||
| 86 | ! | ||
| 87 | |||
| 88 | ! Input arguments | ||
| 89 | !**************************************************************************************** | ||
| 90 | INTEGER, INTENT(IN) :: knon | ||
| 91 | REAL, INTENT(IN) :: dtime | ||
| 92 | REAL, DIMENSION(klon,klev), INTENT(IN) :: coef_in | ||
| 93 | REAL, DIMENSION(klon,klev), INTENT(IN) :: pplay ! pres au milieu de couche (Pa) | ||
| 94 | REAL, DIMENSION(klon,klev+1), INTENT(IN) :: paprs ! pression a inter-couche (Pa) | ||
| 95 | REAL, DIMENSION(klon,klev), INTENT(IN) :: temp ! temperature | ||
| 96 | REAL, DIMENSION(klon,klev), INTENT(IN) :: delp | ||
| 97 | REAL, DIMENSION(klon,klev), INTENT(IN) :: u_old | ||
| 98 | REAL, DIMENSION(klon,klev), INTENT(IN) :: v_old | ||
| 99 | |||
| 100 | ! Output arguments | ||
| 101 | !**************************************************************************************** | ||
| 102 | REAL, DIMENSION(klon), INTENT(OUT) :: Acoef_U_out | ||
| 103 | REAL, DIMENSION(klon), INTENT(OUT) :: Acoef_V_out | ||
| 104 | REAL, DIMENSION(klon), INTENT(OUT) :: Bcoef_U_out | ||
| 105 | REAL, DIMENSION(klon), INTENT(OUT) :: Bcoef_V_out | ||
| 106 | |||
| 107 | !!! nrlmd le 02/05/2011 | ||
| 108 | REAL, DIMENSION(klon,klev), INTENT(OUT) :: Ccoef_U_out | ||
| 109 | REAL, DIMENSION(klon,klev), INTENT(OUT) :: Ccoef_V_out | ||
| 110 | REAL, DIMENSION(klon,klev), INTENT(OUT) :: Dcoef_U_out | ||
| 111 | REAL, DIMENSION(klon,klev), INTENT(OUT) :: Dcoef_V_out | ||
| 112 | REAL, DIMENSION(klon,klev), INTENT(OUT) :: Kcoef_m_out | ||
| 113 | REAL, DIMENSION(klon), INTENT(OUT) :: alf_1_out | ||
| 114 | REAL, DIMENSION(klon), INTENT(OUT) :: alf_2_out | ||
| 115 | !!! | ||
| 116 | |||
| 117 | ! Local variables | ||
| 118 | !**************************************************************************************** | ||
| 119 | REAL, DIMENSION(klon) :: u1lay, v1lay | ||
| 120 | INTEGER :: k, i | ||
| 121 | |||
| 122 | ! Include | ||
| 123 | !**************************************************************************************** | ||
| 124 | INCLUDE "YOMCST.h" | ||
| 125 | INCLUDE "compbl.h" | ||
| 126 | |||
| 127 | !**************************************************************************************** | ||
| 128 | ! Initialize module | ||
| 129 |
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 1919 times.
|
1920 | IF (firstcall) CALL climb_wind_init |
| 130 | |||
| 131 | !**************************************************************************************** | ||
| 132 | ! Calculate the coefficients C and D in : u(k) = C(k) + D(k)*u(k-1) | ||
| 133 | ! | ||
| 134 | !**************************************************************************************** | ||
| 135 | ! - Define alpha (alf1 and alf2) | ||
| 136 |
2/2✓ Branch 0 taken 1908480 times.
✓ Branch 1 taken 1920 times.
|
1910400 | alf1(:) = 1.0 |
| 137 |
2/2✓ Branch 0 taken 1908480 times.
✓ Branch 1 taken 1920 times.
|
1910400 | alf2(:) = 1.0 - alf1(:) |
| 138 | |||
| 139 | ! - Calculate the coefficients K | ||
| 140 |
4/4✓ Branch 0 taken 74880 times.
✓ Branch 1 taken 1920 times.
✓ Branch 2 taken 74430720 times.
✓ Branch 3 taken 74880 times.
|
74507520 | Kcoefm(:,:) = 0.0 |
| 141 |
2/2✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 1920 times.
|
74880 | DO k = 2, klev |
| 142 |
2/2✓ Branch 0 taken 29961176 times.
✓ Branch 1 taken 72960 times.
|
30036056 | DO i=1,knon |
| 143 | Kcoefm(i,k) = coef_in(i,k)*RG*RG*dtime/(pplay(i,k-1)-pplay(i,k)) & | ||
| 144 | 72960 | *(paprs(i,k)*2/(temp(i,k)+temp(i,k-1))/RD)**2 | |
| 145 | END DO | ||
| 146 | END DO | ||
| 147 | |||
| 148 | ! - Calculate the coefficients C and D, component "u" | ||
| 149 | CALL calc_coef(knon, Kcoefm(:,:), delp(:,:), & | ||
| 150 | u_old(:,:), alf1(:), alf2(:), & | ||
| 151 | 1920 | Ccoef_U(:,:), Dcoef_U(:,:), Acoef_U(:), Bcoef_U(:)) | |
| 152 | |||
| 153 | ! - Calculate the coefficients C and D, component "v" | ||
| 154 | CALL calc_coef(knon, Kcoefm(:,:), delp(:,:), & | ||
| 155 | v_old(:,:), alf1(:), alf2(:), & | ||
| 156 | 1920 | Ccoef_V(:,:), Dcoef_V(:,:), Acoef_V(:), Bcoef_V(:)) | |
| 157 | |||
| 158 | !**************************************************************************************** | ||
| 159 | ! 6) | ||
| 160 | ! Return the first layer in output variables | ||
| 161 | ! | ||
| 162 | !**************************************************************************************** | ||
| 163 |
2/2✓ Branch 0 taken 1908480 times.
✓ Branch 1 taken 1920 times.
|
1910400 | Acoef_U_out = Acoef_U |
| 164 |
2/2✓ Branch 0 taken 1908480 times.
✓ Branch 1 taken 1920 times.
|
1910400 | Bcoef_U_out = Bcoef_U |
| 165 |
2/2✓ Branch 0 taken 1908480 times.
✓ Branch 1 taken 1920 times.
|
1910400 | Acoef_V_out = Acoef_V |
| 166 |
2/2✓ Branch 0 taken 1908480 times.
✓ Branch 1 taken 1920 times.
|
1910400 | Bcoef_V_out = Bcoef_V |
| 167 | |||
| 168 | !**************************************************************************************** | ||
| 169 | ! 7) | ||
| 170 | ! If Pbl is split, return also the other layers in output variables | ||
| 171 | ! | ||
| 172 | !**************************************************************************************** | ||
| 173 | !!! jyg le 07/02/2012 | ||
| 174 | !!jyg IF (mod(iflag_pbl_split,2) .eq.1) THEN | ||
| 175 | 1920 | IF (mod(iflag_pbl_split,10) .ge.1) THEN | |
| 176 | !!! nrlmd le 02/05/2011 | ||
| 177 | ✗ | DO k= 1, klev | |
| 178 | ✗ | DO i= 1, klon | |
| 179 | ✗ | Ccoef_U_out(i,k) = Ccoef_U(i,k) | |
| 180 | ✗ | Ccoef_V_out(i,k) = Ccoef_V(i,k) | |
| 181 | ✗ | Dcoef_U_out(i,k) = Dcoef_U(i,k) | |
| 182 | ✗ | Dcoef_V_out(i,k) = Dcoef_V(i,k) | |
| 183 | ✗ | Kcoef_m_out(i,k) = Kcoefm(i,k) | |
| 184 | ENDDO | ||
| 185 | ENDDO | ||
| 186 | ✗ | DO i= 1, klon | |
| 187 | ✗ | alf_1_out(i) = alf1(i) | |
| 188 | ✗ | alf_2_out(i) = alf2(i) | |
| 189 | ENDDO | ||
| 190 | !!! | ||
| 191 | ENDIF ! (mod(iflag_pbl_split,2) .ge.1) | ||
| 192 | !!! | ||
| 193 | |||
| 194 | 1920 | END SUBROUTINE climb_wind_down | |
| 195 | ! | ||
| 196 | !**************************************************************************************** | ||
| 197 | ! | ||
| 198 | 3840 | SUBROUTINE calc_coef(knon, Kcoef, delp, X, alfa1, alfa2, Ccoef, Dcoef, Acoef, Bcoef) | |
| 199 | ! | ||
| 200 | ! Find the coefficients C and D in fonction of alfa, K and delp | ||
| 201 | ! | ||
| 202 | ! Input arguments | ||
| 203 | !**************************************************************************************** | ||
| 204 | INTEGER, INTENT(IN) :: knon | ||
| 205 | REAL, DIMENSION(klon,klev), INTENT(IN) :: Kcoef, delp | ||
| 206 | REAL, DIMENSION(klon,klev), INTENT(IN) :: X | ||
| 207 | REAL, DIMENSION(klon), INTENT(IN) :: alfa1, alfa2 | ||
| 208 | |||
| 209 | ! Output arguments | ||
| 210 | !**************************************************************************************** | ||
| 211 | REAL, DIMENSION(klon), INTENT(OUT) :: Acoef, Bcoef | ||
| 212 | REAL, DIMENSION(klon,klev), INTENT(OUT) :: Ccoef, Dcoef | ||
| 213 | |||
| 214 | ! local variables | ||
| 215 | !**************************************************************************************** | ||
| 216 | INTEGER :: k, i | ||
| 217 | REAL :: buf | ||
| 218 | |||
| 219 | INCLUDE "YOMCST.h" | ||
| 220 | !**************************************************************************************** | ||
| 221 | ! | ||
| 222 | |||
| 223 | ! Calculate coefficients C and D at top level, k=klev | ||
| 224 | ! | ||
| 225 |
4/4✓ Branch 0 taken 3840 times.
✓ Branch 1 taken 149760 times.
✓ Branch 2 taken 148861440 times.
✓ Branch 3 taken 149760 times.
|
149015040 | Ccoef(:,:) = 0.0 |
| 226 |
4/4✓ Branch 0 taken 149760 times.
✓ Branch 1 taken 3840 times.
✓ Branch 2 taken 148861440 times.
✓ Branch 3 taken 149760 times.
|
149015040 | Dcoef(:,:) = 0.0 |
| 227 | |||
| 228 |
2/2✓ Branch 0 taken 1576904 times.
✓ Branch 1 taken 3840 times.
|
1580744 | DO i = 1, knon |
| 229 | 1576904 | buf = delp(i,klev) + Kcoef(i,klev) | |
| 230 | |||
| 231 | 1576904 | Ccoef(i,klev) = X(i,klev)*delp(i,klev)/buf | |
| 232 | 1580744 | Dcoef(i,klev) = Kcoef(i,klev)/buf | |
| 233 | END DO | ||
| 234 | |||
| 235 | ! | ||
| 236 | ! Calculate coefficients C and D at top level (klev-1) <= k <= 2 | ||
| 237 | ! | ||
| 238 |
2/2✓ Branch 0 taken 142080 times.
✓ Branch 1 taken 3840 times.
|
145920 | DO k=(klev-1),2,-1 |
| 239 |
2/2✓ Branch 0 taken 58345448 times.
✓ Branch 1 taken 142080 times.
|
58491368 | DO i = 1, knon |
| 240 | 58345448 | buf = delp(i,k) + Kcoef(i,k) + Kcoef(i,k+1)*(1.-Dcoef(i,k+1)) | |
| 241 | |||
| 242 | 58345448 | Ccoef(i,k) = (X(i,k)*delp(i,k) + Kcoef(i,k+1)*Ccoef(i,k+1))/buf | |
| 243 | 58487528 | Dcoef(i,k) = Kcoef(i,k)/buf | |
| 244 | END DO | ||
| 245 | END DO | ||
| 246 | |||
| 247 | ! | ||
| 248 | ! Calculate coeffiecent A and B at surface | ||
| 249 | ! | ||
| 250 |
2/2✓ Branch 0 taken 1576904 times.
✓ Branch 1 taken 3840 times.
|
1580744 | DO i = 1, knon |
| 251 | 1576904 | buf = delp(i,1) + Kcoef(i,2)*(1-Dcoef(i,2)) | |
| 252 | 1576904 | Acoef(i) = (X(i,1)*delp(i,1) + Kcoef(i,2)*Ccoef(i,2))/buf | |
| 253 | 3840 | Bcoef(i) = -RG/buf | |
| 254 | END DO | ||
| 255 | |||
| 256 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1920 times.
|
29963096 | END SUBROUTINE calc_coef |
| 257 | ! | ||
| 258 | !**************************************************************************************** | ||
| 259 | ! | ||
| 260 | |||
| 261 | 1920 | SUBROUTINE climb_wind_up(knon, dtime, u_old, v_old, flx_u1, flx_v1, & | |
| 262 | !!! nrlmd le 02/05/2011 | ||
| 263 | 1920 | Acoef_U_in, Acoef_V_in, Bcoef_U_in, Bcoef_V_in, & | |
| 264 | 1920 | Ccoef_U_in, Ccoef_V_in, Dcoef_U_in, Dcoef_V_in, & | |
| 265 | Kcoef_m_in, & | ||
| 266 | !!! | ||
| 267 | flx_u_new, flx_v_new, d_u_new, d_v_new) | ||
| 268 | ! | ||
| 269 | ! Diffuse the wind components from the surface layer and up to the top layer. | ||
| 270 | ! Coefficents A, B, C and D are known from before. Start values for the diffusion are the | ||
| 271 | ! momentum fluxes at surface. | ||
| 272 | ! | ||
| 273 | ! u(k=1) = A + B*flx*dtime | ||
| 274 | ! u(k) = C(k) + D(k)*u(k-1) [2 <= k <= klev] | ||
| 275 | ! | ||
| 276 | !**************************************************************************************** | ||
| 277 | |||
| 278 | ! Input arguments | ||
| 279 | !**************************************************************************************** | ||
| 280 | INTEGER, INTENT(IN) :: knon | ||
| 281 | REAL, INTENT(IN) :: dtime | ||
| 282 | REAL, DIMENSION(klon,klev), INTENT(IN) :: u_old | ||
| 283 | REAL, DIMENSION(klon,klev), INTENT(IN) :: v_old | ||
| 284 | REAL, DIMENSION(klon), INTENT(IN) :: flx_u1, flx_v1 ! momentum flux | ||
| 285 | |||
| 286 | !!! nrlmd le 02/05/2011 | ||
| 287 | REAL, DIMENSION(klon), INTENT(IN) :: Acoef_U_in,Acoef_V_in, Bcoef_U_in, Bcoef_V_in | ||
| 288 | REAL, DIMENSION(klon,klev), INTENT(IN) :: Ccoef_U_in, Ccoef_V_in, Dcoef_U_in, Dcoef_V_in | ||
| 289 | REAL, DIMENSION(klon,klev), INTENT(IN) :: Kcoef_m_in | ||
| 290 | !!! | ||
| 291 | |||
| 292 | ! Output arguments | ||
| 293 | !**************************************************************************************** | ||
| 294 | REAL, DIMENSION(klon,klev), INTENT(OUT) :: flx_u_new, flx_v_new | ||
| 295 | REAL, DIMENSION(klon,klev), INTENT(OUT) :: d_u_new, d_v_new | ||
| 296 | |||
| 297 | ! Local variables | ||
| 298 | !**************************************************************************************** | ||
| 299 | 3840 | REAL, DIMENSION(klon,klev) :: u_new, v_new | |
| 300 | INTEGER :: k, i | ||
| 301 | |||
| 302 | ! Include | ||
| 303 | !**************************************************************************************** | ||
| 304 | INCLUDE "YOMCST.h" | ||
| 305 | INCLUDE "compbl.h" | ||
| 306 | |||
| 307 | ! | ||
| 308 | !**************************************************************************************** | ||
| 309 | |||
| 310 | !!! jyg le 07/02/2012 | ||
| 311 | !!jyg IF (mod(iflag_pbl_split,2) .eq.1) THEN | ||
| 312 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1920 times.
|
1920 | IF (mod(iflag_pbl_split,10) .ge.1) THEN |
| 313 | !!! nrlmd le 02/05/2011 | ||
| 314 | ✗ | DO i = 1, knon | |
| 315 | ✗ | Acoef_U(i)=Acoef_U_in(i) | |
| 316 | ✗ | Acoef_V(i)=Acoef_V_in(i) | |
| 317 | ✗ | Bcoef_U(i)=Bcoef_U_in(i) | |
| 318 | ✗ | Bcoef_V(i)=Bcoef_V_in(i) | |
| 319 | ENDDO | ||
| 320 | ✗ | DO k = 1, klev | |
| 321 | ✗ | DO i = 1, knon | |
| 322 | ✗ | Ccoef_U(i,k)=Ccoef_U_in(i,k) | |
| 323 | ✗ | Ccoef_V(i,k)=Ccoef_V_in(i,k) | |
| 324 | ✗ | Dcoef_U(i,k)=Dcoef_U_in(i,k) | |
| 325 | ✗ | Dcoef_V(i,k)=Dcoef_V_in(i,k) | |
| 326 | ✗ | Kcoefm(i,k)=Kcoef_m_in(i,k) | |
| 327 | ENDDO | ||
| 328 | ENDDO | ||
| 329 | !!! | ||
| 330 | ENDIF ! (mod(iflag_pbl_split,2) .ge.1) | ||
| 331 | !!! | ||
| 332 | |||
| 333 | ! Niveau 1 | ||
| 334 |
2/2✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 788452 times.
|
790372 | DO i = 1, knon |
| 335 | 788452 | u_new(i,1) = Acoef_U(i) + Bcoef_U(i)*flx_u1(i)*dtime | |
| 336 | 790372 | v_new(i,1) = Acoef_V(i) + Bcoef_V(i)*flx_v1(i)*dtime | |
| 337 | END DO | ||
| 338 | |||
| 339 | ! Niveau 2 jusqu'au sommet klev | ||
| 340 |
2/2✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 1920 times.
|
74880 | DO k = 2, klev |
| 341 |
2/2✓ Branch 0 taken 29961176 times.
✓ Branch 1 taken 72960 times.
|
30036056 | DO i=1, knon |
| 342 | 29961176 | u_new(i,k) = Ccoef_U(i,k) + Dcoef_U(i,k) * u_new(i,k-1) | |
| 343 | 30034136 | v_new(i,k) = Ccoef_V(i,k) + Dcoef_V(i,k) * v_new(i,k-1) | |
| 344 | END DO | ||
| 345 | END DO | ||
| 346 | |||
| 347 | !**************************************************************************************** | ||
| 348 | ! Calcul flux | ||
| 349 | ! | ||
| 350 | !== flux_u/v est le flux de moment angulaire (positif vers bas) | ||
| 351 | !== dont l'unite est: (kg m/s)/(m**2 s) | ||
| 352 | ! | ||
| 353 | !**************************************************************************************** | ||
| 354 | ! | ||
| 355 |
4/4✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 74880 times.
✓ Branch 2 taken 74430720 times.
✓ Branch 3 taken 74880 times.
|
74507520 | flx_u_new(:,:) = 0.0 |
| 356 |
4/4✓ Branch 0 taken 74880 times.
✓ Branch 1 taken 1920 times.
✓ Branch 2 taken 74430720 times.
✓ Branch 3 taken 74880 times.
|
74507520 | flx_v_new(:,:) = 0.0 |
| 357 | |||
| 358 |
2/2✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 788452 times.
|
790372 | flx_u_new(1:knon,1)=flx_u1(1:knon) |
| 359 |
2/2✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 788452 times.
|
790372 | flx_v_new(1:knon,1)=flx_v1(1:knon) |
| 360 | |||
| 361 | ! Niveau 2->klev | ||
| 362 |
2/2✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 1920 times.
|
74880 | DO k = 2, klev |
| 363 |
2/2✓ Branch 0 taken 29961176 times.
✓ Branch 1 taken 72960 times.
|
30036056 | DO i = 1, knon |
| 364 | flx_u_new(i,k) = Kcoefm(i,k)/RG/dtime * & | ||
| 365 | 29961176 | (u_new(i,k)-u_new(i,k-1)) | |
| 366 | |||
| 367 | flx_v_new(i,k) = Kcoefm(i,k)/RG/dtime * & | ||
| 368 | 30034136 | (v_new(i,k)-v_new(i,k-1)) | |
| 369 | END DO | ||
| 370 | END DO | ||
| 371 | |||
| 372 | !**************************************************************************************** | ||
| 373 | ! Calcul tendances | ||
| 374 | ! | ||
| 375 | !**************************************************************************************** | ||
| 376 |
4/4✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 74880 times.
✓ Branch 2 taken 74430720 times.
✓ Branch 3 taken 74880 times.
|
74507520 | d_u_new(:,:) = 0.0 |
| 377 |
4/4✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 74880 times.
✓ Branch 2 taken 74430720 times.
✓ Branch 3 taken 74880 times.
|
74507520 | d_v_new(:,:) = 0.0 |
| 378 |
2/2✓ Branch 0 taken 74880 times.
✓ Branch 1 taken 1920 times.
|
76800 | DO k = 1, klev |
| 379 |
2/2✓ Branch 0 taken 30749628 times.
✓ Branch 1 taken 74880 times.
|
30826428 | DO i = 1, knon |
| 380 | 30749628 | d_u_new(i,k) = u_new(i,k) - u_old(i,k) | |
| 381 | 30824508 | d_v_new(i,k) = v_new(i,k) - v_old(i,k) | |
| 382 | END DO | ||
| 383 | END DO | ||
| 384 | |||
| 385 | 1576904 | END SUBROUTINE climb_wind_up | |
| 386 | ! | ||
| 387 | !**************************************************************************************** | ||
| 388 | ! | ||
| 389 | END MODULE climb_wind_mod | ||
| 390 |