Changeset 442 for LMDZ.3.3/branches
- Timestamp:
- Jan 27, 2003, 11:45:56 AM (22 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/phyetat0.F
r433 r442 3 3 c 4 4 SUBROUTINE phyetat0 (fichnom,dtime,co2_ppm_etat0,solaire_etat0, 5 . rlat,rlon, pctsrf, tsol,tsoil,deltat,qs ol,snow,5 . rlat,rlon, pctsrf, tsol,tsoil,deltat,qsurf,qsol,snow, 6 6 . albe, evap, rain_fall, snow_fall, solsw, sollw, 7 7 . fder,radsol,frugs,agesno,clesphy0, … … 30 30 REAL tsoil(klon,nsoilmx,nbsrf) 31 31 REAL deltat(klon) 32 REAL qsol(klon,nbsrf) 32 REAL qsurf(klon,nbsrf) 33 REAL qsol(klon) 33 34 REAL snow(klon,nbsrf) 34 35 REAL albe(klon,nbsrf) … … 436 437 PRINT*,'Ecart de la SST deltat:', xmin, xmax 437 438 c 438 c Lecture de l'humidite d u sol:439 c Lecture de l'humidite de l'air juste au dessus du sol: 439 440 c 440 441 ierr = NF_INQ_VARID (nid, "QS", nvarid) … … 454 455 ENDIF 455 456 #ifdef NC_DOUBLE 456 ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qs ol(1,nsrf))457 #else 458 ierr = NF_GET_VAR_REAL(nid, nvarid, qs ol(1,nsrf))457 ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsurf(1,nsrf)) 458 #else 459 ierr = NF_GET_VAR_REAL(nid, nvarid, qsurf(1,nsrf)) 459 460 #endif 460 461 IF (ierr.NE.NF_NOERR) THEN … … 465 466 xmax = -1.0E+20 466 467 DO i = 1, klon 467 xmin = MIN(qs ol(i,nsrf),xmin)468 xmax = MAX(qs ol(i,nsrf),xmax)468 xmin = MIN(qsurf(i,nsrf),xmin) 469 xmax = MAX(qsurf(i,nsrf),xmax) 469 470 ENDDO 470 PRINT*,'Humidite du sol QS**:', nsrf, xmin, xmax471 PRINT*,'Humidite pres du sol QS**:', nsrf, xmin, xmax 471 472 ENDDO 472 473 ELSE … … 474 475 PRINT*, ' J ignore donc les autres humidites QS**' 475 476 #ifdef NC_DOUBLE 476 ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qs ol(1,1))477 #else 478 ierr = NF_GET_VAR_REAL(nid, nvarid, qs ol(1,1))477 ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsurf(1,1)) 478 #else 479 ierr = NF_GET_VAR_REAL(nid, nvarid, qsurf(1,1)) 479 480 #endif 480 481 IF (ierr.NE.NF_NOERR) THEN … … 485 486 xmax = -1.0E+20 486 487 DO i = 1, klon 487 xmin = MIN(qs ol(i,1),xmin)488 xmax = MAX(qs ol(i,1),xmax)489 ENDDO 490 PRINT*,'Humidite du sol <QS>', xmin, xmax488 xmin = MIN(qsurf(i,1),xmin) 489 xmax = MAX(qsurf(i,1),xmax) 490 ENDDO 491 PRINT*,'Humidite pres du sol <QS>', xmin, xmax 491 492 DO nsrf = 2, nbsrf 492 493 DO i = 1, klon 493 qsol(i,nsrf) = qsol(i,1) 494 ENDDO 495 ENDDO 496 ENDIF 494 qsurf(i,nsrf) = qsurf(i,1) 495 ENDDO 496 ENDDO 497 ENDIF 498 C 499 C Eau dans le sol (pour le modele de sol "bucket") 500 C 501 ierr = NF_INQ_VARID (nid, "QSOL", nvarid) 502 IF (ierr .EQ. NF_NOERR) THEN 503 #ifdef NC_DOUBLE 504 ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsol) 505 #else 506 ierr = NF_GET_VAR_REAL(nid, nvarid, qsol) 507 #endif 508 IF (ierr.NE.NF_NOERR) THEN 509 PRINT*, 'phyetat0: Lecture echouee pour <QSOL>' 510 CALL abort 511 ENDIF 512 else 513 PRINT*, 'phyetat0: Le champ <QSOL> est absent' 514 PRINT*, ' Valeur par defaut nulle' 515 qsol(:)=0. 516 c$$$ CALL abort 517 ENDIF 518 xmin = 1.0E+20 519 xmax = -1.0E+20 520 DO i = 1, klon 521 xmin = MIN(qsol(i),xmin) 522 xmax = MAX(qsol(i),xmax) 523 ENDDO 524 PRINT*,'Eau dans le sol (mm) <QSOL>', xmin, xmax 497 525 c 498 526 c Lecture de neige au sol:
Note: See TracChangeset
for help on using the changeset viewer.