Changeset 2455 for LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90

Timestamp:

Feb 28, 2016, 12:30:51 PM (9 years ago)

Author:

jyg

Message:

Implementation of a second order distribution on

sub-surfaces of longwave net radiance (Alain
Lahellec).

Introduction of the flag iflag_order2_sollw:

=1 ==> second order distribution.
=0 ==> linear distribution (Default).

File:

• LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90

@@ -522 +522 @@
 !albedo SB <<<
     REAL, DIMENSION(klon)              :: ztsol
+    REAL, DIMENSION(klon)              :: meansqT ! mean square deviation of subsurface temperatures
     REAL, DIMENSION(klon)              :: alb_m  ! mean albedo for whole SW interval
     REAL, DIMENSION(klon,klev)         :: y_d_t, y_d_q, y_d_t_diss
@@ -1106 +1107 @@
        ENDDO
     ENDDO
+!
+!<al1: second order corrections
+!- net = dwn -up; up=sig( T4 + 4sum%T3T' + 6sum%T2T'2 +...)
+   IF (iflag_order2_sollw == 1) THEN
+    meansqT(:) = 0. ! as working buffer
+    DO nsrf = 1, nbsrf
+     DO i = 1, klon
+      meansqT(i) = meansqT(i)+(ts(i,nsrf)-ztsol(i))**2 *pctsrf(i,nsrf)
+     END DO
+    END DO
+    DO nsrf = 1, nbsrf
+     DO i = 1, klon
+      sollw(i,nsrf) = sollw(i,nsrf) &
+                + 6.0*RSIGMA*ztsol(i)**2 *(meansqT(i)-(ztsol(i)-ts(i,nsrf))**2)
+     ENDDO
+    ENDDO
+   ENDIF   ! iflag_order2_sollw == 1
+!>al1
 
 !****************************************************************************************
@@ -2811 +2830 @@
        END DO
     END DO
+!
+!<al1 order 2 correction to zxtsol, for radiation computations (main atm effect of Ts)
+   IF (iflag_order2_sollw == 1) THEN
+    meansqT(:) = 0. ! as working buffer
+    DO nsrf = 1, nbsrf
+     DO i = 1, klon
+      meansqT(i) = meansqT(i)+(ts(i,nsrf)-zxtsol(i))**2 *pctsrf(i,nsrf)
+     END DO
+    END DO
+    zxtsol(:) = zxtsol(:)+1.5*meansqT(:)/zxtsol(:)
+   ENDIF   ! iflag_order2_sollw == 1
+!>al1
           
 !!! jyg le 07/02/2012