Changeset 4725 for LMDZ6/trunk/libf/phylmd

Timestamp:

Oct 11, 2023, 6:37:33 PM (14 months ago)

Author:

Laurent Fairhead

Message:

Bug that prevented numeric convergence with old physics

File:

• LMDZ6/trunk/libf/phylmd/cdrag_mod.F90 (modified) (4 diffs, 1 prop)

LMDZ6/trunk/libf/phylmd/cdrag_mod.F90

@@ +1 @@
+!
+!$Id$
+!
 !
 MODULE cdrag_mod
@@ -483 +486 @@
            cdmm(i)=cdmn(i)*FM(i)
            cdhh(i)=f_cdrag_ter*cdhn(i)*FH(i)
+           
            IF (choix_bulk == 0) THEN
              cdm(i)=cdmn(i)*FM(i)
@@ -490 +494 @@
            IF (nsrf.EQ.is_oce) THEN
              cdhh(i)=f_cdrag_oce*cdhn(i)*FH(i)
-             cdmm(i)=MIN(cdm(i),cdmmax)
-             cdhh(i)=MIN(cdh(i),cdhmax)
+             cdmm(i)=MIN(cdmm(i),cdmmax)
+             cdhh(i)=MIN(cdhh(i),cdhmax)
            ENDIF
            IF (ok_cdrag_iter) THEN
@@ -622 +626 @@
      ENDIF ! fin du if (nsrf == is_oce)
   END DO   !  Fin de la boucle sur l'horizontal
-     
+
 END SUBROUTINE cdrag