Changeset 1311 for LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell_main.F90

Timestamp:

Feb 18, 2010, 2:14:02 PM (14 years ago)

Author:

Laurent Fairhead

Message:

Modifications to thermals for TKE transport

---

Modifications aux thermiques pour le transport de la TKE

pbl_surface_mode.F90 : ok_flux_surf=.false. seulement pour klon>1
physiq.F : option iflag_pbl=10 pour transporter la TKE avec les thermiques.
calltherm.F90 : passage de iflag_thermals_ed en argument pour thermcell_plume
thermcell_main.F90 : Appel a plusieurs version de thermcell_plume en option
thermcell_plume.F90 : plusieurs versions dans le meme fichier (temporaire)
thermcell_height.F90 : verrue pour les cas ou les thermiques montent tout

en haut

yamada4 : inclusion de la diffusion verticale en option iflag_pbl=9

+ variables anciennement common, puis save/allocatable, remises en local

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell_main.F90 (modified) (6 diffs)

LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell_main.F90

@@ -8 +8 @@
      &                  ,fm0,entr0,detr0,zqta,zqla,lmax  &
      &                  ,ratqscth,ratqsdiff,zqsatth  &
-     &                  ,r_aspect,l_mix,tau_thermals &
+     &                  ,r_aspect,l_mix,tau_thermals,iflag_thermals_ed &
      &                  ,Ale_bl,Alp_bl,lalim_conv,wght_th &
      &                  ,zmax0, f0,zw2,fraca)
@@ -55 +55 @@
       INTEGER ngrid,nlay,w2di
       real tau_thermals
+      integer iflag_thermals_ed
       real ptimestep,l_mix,r_aspect
       REAL pt(ngrid,nlay),pdtadj(ngrid,nlay)
@@ -228 +229 @@
       ENDIF
 !
-     write(lunout,*)'WARNING thermcell_main f0=max(f0,1.e-2)'
+!     write(lunout,*)'WARNING thermcell_main f0=max(f0,1.e-2)'
      do ig=1,klon
       if (prt_level.ge.20) then
@@ -380 +381 @@
       if (prt_level.ge.1) print*,'avant thermcell_plume ',lev_out
 !IM 140508   CALL thermcell_plume(ngrid,nlay,ptimestep,ztv,zthl,po,zl,rhobarz,  &
-      CALL thermcell_plume(itap,ngrid,nlay,ptimestep,ztv,zthl,po,zl,rhobarz,  &
-     &           zlev,pplev,pphi,zpspsk,alim_star,alim_star_tot,  &
-     &           lalim,f0,detr_star,entr_star,f_star,csc,ztva,  &
-     &           ztla,zqla,zqta,zha,zw2,zw_est,ztva_est,zqsatth,lmix,lmix_bis,linter &
-     &            ,lev_out,lunout1,igout)
+
+! Gestion temporaire de plusieurs appels à thermcell_plume au travers
+! de la variable iflag_thermals
+
+!      print*,'THERM thermcell_main iflag_thermals_ed=',iflag_thermals_ed
+      if (iflag_thermals_ed<=9) then
+!         print*,'THERM NOVUELLE/NOUVELLE/ANCIENNE'
+         CALL thermcell_plume(itap,ngrid,nlay,ptimestep,ztv,zthl,po,zl,rhobarz,&
+     &    zlev,pplev,pphi,zpspsk,alim_star,alim_star_tot,  &
+     &    lalim,f0,detr_star,entr_star,f_star,csc,ztva,  &
+     &    ztla,zqla,zqta,zha,zw2,zw_est,ztva_est,zqsatth,lmix,lmix_bis,linter &
+     &    ,lev_out,lunout1,igout)
+
+      elseif (iflag_thermals_ed<=19) then
+! Version d'Arnaud Jam
+!         print*,'THERM RIO et al 2010'
+         CALL thermcellV1_plume(itap,ngrid,nlay,ptimestep,ztv,zthl,po,zl,rhobarz,&
+     &    zlev,pplev,pphi,zpspsk,alim_star,alim_star_tot,  &
+     &    lalim,f0,detr_star,entr_star,f_star,csc,ztva,  &
+     &    ztla,zqla,zqta,zha,zw2,zw_est,ztva_est,zqsatth,lmix,lmix_bis,linter &
+     &    ,lev_out,lunout1,igout)
+
+      endif
+
       if (prt_level.ge.1) print*,'apres thermcell_plume ',lev_out
 
@@ -421 +441 @@
 !
 !
-      write(lunout,*)'THERM NOUVEAU RIO2009 31B'
+!!      write(lunout,*)'THERM NOUVEAU XXXXX'
       CALL thermcell_dry(ngrid,nlay,zlev,pphi,ztv,alim_star,  &
     &                      lalim,lmin,zmax_sec,wmax_sec,lev_out)
@@ -440 +460 @@
 
 
-print*,'THERM 26JJJ'
+!print*,'THERM 26JJJ'
 
 ! Choix de la fonction d'alimentation utilisee pour la fermeture.