Changeset 5082 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/pentes_ini.F
- Timestamp:
- Jul 19, 2024, 5:41:58 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/pentes_ini.F
r2600 r5082 74 74 c if (mode.eq.1.or.mode.eq.3) then 75 75 c if (mode.eq.1) then 76 if (mode .ge.1) then76 if (mode>=1) then 77 77 lati=2 78 78 latf=jjm … … 187 187 c enddo 188 188 CCC 189 if(mode .eq.2) then189 if(mode==2) then 190 190 do l=1,llm 191 191 s0s=0. … … 247 247 endif 248 248 249 if (mode .eq.4) then249 if (mode==4) then 250 250 do l=1,llm 251 251 do i=1,iip1 … … 261 261 call advx( limit,.5*dtvr,pbaru,sm,s0,sx,sy,sz,lati,latf) 262 262 c call minmaxq(zq,1.e33,-1.e33,'avant advy ') 263 if (mode .eq.4) then263 if (mode==4) then 264 264 do l=1,llm 265 265 do i=1,iip1 … … 282 282 call limz(s0,sz,sm,pente_max) 283 283 call advz( limit,dtvr,w,sm,s0,sx,sy,sz ) 284 if (mode .eq.4) then284 if (mode==4) then 285 285 do l=1,llm 286 286 do i=1,iip1 … … 306 306 307 307 c call minmaxq(zq,1.e33,-1.e33,'avant advx ') 308 if (mode .eq.4) then308 if (mode==4) then 309 309 do l=1,llm 310 310 do i=1,iip1 … … 346 346 c Traitements specifiques au pole 347 347 348 if(mode .ge.1) then348 if(mode>=1) then 349 349 DO l=1,llm 350 350 c filtrages aux poles … … 361 361 q( i,jjp1,llm+1-l,0)=qps 362 362 enddo 363 if(mode .eq.3) then363 if(mode==3) then 364 364 dyn1=0. 365 365 dys1=0. … … 382 382 enddo 383 383 endif 384 if(mode .eq.1) then384 if(mode==1) then 385 385 c on filtre les valeurs au bord de la "grande maille pole" 386 386 dyn1=0. … … 440 440 DO j = 1,jjp1 441 441 DO i = 1,iip1 442 IF (q(i,j,l,0) .lt.0.) THEN442 IF (q(i,j,l,0)<0.) THEN 443 443 c PRINT*,'------------ BIP-----------' 444 444 c PRINT*,'Q0(',i,j,l,')=',q(i,j,l,0) … … 459 459 do j=1,jjp1 460 460 do i=1,iip1 461 if(q(i,j,l,0) .lt.qmin)461 if(q(i,j,l,0)<qmin) 462 462 , print*,'apres pentes, s0(',i,',',j,',',l,')=',q(i,j,l,0) 463 463 enddo
Note: See TracChangeset
for help on using the changeset viewer.