Changeset 4171 for LMDZ6/branches/LMDZ-ECRAD/libf/phylmd/thermcell_dq.F90
- Timestamp:
- Jun 17, 2022, 4:24:49 PM (2 years ago)
- Location:
- LMDZ6/branches/LMDZ-ECRAD
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/LMDZ-ECRAD
- Property svn:mergeinfo changed
-
LMDZ6/branches/LMDZ-ECRAD/libf/phylmd/thermcell_dq.F90
r2311 r4171 17 17 !======================================================================= 18 18 19 integer ngrid,nlay,impl 20 21 real ptimestep 22 real masse(ngrid,nlay),fm(ngrid,nlay+1) 23 real entr(ngrid,nlay) 24 real q(ngrid,nlay) 25 real dq(ngrid,nlay) 26 integer lev_out ! niveau pour les print 27 28 real qa(ngrid,nlay),detr(ngrid,nlay),wqd(ngrid,nlay+1) 29 19 ! arguments 20 integer, intent(in) :: ngrid,nlay,impl 21 real, intent(in) :: ptimestep 22 real, intent(in), dimension(ngrid,nlay) :: masse 23 real, intent(inout), dimension(ngrid,nlay) :: entr,q 24 real, intent(in), dimension(ngrid,nlay+1) :: fm 25 real, intent(out), dimension(ngrid,nlay) :: dq,qa 26 integer, intent(in) :: lev_out ! niveau pour les print 27 28 ! Local 29 real, dimension(ngrid,nlay) :: detr,qold 30 real, dimension(ngrid,nlay+1) :: wqd,fqa 30 31 real zzm 31 32 32 integer ig,k 33 33 real cfl 34 34 35 real qold(ngrid,nlay),fqa(ngrid,nlay+1)36 35 integer niter,iter 37 36 CHARACTER (LEN=20) :: modname='thermcell_dq' … … 40 39 41 40 ! Old explicite scheme 42 if (impl<=-1) then 41 if (impl<=-1) then 42 43 43 call thermcell_dq_o(ngrid,nlay,impl,ptimestep,fm,entr, & 44 44 & masse,q,dq,qa,lev_out) 45 return 46 endif 45 46 else 47 47 48 48 49 ! Calcul du critere CFL pour l'advection dans la subsidence … … 143 144 enddo 144 145 145 return 146 endif ! impl=-1 147 RETURN 146 148 end 149 147 150 148 151 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Note: See TracChangeset
for help on using the changeset viewer.