Changeset 2278 for LMDZ5/trunk/libf/phylmd

Timestamp:

May 15, 2015, 6:06:07 PM (9 years ago)

Author:

fhourdin

Message:

Modification sur la prise en compte de la gustinness

Location:

LMDZ5/trunk/libf/phylmd

Files:

• cdrag.F90 (modified) (1 diff)
• clcdrag.F90 (modified) (1 diff)
• physiq.F90 (modified) (1 diff)
• surf_ocean_mod.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/cdrag.F90

@@ -211 +211 @@
            pcfh(i) = zcdn_h(i)* fins(zri(i))
         ENDIF
-! cdrah sur l'ocean cf. Miller et al. (1992)
-        zcr = (0.0016/(zcdn_m(i)*SQRT(zdu2)))*ABS(ztvd-ztsolv)**(1./3.)
-        IF(nsrf.EQ.is_oce) pcfh(i) =f_cdrag_oce* zcdn_h(i)*(1.0+zcr**1.25)**(1./1.25)
+        IF(iflag_gusts==0) THEN
+! cdrah sur l'ocean cf. Miller et al. (1992) - only active when gustiness parameterization is not active
+           zcr = (0.0016/(zcdn_m(i)*SQRT(zdu2)))*ABS(ztvd-ztsolv)**(1./3.)
+           IF(nsrf.EQ.is_oce) pcfh(i) =f_cdrag_oce* zcdn_h(i)*(1.0+zcr**1.25)**(1./1.25)
+        ENDIF
      ENDIF
   END DO

LMDZ5/trunk/libf/phylmd/clcdrag.F90

@@ -129 +129 @@
            pcfh(i) = zcdn(i)* fins(zri(i))
         ENDIF
-        zcr = (0.0016/(zcdn(i)*SQRT(zdu2)))*ABS(ztvd-ztsolv)**(1./3.)
-        IF(nsrf.EQ.is_oce) pcfh(i) =f_cdrag_oce* zcdn(i)*(1.0+zcr**1.25)**(1./1.25)
+        IF(iflag_gusts==0) THEN
+! cdrah sur l'ocean cf. Miller et al. (1992) - only active when gustiness parameterization is not active
+           zcr = (0.0016/(zcdn(i)*SQRT(zdu2)))*ABS(ztvd-ztsolv)**(1./3.)
+           IF(nsrf.EQ.is_oce) pcfh(i) =f_cdrag_oce* zcdn(i)*(1.0+zcr**1.25)**(1./1.25)
+        ENDIF
      ENDIF
   END DO

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -1866 +1866 @@
 !
 !-------gustiness calculation-------!
-     do i = 1, klon
-        gustiness(i)=f_gust_bl*ale_bl(i)+f_gust_wk*ale_wake(i)
-     enddo
+     IF (iflag_gusts==0) THEN
+        gustiness(1:klon)=0
+     ELSE IF (iflag_gusts==1) THEN
+        do i = 1, klon
+           gustiness(i)=f_gust_bl*ale_bl(i)+f_gust_wk*ale_wake(i)
+        enddo
+!     ELSE IF (iflag_gusts==2) THEN
+!        do i = 1, klon
+!           gustiness(i)=f_gust_bl*ale_bl(i)+sigma_wk(i)*f_gust_wk*ale_wake(i) !! need to make sigma_wk accessible here
+!        enddo
+!     ELSE IF (iflag_gusts==3) THEN
+!        do i = 1, klon
+!           gustiness(i)=f_gust_bl*alp_bl(i)+f_gust_wk*alp_wake(i)
+!        enddo
+     ENDIF
+
+
 
      CALL pbl_surface(  &

LMDZ5/trunk/libf/phylmd/surf_ocean_mod.F90

@@ -200 +200 @@
 IF (iflag_z0_oce==0) 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  &
+       tmp = MAX(cepdu2,gustiness(i)+u1(i)**2+v1(i)**2)
+       z0m(i) = 0.018*cdragm(i) * (gustiness(i)+u1(i)**2+v1(i)**2)/RG  &
             +  0.11*14e-6 / SQRT(cdragm(i) * tmp)
        z0m(i) = MAX(1.5e-05,z0m(i))
@@ -209 +209 @@
 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  &
+       tmp = MAX(cepdu2,gustiness(i)+u1(i)**2+v1(i)**2)
+       z0m(i) = 0.018*cdragm(i) * (gustiness(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)
+       z0h(i)=0.4*14e-6 / SQRT(cdragm(i) * tmp)
     ENDDO
 ELSE