Changeset 3587

Timestamp:

Oct 15, 2019, 6:23:45 PM (5 years ago)

Author:

fhourdin

Message:

Suite du bug precedent ...

File:

• LMDZ6/trunk/libf/phylmd/orografi_strato.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/orografi_strato.F90

@@ -1868 +1868 @@
   IF (dist_min/=dist_min_glo) cell=-1
 !ym in future find the point at equator close to (0,0) coordinates.
-  !PRINT *, 'SUGWD distmin dist_min_glo cell=', dist_min,dist_min_glo,cell
+  PRINT *, 'SUGWD distmin dist_min_glo cell=', dist_min,dist_min_glo,cell
 
   nktopg_tmp=nktopg
@@ -1875 +1875 @@
   IF (cell/=-1) THEN
 
+    !print*,'SUGWD shape ',shape(pplay),cell+1
+
     DO jk = 1, nlev
-      zpm1r = pplay(cell+1, jk)/paprs(cell+1, 1)
+      !zpm1r = pplay(cell+1, jk)/paprs(cell+1, 1)
+      zpm1r = pplay(cell, jk)/paprs(cell, 1)
       IF (zpm1r>=zsigt) THEN
         nktopg_tmp = jk