Changeset 1414 for LMDZ4/trunk/libf/cosp/icarus.F
- Timestamp:
- Jul 15, 2010, 5:21:22 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/cosp/icarus.F
r1279 r1414 33 33 &) 34 34 35 ! Id: icarus.f,v 4.0 2009/02/12 13:59:20hadmw Exp $35 !$Id: icarus.f,v 4.1 2010/05/27 16:30:18 hadmw Exp $ 36 36 37 37 ! *****************************COPYRIGHT**************************** … … 72 72 ! *****************************COPYRIGHT******************************* 73 73 ! *****************************COPYRIGHT******************************* 74 ! *****************************COPYRIGHT******************************* 74 75 75 76 implicit none … … 234 235 REAL attropmin (npoints) 235 236 REAL atmax(npoints) 236 REAL atmin(npoints)237 237 REAL btcmin(npoints) 238 238 REAL transmax(npoints) … … 356 356 do j=1,npoints 357 357 ptrop(j)=5000. 358 atmin(j) = 400.359 358 attropmin(j) = 400. 360 359 atmax(j) = 0. … … 373 372 itrop(j)=ilev 374 373 end if 375 if (at(j,ilev) .gt. atmax(j)) atmax(j)=at(j,ilev)376 if (at(j,ilev) .lt. atmin(j)) atmin(j)=at(j,ilev)377 374 enddo 378 375 12 continue 376 377 do 13 ilev=1,nlev 378 do j=1,npoints 379 if (at(j,ilev) .gt. atmax(j) .and. 380 & ilev .ge. itrop(j)) atmax(j)=at(j,ilev) 381 enddo 382 13 continue 379 383 380 384 end if
Note: See TracChangeset
for help on using the changeset viewer.