Changeset 888 for LMDZ4/trunk/libf/phylmd/phyredem.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/phyredem.F (modified) (6 diffs)

LMDZ4/trunk/libf/phylmd/phyredem.F

@@ -5 +5 @@
       SUBROUTINE phyredem (fichnom,dtime,radpas,ocean,
      .           rlat_p,rlon_p, pctsrf_p,tsol_p,
-     .           albedo_p, alblw_p,
+     .           alb1_p, alb2_p,
      .           rain_fall_p, snow_fall_p,solsw_p, sollw_p,
      .           radsol_p,zmea_p,zstd_p,zsig_p,
@@ -46 +46 @@
       REAL qsol_p(klon)
       REAL snow_p(klon,nbsrf)
-      REAL albedo_p(klon,nbsrf)
-cIM BEG
-      REAL alblw_p(klon,nbsrf)
-cIM END
+      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 rain_fall_p(klon)
@@ -80 +78 @@
       REAL qsol(klon_glo)
       REAL snow(klon_glo,nbsrf)
-      REAL albedo(klon_glo,nbsrf)
-      REAL alblw(klon_glo,nbsrf)
+      REAL alb1(klon_glo,nbsrf)
+      REAL alb2(klon_glo,nbsrf)
       REAL evap(klon_glo,nbsrf)
       REAL rain_fall(klon_glo)
@@ -147 +145 @@
       call Gather( qsol_p,qsol)
       call Gather( snow_p,snow)
-      call Gather( albedo_p,albedo)
-      call Gather( alblw_p,alblw)
+      call Gather( alb1_p,alb1)
+      call Gather( alb2_p,alb2)
       call Gather( evap_p,evap)
       call Gather( radsol_p,radsol)
@@ -481 +479 @@
         ENDIF
 #ifdef NC_DOUBLE
-      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,albedo(1,nsrf))
-#else
-      ierr = NF_PUT_VAR_REAL (nid,nvarid,albedo(1,nsrf))
-#endif
-      ENDDO
-
-cIM BEG albedo LW
+      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,alb1(1,nsrf))
+#else
+      ierr = NF_PUT_VAR_REAL (nid,nvarid,alb1(1,nsrf))
+#endif
+      ENDDO
+
         DO nsrf = 1, nbsrf
         IF (nsrf.LE.99) THEN
@@ -505 +502 @@
         ENDIF
 #ifdef NC_DOUBLE
-      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,alblw(1,nsrf))
-#else
-      ierr = NF_PUT_VAR_REAL (nid,nvarid,alblw(1,nsrf))
-#endif
-      ENDDO
-cIM END albedo LW
+      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,alb2(1,nsrf))
+#else
+      ierr = NF_PUT_VAR_REAL (nid,nvarid,alb2(1,nsrf))
+#endif
+      ENDDO
+c
 c
       DO nsrf = 1, nbsrf