Changeset 1047 for LMDZ4/trunk


Ignore:
Timestamp:
Nov 6, 2008, 5:46:57 PM (16 years ago)
Author:
lmdzadmin
Message:

Nettoyage ; print si prt_level >= 10
IM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/phylmd/add_phys_tend.F90

    r987 r1047  
    9292      DO j = 1, jbad
    9393         i=jadrs(j)
    94          print*,'PLANTAGE POUR LE POINT i rlon rlat =',i,rlon(i),rlat(i),text
    95          print*,'l    T     dT       Q     dQ    '
    96          DO k = 1, klev
     94         if(prt_level.ge.debug_level) THEN
     95          print*,'PLANTAGE POUR LE POINT i rlon rlat =',i,rlon(i),rlat(i),text
     96          print*,'l    T     dT       Q     dQ    '
     97          DO k = 1, klev
    9798             write(*,'(i3,2f14.4,2e14.2)') k,t_seri(i,k),zdt(i,k),q_seri(i,k),zdq(i,k)
    98          ENDDO
    99          call print_debug_phys(i,debug_level,text)
     99          ENDDO
     100          call print_debug_phys(i,debug_level,text)
     101         endif
    100102      ENDDO
    101103ENDIF
     
    107109      DO j = 1, jqbad
    108110         i=jqadrs(j)
    109          print*,'WARNING  : EAU POUR LE POINT i rlon rlat =',i,rlon(i),rlat(i),text
    110          print*,'l    T     dT       Q     dQ    '
     111         if(prt_level.ge.debug_level) THEN
     112          print*,'WARNING  : EAU POUR LE POINT i rlon rlat =',i,rlon(i),rlat(i),text
     113          print*,'l    T     dT       Q     dQ    '
     114         endif
    111115         DO k = 1, klev
    112116           zq=q_seri(i,k)+zdq(i,k)
    113117           if (zq.lt.1.e-15) then
    114118              if (q_seri(i,k).lt.1.e-15) then
    115 !              print*,' cas q_seri<1.e-15 i k q_seri zq zdq :',i,k,q_seri(i,k),zq,zdq(i,k)
     119               if(prt_level.ge.debug_level) THEN
     120                print*,' cas q_seri<1.e-15 i k q_seri zq zdq :',i,k,q_seri(i,k),zq,zdq(i,k)
     121               endif
    116122               q_seri(i,k)=1.e-15
    117123               zdq(i,k)=(1.e-15-q_seri(i,k))
     
    135141            jbad = jbad + 1
    136142            jadrs(jbad) = i
    137 !           if(prt_level.ge.10) THEN
    138 !             print*,'cas2 i k t_seri zdt',i,k,t_seri(i,k),zdt(i,k)
    139 !           endif
     143            if(prt_level.ge.debug_level) THEN
     144             print*,'cas2 i k t_seri zdt',i,k,t_seri(i,k),zdt(i,k)
     145            endif
    140146            ENDIF
    141147            IF ( q_seri(i,k)<0. .or. q_seri(i,k)>0.1 .or. abs(zdq(i,k))>1.e-2 ) then
     
    143149            jqadrs(jqbad) = i
    144150            kqadrs(jqbad) = k
    145 !           if(prt_level.ge.10) THEN
    146 !             print*,'cas2 i k q_seri zdq',i,k,q_seri(i,k),zdq(i,k)
    147 !           endif
     151            if(prt_level.ge.debug_level) THEN
     152             print*,'cas2 i k q_seri zdq',i,k,q_seri(i,k),zdq(i,k)
     153            endif
    148154            ENDIF
    149155         ENDDO
     
    153159         i=jadrs(j)
    154160         k=kadrs(j)
    155          print*,'PLANTAGE2 POUR LE POINT i itap rlon rlat txt jbad zdt t',i,itap,rlon(i),rlat(i),text,jbad, &
     161         if(prt_level.ge.debug_level) THEN
     162          print*,'PLANTAGE2 POUR LE POINT i itap rlon rlat txt jbad zdt t',i,itap,rlon(i),rlat(i),text,jbad, &
    156163       &        zdt(i,k),t_seri(i,k)-zdt(i,k)
    157 !        if(prt_level.ge.10) THEN
    158          if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN
     164!!!       if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN
    159165          print*,'l    T     dT       Q     dQ    '
    160166          DO k = 1, klev
     
    170176         i=jqadrs(j)
    171177         k=kqadrs(j)
    172          print*,'WARNING  : EAU2 POUR LE POINT i itap rlon rlat txt jqbad zdq q zdql ql',i,itap,rlon(i),rlat(i),text,jqbad,&
     178         if(prt_level.ge.debug_level) THEN
     179          print*,'WARNING  : EAU2 POUR LE POINT i itap rlon rlat txt jqbad zdq q zdql ql',i,itap,rlon(i),rlat(i),text,jqbad,&
    173180       &        zdq(i,k), q_seri(i,k)-zdq(i,k), zdql(i,k), ql_seri(i,k)-zdql(i,k)
    174 !        if(prt_level.ge.10) THEN
    175          if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN
     181!!!       if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN
    176182          print*,'l    T     dT       Q     dQ    '
    177183          DO k = 1, klev
Note: See TracChangeset for help on using the changeset viewer.