Ignore:
Timestamp:
Jul 7, 2017, 12:02:31 PM (7 years ago)
Author:
fcheruy
Message:

corrected restart for drgtree
Cd_frein is read from physiq.def

Location:
LMDZ5/branches/LMDZ_tree_FC/libf/phylmd
Files:
1 added
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/branches/LMDZ_tree_FC/libf/phylmd/conf_phys_m.F90

    r2934 r2936  
    25782578!FC
    25792579    write(lunout,*)' ifl_pbltree = ', ifl_pbltree
     2580    write(lunout,*)' Cd_frein = ', Cd_frein
    25802581    write(lunout,*)' iflag_pbl_split = ', iflag_pbl_split
    25812582    write(lunout,*)' iflag_order2_sollw = ', iflag_order2_sollw
  • LMDZ5/branches/LMDZ_tree_FC/libf/phylmd/pbl_surface_mod.F90

    r2924 r2936  
    18221822!FC quid qd yveget ylai yheight ne sont pas definit
    18231823!FC  yveget,ylai,yheight, &
    1824             if (ifl_pbltree .eq. 1) then
     1824            if (ifl_pbltree .ge. 1) then
    18251825            CALL   freinage(knon, yu, yv, yt, &
    18261826!              yveget,ylai, yheight,ypaprs,ypplay,y_d_u_frein,y_d_v_frein)
     
    22352235             y_d_v(j,k)  = y_d_v(j,k) * ypct(j)
    22362236!FC
    2237               if  (nsrf .EQ. is_ter .and. ifl_pbltree .eq. 1  ) then
     2237              if  (nsrf .EQ. is_ter .and. ifl_pbltree .ge. 1  ) then
    22382238!            if (y_d_u_frein(j,k).ne.0. ) then
    22392239!        print*, nsrf,'IS_TER ++', y_d_u_frein(j,k)*ypct(j),y_d_u(j,k),j,k
     
    22422242             y_d_v(j,k) =y_d_v(j,k) + y_d_v_frein(j,k)*ypct(j)
    22432243             treedrg(i,k,nsrf)=y_treedrg(j,k)
     2244             else
     2245             treedrg(i,k,nsrf)=0.
    22442246               endif
    22452247!FC
  • LMDZ5/branches/LMDZ_tree_FC/libf/phylmd/phyetat0.F90

    r2924 r2936  
    321321  ENDIF
    322322!FC
     323     IF (ifl_pbltree>1) then
     324  found=phyetat0_srf(1,treedrg,"treedrg","drag from vegetation" , 0.)
    323325  !lecture du profile de freinage des arbres
    324326   IF (.not. found ) then
    325327     treedrg(:,1:klev,1:nbsrf)= 0.0
    326328   else
    327   found=phyetat0_srf(klev,treedrg,"treedrg","freinage arbres",0.)
     329     found=phyetat0_srf(klev,treedrg,"treedrg","freinage arbres",0.)
    328330   ENDIF
     331      endif
    329332
    330333
Note: See TracChangeset for help on using the changeset viewer.