Changeset 2177 for trunk/LMDZ.GENERIC/libf/phystd/thermcell_dq.F90
- Timestamp:
- Nov 12, 2019, 11:55:58 AM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.GENERIC/libf/phystd/thermcell_dq.F90
r2144 r2177 3 3 ! 4 4 SUBROUTINE thermcell_dq(ngrid,nlay,ptimestep,fm,entr,detr,masse, & 5 q,dq,qa,lmin )5 q,dq,qa,lmin,lmax) 6 6 7 7 … … 16 16 ! 17 17 ! Modif 2019/04 (AB alexandre.boissinot@lmd.jussieu.fr) 18 ! dqimpl = 1: implicit scheme19 ! dqimpl = 0: explicit scheme18 ! dqimpl = true : implicit scheme 19 ! dqimpl = false : explicit scheme 20 20 ! 21 21 !=============================================================================== … … 34 34 ! ------- 35 35 36 INTEGER ngrid, nlay 37 INTEGER lmin(ngrid) 36 INTEGER, INTENT(in) :: ngrid 37 INTEGER, INTENT(in) :: nlay 38 INTEGER, INTENT(in) :: lmin(ngrid) 39 INTEGER, INTENT(in) :: lmax(ngrid) 38 40 39 REAL ptimestep 40 REAL masse(ngrid,nlay) 41 REAL fm(ngrid,nlay+1) 42 REAL entr(ngrid,nlay) 43 REAL detr(ngrid,nlay) 44 REAL q(ngrid,nlay) 41 REAL, INTENT(in) :: ptimestep 42 REAL, INTENT(in) :: masse(ngrid,nlay) 43 REAL, INTENT(in) :: fm(ngrid,nlay+1) 44 REAL, INTENT(in) :: entr(ngrid,nlay) 45 REAL, INTENT(in) :: detr(ngrid,nlay) 45 46 46 47 ! Outputs: 47 48 ! -------- 48 49 49 REAL dq(ngrid,nlay) 50 REAL qa(ngrid,nlay) 50 REAL, INTENT(inout) :: q(ngrid,nlay) 51 REAL, INTENT(out) :: dq(ngrid,nlay) 52 REAL, INTENT(out) :: qa(ngrid,nlay) 51 53 52 54 ! Local: … … 82 84 cfl = max(cfl, fm(ig,l) / zzm) 83 85 84 IF (entr(ig,l) .gt.zzm) THEN86 IF (entr(ig,l) > zzm) THEN 85 87 print *, 'ERROR: entrainment is greater than the layer mass!' 86 88 print *, 'ig,l,entr', ig, l, entr(ig,l) 87 print *, 'lmin ', lmin(ig)89 print *, 'lmin,lmax', lmin(ig), lmax(ig) 88 90 print *, '-------------------------------' 89 91 print *, 'entr*dt,mass', entr(ig,l)*ptimestep, masse(ig,l)
Note: See TracChangeset
for help on using the changeset viewer.