Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/phylmd/freinage.F90
- Timestamp:
- Jul 24, 2024, 2:54:37 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/freinage.F90
r5105 r5116 6 6 7 7 !ONLINE: 8 use dimphy, only: klon, klev8 use dimphy, ONLY: klon, klev 9 9 ! USE control, ONLY: nvm 10 10 ! USE indice_sol_mod, ONLY: nvm_orch … … 96 96 sumveg(i,k)= sumveg(i,k)+ veget(i,jv) 97 97 98 ! if ( (height(i,jv) .gt. yzlay(i,k)) .AND. (height(i,jv) .gt. 0.1) .and. LAI(i,jv).gt.0. ) then99 if ( (height(i,jv) > yzlay(i,k)) .AND. (height(i,jv) > 0.1) ) then98 ! 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) > 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==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==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==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==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.