Changeset 2904 for LMDZ5/trunk/libf/phylmd/dyn1d

Timestamp:

Jun 6, 2017, 4:31:18 PM (7 years ago)

Author:

crio

Message:

Modifications pour faire tourner le cas 1D RCE de Caroline dans le format eq_rd_cv (forcing_type=0)

Modifications to run the 1D RCE case-study of Caroline in the eq_rd_cv format (forcing_type=0)

Location:

LMDZ5/trunk/libf/phylmd/dyn1d

Files:

• 1DUTILS.h (modified) (3 diffs)
• lmdz1d.F90 (modified) (1 diff)

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

@@ -5025 +5025 @@
       REAL tau
 !c      DATA tau /3600./
-      DATA tau /5400./
+!      DATA tau /5400./
+       DATA tau /43200./
 !
    INTEGER k,i
@@ -5037 +5038 @@
         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))
@@ -5044 +5046 @@
             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

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

@@ -1077 +1077 @@
        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