Changeset 3821 for LMDZ6

Timestamp:

Feb 4, 2021, 5:25:02 PM (3 years ago)

Author:

musat

Message:

Nouveaux calculs a 10m

Location:

LMDZ6/trunk/libf/phylmd

Files:

• screenc_mod.F90 (modified) (1 diff)
• stdlevvar_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/screenc_mod.F90

@@ -183 +183 @@
                     cdram, cdrah, zri1, pref)
       DO i = 1, knon
-!       IF(ok_prescr_ust) THEN
-! La aussi il faut forcer avec ust (FC + MPL 20160210)
-!       ycdragm(i) = ust*ust/(1.+speed(i))/speed(i)
-!       cdram=ycdragm
-!       delu(i) = ust/sqrt(cdram(i))
-!       ELSE
-!       delu(i) = ustar(i)/sqrt(cdram(i))
-!       ENDIF
-        IF(ok_prescr_ust) THEN
-         print*,'chgmnt calculs à 2m; remplacer par ??? si 1D'
-         STOP
-        ENDIF
         delm(i) = sqrt(cdrm(i))/sqrt(cdram(i))
 !verifier que temp est un temperat potentielle.

LMDZ6/trunk/libf/phylmd/stdlevvar_mod.F90

@@ -417 +417 @@
       LOGICAL, dimension(klon) :: ok_q2m_toosmall, ok_q2m_toobig
       LOGICAL, dimension(klon) :: ok_u2m_toobig
+      LOGICAL, dimension(klon) :: ok_t10m_toosmall, ok_t10m_toobig
+      LOGICAL, dimension(klon) :: ok_q10m_toosmall, ok_q10m_toobig
+      LOGICAL, dimension(klon) :: ok_u10m_toobig
+      INTEGER, dimension(klon, 6) :: n10mout
 
 !------------------------------------------------------------------------- 
@@ -437 +441 @@
  &                   cdram, cdrah, zri1, pref)
 
-! --- special Dice: on force cdragm ( a defaut de forcer ustar) MPL 05082013
-     IF (ok_prescr_ust) then
-      DO i = 1, knon
-       print *,'cdram avant=',cdram(i)
-       cdram(i) = ust*ust/speed(i)/speed(i)
-       print *,'cdram ust speed apres=',cdram(i),ust,speed
-      ENDDO
-     ENDIF
-!
-!---------Star variables----------------------------------------------------
 !
       DO i = 1, knon
@@ -606 +600 @@
 !
       zref = 10.0
-      CALL screenpn(klon, knon, nsrf, speed, tpot, q1, &
- &                 ts1, qsurf, z0m, zri1, &
- &                 zref, &
- &                 delu, delte, delq)
-!
-      DO i = 1, knon
-        u_zref(i) = delu(i)
-        q_zref(i) = max(qsurf(i),0.0) + delq(i)
-        te_zref(i) = ts1(i) + delte(i)
-        temp(i) = te_zref(i) * (psol(i)/pat1(i))**(-RKAPPA)
-        u_zref_p(i) = u_zref(i)
-      ENDDO
-!
-! Iteration of the variables at the reference level zref : corrector ; see Hess & McAvaney, 1995 
+!
+       CALL screencn(klon, knon, nsrf, zxli, &
+ &                   speed, tpot, q1, zref, &
+ &                   ts1, qsurf, z0m, z0h, psol, &           
+ &                   cdram, cdrah,  okri, &
+ &                   ri1, 1, &
+ &                   pref_new, delm_new, delh_new, ri10m)
+!
+       DO i = 1, knon
+         u_zref(i) = delm_new(i)*speed(i)
+         q_zref(i) = delh_new(i)*max(q1(i),0.0) + &
+         &           max(qsurf(i),0.0)*(1-delh_new(i))
+         te_zref(i) = delh_new(i)*tpot(i) + ts1(i)*(1-delh_new(i))
+         temp(i) = te_zref(i) * (psol(i)/pref_new(i))**(-RKAPPA)
+         u_zref_p(i) = u_zref(i)
+!
+! compteurs ici
+!
+         ok_t10m_toosmall(i)=te_zref(i).LT.tpot(i).AND. &
+         & te_zref(i).LT.ts1(i)
+         ok_t10m_toobig(i)=te_zref(i).GT.tpot(i).AND. &
+         & te_zref(i).GT.ts1(i)
+         ok_q10m_toosmall(i)=q_zref(i).LT.q1(i).AND. &
+         & q_zref(i).LT.qsurf(i)
+         ok_q10m_toobig(i)=q_zref(i).GT.q1(i).AND. &
+         & q_zref(i).GT.qsurf(i)
+         ok_u10m_toobig(i)=u_zref(i).GT.speed(i)
+!
+         IF(ok_t10m_toosmall(i).OR.ok_t10m_toobig(i).OR. &
+          & ok_q10m_toosmall(i).OR.ok_q10m_toobig(i).OR. &
+          & ok_u10m_toobig(i)) THEN 
+             delm_new(i)=min(max(delm_new(i),0.),1.)
+             delh_new(i)=min(max(delh_new(i),0.),1.)
+             u_zref(i) = delm_new(i)*speed(i)
+             u_zref_p(i) = u_zref(i)
+             q_zref(i) = delh_new(i)*max(q1(i),0.0) + &
+         &               max(qsurf(i),0.0)*(1-delh_new(i))
+             te_zref(i) = delh_new(i)*tpot(i) + ts1(i)*(1-delh_new(i))
+             temp(i) = te_zref(i) * (psol(i)/pref_new(i))**(-RKAPPA)
+         ENDIF 
+!
+       ENDDO
+!
+! Iteration of the variables at the reference level zref : corrector calculation ; see Hess & McAvaney, 1995 
 !
       DO n = 1, niter
@@ -636 +660 @@
           &           max(qsurf(i),0.0)*(1-delh(i))
           te_zref(i) = delh(i)*tpot(i) + ts1(i)*(1-delh(i))
+!
+! return to normal temperature
           temp(i) = te_zref(i) * (psol(i)/pref(i))**(-RKAPPA)
+!
+! compteurs ici
+!
+          ok_t10m_toosmall(i)=te_zref(i).LT.tpot(i).AND. &
+          & te_zref(i).LT.ts1(i)
+          ok_t10m_toobig(i)=te_zref(i).GT.tpot(i).AND. &
+          & te_zref(i).GT.ts1(i)
+          ok_q10m_toosmall(i)=q_zref(i).LT.q1(i).AND. &
+          & q_zref(i).LT.qsurf(i)
+          ok_q10m_toobig(i)=q_zref(i).GT.q1(i).AND. &
+          & q_zref(i).GT.qsurf(i)
+          ok_u10m_toobig(i)=u_zref(i).GT.speed(i)
+!
+          IF(ok_t10m_toosmall(i).OR.ok_t10m_toobig(i)) THEN
+              n10mout(i,2)=n10mout(i,2)+1
+          ENDIF
+          IF(ok_q10m_toosmall(i).OR.ok_q10m_toobig(i)) THEN
+              n10mout(i,4)=n10mout(i,4)+1
+          ENDIF
+          IF(ok_u10m_toobig(i)) THEN
+              n10mout(i,6)=n10mout(i,6)+1
+          ENDIF
+!
+          IF(ok_t10m_toosmall(i).OR.ok_t10m_toobig(i).OR. &
+           & ok_q10m_toosmall(i).OR.ok_q10m_toobig(i).OR. &
+           & ok_u10m_toobig(i)) THEN 
+              delm(i)=min(max(delm(i),0.),1.)
+              delh(i)=min(max(delh(i),0.),1.)
+              u_zref(i) = delm(i)*speed(i)
+              q_zref(i) = delh(i)*max(q1(i),0.0) + &
+          &           max(qsurf(i),0.0)*(1-delh(i))
+              te_zref(i) = delh(i)*tpot(i) + ts1(i)*(1-delh(i))
+              temp(i) = te_zref(i) * (psol(i)/pref(i))**(-RKAPPA)
+          ENDIF
+!
+!
+          IF(n.EQ.ncon) THEN
+            te_zref_con(i) = te_zref(i)
+            q_zref_con(i) = q_zref(i)
+          ENDIF 
+!
         ENDDO 
 !
-      ENDDO
-!
-      DO i = 1, knon
-        u_zref_c(i) = u_zref(i)
-!
-        u_10m(i) = u_zref_p(i) * ok_pred(i) + u_zref_c(i) * ok_corr(i)
-!
+      ENDDO 
+!
+      DO i = 1, knon
         q_zref_c(i) = q_zref(i)
         temp_c(i) = temp(i)
+!
+        ok_pred(i)=0.
+        ok_corr(i)=1.
+!
         t_10m(i) = temp_p(i) * ok_pred(i) + temp_c(i) * ok_corr(i)
         q_10m(i) = q_zref_p(i) * ok_pred(i) + q_zref_c(i) * ok_corr(i)
-      ENDDO
-! 
+!
+        u_zref_c(i) = u_zref(i)
+        u_10m(i) = u_zref_p(i) * ok_pred(i) + u_zref_c(i) * ok_corr(i)
+      ENDDO
+!
       RETURN
       END subroutine stdlevvarn