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

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

Parent Directory Parent Directory | Revision Log Revision Log

Revision 221 - (hide annotations)
Thu Apr 20 14:44:47 2017 UTC (7 years ago) by guez
File size: 1231 byte(s) 
clcdrag is no longer used in LMDZ. Replaced by cdrag in LMDZ. In cdrag
in LMDZ, zxli is a symbolic constant, false. So removed case zxli true
in LMDZE.

read_sst is called zero (if no ocean point on the whole planet) time or
once per call of physiq. If mod(itap - 1, lmt_pas) == 0 then we have
advanced in time of lmt_pas and deja_lu is necessarily false.

qsat[sl] and dqsat[sl] were never called.

Added output of qsurf in histins, following LMDZ.

Last dummy argument dtime of phystokenc is always the same as first
dummy argument pdtphys, removed dtime.

Removed make rules for nag_xref95, since it does not exist any longer.
1 guez 3 module hgardfou_m
2    
3     IMPLICIT none
4    
5     contains
6    
7 guez 207 SUBROUTINE hgardfou(t, ftsol)
8 guez 3
9 guez 221 ! From phylmd/hgardfou.F, v 1.1.1.1, 2004/05/19 12:53:07
10 guez 3
11     ! This procedure aborts the program if the temperature gets out of range.
12    
13 guez 69 USE indicesol, ONLY: nbsrf
14     USE dimphy, ONLY: klev, klon
15 guez 70 use nr_util, only: ifirstloc
16 guez 3
17 guez 207 REAL, intent(in):: t(klon, klev), ftsol(klon, nbsrf)
18 guez 3
19     ! Variables local to the procedure:
20    
21     real, parameter:: temp_min = 50., temp_max = 370. ! temperature range, in K
22 guez 70 INTEGER k, nsrf, jbad
23 guez 3
24     !----------------------------------------------------------
25    
26     DO k = 1, klev
27 guez 70 jbad = ifirstloc(t(:, k) > temp_max .or. t(:, k) < temp_min)
28     if (jbad <= klon) then
29     PRINT *, 'hgardfou: temperature out of range'
30     print *, "t(", jbad, ", ", k, ") = ", t(jbad, k)
31     stop 1
32     end if
33 guez 3 ENDDO
34    
35     DO nsrf = 1, nbsrf
36 guez 207 jbad = ifirstloc(ftsol(:, nsrf) > temp_max &
37     .or. ftsol(:, nsrf) < temp_min)
38 guez 70 if (jbad <= klon) then
39     PRINT *, 'hgardfou: temperature out of range'
40 guez 207 print *, "ftsol(", jbad, ", ", nsrf, ") = ", ftsol(jbad, nsrf)
41 guez 70 stop 1
42 guez 3 ENDIF
43     ENDDO
44    
45     END SUBROUTINE hgardfou
46    
47     end module hgardfou_m
  ViewVC Help
Powered by ViewVC 1.1.21