Changeset 1146 for LMDZ4/trunk/libf/phylmd/surf_ocean_mod.F90

Timestamp:

Apr 9, 2009, 12:11:35 PM (15 years ago)

Author:

Laurent Fairhead

Message:

Réintegration dans le tronc des modifications issues de la branche LMDZ-dev
comprises entre la révision 1074 et 1145
Validation: une simulation de 1 jour en séquentiel sur PC donne les mêmes
résultats entre la trunk et la dev
LF

Location:

LMDZ4/trunk

Files:

• . (modified) (1 prop)
• libf/phylmd/surf_ocean_mod.F90 (modified) (2 diffs)

LMDZ4/trunk/libf/phylmd/surf_ocean_mod.F90

@@ -79 +79 @@
 !****************************************************************************************
     INTEGER               :: i
+    REAL                  :: tmp
+    REAL, PARAMETER       :: cepdu2=(0.1)**2
     REAL, DIMENSION(klon) :: alb_eau
     REAL, DIMENSION(klon) :: radsol
@@ -153 +155 @@
 !
 !****************************************************************************************
-    z0_new = SQRT(rugos**2 + rugoro**2)
-
-    ! The rugosity is recalculated with another method
-    z0_new(:) = 0.0
     DO i = 1, knon
+       tmp = MAX(cepdu2,u1(i)**2+v1(i)**2)
        z0_new(i) = 0.018*cdragm(i) * (u1(i)**2+v1(i)**2)/RG  &
-            +  0.11*14e-6 / SQRT(cdragm(i) * (u1(i)**2+v1(i)**2))
+            +  0.11*14e-6 / SQRT(cdragm(i) * tmp)
        z0_new(i) = MAX(1.5e-05,z0_new(i))
-    ENDDO
-    
+    ENDDO   
 !
 !****************************************************************************************