Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/dyn3d/iniacademic.F90
- Timestamp:
- Aug 2, 2024, 2:12:03 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d/iniacademic.F90
r5136 r5158 152 152 153 153 PRINT*, 'relief=', minval(relief), maxval(relief), 'g=', g 154 doj = 1, jjp1155 doi = 1, iip1154 DO j = 1, jjp1 155 DO i = 1, iip1 156 156 phis((j - 1) * iip1 + i) = g * relief(i, j) 157 157 enddo … … 264 264 ! CALL writefield('theta_eq',tetajl) 265 265 266 dol = 1, llm267 doj = 1, jjp1268 doi = 1, iip1266 DO l = 1, llm 267 DO j = 1, jjp1 268 DO i = 1, iip1 269 269 ij = (j - 1) * iip1 + i 270 270 tetarappel(ij, l) = tetajl(j, l) … … 303 303 IF (planet_type=="earth") THEN 304 304 ! Earth: first two tracers will be water 305 doiq = 1, nqtot305 DO iq = 1, nqtot 306 306 q(:, :, iq) = 0. 307 307 IF(tracers(iq)%name == addPhase('H2O', 'g')) q(:, :, iq) = 1.e-10 … … 344 344 zz = ran1(idum) 345 345 idum = 0 346 dol = 1, llm347 doij = iip2, ip1jm346 DO l = 1, llm 347 DO ij = iip2, ip1jm 348 348 teta(ij, l) = teta(ij, l) * (1. + tetanoise * ran1(idum)) 349 349 enddo … … 351 351 352 352 ! maintain periodicity in longitude 353 dol = 1, llm354 doij = 1, ip1jmp1, iip1353 DO l = 1, llm 354 DO ij = 1, ip1jmp1, iip1 355 355 teta(ij + iim, l) = teta(ij, l) 356 356 enddo
Note: See TracChangeset
for help on using the changeset viewer.