Changeset 1966

Timestamp:

Feb 11, 2014, 11:52:47 AM (10 years ago)

Author:

Laurent Fairhead

Message:

Correction de bug pour le cas Sandu.
Nudge au dessus de 700hpa à prendre en compte.
MPL

--This line, and those below, will be ignored--

M phy1d/1D_decl_cases.h
A phy1d/1D_nudge_sandu_astex.h
M phy1d/lmdz1d.F

Location:

LMDZ5/trunk/libf/phy1d

Files:

• 1D_decl_cases.h (modified) (1 diff)
• 1D_nudge_sandu_astex.h (added)
• lmdz1d.F (modified) (5 diffs)

LMDZ5/trunk/libf/phy1d/1D_decl_cases.h

@@ -196 +196 @@
         parameter (mth_ini_sandu=7)
         parameter (day_ini_sandu=15)  ! 196 = 15 juillet 2006
-        real dt_sandu
+        real dt_sandu, tau_sandu
         logical  :: trouve_700=.true.
         parameter (dt_sandu=6.*3600.)   ! forcages donnes ttes les 6 heures par ifa_sandu.txt
-!       parameter (tau_sandu=3600.)  ! temps de relaxation u,v,thetal,qt vers profil init et au dessus 700hPa
+        parameter (tau_sandu=3600.)  ! temps de relaxation u,v,thetal,qt vers profil init et au dessus 700hPa
 !!
         real ts_sandu(nt_sandu)

LMDZ5/trunk/libf/phy1d/lmdz1d.F

@@ -133 +133 @@
       real :: pzero=1.e5
       real :: play (llm),zlay (llm),sig_s(llm),plev(llm+1)
+      real :: playd(llm),zlayd(llm),ap_amma(llm+1),bp_amma(llm+1),poub
 
 !---------------------------------------------------------------------
@@ -154 +155 @@
       REAL, ALLOCATABLE, DIMENSION(:,:):: dq_dyn
       REAL, ALLOCATABLE, DIMENSION(:,:):: d_q_adv
+      REAL, ALLOCATABLE, DIMENSION(:):: d_th_adv
 
 !---------------------------------------------------------------------
@@ -424 +426 @@
       allocate(dq_dyn(llm,nqtot)) 
       allocate(d_q_adv(llm,nqtot)) 
+      allocate(d_th_adv(llm)) 
 
 c
@@ -869 +872 @@
 !! Increment state variables
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+! pour les cas sandu et astex, on reclacule u,v,q,temp et teta dans 1D_nudge_sandu_astex.h
+! au dessus de 700hpa, on relaxe vers les profils initiaux
+      if (forcing_sandu .OR. forcing_astex) then
+#include "1D_nudge_sandu_astex.h"
+      else
         u(1:mxcalc)=u(1:mxcalc) + timestep*(
      s              du_phys(1:mxcalc)
@@ -892 +900 @@
      .             +d_th_adv(1:mxcalc)
      .             +dt_cooling(1:mxcalc))  ! Taux de chauffage ou refroid.
+
+      endif  ! forcing_sandu or forcing_astex
 
         teta=temp*(pzero/play)**rkappa