Changeset 1080 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/surf_ocean_mod.F90

Timestamp:

Jan 26, 2009, 11:51:20 AM (15 years ago)

Author:

jghattas

Message:

Modifications pour que la rugosite ne devient jamais zero (ajout des valeurs minimales).

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/surf_ocean_mod.F90 (modified) (2 diffs)

LMDZ4/branches/LMDZ4-dev/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   
 !
 !****************************************************************************************