Changeset 5082 for LMDZ6/branches/Amaury_dev/libf/phylmd/freinage.F90
- Timestamp:
- Jul 19, 2024, 5:41:58 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/freinage.F90
r3319 r5082 97 97 98 98 ! if ( (height(i,jv) .gt. yzlay(i,k)) .AND. (height(i,jv) .gt. 0.1) .and. LAI(i,jv).gt.0. ) then 99 if ( (height(i,jv) .gt. yzlay(i,k)) .AND. (height(i,jv) .gt. 0.1) ) then99 if ( (height(i,jv) > yzlay(i,k)) .AND. (height(i,jv) > 0.1) ) then 100 100 !FC attention veut on le test sur le LAI ? 101 if (ifl_pbltree .eq.1) then101 if (ifl_pbltree==1) then 102 102 drag_pro(i,k)= drag_pro(i,k)+ & 103 103 veget(i,jv) 104 elseif (ifl_pbltree .eq.2) then104 elseif (ifl_pbltree==2) then 105 105 drag_pro(i,k)= drag_pro(i,k)+ & 106 106 6*LAI(i,jv)*veget(i,jv)*( yzlay(i,k)*(height(i,jv)-yzlay(i,k))/(height(i,jv)*height(i,jv)+ 0.01)) 107 elseif (ifl_pbltree .eq.3) then107 elseif (ifl_pbltree==3) then 108 108 drag_pro(i,k)= drag_pro(i,k)+ & 109 109 veget(i,jv)*( yzlay(i,k)*(height(i,jv)-yzlay(i,k))/(height(i,jv)*height(i,jv)+ 0.01)) 110 elseif (ifl_pbltree .eq.0) then110 elseif (ifl_pbltree==0) then 111 111 drag_pro(i,k)=0.0 112 112 endif
Note: See TracChangeset
for help on using the changeset viewer.