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

Diff of /trunk/dyn3d/inter_barxy.f

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

trunk/libf/dyn3d/inter_barxy.f90 revision 36 by guez, Thu Dec 2 17:11:04 2010 UTC trunk/dyn3d/inter_barxy.f revision 82 by guez, Wed Mar 5 14:57:53 2014 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 373  contains Line 380  contains
380      ! reverses their order.      ! reverses their order.
381      ! Finally, the procedure adds 90° as the last value of the array.      ! Finally, the procedure adds 90° as the last value of the array.
382    
383      use nr_util, only: assert_eq      use nr_util, only: assert_eq, pi
     use comconst, only: pi  
384    
385      IMPLICIT NONE      IMPLICIT NONE
386    
# Line 398  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 414  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)
# Line 432  contains Line 440  contains
440      ! This procedure converts to degrees, if necessary, and inverts the      ! This procedure converts to degrees, if necessary, and inverts the
441      ! order.      ! order.
442    
443      use comconst, only: pi      use nr_util, only: pi
444    
445      IMPLICIT NONE      IMPLICIT NONE
446    
Legend:
Removed from v.36  
changed lines
  Added in v.82
  ViewVC Help
Powered by ViewVC 1.1.21