Changeset 1482 for trunk/LMDZ.GENERIC/libf/phystd/condense_co2.F90

Timestamp:

Oct 14, 2015, 5:05:47 PM (9 years ago)

Author:

mturbet

Message:

Implementation of the Spectral Albedo

File:

• trunk/LMDZ.GENERIC/libf/phystd/condense_co2.F90 (modified) (6 diffs)

trunk/LMDZ.GENERIC/libf/phystd/condense_co2.F90

@@ -1 +1 @@
       subroutine condense_co2(ngrid,nlayer,nq,ptimestep, &
-          pcapcal,pplay,pplev,ptsrf,pt, &
-          pphi,pdt,pdu,pdv,pdtsrf,pu,pv,pq,pdq, &
-          piceco2,psolaralb,pemisurf, &
-          pdtc,pdtsrfc,pdpsrf,pduc,pdvc, &
+          pcapcal,pplay,pplev,ptsrf,pt,                  &
+          pphi,pdt,pdu,pdv,pdtsrf,pu,pv,pq,pdq,          &
+          piceco2,albedo,pemisurf,                       &
+          albedo_bareground,albedo_co2_ice_SPECTV,       &
+          pdtc,pdtsrfc,pdpsrf,pduc,pdvc,                 &
           pdqc)
 
-      use radinc_h, only : naerkind
+      use radinc_h, only : L_NSPECTV, naerkind
       use gases_h, only: gfrac, igas_co2
       use radii_mod, only : co2_reffrad
       use aerosol_mod, only : iaero_co2
-      USE surfdat_h, only: albedodat, albedice, emisice, emissiv
+      USE surfdat_h, only: emisice, emissiv
       USE comgeomfi_h, only: lati
       USE tracer_h, only: noms, rho_co2
@@ -32 +33 @@
 !     ptsrf(ngrid)          Surface temperature
 !     
-!     pdt(ngrid,nlayer)   Time derivative before condensation/sublimation of pt
+!     pdt(ngrid,nlayer)     Time derivative before condensation/sublimation of pt
 !     pdtsrf(ngrid)         Time derivative before condensation/sublimation of ptsrf
 !     pqsurf(ngrid,nq)      Sedimentation flux at the surface (kg.m-2.s-1)
@@ -44 +45 @@
 !     Both
 !     ----
-!     piceco2(ngrid)        CO2 ice at the surface (kg/m2)
-!     psolaralb(ngrid)      Albedo at the surface
-!     pemisurf(ngrid)       Emissivity of the surface
+!     piceco2(ngrid)               CO2 ice at the surface (kg/m2)
+!     albedo(ngrid,L_NSPECTV)      Spectral albedo at the surface
+!     pemisurf(ngrid)              Emissivity of the surface
 !
 !     Authors
@@ -77 +78 @@
       REAL,INTENT(IN) :: pq(ngrid,nlayer,nq)
       REAL,INTENT(IN) :: pdq(ngrid,nlayer,nq)
+      REAL,INTENT(IN) :: albedo_bareground(ngrid)
+      REAL,INTENT(IN) :: albedo_co2_ice_SPECTV(L_NSPECTV)
       REAL,INTENT(INOUT) :: piceco2(ngrid)
-      REAL,INTENT(OUT) :: psolaralb(ngrid)
+      REAL,INTENT(OUT) :: albedo(ngrid,L_NSPECTV)
       REAL,INTENT(OUT) :: pemisurf(ngrid)
       REAL,INTENT(OUT) :: pdtc(ngrid,nlayer)
@@ -90 +93 @@
 !     Local variables
 
-      INTEGER l,ig,icap,ilay,it,iq
+      INTEGER l,ig,icap,ilay,it,iq,nw
 
       REAL reffrad(ngrid,nlayer) ! radius (m) of the co2 ice particles
@@ -461 +464 @@
             piceco2(ig)=0.
          endif
-         if (piceco2(ig).gt.0) then
-            psolaralb(ig) = albedice(icap)
+         if (piceco2(ig) .gt. 1.) then  ! CO2 Albedo condition changed to 1 mm coverage. Change by MT2015.
+            DO nw=1,L_NSPECTV
+               albedo(ig,nw) = albedo_co2_ice_SPECTV(nw)
+            ENDDO
             emisref(ig)   = emisice(icap)
          else
-            psolaralb(ig) = albedodat(ig)
+            DO nw=1,L_NSPECTV
+               albedo(ig,nw) = albedo_bareground(ig) ! Note : If you have some water, it will be taken into account in the "hydrol" routine.
+            ENDDO
             emisref(ig)   = emissiv
             pemisurf(ig)  = emissiv