Changeset 23 for codes/icosagcm

Timestamp:

07/18/12 18:54:00 (12 years ago)

Author:

dubos

Message:

Fixed mistakes in transport code - now compiles and runs

Location:

codes/icosagcm/trunk/src

Files:

• advect.f90 (modified) (3 diffs)
• advect_tracer.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/advect.f90

@@ -54 +54 @@
     IMPLICIT NONE
     REAL(rstd),INTENT(IN)  :: qi(iim*jjm,llm)
-    REAL(rstd),INTENT(OUT) :: gradq3d(iim*jm,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,leng
@@ -115 +115 @@
 
   !===================================================================================================   
-  SUBROUTINE compute_advect_horiz(normal,tangent,qi,him,ue,he,bigt)
+  SUBROUTINE compute_advect_horiz(normal,tangent,qi,gradq3d,him,ue,he,bigt)
     USE domain_mod
     USE dimensions
@@ -124 +124 @@
     REAL(rstd),INTENT(IN)    :: tangent(3*iim*jjm,3)
     REAL(rstd),INTENT(INOUT) :: qi(iim*jjm,llm)
-    REAL(rstd),INTENT(IN)    :: gradq3d(iim*jm,llm,3) 
+    REAL(rstd),INTENT(IN)    :: gradq3d(iim*jjm,llm,3) 
     REAL(rstd),INTENT(INOUT) :: him(iim*jjm,llm)
     REAL(rstd),INTENT(IN)    :: ue(iim*3*jjm,llm)

codes/icosagcm/trunk/src/advect_tracer.f90

@@ -19 +19 @@
     REAL(rstd),POINTER :: tangent(:,:)
     REAL(rstd),POINTER :: normal(:,:)
+    INTEGER :: ind
 
     dt=dt_in
@@ -234 +235 @@
        DO k = 1,nqtot
           CALL compute_gradq3d(zq(:,:,k),gradq3d)
-          CALL compute_advect_horiz(zq(:,:,k),zm,u,massflx,bigt) 
+          CALL compute_advect_horiz(tangent,normal,zq(:,:,k),gradq3d,zm,u,massflx,bigt) 
        END DO
     END DO