Changeset 3411 for LMDZ6/branches/DYNAMICO-conv/libf/phylmd/calwake.F90
- Timestamp:
- Nov 5, 2018, 3:24:59 PM (6 years ago)
- Location:
- LMDZ6/branches/DYNAMICO-conv
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/DYNAMICO-conv
- Property svn:mergeinfo changed
/LMDZ6/trunk removed
- Property svn:mergeinfo changed
-
LMDZ6/branches/DYNAMICO-conv/libf/phylmd/calwake.F90
r3356 r3411 4 4 SUBROUTINE calwake(iflag_wake_tend, paprs, pplay, dtime, & 5 5 t, q, omgb, & 6 dt_dwn, dq_dwn, m_dwn, m_up, dt_a, dq_a, wgen,&7 sigd, Cin,&8 wake_deltat, wake_deltaq, wake_s, awake_dens,wake_dens, &6 dt_dwn, dq_dwn, m_dwn, m_up, dt_a, dq_a, & 7 sigd, & 8 wake_deltat, wake_deltaq, wake_s, wake_dens, & 9 9 wake_dth, wake_h, & 10 10 wake_pe, wake_fip, wake_gfl, & … … 14 14 wake_omg, wake_dp_deltomg, & 15 15 wake_spread, wake_cstar, wake_d_deltat_gw, & 16 wake_ddeltat, wake_ddeltaq, wake_ds, awake_ddens,wake_ddens)16 wake_ddeltat, wake_ddeltaq, wake_ds, wake_ddens) 17 17 ! ************************************************************** 18 18 ! * … … 45 45 REAL, DIMENSION(klon, klev), INTENT (IN) :: m_up, m_dwn 46 46 REAL, DIMENSION(klon, klev), INTENT (IN) :: dt_a, dq_a 47 REAL, DIMENSION(klon), INTENT (IN) :: wgen48 47 REAL, DIMENSION(klon), INTENT (IN) :: sigd 49 REAL, DIMENSION(klon), INTENT (IN) :: Cin50 48 ! Input/Output 51 49 ! ------------ 52 50 REAL, DIMENSION(klon, klev), INTENT (INOUT) :: wake_deltat, wake_deltaq 53 51 REAL, DIMENSION(klon), INTENT (INOUT) :: wake_s 54 REAL, DIMENSION(klon), INTENT (INOUT) :: awake_dens,wake_dens52 REAL, DIMENSION(klon), INTENT (INOUT) :: wake_dens 55 53 ! Output 56 54 ! ------ … … 69 67 REAL, DIMENSION(klon), INTENT (OUT) :: wake_cstar 70 68 REAL, DIMENSION(klon, klev), INTENT (OUT) :: wake_ddeltat, wake_ddeltaq 71 REAL, DIMENSION(klon), INTENT (OUT) :: wake_ds, awake_ddens,wake_ddens69 REAL, DIMENSION(klon), INTENT (OUT) :: wake_ds, wake_ddens 72 70 73 71 … … 90 88 REAL, DIMENSION(klon, klev) :: tx, qx 91 89 REAL, DIMENSION(klon) :: hw, wape, fip, gfl 92 REAL, DIMENSION(klon) :: sigmaw, awdens,wdens90 REAL, DIMENSION(klon) :: sigmaw, wdens 93 91 REAL, DIMENSION(klon, klev) :: omgbdth 94 92 REAL, DIMENSION(klon, klev) :: dp_omgb … … 101 99 REAL, DIMENSION(klon, klev) :: d_deltat_gw 102 100 REAL, DIMENSION(klon, klev) :: d_deltatw, d_deltaqw 103 REAL, DIMENSION(klon) :: d_sigmaw, d_ awdens, d_wdens101 REAL, DIMENSION(klon) :: d_sigmaw, d_wdens 104 102 105 103 REAL :: rdcp … … 107 105 108 106 IF (prt_level >= 10) THEN 109 print *, '-> calwake, wake_s , wgen input ', wake_s(1), wgen(1)107 print *, '-> calwake, wake_s input ', wake_s(1) 110 108 ENDIF 111 109 … … 149 147 d_deltaqw(:,:) = 0. 150 148 d_sigmaw(:) = 0. 151 d_awdens(:) = 0.152 149 d_wdens(:) = 0. 153 150 ! … … 182 179 183 180 DO i = 1, klon 184 awdens(i) = max(0., awake_dens(i))185 181 wdens(i) = max(0., wake_dens(i)) 186 182 END DO … … 210 206 CALL wake(znatsurf, p, ph, pi, dtime, & 211 207 te, qe, omgbe, & 212 dtdwn, dqdwn, amdwn, amup, dta, dqa, wgen,&213 sigd0, Cin,&214 dtw, dqw, sigmaw, awdens,wdens, & ! state variables208 dtdwn, dqdwn, amdwn, amup, dta, dqa, & 209 sigd0, & 210 dtw, dqw, sigmaw, wdens, & ! state variables 215 211 dth, hw, wape, fip, gfl, & 216 212 dtls, dqls, ktopw, omgbdth, dp_omgb, tx, qx, & 217 213 dtke, dqke, omg, dp_deltomg, spread, cstar, & 218 214 d_deltat_gw, & 219 d_deltatw, d_deltaqw, d_sigmaw, d_ awdens, d_wdens) ! tendencies215 d_deltatw, d_deltaqw, d_sigmaw, d_wdens) ! tendencies 220 216 221 217 ! … … 278 274 IF (ktopw(i)>0) THEN 279 275 wake_ds(i) = d_sigmaw(i)*dtime 280 awake_ddens(i) = d_awdens(i)*dtime281 276 wake_ddens(i) = d_wdens(i)*dtime 282 277 ELSE … … 303 298 DO i = 1, klon 304 299 wake_s(i) = sigmaw(i) 305 awake_dens(i) = awdens(i)306 300 wake_dens(i) = wdens(i) 307 301 END DO
Note: See TracChangeset
for help on using the changeset viewer.