- Timestamp:
- Jul 19, 2024, 5:41:58 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_down.F90
r4590 r5082 74 74 do ig=1,ngrid 75 75 !if ( lmax(ig) > nlay - 2 ) stop "les thermiques montent trop haut" 76 if (ilay .le.lmax(ig) .and. lmax(ig)>1 ) then76 if (ilay<=lmax(ig) .and. lmax(ig)>1 ) then 77 77 fdn(ig,ilay)=fdn(ig,ilay+1)+edn(ig,ilay)-ddn(ig,ilay) 78 78 if ( fdn(ig,ilay)+ddn(ig,ilay) > 0. ) then … … 90 90 do ilay=1,nlay,1 91 91 do ig=1,ngrid 92 if (ilay .lt.lmax(ig) .and. lmax(ig)>1) then92 if (ilay<lmax(ig) .and. lmax(ig)>1) then 93 93 fup(ig,ilay+1)=fup(ig,ilay)+eup(ig,ilay)-dup(ig,ilay) 94 94 if (fup(ig,ilay+1)+dup(ig,ilay) > 0.) then … … 131 131 !!!!ATTENTION HYPOTHESE de FLUX COMPENSATOIRE DESCENDANT ET DONC comme schema amont on va chercher trac au dessus!!!!! 132 132 !!!! tentative de prise en compte d'un flux compensatoire montant !!!! 133 if (fup(ig,ilay)-fdn(ig,ilay) .lt.0.) then133 if (fup(ig,ilay)-fdn(ig,ilay) < 0.) then 134 134 call abort_physic("thermcell_updown_dq", 'flux compensatoire '& 135 135 // 'montant, cas non traite par thermcell_updown_dq', 1) … … 164 164 do ilay=2,nlay,1 165 165 do ig=1,ngrid 166 if (fup(ig,ilay)-fdn(ig,ilay) .lt.0.) then166 if (fup(ig,ilay)-fdn(ig,ilay) < 0.) then 167 167 call abort_physic("thermcell_updown_dq", 'flux compensatoire ' & 168 168 // 'montant, cas non traite par thermcell_updown_dq', 1) … … 193 193 do ilay=nlay-1,1,-1 194 194 do ig=1,ngrid 195 if((fup(ig,ilay)-fdn(ig,ilay)) .lt.0) then195 if((fup(ig,ilay)-fdn(ig,ilay)) < 0) then 196 196 write(*,*) 'flux compensatoire montant, cas non traite par thermcell_updown_dq dans le cas d une resolution implicite, ilay : ', ilay 197 197 call abort_physic("thermcell_updown_dq", "", 1) … … 273 273 do ilay=nlay,1,-1 274 274 do ig=1,ngrid 275 if (ilay .le.lmax(ig).and.lmax(ig)>1) then275 if (ilay<=lmax(ig).and.lmax(ig)>1) then 276 276 edn(ig,ilay)=fact_thermals_down*dup(ig,ilay) 277 277 ddn(ig,ilay)=fact_thermals_down*eup(ig,ilay)
Note: See TracChangeset
for help on using the changeset viewer.