Changeset 2220 for LMDZ5/branches/testing/libf/phylmd/pbl_surface_mod.F90

Timestamp:

Mar 3, 2015, 2:41:13 PM (10 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes -r2186:2216 into testing branch

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phylmd/pbl_surface_mod.F90 (modified) (12 diffs)

LMDZ5/branches/testing/libf/phylmd/pbl_surface_mod.F90

@@ -33 +33 @@
   REAL, ALLOCATABLE, DIMENSION(:), PRIVATE, SAVE     :: fder   ! flux drift
   !$OMP THREADPRIVATE(fder)
-  REAL, ALLOCATABLE, DIMENSION(:,:), PRIVATE, SAVE   :: snow   ! snow at surface
+  REAL, ALLOCATABLE, DIMENSION(:,:), PUBLIC, SAVE   :: snow   ! snow at surface
   !$OMP THREADPRIVATE(snow)
   REAL, ALLOCATABLE, DIMENSION(:,:), PRIVATE, SAVE   :: qsurf  ! humidity at surface
@@ -172 +172 @@
        debut,     lafin,                              &
        rlon,      rlat,      rugoro,   rmu0,          &
-       zsig,      sollwd_m,  pphi,     cldt,          &
+       zsig,      lwdown_m,  pphi,     cldt,          &
        rain_f,    snow_f,    solsw_m,  sollw_m,       &
        t,         q,         u,        v,             &
@@ -182 +182 @@
        pplay,     paprs,     pctsrf,                  &
        ts,        alb1, alb2,ustar, u10m, v10m,wstar, &
-       lwdown_m,  cdragh,    cdragm,   zu1,    zv1,   &
+       cdragh,    cdragm,   zu1,    zv1,              &
        alb1_m,    alb2_m,    zxsens,   zxevap,        &
        alb3_lic,  runoff,    snowhgt,   qsnow,     to_ice,    sissnow,  &
@@ -327 +327 @@
 ! Martin
     REAL, DIMENSION(klon),        INTENT(IN)        :: zsig    ! slope
-    REAL, DIMENSION(klon),        INTENT(IN)        :: sollwd_m ! net longwave radiation at mean s    
+    REAL, DIMENSION(klon),        INTENT(IN)        :: lwdown_m ! downward longwave radiation at mean s    
     REAL, DIMENSION(klon),        INTENT(IN)        :: cldt    ! total cloud fraction
     REAL, DIMENSION(klon,klev),   INTENT(IN)        :: pphi    ! geopotential (m2/s2)
@@ -367 +367 @@
 ! Output variables
 !****************************************************************************************
-    REAL, DIMENSION(klon),        INTENT(OUT)       :: lwdown_m   ! Downcoming longwave radiation
     REAL, DIMENSION(klon),        INTENT(OUT)       :: cdragh     ! drag coefficient for T and Q
     REAL, DIMENSION(klon),        INTENT(OUT)       :: cdragm     ! drag coefficient for wind
@@ -780 +779 @@
     ! Martin
     REAL, DIMENSION(klon, nbsrf)       :: sollwd ! net longwave radiation at surface
-    REAL, DIMENSION(klon)              :: ysollwd
     REAL, DIMENSION(klon)              :: ytoice
     REAL, DIMENSION(klon)              :: ysnowhgt, yqsnow, ysissnow, yrunoff
@@ -855 +853 @@
 ! 2a) Initialization of all argument variables with INTENT(OUT)
 !****************************************************************************************
- lwdown_m(:)=0.
  cdragh(:)=0. ; cdragm(:)=0.
  zu1(:)=0. ; zv1(:)=0.
@@ -938 +935 @@
     ! Martin
     ysnowhgt = 0.0; yqsnow = 0.0     ; yrunoff = 0.0   ; ytoice =0.0
-    yalb3_new = 0.0  ; ysissnow = 0.0  ; ysollwd = 0.0
+    yalb3_new = 0.0  ; ysissnow = 0.0 
     ypphi = 0.0   ; ycldt = 0.0      ; yrmu0 = 0.0
     ! Martin
@@ -1109 +1106 @@
        DO i = 1, klon
           sollw(i,nsrf) = sollw_m(i) + 4.0*RSIGMA*ztsol(i)**3 * (ztsol(i)-ts(i,nsrf))
-          ! Martin
-          sollwd(i,nsrf)= sollwd_m(i)
-          ! Martin
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+!         ! Martin
+! Apparently introduced for sisvat but not used
+!         sollwd(i,nsrf)= sollwd_m(i)
+!         ! Martin
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
           solsw(i,nsrf) = solsw_m(i) * (1.-alb(i,nsrf)) / (1.-alb_m(i))
        ENDDO
     ENDDO
 
-
-! Downwelling longwave radiation at mean surface
-    lwdown_m(:) = 0.0
-    DO i = 1, klon
-       lwdown_m(i) = sollw_m(i) + RSIGMA*ztsol(i)**4
-    ENDDO
 
 !****************************************************************************************
@@ -1180 +1176 @@
           yagesno(j) = agesno(i,nsrf)
           yfder(j)   = fder(i)
+          ylwdown(j) = lwdown_m(i)
           ysolsw(j)  = solsw(i,nsrf)
           ysollw(j)  = sollw(i,nsrf)
@@ -1703 +1700 @@
      
        CASE(is_ter)
-          ! ylwdown : to be removed, calculation is now done at land surface in surf_land
-          ylwdown(:)=0.0
-          DO i=1,knon
-             ylwdown(i)=lwdown_m(ni(i))
-          END DO
+!          print*,"DEBUGTS",yts(knon/2),ylwdown(knon/2)
           CALL surf_land(itap, dtime, date0, jour, knon, ni,&
                rlon, rlat, &
@@ -1746 +1739 @@
           CALL surf_landice(itap, dtime, knon, ni, &
                rlon, rlat, debut, lafin, &
-               yrmu0, ysollwd, yalb, ypphi(:,1), &
+               yrmu0, ylwdown, yalb, ypphi(:,1), &
                ysolsw, ysollw, yts, ypplay(:,1), &
                ycdragh, ycdragm, yrain_f, ysnow_f, yt(:,1), yq(:,1),&