Ignore:
Timestamp:
Oct 15, 2019, 6:23:45 PM (5 years ago)
Author:
fhourdin
Message:

Suite du bug precedent ...

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/orografi_strato.F90

    r3586 r3587  
    18681868  IF (dist_min/=dist_min_glo) cell=-1
    18691869!ym in future find the point at equator close to (0,0) coordinates.
    1870   !PRINT *, 'SUGWD distmin dist_min_glo cell=', dist_min,dist_min_glo,cell
     1870  PRINT *, 'SUGWD distmin dist_min_glo cell=', dist_min,dist_min_glo,cell
    18711871
    18721872  nktopg_tmp=nktopg
     
    18751875  IF (cell/=-1) THEN
    18761876
     1877    !print*,'SUGWD shape ',shape(pplay),cell+1
     1878
    18771879    DO jk = 1, nlev
    1878       zpm1r = pplay(cell+1, jk)/paprs(cell+1, 1)
     1880      !zpm1r = pplay(cell+1, jk)/paprs(cell+1, 1)
     1881      zpm1r = pplay(cell, jk)/paprs(cell, 1)
    18791882      IF (zpm1r>=zsigt) THEN
    18801883        nktopg_tmp = jk
Note: See TracChangeset for help on using the changeset viewer.