Changeset 6042 for LMDZ6/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

Jan 21, 2026, 1:20:47 PM (8 hours ago)

Author:

evignon

Message:

move diagnostics of relative humidity in call_lscp
loss of convergence in prod mode but not in debug

File:

• LMDZ6/trunk/libf/phylmd/physiq_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -3959 +3959 @@
      abortphy, flag_inhib_tend, itap,                             &
      nqo, pdtphys, missing_val, ok_new_lscp,                      &
-     paprs, pplay, omega, t_seri, q_seri, ql_seri, qs_seri,       &
+     paprs, pplay, omega, t_seri, q_seri,                         &
      zmasse, ptconv, ratqsc, ratqs, ratqs_inter_, sigma_qtherm,   &
      qtc_cv, sigt_cv,detrain_cv,fm_cv,fqd,fqcomp,       &
@@ -3989 +3989 @@
      qraindiag, qsnowdiag, dqreva, dqssub, dqrauto,     &
      dqrcol, dqrmelt, dqrfreez, dqsauto, dqsagg, dqsrim,&
-     dqsmelt, dqsfreez)
+     dqsmelt, dqsfreez, zx_rh, zx_rhl, zx_rhi)
 
     !===============================================================================
@@ -4169 +4169 @@
        ENDDO
     ENDIF
-
-    !
-    ! Calculer l'humidite relative pour diagnostique
-    !
-    ! A inclure dans lmdz_call_lscp via un appel à lmdz_lscp_tools
-    DO k = 1, klev
-       DO i = 1, klon
-          zx_t = t_seri(i,k)
-          IF (thermcep) THEN
-             zdelta = MAX(0.,SIGN(1.,rtt-zx_t))
-             zx_qs  = r2es * FOEEW(zx_t,zdelta)/pplay(i,k)
-             zx_qs  = MIN(0.5,zx_qs)
-             zcor   = 1./(1.-retv*zx_qs)
-             zx_qs  = zx_qs*zcor
-          ELSE
-             IF (zx_t.LT.rtt) THEN                  
-                zx_qs = qsats(zx_t)/pplay(i,k)
-             ELSE
-                zx_qs = qsatl(zx_t)/pplay(i,k)
-             ENDIF
-          ENDIF
-          zx_rh(i,k) = q_seri(i,k)/zx_qs
-          IF (iflag_ice_thermo .GT. 0) THEN
-             zx_rhl(i,k) = MIN(q_seri(i,k)/(qsatl(zx_t)/pplay(i,k)),1.)
-             zx_rhi(i,k) = zx_rhl(i,k)*qsatl(zx_t)/qsats(zx_t)
-          ENDIF
-          zqsat(i,k)=zx_qs
-       ENDDO
-    ENDDO
 
     !===============================================================================