/[lmdze]/trunk/dyn3d/fxhyp_loop_ik.f
ViewVC logotype

Diff of /trunk/dyn3d/fxhyp_loop_ik.f

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

revision 121 by guez, Wed Jan 28 16:10:02 2015 UTC revision 122 by guez, Tue Feb 3 19:30:48 2015 UTC
# Line 109  contains Line 109  contains
109         END IF         END IF
110      END DO      END DO
111    
112      IF (.not. (MINval(xvrai(:iim)) >= - pi_d - 1d-5 &      IF (ik == 1 .and. (MINval(xvrai(:iim)) < - pi_d - 0.1d0 &
113           .and. MAXval(xvrai(:iim)) <= pi_d + 1d-5)) THEN           .or. MAXval(xvrai(:iim)) > pi_d + 0.1d0)) THEN
114         IF (xzoom <= 0.) THEN         IF (xzoom <= 0.) THEN
115            IF (ik == 1) THEN            i = 1
              i = 1  
116    
117               do while (xvrai(i) < - pi_d .and. i < iim)            do while (xvrai(i) < - pi_d .and. i < iim)
118                  i = i + 1               i = i + 1
119               end do            end do
120    
121               if (xvrai(i) < - pi_d) then            if (xvrai(i) < - pi_d) then
122                  print *, 'Xvrai plus petit que - pi !'               print *, 'Xvrai plus petit que - pi !'
123                  STOP 1               STOP 1
124               end if            end if
125    
126               is2 = i            is2 = i
127            END IF         ELSE
128              i = iim
129    
130              do while (xvrai(i) > pi_d .and. i > 1)
131                 i = i - 1
132              end do
133    
134              if (xvrai(i) > pi_d) then
135                 print *, 'Xvrai plus grand que pi !'
136                 STOP 1
137              end if
138    
139              is2 = i
140           END IF
141        END IF
142    
143        if (is2 /= 0) then
144           IF (xzoom <= 0.) THEN
145            IF (is2 /= 1) THEN            IF (is2 /= 1) THEN
146               DO ii = is2, iim               DO ii = is2, iim
147                  xlon(ii-is2 + 1) = xvrai(ii)                  xlon(ii-is2 + 1) = xvrai(ii)
# Line 137  contains Line 152  contains
152                  xprimm(ii + iim-is2 + 1) = xxprim(ii)                  xprimm(ii + iim-is2 + 1) = xxprim(ii)
153               END DO               END DO
154            END IF            END IF
155         ELSE         else
           IF (ik == 1) THEN  
              i = iim  
   
              do while (xvrai(i) > pi_d .and. i > 1)  
                 i = i - 1  
              end do  
   
              if (xvrai(i) > pi_d) then  
                 print *, 'Xvrai plus grand que pi !'  
                 STOP 1  
              end if  
   
              is2 = i  
           END IF  
   
156            idif = iim -is2            idif = iim -is2
157    
158            DO ii = 1, is2            DO ii = 1, is2
# Line 164  contains Line 164  contains
164               xlon(ii) = xvrai(ii + is2) - twopi_d               xlon(ii) = xvrai(ii + is2) - twopi_d
165               xprimm(ii) = xxprim(ii + is2)               xprimm(ii) = xxprim(ii + is2)
166            END DO            END DO
167         END IF         end IF
168      END IF      end if
169    
170      xlon(iim + 1) = xlon(1) + twopi      xlon(iim + 1) = xlon(1) + twopi
171      xprimm(iim + 1) = xprimm(1)      xprimm(iim + 1) = xprimm(1)
Legend:
Removed from v.121  
changed lines
  Added in v.122
  ViewVC Help
Powered by ViewVC 1.1.21