Changeset 192 for codes/icosagcm/trunk/src/etat0.f90

Timestamp:

06/25/14 15:00:24 (10 years ago)

Author:

dubos

Message:

Fixed DCMIP5 physics/etat0

File:

• codes/icosagcm/trunk/src/etat0.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/etat0.f90

@@ -7 +7 @@
   SUBROUTINE etat0(f_ps,f_mass,f_phis,f_theta_rhodz,f_u, f_q)
     USE icosa
+    USE mpipara, ONLY : is_mpi_root
     USE disvert_mod
     USE etat0_williamson_mod, ONLY : etat0_williamson_new
@@ -61 +62 @@
        CALL etat0_dcmip3(f_ps,f_phis,f_theta_rhodz,f_u, f_q)
      CASE ('dcmip4')
+        IF(nqtot<2) THEN
+           IF (is_mpi_root)  THEN
+              PRINT *, "nqtot must be at least 2 for test case DCMIP4"
+           END IF
+           STOP
+        END IF
        CALL etat0_dcmip4(f_ps,f_phis,f_theta_rhodz,f_u, f_q)
      CASE ('dcmip5')