Changeset 2455 for LMDZ5/trunk/libf/phylmd/dyn1d

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/dyn1d

Files:

• 1DUTILS.h (modified) (4 diffs)
• lmdz1d.F90 (modified) (1 diff)

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