Changeset 1624

Timestamp:

May 4, 2012, 12:21:04 PM (13 years ago)

Author:

Laurent Fairhead

Message:

Rajoute la possibilité de ne pas passer par pbl_surface dans la physique FH

---

Added the possibility to not call pbl_surface in the physics FH

File:

• LMDZ5/trunk/libf/phylmd/physiq.F (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/physiq.F

@@ -2044 +2044 @@
 c
 
-      CALL pbl_surface( 
-     e     dtime,     date0,     itap,    days_elapsed+1,
-     e     debut,     lafin,
-     e     rlon,      rlat,      rugoro,  rmu0,     
-     e     rain_fall, snow_fall, solsw,   sollw,    
-     e     t_seri,    q_seri,    u_seri,  v_seri,   
-     e     pplay,     paprs,     pctsrf,            
-     +     ftsol,     falb1,     falb2,   u10m,   v10m,
-     s     sollwdown, cdragh,    cdragm,  u1,    v1,
-     s     albsol1,   albsol2,   sens,    evap,  
-     s     zxtsol,    zxfluxlat, zt2m,    qsat2m, 
-     s     d_t_vdf,   d_q_vdf,   d_u_vdf, d_v_vdf,
-     s     coefh,     coefm,     slab_wfbils,                
-     d     qsol,      zq2m,      s_pblh,  s_lcl,
-     d     s_capCL,   s_oliqCL,  s_cteiCL,s_pblT,
-     d     s_therm,   s_trmb1,   s_trmb2, s_trmb3,
-     d     zxrugs,    zu10m,     zv10m,   fder,
-     d     zxqsurf,   rh2m,      zxfluxu, zxfluxv,
-     d     frugs,     agesno,    fsollw,  fsolsw,
-     d     d_ts,      fevap,     fluxlat, t2m,
-     d     wfbils,    wfbilo,    fluxt,   fluxu,  fluxv,
-     -     dsens,     devap,     zxsnow,
-     -     zxfluxt,   zxfluxq,   q2m,     fluxq, pbl_tke )
+      if (iflag_pbl/=0) then 
+
+        CALL pbl_surface( 
+     e       dtime,     date0,     itap,    days_elapsed+1,
+     e       debut,     lafin,
+     e       rlon,      rlat,      rugoro,  rmu0,     
+     e       rain_fall, snow_fall, solsw,   sollw,    
+     e       t_seri,    q_seri,    u_seri,  v_seri,   
+     e       pplay,     paprs,     pctsrf,            
+     +       ftsol,     falb1,     falb2,   u10m,   v10m,
+     s       sollwdown, cdragh,    cdragm,  u1,    v1,
+     s       albsol1,   albsol2,   sens,    evap,  
+     s       zxtsol,    zxfluxlat, zt2m,    qsat2m, 
+     s       d_t_vdf,   d_q_vdf,   d_u_vdf, d_v_vdf,
+     s       coefh,     coefm,     slab_wfbils,                
+     d       qsol,      zq2m,      s_pblh,  s_lcl,
+     d       s_capCL,   s_oliqCL,  s_cteiCL,s_pblT,
+     d       s_therm,   s_trmb1,   s_trmb2, s_trmb3,
+     d       zxrugs,    zu10m,     zv10m,   fder,
+     d       zxqsurf,   rh2m,      zxfluxu, zxfluxv,
+     d       frugs,     agesno,    fsollw,  fsolsw,
+     d       d_ts,      fevap,     fluxlat, t2m,
+     d       wfbils,    wfbilo,    fluxt,   fluxu,  fluxv,
+     -       dsens,     devap,     zxsnow,
+     -       zxfluxt,   zxfluxq,   q2m,     fluxq, pbl_tke )
 
 
 !-----------------------------------------------------------------------------------------
 ! ajout des tendances de la diffusion turbulente
-      CALL add_phys_tend(d_u_vdf,d_v_vdf,d_t_vdf,d_q_vdf,dql0,'vdf')
+        CALL add_phys_tend(d_u_vdf,d_v_vdf,d_t_vdf,d_q_vdf,dql0,'vdf')
 !-----------------------------------------------------------------------------------------
 
-      if (mydebug) then
-        call writefield_phy('u_seri',u_seri,llm)
-        call writefield_phy('v_seri',v_seri,llm)
-        call writefield_phy('t_seri',t_seri,llm)
-        call writefield_phy('q_seri',q_seri,llm)
-      endif
-
-
-      IF (ip_ebil_phy.ge.2) THEN 
-        ztit='after surface_main'
-        CALL diagetpq(airephy,ztit,ip_ebil_phy,2,2,dtime
+        if (mydebug) then
+          call writefield_phy('u_seri',u_seri,llm)
+          call writefield_phy('v_seri',v_seri,llm)
+          call writefield_phy('t_seri',t_seri,llm)
+          call writefield_phy('q_seri',q_seri,llm)
+        endif
+
+
+        IF (ip_ebil_phy.ge.2) THEN 
+          ztit='after surface_main'
+          CALL diagetpq(airephy,ztit,ip_ebil_phy,2,2,dtime
      e      , t_seri,q_seri,ql_seri,qs_seri,u_seri,v_seri,paprs,pplay
      s      , d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec)
-         call diagphy(airephy,ztit,ip_ebil_phy
+          call diagphy(airephy,ztit,ip_ebil_phy
      e      , zero_v, zero_v, zero_v, zero_v, sens
      e      , evap  , zero_v, zero_v, ztsol
      e      , d_h_vcol, d_qt, d_ec
      s      , fs_bound, fq_bound )
-      END IF 
+        END IF 
+
+      ENDIF
 
 c =================================================================== c