Changeset 2261

Timestamp:

Apr 17, 2015, 4:23:39 PM (9 years ago)

Author:

fhourdin

Message:

Nouvelle option pour les drag océaniques séparant z0m, z0h, et z0q
pour la quantité de mouvement, la température et l'humidité.
La séparation z0m/z0h est activée par la clée iflag_z0_oce==1.
Par ailleurs, on impose le rapport z0q/z0h = f_z0qh_oce
Typiquement, f_z0qh_oce=1.55

New option to separate z0m, z0h and z0q
z0m and z0h differ if iflag_z0_oce==1.
f_z0qh_oce = z0q/z0h is a constant factor

Alina and Frederic

File:

• LMDZ5/trunk/libf/phylmd/surf_ocean_mod.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/surf_ocean_mod.F90

@@ -112 +112 @@
 
     IF ( f_z0qh_oce .ne. 1.) THEN
+! Si on suit les formulations par exemple de Tessel, on 
+! a z0h=0.4*nu/u*, z0q=0.62*nu/u*, d'ou f_z0qh_oce=0.62/0.4=1.55
        cdragq(:)=cdragh(:)*                                      &
        log(z1lay(:)/z0h(:))/log(z1lay(:)/(f_z0qh_oce*z0h(:)))
@@ -205 +207 @@
     z0h(1:knon)=z0m(1:knon) ! En attendant mieux
 
+ELSE IF (iflag_z0_oce==1) THEN
+    DO i = 1, knon
+       tmp = MAX(cepdu2,u1(i)**2+v1(i)**2)
+       z0m(i) = 0.018*cdragm(i) * (u1(i)**2+v1(i)**2)/RG  &
+            +  0.11*14e-6 / SQRT(cdragm(i) * tmp)
+       z0m(i) = MAX(1.5e-05,z0m(i))
+       z0h(i)=0.62*14e-6 / SQRT(cdragm(i) * tmp)
+    ENDDO  
+
 ELSE
-       STOP'Alina, au boulot :)'
+       STOP'version non prevue'
 ENDIF
 !