Changeset 888 for LMDZ4/trunk/libf/phylmd/phyetat0.F

Timestamp:

Feb 4, 2008, 5:24:28 PM (17 years ago)

Author:

Laurent Fairhead

Message:

Modifications sur l'albedo JG
LF

File:

• LMDZ4/trunk/libf/phylmd/phyetat0.F (modified) (12 diffs)

LMDZ4/trunk/libf/phylmd/phyetat0.F

@@ -7 +7 @@
      .           rlat_p, rlon_p, pctsrf_p, tsol_p,
      .           ocean_in, ok_veget_in,
-     .           albe_p, alblw_p,
+     .           alb1_p, alb2_p,
      .           rain_fall_p, snow_fall_p,solsw_p, sollw_p,
      .           radsol_p,clesphy0,
@@ -53 +53 @@
       REAL qsol_p(klon)
       REAL snow_p(klon,nbsrf)
-      REAL albe_p(klon,nbsrf)
-cIM BEG alblw
-      REAL alblw_p(klon,nbsrf)
-cIM END alblw
+      REAL alb1_p(klon,nbsrf)   ! albedo in visible SW interval
+      REAL alb2_p(klon,nbsrf)   ! albedo in near IR interval
       REAL evap_p(klon,nbsrf)
       REAL radsol_p(klon)
@@ -93 +91 @@
       REAL qsol(klon_glo)
       REAL snow(klon_glo,nbsrf)
-      REAL albe(klon_glo,nbsrf)
-      REAL alblw(klon_glo,nbsrf)
+      REAL alb1(klon_glo,nbsrf)
+      REAL alb2(klon_glo,nbsrf)
       REAL evap(klon_glo,nbsrf)
       REAL radsol(klon_glo)
@@ -728 +726 @@
       ENDIF
 c
-c Lecture de albedo au sol:
+c Lecture de albedo de l'interval visible au sol:
 c
       ierr = NF_INQ_VARID (nid, "ALBE", nvarid)
@@ -746 +744 @@
            ENDIF
 #ifdef NC_DOUBLE
-           ierr = NF_GET_VAR_DOUBLE(nid, nvarid, albe(1,nsrf))
-#else
-           ierr = NF_GET_VAR_REAL(nid, nvarid, albe(1,nsrf))
+           ierr = NF_GET_VAR_DOUBLE(nid, nvarid, alb1(1,nsrf))
+#else
+           ierr = NF_GET_VAR_REAL(nid, nvarid, alb1(1,nsrf))
 #endif
            IF (ierr.NE.NF_NOERR) THEN
@@ -757 +755 @@
            xmax = -1.0E+20
            DO i = 1, klon_glo
-              xmin = MIN(albe(i,nsrf),xmin)
-              xmax = MAX(albe(i,nsrf),xmax)
+              xmin = MIN(alb1(i,nsrf),xmin)
+              xmax = MAX(alb1(i,nsrf),xmax)
            ENDDO
            PRINT*,'Albedo du sol ALBE**:', nsrf, xmin, xmax
@@ -766 +764 @@
          PRINT*, '          J ignore donc les autres ALBE**'
 #ifdef NC_DOUBLE
-         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, albe(1,1))
-#else
-         ierr = NF_GET_VAR_REAL(nid, nvarid, albe(1,1))
+         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, alb1(1,1))
+#else
+         ierr = NF_GET_VAR_REAL(nid, nvarid, alb1(1,1))
 #endif
          IF (ierr.NE.NF_NOERR) THEN
@@ -777 +775 @@
          xmax = -1.0E+20
          DO i = 1, klon_glo
-            xmin = MIN(albe(i,1),xmin)
-            xmax = MAX(albe(i,1),xmax)
+            xmin = MIN(alb1(i,1),xmin)
+            xmax = MAX(alb1(i,1),xmax)
          ENDDO
          PRINT*,'Neige du sol <ALBE>', xmin, xmax
          DO nsrf = 2, nbsrf
          DO i = 1, klon_glo
-            albe(i,nsrf) = albe(i,1)
-         ENDDO
-         ENDDO
-      ENDIF
-
-c
-c Lecture de albedo au sol LW:
+            alb1(i,nsrf) = alb1(i,1)
+         ENDDO
+         ENDDO
+      ENDIF
+
+c
+c Lecture de albedo au sol dans l'interval proche infra-rouge:
 c
       ierr = NF_INQ_VARID (nid, "ALBLW", nvarid)
@@ -798 +796 @@
          DO nsrf = 1, nbsrf
            DO i = 1, klon_glo
-             alblw(i,nsrf) = albe(i,nsrf)
+             alb2(i,nsrf) = alb1(i,nsrf)
            ENDDO
          ENDDO
@@ -805 +803 @@
          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))
+         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, alb2(1,1))
+#else
+         ierr = NF_GET_VAR_REAL(nid, nvarid, alb2(1,1))
 #endif
          IF (ierr.NE.NF_NOERR) THEN
@@ -816 +814 @@
          xmax = -1.0E+20
          DO i = 1, klon_glo
-            xmin = MIN(alblw(i,1),xmin)
-            xmax = MAX(alblw(i,1),xmax)
+            xmin = MIN(alb2(i,1),xmin)
+            xmax = MAX(alb2(i,1),xmax)
          ENDDO
          PRINT*,'Neige du sol <ALBLW>', xmin, xmax
          DO nsrf = 2, nbsrf
          DO i = 1, klon_glo
-            alblw(i,nsrf) = alblw(i,1)
+            alb2(i,nsrf) = alb2(i,1)
          ENDDO
          ENDDO
@@ -1583 +1581 @@
       call Scatter( qsol,qsol_p)
       call Scatter( snow,snow_p)
-      call Scatter( albe,albe_p)
-      call Scatter( alblw,alblw_p)
+      call Scatter( alb1,alb1_p)
+      call Scatter( alb2,alb2_p)
       call Scatter( evap,evap_p)
       call Scatter( radsol,radsol_p)