Changeset 2933 for LMDZ5/trunk/libf/phylmd

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)

Location:

LMDZ5/trunk/libf/phylmd/dyn1d

Files:

• 1DUTILS.h (modified) (13 diffs)
• lmdz1d.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/dyn1d/1DUTILS.h

@@ -496 +496 @@
        forc_v =0
        CALL getin('forc_v',forc_v)
-
 !Config  Key  = forc_w
 !Config  Desc = forcage ou non par w
@@ -667 +666 @@
 
       modname = 'dyn1deta0 : '
-      nmq(1)="vap"
-      nmq(2)="cond"
-      do iq=3,nqtot
-        write(nmq(iq),'("tra",i1)') iq-2
-      enddo
+!!      nmq(1)="vap"
+!!      nmq(2)="cond"
+!!      do iq=3,nqtot
+!!        write(nmq(iq),'("tra",i1)') iq-2
+!!      enddo
+      DO iq = 1,nqtot
+        nmq(iq) = trim(tname(iq))
+      ENDDO
       print*,'in dyn1deta0 ',fichnom,klon,klev,nqtot
       CALL open_startphy(fichnom)
@@ -818 +820 @@
       CALL open_restartphy(fichnom)
       print*,'redm1 ',fichnom,klon,klev,nqtot
-      nmq(1)="vap"
-      nmq(2)="cond"
-      nmq(3)="tra1"
-      nmq(4)="tra2"
+!!      nmq(1)="vap"
+!!      nmq(2)="cond"
+!!      nmq(3)="tra1"
+!!      nmq(4)="tra2"
+      DO iq = 1,nqtot
+        nmq(iq) = trim(tname(iq))
+      ENDDO
 
       modname = 'dyn1dredem'
@@ -1404 +1409 @@
       cor(:) = rkappa*temp*(1.+q(:,1)*rv/rd)/(play*(1.+q(:,1)))
 
+
       do k=2,llm-1
 
@@ -2796 +2802 @@
          hq_mod_cas(l)= hq_prof_cas(k2) - frac*(hq_prof_cas(k2)-hq_prof_cas(k1))
          vq_mod_cas(l)= vq_prof_cas(k2) - frac*(vq_prof_cas(k2)-vq_prof_cas(k1))
+         dtrad_mod_cas(l)= dtrad_prof_cas(k2) - frac*(dtrad_prof_cas(k2)-dtrad_prof_cas(k1))
      
          else !play>plev_prof_cas(1)
@@ -2822 +2829 @@
          hq_mod_cas(l)= frac1*hq_prof_cas(k1) - frac2*hq_prof_cas(k2)
          vq_mod_cas(l)= frac1*vq_prof_cas(k1) - frac2*vq_prof_cas(k2)
+         dtrad_mod_cas(l)= frac1*dtrad_prof_cas(k1) - frac2*dtrad_prof_cas(k2)
 
          endif ! play.le.plev_prof_cas(1)
@@ -2850 +2858 @@
          hq_mod_cas(l)= hq_prof_cas(nlev_cas)*fact                            !jyg
          vq_mod_cas(l)= vq_prof_cas(nlev_cas)*fact                            !jyg
+         dtrad_mod_cas(l)= dtrad_prof_cas(nlev_cas)*fact                      !jyg
  
         endif ! play
@@ -5029 +5038 @@
       REAL tau
 !c      DATA tau /3600./
-      DATA tau /5400./
+!      DATA tau /5400./
+       DATA tau /43200./
 !
    INTEGER k,i
@@ -5041 +5051 @@
         DO k = 1,klev
          DO i = 1,klon
-           IF (paprs(i,1)-pplay(i,k) .GT. 10000.) THEN
+!CR: nudging everywhere
+!           IF (paprs(i,1)-pplay(i,k) .GT. 10000.) THEN
 !
             d_u(i,k) = d_u(i,k) + 1./tau*(u_targ(i,k)-u(i,k))
@@ -5048 +5059 @@
             print *,' k,u,d_u,v,d_v ',    &
                       k,u(i,k),d_u(i,k),v(i,k),d_v(i,k)
-           ENDIF
+!           ENDIF
 !
          ENDDO
@@ -5175 +5186 @@
          hq_mod_cas(l)= hq_prof_cas(k2) - frac*(hq_prof_cas(k2)-hq_prof_cas(k1))
          vq_mod_cas(l)= vq_prof_cas(k2) - frac*(vq_prof_cas(k2)-vq_prof_cas(k1))
+         dtrad_mod_cas(l)= dtrad_prof_cas(k2) - frac*(dtrad_prof_cas(k2)-dtrad_prof_cas(k1))
      
          else !play>plev_prof_cas(1)
@@ -5211 +5223 @@
          hq_mod_cas(l)= frac1*hq_prof_cas(k1) - frac2*hq_prof_cas(k2)
          vq_mod_cas(l)= frac1*vq_prof_cas(k1) - frac2*vq_prof_cas(k2)
+         dtrad_mod_cas(l)= frac1*dtrad_prof_cas(k1) - frac2*dtrad_prof_cas(k2)
 
          endif ! play.le.plev_prof_cas(1)
@@ -5247 +5260 @@
          hq_mod_cas(l)= hq_prof_cas(nlev_cas)*fact                     !jyg
          vq_mod_cas(l)= vq_prof_cas(nlev_cas)*fact                     !jyg
+         dtrad_mod_cas(l)= dtrad_prof_cas(nlev_cas)*fact               !jyg
  
         endif ! play

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