- Timestamp:
- 06/13/19 16:45:40 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/transport/advect.F90
r894 r896 49 49 !======================================================================================= 50 50 51 SUBROUTINE compute_gradq3d(qi _in,sqrt_leng_in,gradq3d_out,xyz_i,xyz_v)51 SUBROUTINE compute_gradq3d(qi,sqrt_leng,gradq3d,xyz_i,xyz_v) 52 52 USE trace 53 53 USE omp_para 54 54 IMPLICIT NONE 55 REAL(rstd),INTENT(IN) :: qi _in(iim*jjm,llm)56 REAL(rstd),INTENT(IN) :: sqrt_leng _in(iim*jjm)55 REAL(rstd),INTENT(IN) :: qi(iim*jjm,llm) 56 REAL(rstd),INTENT(IN) :: sqrt_leng(iim*jjm) 57 57 REAL(rstd),INTENT(IN) :: xyz_i(iim*jjm,3) 58 58 REAL(rstd),INTENT(IN) :: xyz_v(2*iim*jjm,3) 59 REAL(rstd),INTENT(OUT) :: gradq3d _out(iim*jjm,llm,3)59 REAL(rstd),INTENT(OUT) :: gradq3d(iim*jjm,llm,3) 60 60 REAL(rstd) :: maxq,minq,minq_c,maxq_c 61 61 REAL(rstd) :: alphamx,alphami,alpha ,maggrd … … 65 65 REAL(rstd) :: gradtri(2*iim*jjm,llm,3) 66 66 INTEGER :: ij,k,ind,l 67 REAL(rstd) :: qi(iim*jjm,llm)68 REAL(rstd) :: sqrt_leng(iim*jjm)69 REAL(rstd) :: gradq3d(iim*jjm,llm,3)70 67 REAL(rstd) :: detx,dety,detz,det 71 68 REAL(rstd) :: A(3,3), a11,a12,a13,a21,a22,a23,a31,a32,a33 … … 73 70 REAL(rstd) :: dq(3) 74 71 75 qi=qi_in76 sqrt_leng=sqrt_leng_in77 78 72 CALL trace_start("compute_gradq3d1") 79 73 … … 272 266 CALL trace_end("compute_gradq3d3") 273 267 274 275 DO k=1,3276 DO l = ll_begin,ll_end277 DO ij=ij_begin,ij_end278 gradq3d_out(ij,l,k)=gradq3d(ij,l,k)279 ENDDO280 ENDDO281 ENDDO282 283 268 CONTAINS 284 269
Note: See TracChangeset
for help on using the changeset viewer.