- Timestamp:
- Jul 23, 2024, 3:29:36 PM (8 weeks ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_alp.F90
r5101 r5103 136 136 !------------Test sur le LCL des thermiques 137 137 do ig=1,ngrid 138 ok_lcl(ig)=. false.139 if ( (pcon(ig) > pplay(ig,nlay-1)) .and. (pcon(ig) < pplay(ig,1)) ) ok_lcl(ig)=. true.138 ok_lcl(ig)=.FALSE. 139 if ( (pcon(ig) > pplay(ig,nlay-1)) .and. (pcon(ig) < pplay(ig,1)) ) ok_lcl(ig)=.TRUE. 140 140 enddo 141 141 … … 199 199 200 200 !-----Calcul de la TKE transportée par les thermiques : therm_tke_max 201 callthermcell_tke_transport(ngrid,nlay,ptimestep,fm0,entr0, & ! in201 CALL thermcell_tke_transport(ngrid,nlay,ptimestep,fm0,entr0, & ! in 202 202 rg,pplev,therm_tke_max) ! out 203 203 ! print *,' thermcell_tke_transport -> ' !!jyg … … 337 337 !!! fin nrlmd le 10/04/2012 338 338 339 ! print*,'avant calcul ale et alp'339 ! PRINT*,'avant calcul ale et alp' 340 340 !calcul de ALE et ALP pour la convection 341 341 alp_bl(:)=0. 342 342 ale_bl(:)=0. 343 ! print*,'ALE,ALP ,l,zw2(ig,l),ale_bl(ig),alp_bl(ig)'343 ! PRINT*,'ALE,ALP ,l,zw2(ig,l),ale_bl(ig),alp_bl(ig)' 344 344 do l=1,nlay 345 345 do ig=1,ngrid 346 346 alp_bl(ig)=max(alp_bl(ig),0.5*rhobarz(ig,l)*wth3(ig,l) ) 347 347 ale_bl(ig)=max(ale_bl(ig),0.5*zw2(ig,l)**2) 348 ! print*,'ALE,ALP',l,zw2(ig,l),ale_bl(ig),alp_bl(ig)348 ! PRINT*,'ALE,ALP',l,zw2(ig,l),ale_bl(ig),alp_bl(ig) 349 349 enddo 350 350 enddo … … 377 377 enddo 378 378 379 ! print*,'apres wght_th'379 ! PRINT*,'apres wght_th' 380 380 !test pour prolonger la convection 381 381 do ig=1,ngrid … … 384 384 lalim_conv(ig)=1 385 385 wght_th(ig,1)=1. 386 ! print*,'lalim_conv ok',lalim_conv(ig),wght_th(ig,1)386 ! PRINT*,'lalim_conv ok',lalim_conv(ig),wght_th(ig,1) 387 387 endif 388 388 enddo
Note: See TracChangeset
for help on using the changeset viewer.