Changeset 467 for LMDZ.3.3/branches/rel-LF/libf/phylmd/phyetat0.F

Timestamp:

Aug 6, 2003, 4:50:49 PM (21 years ago)

Author:

lmdzadmin

Message:

Modifs sur les seuils (cdrag etc...), inclusion des diagnostics ISCCP par Ionela
LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/phyetat0.F (modified) (3 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/phyetat0.F

@@ -4 +4 @@
       SUBROUTINE phyetat0 (fichnom,dtime,co2_ppm_etat0,solaire_etat0,
      .            rlat,rlon, pctsrf, tsol,tsoil,deltat,qsurf,qsol,snow,
-     .           albe, evap, rain_fall, snow_fall, solsw, sollw,
+     .           albe, alblw, evap, rain_fall, snow_fall, solsw, sollw,
      .           fder,radsol,frugs,agesno,clesphy0,
      .           zmea,zstd,zsig,zgam,zthe,zpic,zval,rugsrel,tabcntr0,
@@ -34 +34 @@
       REAL snow(klon,nbsrf)
       REAL albe(klon,nbsrf)
+cIM BEG alblw
+      REAL alblw(klon,nbsrf)
+cIM END alblw
       REAL evap(klon,nbsrf)
       REAL radsol(klon)
@@ -645 +648 @@
 
 c
+cIM BEG alblw
+c Lecture de albedo au sol LW:
+c
+      ierr = NF_INQ_VARID (nid, "ALBLW", nvarid)
+      IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, 'phyetat0: Le champ <ALBLW> est absent'
+c        PRINT*, '          Mais je vais essayer de lire ALBLW**'
+         PRINT*, '          Mais je vais prendre ALBE**'
+         DO nsrf = 1, nbsrf
+           DO i = 1, klon
+             alblw(i,nsrf) = albe(i,nsrf)
+           ENDDO
+         ENDDO
+c          IF (nsrf.GT.99) THEN
+c            PRINT*, "Trop de sous-mailles"
+c            CALL abort
+c          ENDIF
+c          WRITE(str2,'(i2.2)') nsrf
+c          ierr = NF_INQ_VARID (nid, "ALBLW"//str2, nvarid)
+c           IF (ierr.NE.NF_NOERR) THEN
+c             PRINT*, "phyetat0: Le champ <ALBLW"//str2//"> est absent"
+c             CALL abort
+c          ENDIF
+c#ifdef NC_DOUBLE
+c           ierr = NF_GET_VAR_DOUBLE(nid, nvarid, alblw(1,nsrf))
+c#else
+c           ierr = NF_GET_VAR_REAL(nid, nvarid, alblw(1,nsrf))
+c#endif
+c          IF (ierr.NE.NF_NOERR) THEN
+c            PRINT*, "phyetat0: Lecture echouee pour <ALBLW"//str2//">"
+c            CALL abort
+c          ENDIF
+c          xmin = 1.0E+20
+c          xmax = -1.0E+20
+c          DO i = 1, klon
+c             xmin = MIN(alblw(i,nsrf),xmin)
+c             xmax = MAX(alblw(i,nsrf),xmax)
+c          ENDDO
+c          PRINT*,'Albedo du sol ALBLW**:', nsrf, xmin, xmax
+c        ENDDO
+      ELSE
+         PRINT*, 'phyetat0: Le champ <ALBLW> est present'
+         PRINT*, '          J ignore donc les autres ALBLW**'
+#ifdef NC_DOUBLE
+         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, alblw(1,1))
+#else
+         ierr = NF_GET_VAR_REAL(nid, nvarid, alblw(1,1))
+#endif
+         IF (ierr.NE.NF_NOERR) THEN
+            PRINT*, "phyetat0: Lecture echouee pour <ALBLW>"
+            CALL abort
+         ENDIF
+         xmin = 1.0E+20
+         xmax = -1.0E+20
+         DO i = 1, klon
+            xmin = MIN(alblw(i,1),xmin)
+            xmax = MAX(alblw(i,1),xmax)
+         ENDDO
+         PRINT*,'Neige du sol <ALBLW>', xmin, xmax
+         DO nsrf = 2, nbsrf
+         DO i = 1, klon
+            alblw(i,nsrf) = alblw(i,1)
+         ENDDO
+         ENDDO
+      ENDIF
+
+cIM END alblw
+
+c
 c Lecture de evaporation:  
 c