Changeset 580 for codes/icosagcm/trunk/src/dynamics/caldyn_hevi.f90

Timestamp:

10/13/17 15:55:07 (7 years ago)

Author:

dubos

Message:

trunk : upgrading to devel

File:

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

codes/icosagcm/trunk/src/dynamics/caldyn_hevi.f90

@@ -48 +48 @@
     TYPE(t_field) :: f_dPhi_fast(:)
     
-    REAL(rstd),POINTER :: ps(:), dps(:)
+    REAL(rstd),POINTER :: ps(:), dps(:), phis(:)
     REAL(rstd),POINTER :: mass(:,:), theta_rhodz(:,:,:), dtheta_rhodz(:,:,:)
     REAL(rstd),POINTER :: du(:,:), dW(:,:), dPhi(:,:), hflux(:,:), wflux(:,:)
@@ -115 +115 @@
           CALL compute_geopot(mass,theta, ps,pk,geopot)
        ELSE
+          phis = f_phis(ind)
           W = f_W(ind)
           dW = f_dW_fast(ind)
@@ -121 +122 @@
           m_il = f_wil(ind)
           pres = f_gradPhi2(ind)
-          CALL compute_caldyn_solver(tau,mass,theta,pk,geopot,W, m_il,pres, dPhi,dW,du) ! computes d(Phi,W,du)_fast and updates Phi,W
+          CALL compute_caldyn_solver(tau,phis, mass,theta,pk,geopot,W, m_il,pres, dPhi,dW,du) ! computes d(Phi,W,du)_fast and updates Phi,W
        END IF
        u=f_u(ind)