Changeset 1966 for LMDZ5/trunk/libf/phy1d/lmdz1d.F

Timestamp:

Feb 11, 2014, 11:52:47 AM (11 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

File:

• LMDZ5/trunk/libf/phy1d/lmdz1d.F (modified) (5 diffs)

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