source: LMDZ6/branches/LMDZ_DECOUPLE/libf/phylmd/minmaxqfi.F90 @ 5474

Last change on this file since 5474 was 4795, checked in by nfevrier, 12 months ago

First save of N. Février's modifications
  • 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: 1.1 KB
RevLine 
[524]1!
[1191]2! $Id: minmaxqfi.F90 4795 2024-02-01 15:53:24Z dcugnet $
[524]3!
[4795]4SUBROUTINE minmaxqfi(zq,qmin,qmax,comment,it)
[1191]5  USE dimphy
[2311]6  USE print_control_mod, ONLY: prt_level
[4795]7  USE infotrac_phy, ONLY: nbtr, type_trac, conv_flg, solsym
[1191]8  IMPLICIT NONE
[524]9
[1191]10! Entrees
11  REAL,DIMENSION(klon,klev), INTENT(IN)   :: zq
12  REAL,INTENT(IN)                         :: qmin,qmax
13  CHARACTER(LEN=*),INTENT(IN)             :: comment
[4795]14  INTEGER,INTENT(IN)                      :: it
[524]15
[1191]16! Local 
17  INTEGER,DIMENSION(klon)     :: jadrs
18  INTEGER                     :: i, jbad, k
[2098]19
[1191]20  DO k = 1, klev
21     jbad = 0
22     DO i = 1, klon
23        IF (zq(i,k).GT.qmax .OR. zq(i,k).LT.qmin) THEN
24           jbad = jbad + 1
25           jadrs(jbad) = i
26        ENDIF
27     ENDDO
[4795]28     IF ((jbad.GT.0) .AND. ((TRIM(solsym(it))=='AIBCM') .OR. (TRIM(solsym(it))=='ASBCM'))) THEN ! Nicolas Février : alerte seulement si le BC est concerné
[2098]29        WRITE(*,*)comment, "k = ", k, "jbad > 0"
[4795]30        if (prt_level >= 1) then ! Nicolas Février on print tout
[2098]31           DO i = 1, jbad
32              WRITE(*,*) "i,q=", jadrs(i), zq(jadrs(i),k)
33           ENDDO
34        end if
[1191]35     ENDIF
36  ENDDO
[4795]37
[1191]38END SUBROUTINE minmaxqfi
Note: See TracBrowser for help on using the repository browser.