Changeset 2455 for LMDZ5

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).

Location:

LMDZ5/trunk/libf/phylmd

Files:

• compbl.h (modified) (1 diff)
• conf_phys_m.F90 (modified) (6 diffs)
• dyn1d/1DUTILS.h (modified) (4 diffs)
• dyn1d/lmdz1d.F90 (modified) (1 diff)
• pbl_surface_mod.F90 (modified) (3 diffs)
• surf_ocean_mod.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/compbl.h

@@ -2 +2 @@
       ! $Header$
       !
-!jyg+nrlmd<
-!!!      integer iflag_pbl
-!!!      common/compbl/iflag_pbl
-      integer iflag_pbl,iflag_pbl_split
-      common/compbl/iflag_pbl,iflag_pbl_split
-!>jyg+nrlmd
+!jyg+al1<
+!!      integer iflag_pbl,iflag_pbl_split
+!!      common/compbl/iflag_pbl,iflag_pbl_split
+      integer iflag_pbl, iflag_pbl_split, iflag_order2_sollw
+      common/compbl/iflag_pbl, iflag_pbl_split, iflag_order2_sollw
+!>jyg+al1
 !$OMP THREADPRIVATE(/compbl/)

LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

@@ -179 +179 @@
     INTEGER,SAVE :: iflag_pbl_omp,lev_histhf_omp,lev_histday_omp,lev_histmth_omp
     INTEGER,SAVE :: iflag_pbl_split_omp
+    INTEGER,SAVE :: iflag_order2_sollw_omp
     Integer, save :: lev_histins_omp, lev_histLES_omp 
     INTEGER, SAVE :: lev_histdayNMC_omp
@@ -1287 +1288 @@
     call getin('iflag_pbl_split',iflag_pbl_split_omp)
     !
+    !Config Key  = iflag_order2_sollw
+    !Config Desc =
+    !Config Def  = 0
+    !Config Help =
+    !
+    iflag_order2_sollw_omp = 0
+    call getin('iflag_order2_sollw',iflag_order2_sollw_omp)
+    !
     !Config Key  = iflag_thermals
     !Config Desc =
@@ -1703 +1712 @@
     f_gust_wk_omp = 0.
     call getin('f_gust_wk',f_gust_wk_omp)
+    !
+    !Config Key  = iflag_z0_oce
+    !Config Desc = 0 (z0h=z0m), 1 (diff. equ. for z0h and z0m), -1 (z0m=z0h=z0min)
+    !Config Def  = 0   ! z0h = z0m
+    !Config Help =
     !
     iflag_z0_oce_omp=0
@@ -2025 +2039 @@
     iflag_pbl = iflag_pbl_omp
     iflag_pbl_split = iflag_pbl_split_omp
+    iflag_order2_sollw = iflag_order2_sollw_omp
     lev_histhf = lev_histhf_omp
     lev_histday = lev_histday_omp
@@ -2224 +2239 @@
        END IF
     END IF
-
-    ! Flag_aerosol cannot be to zero if we are in coupled mode for aerosol 
-    IF (aerosol_couple .AND. flag_aerosol .eq. 0 ) THEN 
-       CALL abort_physic('conf_phys', 'flag_aerosol cannot be to zero if aerosol_couple=y ', 1) 
-    ENDIF
-
-    ! flag_aerosol need to be different to zero if ok_cdnc is activated
-    IF (ok_cdnc .AND. flag_aerosol .eq. 0) THEN 
-       CALL abort_physic('conf_phys', 'flag_aerosol cannot be to zero if ok_cdnc is activated ', 1) 
-    ENDIF
 
     ! ok_cdnc must be set to y if ok_aie is activated
@@ -2357 +2362 @@
     write(lunout,*)' iflag_pbl = ', iflag_pbl
     write(lunout,*)' iflag_pbl_split = ', iflag_pbl_split
+    write(lunout,*)' iflag_order2_sollw = ', iflag_order2_sollw
     write(lunout,*)' iflag_thermals = ', iflag_thermals
     write(lunout,*)' iflag_thermals_ed = ', iflag_thermals_ed

LMDZ5/trunk/libf/phylmd/dyn1d/1DUTILS.h

@@ -4344 +4344 @@
 !
    INTEGER k,i
-   REAL zx_qs, rh, tnew, d_rh
+   REAL zx_qs, rh, tnew, d_rh, rhnew
 
 ! Declaration des constantes et des fonctions thermodynamiques
@@ -4361 +4361 @@
         print *,'temp ',t
         print *,'hum ',q
+!
         DO k = 1,klev
          DO i = 1,klon
-!!           IF (paprs(i,1)-pplay(i,k) .GT. 10000.) THEN
+           IF (paprs(i,1)-pplay(i,k) .GT. 10000.) THEN
             IF (t(i,k).LT.RTT) THEN
                zx_qs = qsats(t(i,k))/(pplay(i,k))
@@ -4374 +4375 @@
             d_rh = 1./tau*(rh_targ(i,k)-rh)
 !
-            tnew = t(i,k)+d_t(i,k)
+            tnew = t(i,k)+d_t(i,k)*dtime
+!jyg<
+!   Formule pour q :
+!                         d_q = (1/tau) [rh_targ*qsat(T_new) - q] 
+!
+!  Cette formule remplace d_q = (1/tau) [rh_targ - rh] qsat(T_new)
+!   qui n'était pas correcte.
+!
             IF (tnew.LT.RTT) THEN
                zx_qs = qsats(tnew)/(pplay(i,k))
@@ -4380 +4388 @@
                zx_qs = qsatl(tnew)/(pplay(i,k))
             ENDIF
-            d_q(i,k) = d_q(i,k) + d_rh*zx_qs
-!
-            print *,' k,d_t,rh,d_rh,d_q ',    &
-                      k,d_t(i,k),rh,d_rh,d_q(i,k)
-!!           ENDIF
+!!            d_q(i,k) = d_q(i,k) + d_rh*zx_qs
+            d_q(i,k) = d_q(i,k) + (1./tau)*(rh_targ(i,k)*zx_qs - q(i,k))
+            rhnew = (q(i,k)+d_q(i,k)*dtime)/zx_qs
+!
+            print *,' k,d_t,rh,d_rh,rhnew,d_q ',    &
+                      k,d_t(i,k),rh,d_rh,rhnew,d_q(i,k)
+           ENDIF
 !
          ENDDO

LMDZ5/trunk/libf/phylmd/dyn1d/lmdz1d.F90

@@ -932 +932 @@
 !   Call physiq :
 !---------------------------------------------------------------------
-
        call physiq(ngrid,llm, &
-            firstcall,lastcall,timestep, &
-            plev,play,phi,phis,presnivs, &
-            u,v, rot, temp,q,omega2, &
-            du_phys,dv_phys,dt_phys,dq,dpsrf)
-        firstcall=.false.
+                    firstcall,lastcall,timestep, &
+                    plev,play,phi,phis,presnivs, &
+                    u,v, rot, temp,q,omega2, &
+                    du_phys,dv_phys,dt_phys,dq,dpsrf)
+                firstcall=.false.
 
 !---------------------------------------------------------------------

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

LMDZ5/trunk/libf/phylmd/surf_ocean_mod.F90

@@ -36 +36 @@
 
     include "clesphys.h"
-    ! for cycle_diurne
+    ! for cycle_diurne and for iflag_z0_oce==-1 (prescribed z0)
 
 ! Input variables
@@ -244 +244 @@
        z0h(i)=0.4*14e-6 / SQRT(cdragm(i) * tmp)
     ENDDO
+ELSE IF (iflag_z0_oce==-1) THEN
+    DO i = 1, knon
+       z0m(i) = z0min
+       z0h(i) = z0min
+    ENDDO
 ELSE
        CALL abort_physic(modname,'version non prevue',1)