Changeset 485 for codes

Timestamp:

10/12/16 10:07:58 (8 years ago)

Author:

ymipsl

Message:

Fix extra-dimension to theta_rhodz for etat0 h&Z

YM

File:

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

codes/icosagcm/trunk/src/etat0_heldsz.f90

@@ -36 +36 @@
     CALL allocate_field(f_ps,field_t,type_real)
     CALL allocate_field(f_phis,field_t,type_real)
-    CALL allocate_field(f_theta_rhodz,field_t,type_real,llm)
+    CALL allocate_field(f_theta_rhodz,field_t,type_real,llm,2)
     CALL allocate_field(f_u,field_u,type_real,llm)
     CALL allocate_field(f_Ki,field_t,type_real,llm)
@@ -58 +58 @@
     REAL(rstd),POINTER :: ps(:)
     REAL(rstd),POINTER :: phis(:)
-    REAL(rstd),POINTER :: theta_rhodz(:,:)
+    REAL(rstd),POINTER :: theta_rhodz(:,:,:)
     REAL(rstd),POINTER :: u(:,:)
     REAL(rstd),POINTER :: q(:,:,:)
@@ -85 +85 @@
        theta=f_theta(ind)
        CALL compute_etat0_heldsz(theta_eq,theta)
-       CALL compute_theta2theta_rhodz(ps,theta,theta_rhodz,1)
+       CALL compute_theta2theta_rhodz(ps,theta,theta_rhodz(:,:,1),1)
        IF(nqtot>0) THEN
           q=f_q(ind)