Changeset 2565

Timestamp:

Jun 10, 2016, 4:01:28 PM (8 years ago)

Author:

jyg

Message:

Some small improvements to the 1D model:
(1) initialization of pbl_tke (-> 1+1=2);
(2) bug fix concerning nudge_tsoil=y;
(3) possibility of a constant moisture convergence
when forcing_type=0.

Location:

LMDZ5/trunk/libf/phylmd

Files:

• dyn1d/1DUTILS.h (modified) (3 diffs)
• dyn1d/1D_interp_cases.h (modified) (1 diff)
• dyn1d/lmdz1d.F90 (modified) (7 diffs)
• physiq_mod.F90 (modified) (2 diffs)

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

@@ -3236 +3236 @@
      &           thlpcar(nlev_max),tracer(nlev_max,ntrac)
 
+        real height1(nlev_max)
+
         integer, parameter :: ilesfile=1
         integer :: ierr,k,itrac,nt1,nt2
@@ -3245 +3247 @@
         read (ilesfile,*) kmax
         do k=1,kmax
-          read (ilesfile,*) height(k),thlprof(k),qtprof (k),               &
+          read (ilesfile,*) height1(k),thlprof(k),qtprof (k),               &
      &                      uprof (k),vprof  (k),e12prof(k)
         enddo
@@ -3261 +3263 @@
           read (ilesfile,*) height(k),ugprof(k),vgprof(k),wfls(k),         &
      &                      dqtdxls(k),dqtdyls(k),dqtdtls(k),thlpcar(k)
+        end do
+        do k=1,kmax
+          if (height(k) .ne. height1(k)) then
+            print *, 'fichiers prof.inp et lscale.inp incompatibles :'
+            print *, 'les niveaux different : ',k,height1(k), height(k)
+            stop
+          endif
         end do
         close(ilesfile)

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

@@ -2 +2 @@
 ! $Id$
 !
+!---------------------------------------------------------------------
+! Forcing_LES case: constant dq_dyn
+!---------------------------------------------------------------------
+      if (forcing_LES) then
+        DO l = 1,llm
+          d_q_adv(l,1) = dq_dyn(l,1)
+        ENDDO
+      endif ! forcing_LES
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !---------------------------------------------------------------------
 ! Interpolation forcing in time and onto model levels

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

@@ -33 +33 @@
    USE mod_1D_cases_read
    USE mod_1D_amma_read
-   USE print_control_mod, ONLY: prt_level
+   USE print_control_mod, ONLY: lunout, prt_level
    USE iniphysiq_mod, ONLY: iniphysiq
    USE mod_const_mpi, ONLY: comm_lmdz
@@ -54 +54 @@
 #include "fcg_gcssold.h"
 !!!#include "fbforcing.h"
+#include "compbl.h"
 
 !=====================================================================
@@ -242 +243 @@
 !---------------------------------------------------------------------
       integer :: k,l,i,it=1,mxcalc
+      integer :: nsrf
       integer jcode
       INTEGER read_climoz 
@@ -738 +740 @@
         t_ancien(1,:)=temp(:)
         q_ancien(1,:)=q(:,1)
-        pbl_tke(:,:,:)=1.e-8
-        wake_delta_pbl_tke(:,:,:)=0.
+!jyg<
+!!        pbl_tke(:,:,:)=1.e-8
+        pbl_tke(:,:,:)=0.
+        pbl_tke(:,2,:)=1.e-2
+        PRINT *, ' pbl_tke dans lmdz1d '
+       DO nsrf = 1,4
+         PRINT *,'pbl_tke(1,:,',nsrf,') ',pbl_tke(1,:,nsrf)
+       ENDDO
+
+!>jyg
 
         rain_fall=0.
@@ -762 +772 @@
         wake_deltaq = 0.
         wake_deltat = 0.
-        wake_delta_pbl_TKE = 0.
+        wake_delta_pbl_TKE(:,:,:) = 0.
         delta_tsurf = 0.
         wake_fip = 0.
@@ -791 +801 @@
 ! run_off_lic_0,pbl_tke(:,1:klev,nsrf), zmax0,f0,sig1,w01
 ! wake_deltat,wake_deltaq,wake_s,wake_cstar,wake_fip,wake_delta_pbl_tke(:,1:klev,nsrf)
+!
+! NB2: The content of the startphy.nc file depends on some flags defined in
+! the ".def" files. However, since conf_phys is not called in lmdz1d.F90, these flags have 
+! to be set at some arbitratry convenient values.
 !------------------------------------------------------------------------
 !Al1 =============== restart option ==========================
         if (.not.restart) then
+          iflag_pbl = 5
           call phyredem ("startphy.nc")
         else
@@ -1075 +1090 @@
 !---------------------------------------------------------------------
 
-      IF (nudge_tsoil) THEN
+      IF (nudge_tsoil .AND. .NOT. lastcall) THEN
        ftsoil(1,isoil_nudge,:) = ftsoil(1,isoil_nudge,:)                     &
      &  -timestep/tau_soil_nudge*(ftsoil(1,isoil_nudge,:)-Tsoil_nudge)

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -1165 +1165 @@
 
        CALL phyetat0 ("startphy.nc",clesphy0,tabcntr0)
+!jyg<
        IF (klon_glo==1) THEN
-          coefh=0. ; coefm=0. ; pbl_tke=0.
-          coefh(:,2,:)=1.e-2 ; coefm(:,2,:)=1.e-2 ; pbl_tke(:,2,:)=1.e-2
-          PRINT*,'FH WARNING : lignes a supprimer'
+          pbl_tke(:,:,is_ave) = 0.
+          DO nsrf=1,nbsrf
+            DO k = 1,klev+1
+                 pbl_tke(:,k,is_ave) = pbl_tke(:,k,is_ave) &
+                     +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf)
+            ENDDO
+          ENDDO
+!>jyg
        ENDIF
        !IM begin
@@ -2776 +2782 @@
 !jyg<
 !
-          CALL alpale_th( dtime, lmax_th, t_seri,  &
+          CALL alpale_th( dtime, lmax_th, t_seri, cell_area,  &
                           cin, s2, n2,  &
                           ale_bl_trig, ale_bl_stat, ale_bl,  &