Changeset 2160 for LMDZ5/branches/testing/libf/phylmd/pbl_surface_mod.F90

Timestamp:

Nov 28, 2014, 4:36:29 PM (10 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes -r2070:2158 into testing branch. Compilation problems introduced by revision r2155 have been corrected by hand

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phylmd/pbl_surface_mod.F90 (modified) (4 diffs)

LMDZ5/branches/testing/libf/phylmd/pbl_surface_mod.F90

@@ -319 +319 @@
     REAL, DIMENSION(klon),        INTENT(OUT)       :: alb2_m     ! mean albedo in near IR SW interval
     ! Martin
-        REAL, DIMENSION(klon),        INTENT(OUT)       :: alb3_lic
+    REAL, DIMENSION(klon),        INTENT(OUT)       :: alb3_lic
     ! Martin
     REAL, DIMENSION(klon),        INTENT(OUT)       :: zxsens     ! sensible heat flux at surface with inversed sign 
@@ -475 +475 @@
     LOGICAL, PARAMETER                 :: check=.FALSE.
     REAL, DIMENSION(klon)              :: Kech_h       ! Coefficient d'echange pour l'energie
+    REAL                               :: vent
 
 ! For debugging with IOIPSL
@@ -797 +798 @@
             yts, yqsurf, yrugos, &
             ycdragm, ycdragh )
+! --- special Dice: on force cdragm ( a defaut de forcer ustar) MPL 05082013
+     IF (ok_prescr_ust) then
+      DO i = 1, knon
+       print *,'ycdragm avant=',ycdragm(i)
+       vent= sqrt(yu(i,1)*yu(i,1)+yv(i,1)*yv(i,1))
+!      ycdragm(i) = ust*ust/(1.+(yu(i,1)*yu(i,1)+yv(i,1)*yv(i,1)))
+!      ycdragm(i) = ust*ust/((1.+sqrt(yu(i,1)*yu(i,1)+yv(i,1)*yv(i,1))) &
+!     *sqrt(yu(i,1)*yu(i,1)+yv(i,1)*yv(i,1)))
+       ycdragm(i) = ust*ust/(1.+vent)/vent
+       print *,'ycdragm ust yu yv apres=',ycdragm(i),ust,yu(i,1),yv(i,1)
+      ENDDO
+     ENDIF
+
 
 !****************************************************************************************
@@ -905 +919 @@
                y_flux_u1, y_flux_v1 )
                
+! Special DICE MPL 05082013
+       IF (ok_prescr_ust) THEN
+!         ysnow(:)=0.
+!         yqsol(:)=0.
+!         yagesno(:)=50.
+!         ytsoil(:,:)=300.
+!         yz0_new(:)=0.001
+!         yalb1_new(:)=0.22
+!         yalb2_new(:)=0.22
+!         yevap(:)=flat/RLVTT
+!         yfluxlat(:)=-flat
+!         yfluxsens(:)=-fsens
+!         yqsurf(:)=0.
+!         ytsurf_new(:)=tg
+!         y_dflux_t(:)=0.
+!         y_dflux_q(:)=0.
+          y_flux_u1(:)=ycdragm(:)*(1.+sqrt(yu(:,1)*yu(:,1)+yv(:,1)*yv(:,1)))*yu(:,1)*ypplay(:,1)/RD/yt(:,1)
+          y_flux_v1(:)=ycdragm(:)*(1.+sqrt(yu(:,1)*yu(:,1)+yv(:,1)*yv(:,1)))*yv(:,1)*ypplay(:,1)/RD/yt(:,1)
+      ENDIF
+
      
        CASE(is_lic)