Ignore:
Timestamp:
Apr 17, 2015, 10:25:14 PM (10 years ago)
Author:
fhourdin
Message:

Bug fixing

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/surf_ocean_mod.F90

    r2261 r2264  
    208208
    209209ELSE IF (iflag_z0_oce==1) THEN
    210     DO i = 1, knon
    211        tmp = MAX(cepdu2,u1(i)**2+v1(i)**2)
    212        z0m(i) = 0.018*cdragm(i) * (u1(i)**2+v1(i)**2)/RG  &
    213             +  0.11*14e-6 / SQRT(cdragm(i) * tmp)
    214        z0m(i) = MAX(1.5e-05,z0m(i))
    215        z0h(i)=0.62*14e-6 / SQRT(cdragm(i) * tmp)
    216     ENDDO  
    217 
     210    DO i = 1, knon
     211       tmp = MAX(cepdu2,u1(i)**2+v1(i)**2)
     212       z0m(i) = 0.018*cdragm(i) * (u1(i)**2+v1(i)**2)/RG  &
     213            + 0.11*14e-6 / SQRT(cdragm(i) * tmp)
     214       z0m(i) = MAX(1.5e-05,z0m(i))
     215       z0h(i)=0.62*14e-6 / SQRT(cdragm(i) * tmp)
     216    ENDDO
    218217ELSE
    219218       STOP'version non prevue'
Note: See TracChangeset for help on using the changeset viewer.