Changeset 402 for codes/icosagcm/trunk/src/dcmip2016_simple_physics.f90
- Timestamp:
- 06/08/16 02:50:01 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/dcmip2016_simple_physics.f90
r381 r402 484 484 u(i,k) = CEm(i,k)*u(i,k-1)+CFu(i,k) 485 485 v(i,k) = CEm(i,k)*v(i,k-1)+CFv(i,k) 486 dtdt(i,k) = dtdt(i,k)+((CE(i,k)*t(i,k-1)*(p0/pmid(i,k-1))**(rair/cpair)+CFt(i,k))*(pmid(i,k)/p0)**(rair/cpair)-t(i,k))/dtime 486 dtdt(i,k) = dtdt(i,k)+ & 487 ((CE(i,k)*t(i,k-1)*(p0/pmid(i,k-1))**(rair/cpair)+CFt(i,k))*(pmid(i,k)/p0)**(rair/cpair)-t(i,k))/dtime 487 488 t(i,k) = (CE(i,k)*t(i,k-1)*(p0/pmid(i,k-1))**(rair/cpair)+CFt(i,k))*(pmid(i,k)/p0)**(rair/cpair) 488 489 dqdt(i,k) = dqdt(i,k)+(CE(i,k)*q(i,k-1)+CFq(i,k)-q(i,k))/dtime
Note: See TracChangeset
for help on using the changeset viewer.