Changeset 824 for codes/icosagcm/devel/src/dcmip
- Timestamp:
- 04/30/19 16:59:07 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/dcmip/dcmip2016_kessler_physic.f90
r531 r824 158 158 159 159 ! Autoconversion and accretion rates following KW eq. 2.13a,b 160 qrprod = qc(k) - (qc(k)-dt0* amax1(.001*(qc(k)-.001d0),0.))/(1.d0+dt0*2.2d0*qr(k)**.875)161 qc(k) = amax1(qc(k)-qrprod,0.)162 qr(k) = amax1(qr(k)+qrprod+sed(k),0.)160 qrprod = qc(k) - (qc(k)-dt0*max(.001*(qc(k)-.001d0),0.))/(1.d0+dt0*2.2d0*qr(k)**.875) 161 qc(k) = max(qc(k)-qrprod,0.) 162 qr(k) = max(qr(k)+qrprod+sed(k),0.) 163 163 164 164 ! Saturation vapor mixing ratio (gm/gm) following KW eq. 2.11 … … 168 168 169 169 ! Evaporation rate following KW eq. 2.14a,b 170 ern = amin1(dt0*(((1.6d0+124.9d0*(r(k)*qr(k))**.2046) &170 ern = min(dt0*(((1.6d0+124.9d0*(r(k)*qr(k))**.2046) & 171 171 *(r(k)*qr(k))**.525)/(2550000d0*pc(k) & 172 172 /(3.8d0 *qvs)+540000d0))*(dim(qvs,qv(k)) & 173 /(r(k)*qvs)), amax1(-prod-qc(k),0.),qr(k))173 /(r(k)*qvs)),max(-prod-qc(k),0.),qr(k)) 174 174 175 175 ! Saturation adjustment following KW eq. 3.10 176 theta(k)= theta(k) + 2500000d0/(1003.d0*pk(k))*( amax1( prod,-qc(k))-ern)177 qv(k) = amax1(qv(k)-max(prod,-qc(k))+ern,0.)176 theta(k)= theta(k) + 2500000d0/(1003.d0*pk(k))*(max( prod,-qc(k))-ern) 177 qv(k) = max(qv(k)-max(prod,-qc(k))+ern,0.) 178 178 qc(k) = qc(k)+max(prod,-qc(k)) 179 179 qr(k) = qr(k)-ern
Note: See TracChangeset
for help on using the changeset viewer.