Changeset 2952 for LMDZ5/trunk/libf/phylmd/phyetat0.F90

Timestamp:

Jul 14, 2017, 10:23:40 PM (7 years ago)

Author:

Laurent Fairhead

Message:

Parametrization of drag by copses
Need version 4465 of ORCHIDEE at least

6. Cheruy

Location:

LMDZ5/trunk

Files:

• . (modified) (1 prop)
• libf/phylmd/phyetat0.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/phyetat0.F90

@@ -18 +18 @@
        wake_deltat, wake_delta_pbl_TKE, delta_tsurf, wake_fip, wake_pe, &
        wake_s, wake_dens, zgam, zmax0, zmea, zpic, zsig, &
-       zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl, u10m, v10m
+       zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl, u10m, v10m, treedrg
+!FC
   USE geometry_mod, ONLY : longitude_deg, latitude_deg
   USE iostart, ONLY : close_startphy, get_field, get_var, open_startphy
@@ -52 +53 @@
   REAL trs(klon, nbtr)
   REAL zts(klon)
+  ! pour drag arbres FC
+  REAL drg_ter(klon,klev)
 
   CHARACTER*6 ocean_in
@@ -319 +322 @@
      found=phyetat0_srf(1,z0h,"Z0h","Roughness length, enthalpy ",0.001)
   ENDIF
+!FC
+     IF (ifl_pbltree>0) then
+!CALL get_field("FTER", pctsrf(:, is_ter), found)
+     treedrg(:,1:klev,1:nbsrf)= 0.0
+   CALL get_field("treedrg_ter", drg_ter(:,:), found)
+!  found=phyetat0_srf(1,treedrg,"treedrg","drag from vegetation" , 0.)
+  !lecture du profile de freinage des arbres
+   IF (.not. found ) then
+     treedrg(:,1:klev,1:nbsrf)= 0.0
+   else
+     treedrg(:,1:klev,is_ter)= drg_ter(:,:)
+!     found=phyetat0_srf(klev,treedrg,"treedrg","freinage arbres",0.)
+   ENDIF
+      endif
+
 
   ! Lecture de l'age de la neige: