Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/phylmdiso/lmdz_wake.F90
- Timestamp:
- Jul 22, 2024, 9:29:09 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmdiso/lmdz_wake.F90
r5093 r5099 159 159 #endif 160 160 161 !162 161 ! Input/Output 163 162 ! State variables … … 425 424 CALL getin_p('wdens_ref_l',wdens_ref(2)) !wake number per unit area ; land 426 425 !>jyg 427 ! 426 428 427 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 429 428 !!!!!!!!! Population dynamics parameters !!!!!!!!!!!!!!!!!!!!!!!!!!!! … … 441 440 CALL getin_p('rzero_wk', rzero) 442 441 aa0 = 3.14*rzero*rzero 443 ! 442 444 443 tau_cv = 4000. 445 444 CALL getin_p('tau_cv', tau_cv) … … 574 573 END DO 575 574 ENDIF ! (iflag_wk_pop_dyn >=1) 576 ! 575 577 576 DO i = 1, klon 578 577 ! c sigmaw(i) = amax1(sigmaw(i),sigd_con(i)) … … 586 585 END DO 587 586 588 !589 587 IF (iflag_wk_pop_dyn >= 1) THEN 590 588 awdens_in(:) = awdens(:) … … 599 597 d_sigmaw(:) = 0. 600 598 ktopw(:) = 0 601 ! 599 602 600 !<jyg 603 601 dth(:,:) = 0. … … 617 615 cstar(:) = 0. 618 616 ktopw(:) = 0 619 ! 617 620 618 ! Vertical advection local variables 621 619 omgbw(:,:) = 0. … … 629 627 #endif 630 628 !>jyg 631 ! 629 632 630 IF (prt_level>=10) THEN 633 631 PRINT *, 'wake-1, sigmaw(igout) ', sigmaw(igout) … … 1644 1642 ! print*,'dtKE= ',dtKE(i,k),' dqKE= ',dqKE(i,k) 1645 1643 1646 !1647 1648 1644 ! cc nrlmd Prise en compte du taux de mortalité 1649 1645 ! cc Définitions de entr, detr … … 1834 1830 END DO 1835 1831 END DO 1836 ! 1832 1837 1833 DO i = 1, klon 1838 1834 IF (wk_adv(i)) THEN … … 1879 1875 d_wdens2(i) = d_wdens2(i) + wdens_targ - wdens(i) 1880 1876 wdens(i) = wdens_targ 1881 ! 1877 1882 1878 wdens_targ = min( max(awdens(i),0.), wdens(i) ) 1883 1879 d_awdens2(i) = d_awdens2(i) + wdens_targ - awdens(i) … … 2279 2275 sum_dtdwn(i) = sum_dtdwn(i) + dtdwn(i, k)*dz(i) 2280 2276 sum_dqdwn(i) = sum_dqdwn(i) + dqdwn(i, k)*dz(i) 2281 ! 2277 2282 2278 dthmin(i) = min(dthmin(i), dth(i,k)) 2283 2279 END IF
Note: See TracChangeset
for help on using the changeset viewer.