Changeset 5082 for LMDZ6/branches/Amaury_dev/libf/phylmd/yamada4.F90
- Timestamp:
- Jul 19, 2024, 5:41:58 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/yamada4.F90
r4884 r5082 1038 1038 hlim(1:ngrid)=0. 1039 1039 1040 IF (nsrf .EQ.1) THEN1040 IF (nsrf == 1) THEN 1041 1041 1042 1042 ! coefficients … … 1076 1076 DO ig=1,ngrid 1077 1077 winds(ig,k)=sqrt(u(ig,k)**2+v(ig,k)**2) 1078 IF (zlev(ig,k) .LE.h_oro(ig)) THEN ! sous l'orographie1078 IF (zlev(ig,k) <= h_oro(ig)) THEN ! sous l'orographie 1079 1079 l2strat= kapb*pbl_lmixmin_alpha*winds(ig,k)/sqrt(max(n2(ig,k),1.E-10)) ! si stratifi??, amplitude d'oscillation * kappab (voir Van de Wiel et al 2008) 1080 1080 l2neutre=kap*zlev(ig,k)*h_oro(ig)/(kap*zlev(ig,k)+h_oro(ig)) ! Dans le cas neutre, formule de blackadar. tend asymptotiquement vers h … … 1083 1083 l2(ig,k)=l2limit(ig) 1084 1084 1085 ELSE IF (zlev(ig,k) .LE.hlim(ig)) THEN ! Si on est au dessus des montagnes, mais affect?? encore par elles1085 ELSE IF (zlev(ig,k) <= hlim(ig)) THEN ! Si on est au dessus des montagnes, mais affect?? encore par elles 1086 1086 1087 1087 ! Au dessus des montagnes, on prend la l2limit au sommet des montagnes
Note: See TracChangeset
for help on using the changeset viewer.