Changeset 2897 for LMDZ5/trunk/libf/phylmd/orografi_strato.F90
- Timestamp:
- May 31, 2017, 12:34:09 AM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/orografi_strato.F90
r2357 r2897 1 SUBROUTINE drag_noro_strato( nlon, nlev, dtime, paprs, pplay, pmea, pstd, &1 SUBROUTINE drag_noro_strato(partdrag, nlon, nlev, dtime, paprs, pplay, pmea, pstd, & 2 2 psig, pgam, pthe, ppic, pval, kgwd, kdx, ktest, t, u, v, pulow, pvlow, & 3 3 pustr, pvstr, d_t, d_u, d_v) … … 15 15 ! Explicit Arguments: 16 16 ! ================== 17 ! partdrag-input-I-control which part of the drag we consider (total part or GW part) 17 18 ! nlon----input-I-Total number of horizontal points that get into physics 18 19 ! nlev----input-I-Number of vertical levels … … 66 67 ! ARGUMENTS 67 68 68 INTEGER nlon, nlev69 INTEGER partdrag,nlon, nlev 69 70 REAL dtime 70 71 REAL paprs(nlon, nlev+1) … … 134 135 ! CALL SSO DRAG ROUTINES 135 136 136 CALL orodrag_strato( klon, klev, kgwd, kdx, ktest, dtime, papmh, papmf, &137 CALL orodrag_strato(partdrag,klon, klev, kgwd, kdx, ktest, dtime, papmh, papmf, & 137 138 zgeom, pt, pu, pv, pmea, pstd, psig, pgam, pthe, ppic, pval, pulow, & 138 139 pvlow, pdudt, pdvdt, pdtdt) … … 153 154 END SUBROUTINE drag_noro_strato 154 155 155 SUBROUTINE orodrag_strato( nlon, nlev, kgwd, kdx, ktest, ptsphy, paphm1, &156 SUBROUTINE orodrag_strato(partdrag,nlon, nlev, kgwd, kdx, ktest, ptsphy, paphm1, & 156 157 papm1, pgeom1, ptm1, pum1, pvm1, pmea, pstd, psig, pgam, pthe, ppic, pval & 157 158 ! outputs … … 182 183 ! -------------------- 183 184 ! ==== inputs === 185 ! partdrag-input-I-control which part of the drag we consider (total part or GW part) 184 186 ! nlon----input-I-Total number of horizontal points that get into physics 185 187 ! nlev----input-I-Number of vertical levels … … 201 203 ! pval----input-R-SSO Valleys elevation (m) 202 204 203 INTEGER nlon, nlev, kgwd205 INTEGER nlon, nlev, kgwd 204 206 REAL ptsphy 205 207 … … 239 241 include "YOMCST.h" 240 242 include "YOEGWD.h" 243 241 244 ! ----------------------------------------------------------------------- 242 245 … … 244 247 ! --------- 245 248 246 249 INTEGER partdrag 247 250 REAL pte(nlon, nlev), pvol(nlon, nlev), pvom(nlon, nlev), pulow(nlon), & 248 251 pvlow(nlon) … … 269 272 270 273 INTEGER jl, jk, ji 271 REAL ztmst, zdelp, ztemp, zforc, ztend, rover 274 REAL ztmst, zdelp, ztemp, zforc, ztend, rover, facpart 272 275 REAL zb, zc, zconb, zabsv, zzd1, ratio, zbet, zust, zvst, zdis 273 276 … … 393 396 ! ----------------- 394 397 398 IF (partdrag .GE. 2) THEN 399 facpart=0. 400 ELSE 401 facpart=gkwake 402 ENDIF 403 404 395 405 IF (jk>ikenvh(ji)) THEN 396 406 zb = 1.0 - 0.18*pgam(ji) - 0.04*pgam(ji)**2 397 407 zc = 0.48*pgam(ji) + 0.3*pgam(ji)**2 398 zconb = 2.*ztmst* gkwake*psig(ji)/(4.*pstd(ji))408 zconb = 2.*ztmst*facpart*psig(ji)/(4.*pstd(ji)) 399 409 zabsv = sqrt(pum1(ji,jk)**2+pvm1(ji,jk)**2)/2. 400 410 zzd1 = zb*cos(zpsi(ji,jk))**2 + zc*sin(zpsi(ji,jk))**2 … … 1894 1904 RETURN 1895 1905 END SUBROUTINE sugwd_strato 1896
Note: See TracChangeset
for help on using the changeset viewer.