Changeset 565 for codes/icosagcm/devel/src
- Timestamp:
- 10/09/17 01:22:21 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/dynamics/caldyn_kernels_hevi.F90
r562 r565 171 171 B_il(ij,1) = A_ik(ij,1) + ml_g2 + tau2_g*rho_bot 172 172 R_il(ij,1) = ml_g2*( Phi_il(ij,1)-Phi_star_il(ij,1)) & 173 + tau2_g*( p_ik(ij,1)-pbot+rho_bot*(Phi_il(ij,1)- Phi_bot) )173 + tau2_g*( p_ik(ij,1)-pbot+rho_bot*(Phi_il(ij,1)-phis(ij)) ) 174 174 ENDDO 175 175 ! inner interfaces … … 320 320 321 321 IF(tau>0) THEN ! solve implicit problem for geopotential 322 CALL compute_NH_geopot(tau, rhodz, m_il, theta, W, geopot)322 CALL compute_NH_geopot(tau, phis, rhodz, m_il, theta, W, geopot) 323 323 END IF 324 324
Note: See TracChangeset
for help on using the changeset viewer.