Changeset 2761 for LMDZ5/trunk/libf/phylmd/calwake.F90

Timestamp:

Jan 11, 2017, 12:39:50 PM (7 years ago)

Author:

jyg

Message:

Two options introduced concerning deep convection
and wakes:

+ deep convection is allowed only when top
temperature is lower than t_top_max.
(Default: t_top_max=1000.).
+ wake number per unit area may now be
different over ocean and over land:
wdens_ref_o and wdens_ref_l.
(Default is:wdens_ref_o=8.e-12,
wdens_ref_l=8.e-12)

File:

• LMDZ5/trunk/libf/phylmd/calwake.F90 (modified) (8 diffs)

LMDZ5/trunk/libf/phylmd/calwake.F90

@@ -26 +26 @@
 
   USE dimphy
+  USE phys_state_var_mod, ONLY: pctsrf
+  USE indice_sol_mod, ONLY: is_oce
+  USE print_control_mod, ONLY: mydebug=>debug , lunout, prt_level
   IMPLICIT NONE
   ! ======================================================================
@@ -51 +54 @@
   ! ------
   REAL, DIMENSION(klon, klev),   INTENT (OUT)        :: dt_wake, dq_wake
-  INTEGER, DIMENSION(klon),      INTENT (OUT)        :: wake_k
+  REAL, DIMENSION(klon),         INTENT (OUT)        :: wake_k
   REAL, DIMENSION(klon, klev),   INTENT (OUT)        :: wake_d_deltat_gw
   REAL, DIMENSION(klon),         INTENT (OUT)        :: wake_h
@@ -69 +72 @@
   ! -----------------
   INTEGER                                            :: i, l
+  INTEGER, DIMENSION(klon)                           :: znatsurf    ! 0 if pctsrf(is_oce)>0.1; 1 else.
   REAL                                               :: aire
   REAL, DIMENSION(klon, klev)                        :: p,  pi
@@ -97 +101 @@
 
 
-  ! print *, '-> calwake, wake_s ', wake_s(1)
+  IF (prt_level >= 10) THEN
+    print *, '-> calwake, wake_s input ', wake_s(1)
+  ENDIF
 
   rdcp = 1./3.5
+
+  znatsurf(:) = 0
+  DO i = 1,klon
+    IF (pctsrf(i,is_oce) < 0.1) znatsurf(i) = 1
+  ENDDO
 
 
@@ -124 +135 @@
   END DO
 
+!----------------------------------------------------------------
+!         Initialize tendencies to zero
+!----------------------------------------------------------------
+dtls(:,:) = 0.
+dqls(:,:) = 0.
+d_deltat_gw(:,:) = 0.
+d_deltatw(:,:) = 0.
+d_deltaqw(:,:) = 0.
+d_sigmaw(:) = 0.
+d_wdens(:) = 0.
+!
+
   DO i = 1, klon
     sigd0(i) = sigd(i)
@@ -132 +155 @@
   END DO
 
-  DO i = 1, klon
-    ktopw(i) = wake_k(i)
-  END DO
+!!jyg!  DO i = 1, klon                  
+!!jyg!    ktopw(i) = NINT(wake_k(i))    
+!!jyg!  END DO                          
 
   DO i = 1, klon
@@ -178 +201 @@
   END DO
 
-  CALL wake(p, ph, pi, dtime, &
+  CALL wake(znatsurf, p, ph, pi, dtime, &
     te, qe, omgbe, &
     dtdwn, dqdwn, amdwn, amup, dta, dqa, &
@@ -254 +277 @@
     END IF
   END DO
+!
 
 !jyg<