Changeset 2979 for LMDZ5/trunk


Ignore:
Timestamp:
Aug 3, 2017, 10:25:17 AM (7 years ago)
Author:
Ehouarn Millour
Message:

Missing initializations that caused problems when running in aquaplanet mode.
EM

Location:
LMDZ5/trunk/libf/phylmd
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/phyaqua_mod.F90

    r2827 r2979  
    349349    wake_deltaq = 0.
    350350    wake_s = 0.
     351    wake_dens = 0.
    351352    wake_cstar = 0.
    352353    wake_pe = 0.
     
    358359    ale_bl_trig =0.
    359360    alp_bl =0.
    360 
     361    treedrg(:,:,:)=0.
    361362
    362363
  • LMDZ5/trunk/libf/phylmd/phyetat0.F90

    r2952 r2979  
    323323  ENDIF
    324324!FC
    325      IF (ifl_pbltree>0) then
     325  IF (ifl_pbltree>0) THEN
    326326!CALL get_field("FTER", pctsrf(:, is_ter), found)
    327      treedrg(:,1:klev,1:nbsrf)= 0.0
    328    CALL get_field("treedrg_ter", drg_ter(:,:), found)
     327    treedrg(:,1:klev,1:nbsrf)= 0.0
     328    CALL get_field("treedrg_ter", drg_ter(:,:), found)
    329329!  found=phyetat0_srf(1,treedrg,"treedrg","drag from vegetation" , 0.)
    330   !lecture du profile de freinage des arbres
    331    IF (.not. found ) then
    332      treedrg(:,1:klev,1:nbsrf)= 0.0
    333    else
    334      treedrg(:,1:klev,is_ter)= drg_ter(:,:)
     330    !lecture du profile de freinage des arbres
     331    IF (.not. found ) THEN
     332      treedrg(:,1:klev,1:nbsrf)= 0.0
     333    ELSE
     334      treedrg(:,1:klev,is_ter)= drg_ter(:,:)
    335335!     found=phyetat0_srf(klev,treedrg,"treedrg","freinage arbres",0.)
    336    ENDIF
    337       endif
    338 
     336    ENDIF
     337  ELSE
     338    ! initialize treedrg(), because it will be written in restartphy.nc
     339    treedrg(:,:,:) = 0.0
     340  ENDIF
    339341
    340342  ! Lecture de l'age de la neige:
  • LMDZ5/trunk/libf/phylmd/physiq_mod.F90

    r2973 r2979  
    20822082       CALL zenang_an(iflag_cycle_diurne.GE.1,jH_cur, &
    20832083            latitude_deg,longitude_deg,rmu0,fract)
    2084        JrNt = 1.0
     2084       swradcorr(:) = 1.0
     2085       JrNt(:) = 1.0
     2086       zrmu0(:) = rmu0(:)
    20852087    ELSE
    20862088       ! recode par Olivier Boucher en sept 2015
Note: See TracChangeset for help on using the changeset viewer.