Changeset 1894 for palm/trunk
- Timestamp:
- Apr 27, 2016 9:01:48 AM (9 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/pmc_interface_mod.f90 (modified) (13 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/pmc_interface_mod.f90
r1893 r1894 20 20 ! Current revisions: 21 21 ! ------------------ 22 ! 22 ! bugfix: pt interpolations are omitted in case that the temperature equation is 23 ! switched off 23 24 ! 24 25 ! Former revisions: … … 2618 2619 CALL pmci_interp_tril_all ( e, ec, ico, jco, kco, r1xo, r2xo, r1yo, & 2619 2620 r2yo, r1zo, r2zo, nzb_s_inner, 'e' ) 2620 CALL pmci_interp_tril_all ( pt, ptc, ico, jco, kco, r1xo, r2xo, r1yo, & 2621 r2yo, r1zo, r2zo, nzb_s_inner, 's' ) 2621 IF ( .NOT. neutral ) THEN 2622 CALL pmci_interp_tril_all ( pt, ptc, ico, jco, kco, r1xo, r2xo, & 2623 r1yo, r2yo, r1zo, r2zo, nzb_s_inner, 's' ) 2624 ENDIF 2622 2625 IF ( humidity .OR. passive_scalar ) THEN 2623 2626 CALL pmci_interp_tril_all ( q, qc, ico, jco, kco, r1xo, r2xo, r1yo, & … … 3255 3258 logc_ratio_u_l, nzt_topo_nestbc_l, 'l', & 3256 3259 'e' ) 3257 CALL pmci_interp_tril_lr( pt, ptc, ico, jco, kco, r1xo, r2xo, r1yo, & 3258 r2yo, r1zo, r2zo, nzb_s_inner, logc_u_l, & 3259 logc_ratio_u_l, nzt_topo_nestbc_l, 'l', & 3260 's' ) 3260 IF ( .NOT. neutral ) THEN 3261 CALL pmci_interp_tril_lr( pt, ptc, ico, jco, kco, r1xo, r2xo, & 3262 r1yo, r2yo, r1zo, r2zo, nzb_s_inner, & 3263 logc_u_l, logc_ratio_u_l, & 3264 nzt_topo_nestbc_l, 'l', 's' ) 3265 ENDIF 3261 3266 IF ( humidity .OR. passive_scalar ) THEN 3262 3267 CALL pmci_interp_tril_lr( q, qc, ico, jco, kco, r1xo, r2xo, r1yo, & … … 3271 3276 CALL pmci_extrap_ifoutflow_lr( w, nzb_w_inner, 'l', 'w' ) 3272 3277 CALL pmci_extrap_ifoutflow_lr( e, nzb_s_inner, 'l', 'e' ) 3273 CALL pmci_extrap_ifoutflow_lr( pt,nzb_s_inner, 'l', 's' ) 3278 IF ( .NOT. neutral ) THEN 3279 CALL pmci_extrap_ifoutflow_lr( pt,nzb_s_inner, 'l', 's' ) 3280 ENDIF 3274 3281 IF ( humidity .OR. passive_scalar ) THEN 3275 3282 CALL pmci_extrap_ifoutflow_lr( q, nzb_s_inner, 'l', 's' ) … … 3297 3304 logc_ratio_u_r, nzt_topo_nestbc_r, 'r', & 3298 3305 'e' ) 3299 CALL pmci_interp_tril_lr( pt, ptc, ico, jco, kco, r1xo, r2xo, r1yo, & 3300 r2yo, r1zo, r2zo, nzb_s_inner, logc_u_r, & 3301 logc_ratio_u_r, nzt_topo_nestbc_r, 'r', & 3302 's' ) 3306 IF ( .NOT. neutral ) THEN 3307 CALL pmci_interp_tril_lr( pt, ptc, ico, jco, kco, r1xo, r2xo, & 3308 r1yo, r2yo, r1zo, r2zo, nzb_s_inner, & 3309 logc_u_r, logc_ratio_u_r, & 3310 nzt_topo_nestbc_r, 'r', 's' ) 3311 ENDIF 3303 3312 IF ( humidity .OR. passive_scalar ) THEN 3304 3313 CALL pmci_interp_tril_lr( q, qc, ico, jco, kco, r1xo, r2xo, r1yo, & … … 3313 3322 CALL pmci_extrap_ifoutflow_lr( w, nzb_w_inner, 'r', 'w' ) 3314 3323 CALL pmci_extrap_ifoutflow_lr( e, nzb_s_inner, 'r', 'e' ) 3315 CALL pmci_extrap_ifoutflow_lr( pt,nzb_s_inner, 'r', 's' ) 3324 IF ( .NOT. neutral ) THEN 3325 CALL pmci_extrap_ifoutflow_lr( pt,nzb_s_inner, 'r', 's' ) 3326 ENDIF 3316 3327 IF ( humidity .OR. passive_scalar ) THEN 3317 3328 CALL pmci_extrap_ifoutflow_lr( q, nzb_s_inner, 'r', 's' ) … … 3339 3350 logc_ratio_u_s, nzt_topo_nestbc_s, 's', & 3340 3351 'e' ) 3341 CALL pmci_interp_tril_sn( pt, ptc, ico, jco, kco, r1xo, r2xo, r1yo, & 3342 r2yo, r1zo, r2zo, nzb_s_inner, logc_u_s, & 3343 logc_ratio_u_s, nzt_topo_nestbc_s, 's', & 3344 's' ) 3352 IF ( .NOT. neutral ) THEN 3353 CALL pmci_interp_tril_sn( pt, ptc, ico, jco, kco, r1xo, r2xo, & 3354 r1yo, r2yo, r1zo, r2zo, nzb_s_inner, & 3355 logc_u_s, logc_ratio_u_s, & 3356 nzt_topo_nestbc_s, 's', 's' ) 3357 ENDIF 3345 3358 IF ( humidity .OR. passive_scalar ) THEN 3346 3359 CALL pmci_interp_tril_sn( q, qc, ico, jco, kco, r1xo, r2xo, r1yo, & … … 3355 3368 CALL pmci_extrap_ifoutflow_sn( w, nzb_w_inner, 's', 'w' ) 3356 3369 CALL pmci_extrap_ifoutflow_sn( e, nzb_s_inner, 's', 'e' ) 3357 CALL pmci_extrap_ifoutflow_sn( pt,nzb_s_inner, 's', 's' ) 3370 IF ( .NOT. neutral ) THEN 3371 CALL pmci_extrap_ifoutflow_sn( pt,nzb_s_inner, 's', 's' ) 3372 ENDIF 3358 3373 IF ( humidity .OR. passive_scalar ) THEN 3359 3374 CALL pmci_extrap_ifoutflow_sn( q, nzb_s_inner, 's', 's' ) … … 3381 3396 logc_ratio_u_n, nzt_topo_nestbc_n, 'n', & 3382 3397 'e' ) 3383 CALL pmci_interp_tril_sn( pt, ptc, ico, jco, kco, r1xo, r2xo, r1yo, & 3384 r2yo, r1zo, r2zo, nzb_s_inner, logc_u_n, & 3385 logc_ratio_u_n, nzt_topo_nestbc_n, 'n', & 3386 's' ) 3398 IF ( .NOT. neutral ) THEN 3399 CALL pmci_interp_tril_sn( pt, ptc, ico, jco, kco, r1xo, r2xo, & 3400 r1yo, r2yo, r1zo, r2zo, nzb_s_inner, & 3401 logc_u_n, logc_ratio_u_n, & 3402 nzt_topo_nestbc_n, 'n', 's' ) 3403 ENDIF 3387 3404 IF ( humidity .OR. passive_scalar ) THEN 3388 3405 CALL pmci_interp_tril_sn( q, qc, ico, jco, kco, r1xo, r2xo, r1yo, & … … 3397 3414 CALL pmci_extrap_ifoutflow_sn( w, nzb_w_inner, 'n', 'w' ) 3398 3415 CALL pmci_extrap_ifoutflow_sn( e, nzb_s_inner, 'n', 'e' ) 3399 CALL pmci_extrap_ifoutflow_sn( pt,nzb_s_inner, 'n', 's' ) 3416 IF ( .NOT. neutral ) THEN 3417 CALL pmci_extrap_ifoutflow_sn( pt,nzb_s_inner, 'n', 's' ) 3418 ENDIF 3400 3419 IF ( humidity .OR. passive_scalar ) THEN 3401 3420 CALL pmci_extrap_ifoutflow_sn( q, nzb_s_inner, 'n', 's' ) … … 3416 3435 CALL pmci_interp_tril_t( e, ec, ico, jco, kco, r1xo, r2xo, r1yo, & 3417 3436 r2yo, r1zo, r2zo, 'e' ) 3418 CALL pmci_interp_tril_t( pt, ptc, ico, jco, kco, r1xo, r2xo, r1yo, & 3419 r2yo, r1zo, r2zo, 's' ) 3437 IF ( .NOT. neutral ) THEN 3438 CALL pmci_interp_tril_t( pt, ptc, ico, jco, kco, r1xo, r2xo, r1yo, & 3439 r2yo, r1zo, r2zo, 's' ) 3440 ENDIF 3420 3441 IF ( humidity .OR. passive_scalar ) THEN 3421 3442 CALL pmci_interp_tril_t( q, qc, ico, jco, kco, r1xo, r2xo, r1yo, & … … 3428 3449 CALL pmci_extrap_ifoutflow_t( w, 'w' ) 3429 3450 CALL pmci_extrap_ifoutflow_t( e, 'e' ) 3430 CALL pmci_extrap_ifoutflow_t( pt, 's' ) 3451 IF ( .NOT. neutral ) THEN 3452 CALL pmci_extrap_ifoutflow_t( pt, 's' ) 3453 ENDIF 3431 3454 IF ( humidity .OR. passive_scalar ) THEN 3432 3455 CALL pmci_extrap_ifoutflow_t( q, 's' ) … … 3451 3474 CALL pmci_anterp_tophat( w, wc, kctw, iflo, ifuo, jflo, jfuo, kflw, & 3452 3475 kfuw, ijfc_s, 'w' ) 3453 CALL pmci_anterp_tophat( pt, ptc, kctu, iflo, ifuo, jflo, jfuo, kflo, & 3454 kfuo, ijfc_s, 's' ) 3476 IF ( .NOT. neutral ) THEN 3477 CALL pmci_anterp_tophat( pt, ptc, kctu, iflo, ifuo, jflo, jfuo, kflo, & 3478 kfuo, ijfc_s, 's' ) 3479 ENDIF 3455 3480 IF ( humidity .OR. passive_scalar ) THEN 3456 3481 CALL pmci_anterp_tophat( q, qc, kctu, iflo, ifuo, jflo, jfuo, kflo, &
Note: See TracChangeset
for help on using the changeset viewer.
