Changeset 3162 for trunk/LMDZ.MARS/libf
- Timestamp:
- Dec 19, 2023, 4:17:44 PM (12 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F
r3153 r3162 592 592 593 593 c----------------------------------------------------------------------- 594 c Using the wind modified by friction for lifting and sublimation 595 c ---------------------------------------------------------------- 596 597 ! This is computed above and takes into account surface-atmosphere flux 598 ! enhancement by subgrid gustiness and atmospheric-stability related 599 ! variations of transfer coefficients. 600 ! Calculate Cd again with wind slowed by friction 601 c ------------------------------------------- 602 603 CALL vdif_cd(ngrid,nlay,pz0,g,pzlay,zu,zv,wstar,ptsrf_tmp 604 & ,zh,zcdv_true,zcdh_true) 605 606 zu2(:)=zu(:,1)*zu(:,1)+zv(:,1)*zv(:,1) 607 608 IF (callrichsl) THEN 609 zcdv(:)=zcdv_true(:)*sqrt(zu2(:)+ 610 & (log(1.+0.7*wstar(:) + 2.3*wstar(:)**2))**2) 611 zcdh(:)=zcdh_true(:)*sqrt(zu2(:)+ 612 & (log(1.+0.7*wstar(:) + 2.3*wstar(:)**2))**2) 613 614 ustar(:)=sqrt(zcdv_true(:))*sqrt(zu2(:)+ 615 & (log(1.+0.7*wstar(:) + 2.3*wstar(:)**2))**2) 616 617 ELSE 618 zcdv(:)=zcdv_true(:)*sqrt(zu2(:)) ! 1 / bulk aerodynamic momentum conductance 619 zcdh(:)=zcdh_true(:)*sqrt(zu2(:)) ! 1 / bulk aerodynamic heat conductance 620 ustar(:)=sqrt(zcdv_true(:))*sqrt(zu2(:)) 621 ENDIF 622 623 624 c----------------------------------------------------------------------- 594 625 c 6. inversion pour l'implicite sur h sans oublier le couplage 595 626 c avec le sol (conduction) … … 776 807 c TRACERS 777 808 c ------- 778 779 c Using the wind modified by friction for lifting and sublimation780 c ----------------------------------------------------------------781 782 ! This is computed above and takes into account surface-atmosphere flux783 ! enhancement by subgrid gustiness and atmospheric-stability related784 ! variations of transfer coefficients.785 786 ! DO ig=1,ngrid787 ! zu2(ig)=zu(ig,1)*zu(ig,1)+zv(ig,1)*zv(ig,1)788 ! zcdv(ig)=zcdv_true(ig)*sqrt(zu2(ig))789 ! zcdh(ig)=zcdh_true(ig)*sqrt(zu2(ig))790 ! ENDDO791 792 809 c Calcul du flux vertical au bas de la premiere couche (dust) : 793 810 c -----------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.