source: LMDZ6/branches/LMDZ_cdrag_LSCE/libf/phylmd/minmaxqfi.F90 @ 5423

Last change on this file since 5423 was 2311, checked in by Ehouarn Millour, 10 years ago

Further modifications to enforce physics/dynamics separation:

  • moved iniprint.h and misc_mod back to dyn3d_common, as these should only be used by dynamics.
  • created print_control_mod in the physics to store flags prt_level, lunout, debug to be local to physics (should be used rather than iniprint.h)
  • created abort_physic.F90 , which does the same job as abort_gcm() did, but should be used instead when in physics.
  • reactivated inifis (turned it into a module, inifis_mod.F90) to initialize physical constants and print_control_mod flags.

EM

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 842 bytes
Line 
1!
2! $Id: minmaxqfi.F90 2311 2015-06-25 07:45:24Z evignon $
3!
4SUBROUTINE minmaxqfi(zq,qmin,qmax,comment)
5  USE dimphy
6  USE print_control_mod, ONLY: prt_level
7  IMPLICIT NONE
8
9! Entrees
10  REAL,DIMENSION(klon,klev), INTENT(IN)   :: zq
11  REAL,INTENT(IN)                         :: qmin,qmax
12  CHARACTER(LEN=*),INTENT(IN)             :: comment
13
14! Local 
15  INTEGER,DIMENSION(klon)     :: jadrs
16  INTEGER                     :: i, jbad, k
17
18  DO k = 1, klev
19     jbad = 0
20     DO i = 1, klon
21        IF (zq(i,k).GT.qmax .OR. zq(i,k).LT.qmin) THEN
22           jbad = jbad + 1
23           jadrs(jbad) = i
24        ENDIF
25     ENDDO
26     IF (jbad.GT.0) THEN
27        WRITE(*,*)comment, "k = ", k, "jbad > 0"
28        if (prt_level >= 1) then
29           DO i = 1, jbad
30              WRITE(*,*) "i,q=", jadrs(i), zq(jadrs(i),k)
31           ENDDO
32        end if
33     ENDIF
34  ENDDO
35 
36END SUBROUTINE minmaxqfi
Note: See TracBrowser for help on using the repository browser.