Changeset 2569

Timestamp:

Jun 14, 2016, 9:16:45 AM (9 years ago)

Author:

Laurent Fairhead

Message:

Pour retrouver 1+1=2 avec iflag_albedo=1
LF

Location:

LMDZ5/trunk/libf/phylmd

Files:

• phyetat0.F90 (modified) (2 diffs)
• phyredem.F90 (modified) (2 diffs)
• physiq_mod.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/phyetat0.F90

@@ -18 +18 @@
        wake_deltat, wake_delta_pbl_TKE, delta_tsurf, wake_fip, wake_pe, &
        wake_s, zgam, zmax0, zmea, zpic, zsig, &
-       zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl
+       zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl, u10m, v10m
   USE geometry_mod, ONLY : longitude_deg, latitude_deg
   USE iostart, ONLY : close_startphy, get_field, get_var, open_startphy
@@ -260 +260 @@
   ENDDO
 
+  found=phyetat0_srf(1,u10m,"U10M","u a 10m",0.)
+  found=phyetat0_srf(1,v10m,"V10M","v a 10m",0.)
+
 !===================================================================
   ! Lecture des temperatures du sol profond:

LMDZ5/trunk/libf/phylmd/phyredem.F90

@@ -21 +21 @@
                                 wake_pe, wake_fip, fm_therm, entr_therm,     &
                                 detr_therm, Ale_bl, Ale_bl_trig, Alp_bl,     &
-                                du_gwd_rando, du_gwd_front
+                                du_gwd_rando, du_gwd_front, u10m, v10m
   USE geometry_mod, ONLY : longitude_deg, latitude_deg
   USE iostart, ONLY: open_restartphy, close_restartphy, put_field, put_var
@@ -151 +151 @@
   CALL put_field_srf2("A_dir_SW","Albedo direct",falb_dir(:,:,:))
   CALL put_field_srf2("A_dif_SW","Albedo diffus",falb_dif(:,:,:))
+
+  CALL put_field_srf1("U10M", "u a 10m", u10m)
+
+  CALL put_field_srf1("V10M", "v a 10m", v10m)
+
 
 ! ================== Tsoil =========================================

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -1077 +1077 @@
        !RC
        ustar(:,:)=0.
-       u10m(:,:)=0.
-       v10m(:,:)=0.
+!       u10m(:,:)=0.
+!       v10m(:,:)=0.
        rain_con(:)=0.
        snow_con(:)=0.