Changeset 387 for codes/icosagcm/trunk/src/caldyn_hevi.f90

Timestamp:

05/28/16 00:32:21 (8 years ago)

Author:

dubos

Message:

Infrastructure for multiple dynamical tracers - tested with JW06 and moist baroclinic wave

File:

• codes/icosagcm/trunk/src/caldyn_hevi.f90 (modified) (3 diffs)

codes/icosagcm/trunk/src/caldyn_hevi.f90

@@ -49 +49 @@
     
     REAL(rstd),POINTER :: ps(:), dps(:)
-    REAL(rstd),POINTER :: mass(:,:), theta_rhodz(:,:), dtheta_rhodz(:,:)
+    REAL(rstd),POINTER :: mass(:,:), theta_rhodz(:,:,:), dtheta_rhodz(:,:,:)
     REAL(rstd),POINTER :: du(:,:), dW(:,:), dPhi(:,:), hflux(:,:), wflux(:,:)
     REAL(rstd),POINTER :: u(:,:), w(:,:), qu(:,:), qv(:,:)
@@ -107 +107 @@
        mass=f_mass(ind)
        theta = f_theta(ind)
-       CALL compute_theta(ps,theta_rhodz, mass,theta)
+       CALL compute_theta(ps,theta_rhodz(:,:,1), mass,theta)
        pk = f_pk(ind)
        geopot = f_geopot(ind)
@@ -163 +163 @@
           CALL compute_caldyn_slow_NH(u,mass,geopot,W, hflux,du,dPhi,dW)
        END IF
-       CALL compute_caldyn_Coriolis(hflux,theta,qu, convm,dtheta_rhodz,du)
+       CALL compute_caldyn_Coriolis(hflux,theta,qu, convm,dtheta_rhodz(:,:,1),du)
        IF(caldyn_eta==eta_mass) THEN
           wflux=f_wflux(ind)