Changeset 2933 for LMDZ5/trunk/libf/phylmd/dyn1d/lmdz1d.F90

Timestamp:

Jul 3, 2017, 6:15:15 PM (7 years ago)

Author:

fhourdin

Message:

Les insectes ont parfois une carapace solide.
Correction d'une correction qui en avait oublié un précédente.
(Frédéric pour Marie-Pierre)

File:

• LMDZ5/trunk/libf/phylmd/dyn1d/lmdz1d.F90 (modified) (5 diffs)

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

@@ -20 +20 @@
        wake_deltaq, wake_deltat, wake_s, wake_dens, &
        zgam, zmax0, zmea, zpic, zsig, &
-       zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl
+       zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl, ql_ancien, qs_ancien, &
+       prlw_ancien, prsw_ancien, prw_ancien
   
    USE dimphy
@@ -206 +207 @@
       REAL, ALLOCATABLE, DIMENSION(:,:):: d_q_adv
       REAL, ALLOCATABLE, DIMENSION(:,:):: d_q_nudge
-!      REAL, ALLOCATABLE, DIMENSION(:):: d_t_adv
+!      REAL, ALLOCATABLE, DIMENSION(:):: d_th_adv
 
 !---------------------------------------------------------------------
@@ -580 +581 @@
       allocate(d_q_adv(llm,nqtot)) 
       allocate(d_q_nudge(llm,nqtot)) 
-!      allocate(d_t_adv(llm)) 
+!      allocate(d_th_adv(llm)) 
 
       q(:,:) = 0.
@@ -813 +814 @@
         t_ancien(1,:)=temp(:)
         q_ancien(1,:)=q(:,1)
+        ql_ancien = 0.
+        qs_ancien = 0.
+        prlw_ancien = 0.
+        prsw_ancien = 0.
+        prw_ancien = 0.
 !jyg<
 !!        pbl_tke(:,:,:)=1.e-8
@@ -1078 +1084 @@
        endif
 
+!CRio:Attention modif spécifique cas de Caroline
+      if (forcing_les) then
+         fcoriolis=0.
+!Nudging
+        
+!on calcule dt_cooling
+        do l=1,llm
+        if (play(l).ge.20000.) then
+            dt_cooling(l)=-1.5/86400.
+        elseif ((play(l).ge.10000.).and.((play(l).lt.20000.))) then
+            dt_cooling(l)=-1.5/86400.*(play(l)-10000.)/(10000.)-1./86400.*(20000.-play(l))/10000.*(temp(l)-200.)
+        else
+            dt_cooling(l)=-1.*(temp(l)-200.)/86400.
+        endif
+        enddo
+
+      endif     
+!RC
+
       IF (prt_level >= 5) print*, 'fcoriolis, xlat,mxcalc ', &
                                    fcoriolis, xlat,mxcalc