Changeset 899 for codes/icosagcm/trunk/src/dynamics/caldyn_kernels_hevi.F90
- Timestamp:
- 06/13/19 16:45:41 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/dynamics/caldyn_kernels_hevi.F90
r604 r899 117 117 REAL(rstd) :: C_ik(iim*jjm,llm) ! Thomas algorithm 118 118 REAL(rstd) :: D_il(iim*jjm,llm+1) ! Thomas algorithm 119 REAL(rstd) :: gamma, rho_ij, X_ij , Y_ij120 REAL(rstd) :: wil,tau2_g, g2, gm2, ml_g2, c2_mik119 REAL(rstd) :: gamma, rho_ij, X_ij 120 REAL(rstd) :: tau2_g, g2, gm2, ml_g2, c2_mik 121 121 122 122 INTEGER :: iter, ij, l, ij_omp_begin_ext, ij_omp_end_ext … … 133 133 tau2_g=tau*tau/g 134 134 g2=g*g 135 gm2 = g** -2135 gm2 = g**(-2) 136 136 gamma = 1./(1.-kappa) 137 137 … … 256 256 PRINT *, '[hevi_solver] C,D', iter, MAXVAL(ABS(C_ik)),MAXVAL(ABS(D_il)) 257 257 DO l=1,llm+1 258 WRITE(*,'(A,I2.1,I3.2,E9.2)') ,'[hevi_solver] x', iter,l, MAXVAL(ABS(x_il(:,l)))258 WRITE(*,'(A,I2.1,I3.2,E9.2)') '[hevi_solver] x', iter,l, MAXVAL(ABS(x_il(:,l))) 259 259 END DO 260 260 END IF … … 282 282 REAL(rstd) :: berni(iim*jjm,llm) ! (W/m_il)^2 283 283 REAL(rstd) :: berni1(iim*jjm) ! (W/m_il)^2 284 REAL(rstd) :: gamma, rho_ij, T_ij, X_ij, Y_ij,vreff, Rd, Cvd284 REAL(rstd) :: gamma, rho_ij, T_ij, X_ij, vreff, Rd, Cvd 285 285 INTEGER :: ij, l 286 286 … … 401 401 REAL(rstd) :: berniv(iim*jjm,llm) ! moist Bernoulli function 402 402 403 INTEGER :: i ,j,ij,l403 INTEGER :: ij,l 404 404 REAL(rstd) :: Rd, qv, temp, chi, nu, due, due_right, due_lup, due_ldown 405 405 … … 525 525 REAL(rstd) :: Ftheta(3*iim*jjm,llm) ! potential temperature flux 526 526 REAL(rstd) :: uu_right, uu_lup, uu_ldown, du_trisk, divF 527 INTEGER :: ij,iq,l ,kdown527 INTEGER :: ij,iq,l 528 528 529 529 CALL trace_start("compute_caldyn_Coriolis")
Note: See TracChangeset
for help on using the changeset viewer.