/[lmdze]/trunk/dyn3d/Inter_barxy/inter_barxy.f
ViewVC logotype

Diff of /trunk/dyn3d/Inter_barxy/inter_barxy.f

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

revision 39 by guez, Tue Jan 25 15:11:05 2011 UTC revision 73 by guez, Fri Nov 15 17:48:30 2013 UTC
# Line 55  contains Line 55  contains
55    
56      ! Check decreasing order for "rlatimod":      ! Check decreasing order for "rlatimod":
57      DO i = 2, jjm      DO i = 2, jjm
58         IF (rlatimod(i) >= rlatimod(i-1)) stop &         IF (rlatimod(i) >= rlatimod(i-1)) then
59              '"inter_barxy": "rlatimod" should be strictly decreasing'            print *, '"inter_barxy": "rlatimod" should be strictly decreasing'
60              stop 1
61           end IF
62      ENDDO      ENDDO
63    
64      yjmod(:jjm) = ord_coordm(rlatimod)      yjmod(:jjm) = ord_coordm(rlatimod)
65      IF (jmods == jjm + 1) THEN      IF (jmods == jjm + 1) THEN
66         IF (90. - yjmod(jjm) < 0.01) stop &         IF (90. - yjmod(jjm) < 0.01) then
67              '"inter_barxy": with jmods = jjm + 1, yjmod(jjm) should be < 90.'            print *, '"inter_barxy": with jmods = jjm + 1, ' &
68                   // 'yjmod(jjm) should be < 90.'
69              stop 1
70           end IF
71      ELSE      ELSE
72         ! jmods = jjm         ! jmods = jjm
73         IF (ABS(yjmod(jjm) - 90.) > 0.01) stop &         IF (ABS(yjmod(jjm) - 90.) > 0.01) then
74              '"inter_barxy": with jmods = jjm, yjmod(jjm) should be 90.'            print *, '"inter_barxy": with jmods = jjm, yjmod(jjm) should be 90.'
75              stop 1
76           end IF
77      ENDIF      ENDIF
78    
79      if (jmods == jjm + 1) yjmod(jjm + 1) = 90.      if (jmods == jjm + 1) yjmod(jjm + 1) = 90.
# Line 397  contains Line 404  contains
404      ! Check monotonicity:      ! Check monotonicity:
405      decrois = xi(2) < xi(1)      decrois = xi(2) < xi(1)
406      DO i = 3, nmax      DO i = 3, nmax
407         IF (decrois .neqv. xi(i) < xi(i-1)) stop &         IF (decrois .neqv. xi(i) < xi(i-1)) then
408              '"ord_coord":  latitudes are not monotonic'            print *, '"ord_coord":  latitudes are not monotonic'
409              stop 1
410           end IF
411      ENDDO      ENDDO
412    
413      IF (abs(xi(1)) < pi) then      IF (abs(xi(1)) < pi) then
# Line 413  contains Line 422  contains
422         print *, "ord_coord"         print *, "ord_coord"
423         PRINT *, '"xi" should contain the latitudes of the boundaries of ' &         PRINT *, '"xi" should contain the latitudes of the boundaries of ' &
424              // 'grid cells, not the centers of grid cells.'              // 'grid cells, not the centers of grid cells.'
425         STOP         STOP 1
426      ENDIF      ENDIF
427    
428      IF (decrois) xo(:nmax) = xo(nmax:1:- 1)      IF (decrois) xo(:nmax) = xo(nmax:1:- 1)
Legend:
Removed from v.39  
changed lines
  Added in v.73
  ViewVC Help
Powered by ViewVC 1.1.21