Changeset 340

Timestamp:

07/30/15 18:00:12 (9 years ago)

Author:

dubos

Message:

etat0_dcmip5 robustness fix

File:

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

codes/icosagcm/trunk/src/etat0_dcmip5.f90

@@ -1 +1 @@
 MODULE etat0_dcmip5_mod
   USE icosa
+  IMPLICIT NONE
   PRIVATE
   REAL(rstd),PARAMETER :: zt=15000
@@ -27 +28 @@
   
   SUBROUTINE getin_etat0
+    USE mpipara, ONLY : is_mpi_root
     dcmip5_testcase=1
     CALL getin("dcmip5_testcase",dcmip5_testcase)
+    IF(nqtot<1) THEN
+       IF (is_mpi_root)  THEN
+          PRINT *, "nqtot must be at least 1 for test case DCMIP5"
+       END IF
+       STOP
+    END IF
   END SUBROUTINE getin_etat0
 
   SUBROUTINE compute_etat0(ngrid,lon,lat, phis, ps, Temp, ulon, ulat, q)
     USE disvert_mod
-    IMPLICIT NONE
     INTEGER, INTENT(IN)    :: ngrid
     REAL(rstd),INTENT(IN)  :: lon(ngrid)
@@ -79 +86 @@
           d2=cos(latc)*sin(lon(ij)-lonc)
           d=MAX(1e-25,sqrt(d1**2+d2**2))
-          ulon(ij+u_right,l)=vt*d1/d
-          ulat(ij+u_right,l)=vt*d2/d
+          ulon(ij,l)=vt*d1/d
+          ulat(ij,l)=vt*d2/d
        END DO
     ENDDO