/[lmdze]/trunk/phylmd/tetalevel.f
ViewVC logotype

Diff of /trunk/phylmd/tetalevel.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 87 by guez, Wed Mar 5 14:57:53 2014 UTC revision 88 by guez, Tue Mar 11 15:09:02 2014 UTC
# Line 28  SUBROUTINE tetalevel(ilon, ilev, lnew, p Line 28  SUBROUTINE tetalevel(ilon, ilev, lnew, p
28    ! arguments :    ! arguments :
29    ! -----------    ! -----------
30    
31    INTEGER ilon, ilev    INTEGER, intent(in):: ilon, ilev
32    LOGICAL lnew    LOGICAL lnew
33    
34    REAL pgcm(ilon, ilev)    REAL pgcm(ilon, ilev)
35    REAL, INTENT (IN) :: qgcm(ilon, ilev)    REAL, INTENT (IN) :: qgcm(ilon, ilev)
36    REAL pres    REAL, INTENT (IN) :: pres
37    REAL qpres(ilon)    REAL qpres(ilon)
38    
39    ! local :    ! local :
# Line 58  SUBROUTINE tetalevel(ilon, ilev, lnew, p Line 58  SUBROUTINE tetalevel(ilon, ilev, lnew, p
58      ! Eventuellement, faire l'extrapolation a partir des deux couches      ! Eventuellement, faire l'extrapolation a partir des deux couches
59      ! les plus basses ou les deux couches les plus hautes:      ! les plus basses ou les deux couches les plus hautes:
60      DO i = 1, ilon      DO i = 1, ilon
       ! IM      IF ( ABS(pres-pgcm(i,ilev) ) .LT.  
61        IF (abs(pres-pgcm(i,ilev))>abs(pres-pgcm(i,1))) THEN        IF (abs(pres-pgcm(i,ilev))>abs(pres-pgcm(i,1))) THEN
62          lt(i) = ilev ! 2          lt(i) = ilev ! 2
63          lb(i) = ilev - 1 ! 1          lb(i) = ilev - 1 ! 1
# Line 66  SUBROUTINE tetalevel(ilon, ilev, lnew, p Line 65  SUBROUTINE tetalevel(ilon, ilev, lnew, p
65          lt(i) = 2          lt(i) = 2
66          lb(i) = 1          lb(i) = 1
67        END IF        END IF
       ! IM   PRINT*,'i, ABS(pres-pgcm),ABS(pres-pgcm)',  
       ! IM  .i, ABS(pres-pgcm(i,ilev)),ABS(pres-pgcm(i,1))  
68      END DO      END DO
69      DO k = 1, ilev - 1      DO k = 1, ilev - 1
70        DO i = 1, ilon        DO i = 1, ilon
Legend:
Removed from v.87  
changed lines
  Added in v.88
  ViewVC Help
Powered by ViewVC 1.1.21