Changeset 513

Timestamp:

12/16/16 18:46:59 (7 years ago)

Author:

dubos

Message:

Trunk bugfix : sponge layer

File:

• codes/icosagcm/trunk/src/sponge.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/sponge.f90

@@ -103 +103 @@
     REAL(rstd),POINTER         :: due(:,:)!,due_sponge(:,:)
     REAL(rstd),POINTER         :: ue(:,:)
-    REAL(rstd),POINTER         :: theta_rhodz(:,:)
-    REAL(rstd),POINTER         :: dtheta_rhodz(:,:)!,dtheta_sponge(:,:)
+    REAL(rstd),POINTER         :: theta_rhodz(:,:,:)
+    REAL(rstd),POINTER         :: dtheta_rhodz(:,:,:)!,dtheta_sponge(:,:)
 
     INTEGER :: ind
@@ -131 +131 @@
             due(ij+u_ldown,l) = -rdamp(l)*ue(ij+u_ldown,l)
 
-            dtheta_rhodz(ij,l) = 0.0
+            dtheta_rhodz(ij,l,:) = 0.0
           ENDDO
         ENDDO