- 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_closure.F90
r5099 r5103 36 36 37 37 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 38 ! print*,'THERMCELL CLOSURE 26E'38 !PRINT*,'THERMCELL CLOSURE 26E' 39 39 40 40 alim_star2(:)=0. … … 44 44 ! Indice vertical max (max de lalim) atteint par les thermiques sur le domaine 45 45 llmax=1 46 doig=1,ngrid46 DO ig=1,ngrid 47 47 if (lalim(ig)>llmax) llmax=lalim(ig) 48 enddo 48 END DO 49 49 50 50 51 51 ! Calcul des integrales sur la verticale de alim_star et de 52 52 ! alim_star^2/(rho dz) 53 dok=1,llmax-153 DO k=1,llmax-1 54 54 do ig=1,ngrid 55 55 if (k<lalim(ig)) then … … 59 59 endif 60 60 enddo 61 enddo 61 END DO 62 62 63 63 64 doig=1,ngrid64 DO ig=1,ngrid 65 65 if (alim_star2(ig)>1.e-10) then 66 66 f(ig)=wmax(ig)*alim_star_tot(ig)/ & 67 67 & (max(500.,zmax(ig))*r_aspect*alim_star2(ig)) 68 68 endif 69 enddo 69 END DO 70 70 71 71
Note: See TracChangeset
for help on using the changeset viewer.