Changeset 1080 for LMDZ4/branches/LMDZ4-dev/libf

Timestamp:

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

Author:

jghattas

Message:

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

Location:

LMDZ4/branches/LMDZ4-dev/libf/phylmd

Files:

• surf_land_bucket_mod.F90 (modified) (1 diff)
• surf_land_mod.F90 (modified) (1 diff)
• surf_landice_mod.F90 (modified) (1 diff)
• surf_ocean_mod.F90 (modified) (2 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/surf_land_bucket_mod.F90

@@ -153 +153 @@
 !
     DO i = 1, knon
-       z0_new(i) = SQRT(z0_new(i)**2+rugoro(i)**2)
+       z0_new(i) = MAX(1.5e-05,SQRT(z0_new(i)**2 + rugoro(i)**2))
     END DO
 

LMDZ4/branches/LMDZ4-dev/libf/phylmd/surf_land_mod.F90

@@ -133 +133 @@
 !  
        DO i=1,knon
-          z0_new(i) = SQRT(z0_new(i)**2 + rugoro(i)**2)
+          z0_new(i) = MAX(1.5e-05,SQRT(z0_new(i)**2 + rugoro(i)**2))
        ENDDO
 

LMDZ4/branches/LMDZ4-dev/libf/phylmd/surf_landice_mod.F90

@@ -164 +164 @@
 !
 !****************************************************************************************
-    z0_new(:) = rugoro(:)
+    z0_new(:) = MAX(1.E-3,rugoro(:))
 
 !****************************************************************************************

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   
 !
 !****************************************************************************************