/[lmdze]/trunk/Sources/phylmd/hgardfou.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/hgardfou.f

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

revision 222 by guez, Tue Apr 25 15:31:48 2017 UTC revision 227 by guez, Thu Nov 2 15:47:03 2017 UTC
# Line 10  contains Line 10  contains
10    
11      ! This procedure aborts the program if the temperature gets out of range.      ! This procedure aborts the program if the temperature gets out of range.
12    
13      USE indicesol, ONLY: nbsrf      use abort_gcm_m, only: abort_gcm
14        USE indicesol, ONLY: nbsrf, clnsurf
15      USE dimphy, ONLY: klev, klon      USE dimphy, ONLY: klev, klon
16      use nr_util, only: ifirstloc      use nr_util, only: ifirstloc
17        use phyetat0_m, only: rlon, rlat
18    
19      REAL, intent(in):: t_seri(:, :) ! (klon, klev)      REAL, intent(in):: t_seri(:, :) ! (klon, klev)
20      REAL, intent(in):: ftsol(:, :) ! (klon, nbsrf)      REAL, intent(in):: ftsol(:, :) ! (klon, nbsrf)
# Line 27  contains Line 29  contains
29      DO k = 1, klev      DO k = 1, klev
30         jbad = ifirstloc(t_seri(:, k) > temp_max .or. t_seri(:, k) < temp_min)         jbad = ifirstloc(t_seri(:, k) > temp_max .or. t_seri(:, k) < temp_min)
31         if (jbad <= klon) then         if (jbad <= klon) then
           PRINT *, 'hgardfou: temperature out of range'  
32            print *, "t_seri(", jbad, ", ", k, ") = ", t_seri(jbad, k)            print *, "t_seri(", jbad, ", ", k, ") = ", t_seri(jbad, k)
33            stop 1            call abort_gcm('hgardfou', 'temperature out of range')
34         end if         end if
35      ENDDO      ENDDO
36    
# Line 37  contains Line 38  contains
38         jbad = ifirstloc(ftsol(:, nsrf) > temp_max &         jbad = ifirstloc(ftsol(:, nsrf) > temp_max &
39              .or. ftsol(:, nsrf) < temp_min)              .or. ftsol(:, nsrf) < temp_min)
40         if (jbad <= klon) then         if (jbad <= klon) then
41            PRINT *, 'hgardfou: temperature out of range'            print *, "ftsol(position index =", jbad, ", sub-surface index =", &
42            print *, "ftsol(", jbad, ", ", nsrf, ") = ", ftsol(jbad, nsrf)                 nsrf, ") =", ftsol(jbad, nsrf)
43            stop 1            print *, "sub-surface name: ", clnsurf(nsrf)
44              print *, "longitude:", rlon(jbad), "degrees east"
45              print *, "latitude:", rlat(jbad), "degrees north"
46              call abort_gcm('hgardfou', 'temperature out of range')
47         ENDIF         ENDIF
48      ENDDO      ENDDO
49    
Legend:
Removed from v.222  
changed lines
  Added in v.227
  ViewVC Help
Powered by ViewVC 1.1.21