Ignore:
Timestamp:
Jun 13, 2016, 4:55:29 PM (8 years ago)
Author:
Laurent Fairhead
Message:

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

Location:
LMDZ5/branches/testing/libf/phylmd
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/branches/testing/libf/phylmd/phyetat0.F90

    r2542 r2568  
    1818       wake_deltat, wake_delta_pbl_TKE, delta_tsurf, wake_fip, wake_pe, &
    1919       wake_s, zgam, zmax0, zmea, zpic, zsig, &
    20        zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl
     20       zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl, u10m, v10m
    2121  USE geometry_mod, ONLY : longitude_deg, latitude_deg
    2222  USE iostart, ONLY : close_startphy, get_field, get_var, open_startphy
    … …  
    260260  ENDDO
    261261
     262  found=phyetat0_srf(1,u10m,"U10M","u a 10m",0.)
     263  found=phyetat0_srf(1,v10m,"V10M","v a 10m",0.)
     264
     265 
    262266!===================================================================
    263267  ! Lecture des temperatures du sol profond:
    … …  
    347351  found=phyetat0_get(klev,rnebcon,"RNEBCON","RNEBCON",0.)
    348352  found=phyetat0_get(klev,ratqs,"RATQS","RATQS",0.)
     353
    349354
    350355  found=phyetat0_get(1,run_off_lic_0,"RUNOFFLIC0","RUNOFFLIC0",0.)

  • LMDZ5/branches/testing/libf/phylmd/phyredem.F90

    r2542 r2568  
    2121                                wake_pe, wake_fip, fm_therm, entr_therm,     &
    2222                                detr_therm, Ale_bl, Ale_bl_trig, Alp_bl,     &
    23                                 du_gwd_rando, du_gwd_front
     23                                du_gwd_rando, du_gwd_front, u10m, v10m
    2424  USE geometry_mod, ONLY : longitude_deg, latitude_deg
    2525  USE iostart, ONLY: open_restartphy, close_restartphy, put_field, put_var
    … …  
    152152  CALL put_field_srf2("A_dif_SW","Albedo diffus",falb_dif(:,:,:))
    153153
     154 
     155  CALL put_field_srf1("U10M", "u a 10m", u10m)
     156
     157  CALL put_field_srf1("V10M", "v a 10m", v10m)
     158 
    154159! ================== Tsoil =========================================
    155160  CALL put_field_srf2("Tsoil","Temperature",tsoil(:,:,:))

  • LMDZ5/branches/testing/libf/phylmd/physiq_mod.F90

    r2542 r2568  
    10771077       !RC
    10781078       ustar(:,:)=0.
    1079        u10m(:,:)=0.
    1080        v10m(:,:)=0.
     1079!       u10m(:,:)=0.
     1080!       v10m(:,:)=0.
    10811081       rain_con(:)=0.
    10821082       snow_con(:)=0.
Note: See TracChangeset for help on using the changeset viewer.