- Timestamp:
- Nov 6, 2008, 5:46:57 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/phylmd/add_phys_tend.F90
r987 r1047 92 92 DO j = 1, jbad 93 93 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 97 98 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 100 102 ENDDO 101 103 ENDIF … … 107 109 DO j = 1, jqbad 108 110 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 111 115 DO k = 1, klev 112 116 zq=q_seri(i,k)+zdq(i,k) 113 117 if (zq.lt.1.e-15) then 114 118 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 116 122 q_seri(i,k)=1.e-15 117 123 zdq(i,k)=(1.e-15-q_seri(i,k)) … … 135 141 jbad = jbad + 1 136 142 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 !endif143 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 140 146 ENDIF 141 147 IF ( q_seri(i,k)<0. .or. q_seri(i,k)>0.1 .or. abs(zdq(i,k))>1.e-2 ) then … … 143 149 jqadrs(jqbad) = i 144 150 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 !endif151 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 148 154 ENDIF 149 155 ENDDO … … 153 159 i=jadrs(j) 154 160 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, & 156 163 & 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 159 165 print*,'l T dT Q dQ ' 160 166 DO k = 1, klev … … 170 176 i=jqadrs(j) 171 177 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,& 173 180 & 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 176 182 print*,'l T dT Q dQ ' 177 183 DO k = 1, klev
Note: See TracChangeset
for help on using the changeset viewer.