- Timestamp:
- Jul 19, 2024, 5:41:58 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_plume.F90
r4590 r5082 255 255 zdz2=zlev(ig,lt)-zlev(ig,l) 256 256 257 do while (lmel .gt.zdz2)257 do while (lmel>zdz2) 258 258 lt=lt+1 259 259 zlt=zlev(ig,lt) … … 321 321 ! En dessous de lalim, on prend le max de alim_star et entr_star pour 322 322 ! alim_star et 0 sinon 323 if (l .lt.lalim(ig)) then323 if (l<lalim(ig)) then 324 324 alim_star(ig,l)=max(alim_star(ig,l),entr_star(ig,l)) 325 325 entr_star(ig,l)=0. … … 375 375 enddo 376 376 377 if (prt_level .ge.20) print*,'coucou calcul detr 460: ig, l',ig, l377 if (prt_level>=20) print*,'coucou calcul detr 460: ig, l',ig, l 378 378 ! 379 379 !=========================================================================== … … 383 383 nbpb=0 384 384 do ig=1,ngrid 385 if (zw2(ig,l+1)>0. .and. zw2(ig,l+1) .lt.1.e-10) then385 if (zw2(ig,l+1)>0. .and. zw2(ig,l+1)<1.e-10) then 386 386 ! stop'On tombe sur le cas particulier de thermcell_dry' 387 387 ! print*,'On tombe sur le cas particulier de thermcell_plume' … … 391 391 endif 392 392 393 if (zw2(ig,l+1) .lt.0.) then393 if (zw2(ig,l+1)<0.) then 394 394 linter(ig)=(l*(zw2(ig,l+1)-zw2(ig,l)) & 395 395 & -zw2(ig,l))/(zw2(ig,l+1)-zw2(ig,l)) 396 396 zw2(ig,l+1)=0. 397 397 !+CR:04/05/12:correction calcul linter pour calcul de zmax continu 398 elseif (f_star(ig,l+1) .lt.0.) then398 elseif (f_star(ig,l+1)<0.) then 399 399 linter(ig)=(l*(f_star(ig,l+1)-f_star(ig,l)) & 400 400 & -f_star(ig,l))/(f_star(ig,l+1)-f_star(ig,l)) … … 405 405 wa_moy(ig,l+1)=sqrt(zw2(ig,l+1)) 406 406 407 if (wa_moy(ig,l+1) .gt.wmaxa(ig)) then407 if (wa_moy(ig,l+1)>wmaxa(ig)) then 408 408 ! lmix est le niveau de la couche ou w (wa_moy) est maximum 409 409 !on rajoute le calcul de lmix_bis 410 if (zqla(ig,l) .lt.1.e-10) then410 if (zqla(ig,l)<1.e-10) then 411 411 lmix_bis(ig)=l+1 412 412 endif … … 436 436 437 437 438 if (prt_level .ge.20) print*,'coucou calcul detr 470: ig, l', ig, l438 if (prt_level>=20) print*,'coucou calcul detr 470: ig, l', ig, l 439 439 440 440 #undef wrgrads_thermcell
Note: See TracChangeset
for help on using the changeset viewer.