Ignore:
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:
4 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/surf_land_bucket_mod.F90

    r1072 r1080  
    153153!
    154154    DO i = 1, knon
    155        z0_new(i) = SQRT(z0_new(i)**2+rugoro(i)**2)
     155       z0_new(i) = MAX(1.5e-05,SQRT(z0_new(i)**2 + rugoro(i)**2))
    156156    END DO
    157157
  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/surf_land_mod.F90

    r1067 r1080  
    133133
    134134       DO i=1,knon
    135           z0_new(i) = SQRT(z0_new(i)**2 + rugoro(i)**2)
     135          z0_new(i) = MAX(1.5e-05,SQRT(z0_new(i)**2 + rugoro(i)**2))
    136136       ENDDO
    137137
  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/surf_landice_mod.F90

    r1067 r1080  
    164164!
    165165!****************************************************************************************
    166     z0_new(:) = rugoro(:)
     166    z0_new(:) = MAX(1.E-3,rugoro(:))
    167167
    168168!****************************************************************************************
  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/surf_ocean_mod.F90

    r1067 r1080  
    7979!****************************************************************************************
    8080    INTEGER               :: i
     81    REAL                  :: tmp
     82    REAL, PARAMETER       :: cepdu2=(0.1)**2
    8183    REAL, DIMENSION(klon) :: alb_eau
    8284    REAL, DIMENSION(klon) :: radsol
     
    153155!
    154156!****************************************************************************************
    155     z0_new = SQRT(rugos**2 + rugoro**2)
    156 
    157     ! The rugosity is recalculated with another method
    158     z0_new(:) = 0.0
    159157    DO i = 1, knon
     158       tmp = MAX(cepdu2,u1(i)**2+v1(i)**2)
    160159       z0_new(i) = 0.018*cdragm(i) * (u1(i)**2+v1(i)**2)/RG  &
    161             +  0.11*14e-6 / SQRT(cdragm(i) * (u1(i)**2+v1(i)**2))
     160            +  0.11*14e-6 / SQRT(cdragm(i) * tmp)
    162161       z0_new(i) = MAX(1.5e-05,z0_new(i))
    163     ENDDO
    164    
     162    ENDDO   
    165163!
    166164!****************************************************************************************
Note: See TracChangeset for help on using the changeset viewer.