- Timestamp:
- Jul 24, 2024, 2:54:37 PM (8 weeks ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_closure.F90
r5103 r5116 53 53 DO k=1,llmax-1 54 54 do ig=1,ngrid 55 if (k<lalim(ig)) then55 if (k<lalim(ig)) THEN 56 56 alim_star2(ig)=alim_star2(ig)+alim_star(ig,k)**2 & 57 57 & /(rho(ig,k)*(zlev(ig,k+1)-zlev(ig,k))) … … 63 63 64 64 DO ig=1,ngrid 65 if (alim_star2(ig)>1.e-10) then65 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))
Note: See TracChangeset
for help on using the changeset viewer.