Changeset 359

Timestamp:

09/08/15 15:22:15 (9 years ago)

Author:

dubos

Message:

Fix Held&Suarez when nqtot<3

File:

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

codes/icosagcm/trunk/src/etat0_heldsz.f90

@@ -89 +89 @@
        CALL compute_etat0_heldsz(theta_eq,theta)
        CALL compute_theta2theta_rhodz(ps,theta,theta_rhodz,1)
-       q=f_q(ind)
-       q(:,:,1)=1e-2
-       q(:,:,2)=0
-       q(:,:,3:)=1e-2
+       IF(nqtot>0) THEN
+          q=f_q(ind)
+          q(:,:,1)=1e-2
+          IF(nqtot>1) q(:,:,2)=0
+          IF(nqtot>2) q(:,:,3:)=1e-2
+       END IF
     ENDDO
   END SUBROUTINE etat0