Changeset 377 for codes/icosagcm/trunk/src/caldyn_kernels_base.f90

Timestamp:

04/14/16 23:12:41 (8 years ago)

Author:

dubos

Message:

New : positive advection option for theta

File:

• codes/icosagcm/trunk/src/caldyn_kernels_base.f90 (modified) (1 diff)

codes/icosagcm/trunk/src/caldyn_kernels_base.f90

@@ -302 +302 @@
        !DIR$ SIMD
        DO ij=ij_begin,ij_end
-          dW(ij,l) = dW(ij,l) + W_etadot(ij,l-1) - W_etadot(ij,l)
           dPhi(ij,l) = dPhi(ij,l) - wflux(ij,l) &
                * (geopot(ij,l+1)-geopot(ij,l-1))/(mass(ij,l-1)+mass(ij,l))
        END DO
     END DO
+    DO l=ll_begin,ll_end
+       !DIR$ SIMD
+       DO ij=ij_begin,ij_end
+          dW(ij,l+1) = dW(ij,l+1) + W_etadot(ij,l) ! update inner+top interfaces
+          dW(ij,l)   = dW(ij,l)   - W_etadot(ij,l) ! update bottom+inner interfaces
+       END DO
+    END DO
     CALL trace_end("compute_caldyn_vert_nh")