Changeset 686 for codes/icosagcm/devel/src/kernels_unst/compute_geopot.k90

Timestamp:

03/08/18 12:48:37 (6 years ago)

Author:

dubos

Message:

devel/unstructured : piecewise-constant vertical remapping

File:

• codes/icosagcm/devel/src/kernels_unst/compute_geopot.k90 (modified) (4 diffs)

codes/icosagcm/devel/src/kernels_unst/compute_geopot.k90

@@ -6 +6 @@
       !$OMP DO SCHEDULE(STATIC)
       DO ij=1,primal_num
-         l=llm
-         pk(l,ij) = ptop + .5*g* theta(l,ij,1)*rhodz(l,ij)
+         pk(llm,ij) = ptop + .5*g* theta(llm,ij,1)*rhodz(llm,ij)
          DO l = llm-1,1,-1
             pk(l,ij) = pk(l+1,ij) + (.5*g)*( theta(l,ij,1)*rhodz(l,ij) + theta(l+1,ij,1)*rhodz(l+1,ij) )
          END DO
          IF(caldyn_eta == eta_lag) THEN
-            l=1
-            ps(ij) = pk(l,ij) + .5*g* theta(l,ij,1)*rhodz(l,ij)
+            ps(ij) = pk(1,ij) + .5*g* theta(1,ij,1)*rhodz(1,ij)
          END IF
       END DO
@@ -29 +27 @@
       !$OMP DO SCHEDULE(STATIC)
       DO ij=1,primal_num
-         l=llm
-         pk(l,ij) = ptop + .5*g* rhodz(l,ij)
+         pk(llm,ij) = ptop + .5*g* rhodz(llm,ij)
          DO l = llm-1,1,-1
             pk(l,ij) = pk(l+1,ij) + (.5*g)*( rhodz(l,ij) + rhodz(l+1,ij) )
          END DO
          IF(caldyn_eta == eta_lag) THEN
-            l=1
-            ps(ij) = pk(l,ij) + .5*g* rhodz(l,ij)
+            ps(ij) = pk(1,ij) + .5*g* rhodz(1,ij)
          END IF
       END DO
@@ -54 +50 @@
       !$OMP DO SCHEDULE(STATIC)
       DO ij=1,primal_num
-         l=llm
-         pk(l,ij) = ptop + .5*g* rhodz(l,ij)
+         pk(llm,ij) = ptop + .5*g* rhodz(llm,ij)
          DO l = llm-1,1,-1
             pk(l,ij) = pk(l+1,ij) + (.5*g)*( rhodz(l,ij) + rhodz(l+1,ij) )
          END DO
          IF(caldyn_eta == eta_lag) THEN
-            l=1
-            ps(ij) = pk(l,ij) + .5*g* rhodz(l,ij)
+            ps(ij) = pk(1,ij) + .5*g* rhodz(1,ij)
          END IF
       END DO
@@ -79 +73 @@
       !$OMP DO SCHEDULE(STATIC)
       DO ij=1,primal_num
-         l=llm
-         pk(l,ij) = ptop + .5*g* rhodz(l,ij)*(1.+theta(l,ij,2))
+         pk(llm,ij) = ptop + .5*g* rhodz(llm,ij)*(1.+theta(llm,ij,2))
          DO l = llm-1,1,-1
             pk(l,ij) = pk(l+1,ij) + (.5*g)*( rhodz(l,ij)*(1.+theta(l,ij,2)) + rhodz(l+1,ij)*(1.+theta(l+1,ij,2)) )
          END DO
          IF(caldyn_eta == eta_lag) THEN
-            l=1
-            ps(ij) = pk(l,ij) + .5*g* rhodz(l,ij)*(1.+theta(l,ij,2))
+            ps(ij) = pk(1,ij) + .5*g* rhodz(1,ij)*(1.+theta(1,ij,2))
          END IF
       END DO