Changeset 896 for codes/icosagcm/trunk/src/transport

Timestamp:

06/13/19 16:45:40 (5 years ago)

Author:

adurocher

Message:

trunk : compute_gradq3d: Remove all useless temporary arrays.

File:

• codes/icosagcm/trunk/src/transport/advect.F90 (modified) (4 diffs)

codes/icosagcm/trunk/src/transport/advect.F90

@@ -49 +49 @@
   !=======================================================================================
 
-  SUBROUTINE compute_gradq3d(qi_in,sqrt_leng_in,gradq3d_out,xyz_i,xyz_v)
+  SUBROUTINE compute_gradq3d(qi,sqrt_leng,gradq3d,xyz_i,xyz_v)
     USE trace
     USE omp_para
     IMPLICIT NONE
-    REAL(rstd),INTENT(IN)  :: qi_in(iim*jjm,llm)
-    REAL(rstd),INTENT(IN)  :: sqrt_leng_in(iim*jjm)
+    REAL(rstd),INTENT(IN)  :: qi(iim*jjm,llm)
+    REAL(rstd),INTENT(IN)  :: sqrt_leng(iim*jjm)
     REAL(rstd),INTENT(IN)  :: xyz_i(iim*jjm,3)
     REAL(rstd),INTENT(IN)  :: xyz_v(2*iim*jjm,3)
-    REAL(rstd),INTENT(OUT) :: gradq3d_out(iim*jjm,llm,3) 
+    REAL(rstd),INTENT(OUT) :: gradq3d(iim*jjm,llm,3)
     REAL(rstd) :: maxq,minq,minq_c,maxq_c 
     REAL(rstd) :: alphamx,alphami,alpha ,maggrd
@@ -65 +65 @@
     REAL(rstd) :: gradtri(2*iim*jjm,llm,3) 
     INTEGER :: ij,k,ind,l
-    REAL(rstd)  :: qi(iim*jjm,llm)
-    REAL(rstd)  :: sqrt_leng(iim*jjm)
-    REAL(rstd) :: gradq3d(iim*jjm,llm,3) 
     REAL(rstd) :: detx,dety,detz,det
     REAL(rstd) :: A(3,3), a11,a12,a13,a21,a22,a23,a31,a32,a33
@@ -73 +70 @@
     REAL(rstd) :: dq(3)
 
-    qi=qi_in
-    sqrt_leng=sqrt_leng_in
-    
     CALL trace_start("compute_gradq3d1")
 
@@ -272 +266 @@
   CALL trace_end("compute_gradq3d3")
   
-  
-  DO k=1,3
-    DO l = ll_begin,ll_end
-      DO ij=ij_begin,ij_end 
-        gradq3d_out(ij,l,k)=gradq3d(ij,l,k)
-      ENDDO
-    ENDDO
-  ENDDO
-  
   CONTAINS