Ignore:
Timestamp:
Feb 5, 2021, 10:00:13 AM (4 years ago)
Author:
musat
Message:

Nouveaux calculs a 2m et 10m

Location:
LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd
Files:
4 added
5 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd/pbl_surface_mod.F90

    r3435 r3823  
    2323  USE climb_wind_mod,      ONLY : climb_wind_down, climb_wind_up
    2424  USE coef_diff_turb_mod,  ONLY : coef_diff_turb
     25  USE ioipsl_getin_p_mod,  ONLY : getin_p
     26  USE cdrag_mod
     27  USE stdlevvar_mod
    2528  USE wx_pbl_mod,          ONLY : wx_pbl_init, wx_pbl_final, &
    2629!!                                  wx_pbl_fuse_no_dts, wx_pbl_split_no_dts, &
    … …  
    4245  INTEGER, SAVE :: iflag_pbl_surface_t2m_bug
    4346  !$OMP THREADPRIVATE(iflag_pbl_surface_t2m_bug)
     47  integer, save :: iflag_new_t2mq2m
     48  !$OMP THREADPRIVATE(iflag_new_t2mq2m)
    4449!FC
    4550!  integer, save :: iflag_frein
    … …  
    177182       alb_dir_m,    alb_dif_m,  zxsens,   zxevap,    &
    178183       alb3_lic,  runoff,    snowhgt,   qsnow,     to_ice,    sissnow,  &
    179        zxtsol,    zxfluxlat, zt2m,     qsat2m,       &
     184       zxtsol,    zxfluxlat, zt2m,     qsat2m, zn2mout, &
    180185       d_t,       d_q,       d_u,      d_v, d_t_diss, &
    181186!!! nrlmd+jyg le 02/05/2011 et le 20/02/2012
    … …  
    394399    REAL, DIMENSION(klon),        INTENT(OUT)       :: zxfluxlat  ! latent flux, mean for each grid point
    395400    REAL, DIMENSION(klon),        INTENT(OUT)       :: zt2m       ! temperature at 2m, mean for each grid point
     401    INTEGER, DIMENSION(klon, 6),  INTENT(OUT)       :: zn2mout ! nb of times t2m is out of the [tsol,temp]
    396402    REAL, DIMENSION(klon),        INTENT(OUT)       :: qsat2m
    397403    REAL, DIMENSION(klon, klev),  INTENT(OUT)       :: d_t        ! change in temperature
    … …  
    537543    REAL, DIMENSION(klon)              :: y_flux_u1, y_flux_v1
    538544    REAL, DIMENSION(klon)              :: yt2m, yq2m, yu10m
     545    INTEGER, DIMENSION(klon, nbsrf, 6) :: yn2mout, yn2mout_x, yn2mout_w
     546    INTEGER, DIMENSION(klon, nbsrf, 6) :: n2mout, n2mout_x, n2mout_w
    539547    REAL, DIMENSION(klon)              :: yustar
    540548    REAL, DIMENSION(klon)              :: ywstar
    … …  
    838846
    839847    IF (first_call) THEN
     848
     849       iflag_new_t2mq2m=0
     850       CALL getin_p('iflag_new_t2mq2m',iflag_new_t2mq2m)
     851       print*,'pbl_iflag_new_t2mq2m=',iflag_new_t2mq2m
     852
    840853       print*,'PBL SURFACE AVEC GUSTINESS'
    841854       first_call=.FALSE.
    … …  
    903916 zxfluxlat(:)=0.
    904917 zt2m(:)=0. ; zq2m(:)=0. ; qsat2m(:)=0. ; rh2m(:)=0.
     918 zn2mout(:,:)=0 ;
    905919 d_t(:,:)=0. ; d_t_diss(:,:)=0. ; d_q(:,:)=0. ; d_u(:,:)=0. ; d_v(:,:)=0.
    906920 zcoefh(:,:,:)=0. ; zcoefm(:,:,:)=0.
    … …  
    18211835
    18221836          ! Calculate the temperature et relative humidity at 2m and the wind at 10m
     1837          IF (iflag_new_t2mq2m==1) THEN
     1838           CALL stdlevvarn(klon, knon, is_ter, zxli, &
     1839               yu(:,1), yv(:,1), yt(:,1), yq(:,1), zgeo1, &
     1840               yts, yqsurf, yz0m, yz0h, ypaprs(:,1), ypplay(:,1), &
     1841               yt2m, yq2m, yt10m, yq10m, yu10m, yn2mout(:, nsrf, :))
     1842          ELSE
    18231843          CALL stdlevvar(klon, knon, is_ter, zxli, &
    18241844               yu(:,1), yv(:,1), yt(:,1), yq(:,1), zgeo1, &
    18251845               yts, yqsurf, yz0m, yz0h, ypaprs(:,1), ypplay(:,1), &
    18261846               yt2m, yq2m, yt10m, yq10m, yu10m, yustar)
    1827          
     1847          ENDIF
    18281848       ENDIF
    18291849
    … …  
    26142634!!! jyg le 07/02/2012
    26152635       IF (iflag_split .eq.0) THEN
     2636        IF (iflag_new_t2mq2m==1) THEN
     2637         CALL stdlevvarn(klon, knon, nsrf, zxli, &
     2638            uzon, vmer, tair1, qair1, zgeo1, &
     2639            tairsol, qairsol, yz0m, yz0h_oupas, psfce, patm, &
     2640            yt2m, yq2m, yt10m, yq10m, yu10m, yn2mout(:, nsrf, :))
     2641        ELSE
    26162642        CALL stdlevvar(klon, knon, nsrf, zxli, &
    26172643            uzon, vmer, tair1, qair1, zgeo1, &
    26182644            tairsol, qairsol, yz0m, yz0h_oupas, psfce, patm, &
    26192645            yt2m, yq2m, yt10m, yq10m, yu10m, yustar)
     2646        ENDIF
    26202647       ELSE  !(iflag_split .eq.0)
     2648        IF (iflag_new_t2mq2m==1) THEN
     2649         CALL stdlevvarn(klon, knon, nsrf, zxli, &
     2650            uzon_x, vmer_x, tair1_x, qair1_x, zgeo1_x, &
     2651            tairsol_x, qairsol, yz0m, yz0h_oupas, psfce, patm, &
     2652            yt2m_x, yq2m_x, yt10m_x, yq10m_x, yu10m_x, yn2mout_x(:, nsrf, :))
     2653         CALL stdlevvarn(klon, knon, nsrf, zxli, &
     2654            uzon_w, vmer_w, tair1_w, qair1_w, zgeo1_w, &
     2655            tairsol_w, qairsol, yz0m, yz0h_oupas, psfce, patm, &
     2656            yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yn2mout_w(:, nsrf, :))
     2657        ELSE
    26212658        CALL stdlevvar(klon, knon, nsrf, zxli, &
    26222659            uzon_x, vmer_x, tair1_x, qair1_x, zgeo1_x, &
    … …  
    26272664            tairsol_w, qairsol, yz0m, yz0h_oupas, psfce, patm, &
    26282665            yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yustar_w)
     2666        ENDIF
    26292667!!!
    26302668       ENDIF  ! (iflag_split .eq.0)
    … …  
    26402678          u10m(i,nsrf)=(yu10m(j) * uzon(j))/SQRT(uzon(j)**2+vmer(j)**2)
    26412679          v10m(i,nsrf)=(yu10m(j) * vmer(j))/SQRT(uzon(j)**2+vmer(j)**2)
     2680!
     2681          DO k = 1, 6
     2682           n2mout(i,nsrf,k) = yn2mout(j,nsrf,k)
     2683          END DO
     2684!
    26422685        ENDDO
    26432686       ELSE  !(iflag_split .eq.0)
    … …  
    26502693          u10m_x(i,nsrf)=(yu10m_x(j) * uzon_x(j))/SQRT(uzon_x(j)**2+vmer_x(j)**2)
    26512694          v10m_x(i,nsrf)=(yu10m_x(j) * vmer_x(j))/SQRT(uzon_x(j)**2+vmer_x(j)**2)
     2695!
     2696          DO k = 1, 6
     2697           n2mout_x(i,nsrf,k) = yn2mout_x(j,nsrf,k)
     2698          END DO
     2699!
    26522700        ENDDO
    26532701        DO j=1, knon
    … …  
    26632711          u10m(i,nsrf) = u10m_x(i,nsrf) + wake_s(i)*(u10m_w(i,nsrf)-u10m_x(i,nsrf))
    26642712          v10m(i,nsrf) = v10m_x(i,nsrf) + wake_s(i)*(v10m_w(i,nsrf)-v10m_x(i,nsrf))
     2713!
     2714          DO k = 1, 6
     2715           n2mout_w(i,nsrf,k) = yn2mout_w(j,nsrf,k)
     2716          END DO
     2717!
    26652718        ENDDO
    26662719!!!
    … …  
    29192972!
    29202973    zxtsol(:) = 0.0  ; zxfluxlat(:) = 0.0
    2921     zt2m(:) = 0.0    ; zq2m(:) = 0.0
     2974    zt2m(:) = 0.0    ; zq2m(:) = 0.0 ; zn2mout(:,:) = 0
    29222975    zustar(:)=0.0 ; zu10m(:) = 0.0   ; zv10m(:) = 0.0
    29232976    s_pblh(:) = 0.0  ; s_plcl(:) = 0.0
    … …  
    29723025          zt2m(i)  = zt2m(i)  + t2m(i,nsrf)  * pctsrf(i,nsrf)
    29733026          zq2m(i)  = zq2m(i)  + q2m(i,nsrf)  * pctsrf(i,nsrf)
     3027!
     3028          DO k = 1, 6
     3029           zn2mout(i,k)  = zn2mout(i,k)  + n2mout(i,nsrf,k)  * pctsrf(i,nsrf)
     3030          ENDDO
     3031!
    29743032          zustar(i) = zustar(i) + ustar(i,nsrf) * pctsrf(i,nsrf)
    29753033          wstar(i,is_ave)=wstar(i,is_ave)+wstar(i,nsrf)*pctsrf(i,nsrf)
    … …  
    30033061          zt2m(i)  = zt2m(i)  + (t2m_x(i,nsrf)+wake_s(i)*(t2m_w(i,nsrf)-t2m_x(i,nsrf))) * pctsrf(i,nsrf)
    30043062          zq2m(i)  = zq2m(i)  + q2m_x(i,nsrf)  * pctsrf(i,nsrf)
     3063!
     3064          DO k = 1, 6
     3065           zn2mout(i,k)  = zn2mout(i,k)  + n2mout_x(i,nsrf,k)  * pctsrf(i,nsrf)
     3066          ENDDO
     3067!
    30053068          zustar(i) = zustar(i) + ustar_x(i,nsrf) * pctsrf(i,nsrf)
    30063069          wstar(i,is_ave)=wstar(i,is_ave)+wstar_x(i,nsrf)*pctsrf(i,nsrf)
    … …  
    30813144    DO nsrf = 1, nbsrf
    30823145       DO i = 1, klon
    3083           zxqsurf(i) = zxqsurf(i) + qsurf(i,nsrf) * pctsrf(i,nsrf)
     3146          zxqsurf(i) = zxqsurf(i) + MAX(qsurf(i,nsrf),0.0) * pctsrf(i,nsrf)
    30843147          zxsnow(i)  = zxsnow(i)  + snow(i,nsrf)  * pctsrf(i,nsrf)
    30853148       ENDDO

  • LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd/phys_local_var_mod.F90

    r3662 r3823  
    380380      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: t2m_min_mon, t2m_max_mon
    381381!$OMP THREADPRIVATE(t2m_min_mon, t2m_max_mon)
    382       REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zq2m_cor, zt2m_cor
    383 !$OMP THREADPRIVATE(zq2m_cor, zt2m_cor)
    384       REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zu10m_cor, zv10m_cor
    385 !$OMP THREADPRIVATE(zu10m_cor, zv10m_cor)
    386       REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zrh2m_cor, zqsat2m_cor
    387 !$OMP THREADPRIVATE(zrh2m_cor, zqsat2m_cor)
    388382      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: weak_inversion
    389383!$OMP THREADPRIVATE(weak_inversion)
    … …  
    475469      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: p_tropopause, z_tropopause, t_tropopause
    476470!$OMP THREADPRIVATE(p_tropopause, z_tropopause, t_tropopause)
     471
     472      INTEGER,ALLOCATABLE,SAVE,DIMENSION(:,:) :: zn2mout
     473!$OMP THREADPRIVATE(zn2mout)
    477474
    478475#ifdef CPP_StratAer
    … …  
    760757      ALLOCATE(zt2m_min_mon(klon), zt2m_max_mon(klon))
    761758      ALLOCATE(t2m_min_mon(klon), t2m_max_mon(klon))
    762       ALLOCATE(zq2m_cor(klon), zt2m_cor(klon), zu10m_cor(klon), zv10m_cor(klon))
    763       ALLOCATE(zrh2m_cor(klon), zqsat2m_cor(klon))
    764759      ALLOCATE(sens(klon), flwp(klon), fiwp(klon))
    765760      ALLOCATE(alp_bl_conv(klon), alp_bl_det(klon))
    … …  
    831826      ALLOCATE (z_tropopause(klon))
    832827      ALLOCATE (t_tropopause(klon))
     828
     829      ALLOCATE(zn2mout(klon,6))
    833830
    834831#ifdef CPP_StratAer
    … …  
    10571054      DEALLOCATE(zt2m_min_mon, zt2m_max_mon)
    10581055      DEALLOCATE(t2m_min_mon, t2m_max_mon)
    1059       DEALLOCATE(zq2m_cor, zt2m_cor, zu10m_cor, zv10m_cor)
    1060       DEALLOCATE(zrh2m_cor, zqsat2m_cor)
    10611056      DEALLOCATE(sens, flwp, fiwp)
    10621057      DEALLOCATE(alp_bl_conv,alp_bl_det)
    … …  
    11191114      DEALLOCATE (z_tropopause)
    11201115      DEALLOCATE (t_tropopause)
     1116
     1117      DEALLOCATE(zn2mout)
    11211118
    11221119#ifdef CPP_StratAer

  • LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd/phys_output_ctrlout_mod.F90

    r3702 r3823  
    271271      ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11, 11/), &
    272272    't2m_sic', "Temp 2m "//clnsurf(4), "K", (/ ('', i=1, 10) /)) /)
     273
     274  TYPE(ctrl_out), SAVE :: o_nt2mout = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), &
     275    'nt2mout', 'Nbt2m out of range complete computation', '-', (/ ('', i=1, 10) /))
     276  TYPE(ctrl_out), SAVE :: o_nq2mout = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), &
     277    'nq2mout', 'Nbq2m out of range complete computation', '-', (/ ('', i=1, 10) /))
     278  TYPE(ctrl_out), SAVE :: o_nu2mout = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), &
     279    'nu2mout', 'Nbu2m out of range complete computation', '-', (/ ('', i=1, 10) /))
     280
     281  TYPE(ctrl_out), SAVE :: o_nt2moutfg = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), &
     282    'nt2moutfg', 'Nbt2m out of range complete/fgRi1 computation', '-', (/ ('', i=1, 10) /))
     283  TYPE(ctrl_out), SAVE :: o_nq2moutfg = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), &
     284    'nq2moutfg', 'Nbq2m out of range complete/fgRi1 computation', '-', (/ ('', i=1, 10) /))
     285  TYPE(ctrl_out), SAVE :: o_nu2moutfg = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), &
     286    'nu2moutfg', 'Nbu2m out of range complete/fgRi1 computation', '-', (/ ('', i=1, 10) /))
    273287
    274288  TYPE(ctrl_out), SAVE :: o_gusts = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11, 11/), &

  • LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd/phys_output_write_mod.F90

    r3702 r3823  
    3838         o_t2m, o_t2m_min, o_t2m_max, &
    3939         o_t2m_min_mon, o_t2m_max_mon, &
     40         o_nt2mout, o_nt2moutfg, &
     41         o_nq2mout, o_nq2moutfg, &
     42         o_nu2mout, o_nu2moutfg, &
    4043         o_q2m, o_ustar, o_u10m, o_v10m, &
    4144         o_wind10m, o_wind10max, o_wind100m, o_gusts, o_sicf, &
    … …  
    249252
    250253    USE phys_local_var_mod, ONLY: zxfluxlat, slp, ptstar, pt0, zxtsol, zt2m, &
    251          zt2m_cor,zq2m_cor,zu10m_cor,zv10m_cor, zrh2m_cor, zqsat2m_cor, &
    252          t2m_min_mon, t2m_max_mon, evap, &
     254         zn2mout, t2m_min_mon, t2m_max_mon, evap, &
    253255         l_mixmin,l_mix, &
    254256         zu10m, zv10m, zq2m, zustar, zxqsurf, &
    … …  
    676678       CALL histwrite_phy(o_slp, slp)
    677679       CALL histwrite_phy(o_tsol, zxtsol)
    678        CALL histwrite_phy(o_t2m, zt2m_cor)
    679        CALL histwrite_phy(o_t2m_min, zt2m_cor)
    680        CALL histwrite_phy(o_t2m_max, zt2m_cor)
     680       CALL histwrite_phy(o_t2m, zt2m)
     681       CALL histwrite_phy(o_t2m_min, zt2m)
     682       CALL histwrite_phy(o_t2m_max, zt2m)
    681683       CALL histwrite_phy(o_t2m_max_mon, t2m_max_mon)
    682684       CALL histwrite_phy(o_t2m_min_mon, t2m_min_mon)
    … …  
    684686       IF (vars_defined) THEN
    685687          DO i=1, klon
    686              zx_tmp_fi2d(i)=SQRT(zu10m_cor(i)*zu10m_cor(i)+zv10m_cor(i)*zv10m_cor(i))
     688             zx_tmp_fi2d(i)=real(zn2mout(i,1))
     689          ENDDO
     690       ENDIF
     691       CALL histwrite_phy(o_nt2mout, zx_tmp_fi2d)
     692
     693       IF (vars_defined) THEN
     694          DO i=1, klon
     695             zx_tmp_fi2d(i)=real(zn2mout(i,2))
     696          ENDDO
     697       ENDIF
     698       CALL histwrite_phy(o_nt2moutfg, zx_tmp_fi2d)
     699
     700       IF (vars_defined) THEN
     701          DO i=1, klon
     702             zx_tmp_fi2d(i)=real(zn2mout(i,3))
     703          ENDDO
     704       ENDIF
     705       CALL histwrite_phy(o_nq2mout, zx_tmp_fi2d)
     706
     707       IF (vars_defined) THEN
     708          DO i=1, klon
     709             zx_tmp_fi2d(i)=real(zn2mout(i,4))
     710          ENDDO
     711       ENDIF
     712       CALL histwrite_phy(o_nq2moutfg, zx_tmp_fi2d)
     713
     714       IF (vars_defined) THEN
     715          DO i=1, klon
     716             zx_tmp_fi2d(i)=real(zn2mout(i,5))
     717          ENDDO
     718       ENDIF
     719       CALL histwrite_phy(o_nu2mout, zx_tmp_fi2d)
     720
     721       IF (vars_defined) THEN
     722          DO i=1, klon
     723             zx_tmp_fi2d(i)=real(zn2mout(i,6))
     724          ENDDO
     725       ENDIF
     726       CALL histwrite_phy(o_nu2moutfg, zx_tmp_fi2d)
     727
     728       IF (vars_defined) THEN
     729          DO i=1, klon
     730             zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
    687731          ENDDO
    688732       ENDIF
    … …  
    691735       IF (vars_defined) THEN
    692736          DO i=1, klon
    693              zx_tmp_fi2d(i)=SQRT(zu10m_cor(i)*zu10m_cor(i)+zv10m_cor(i)*zv10m_cor(i))
     737             zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
    694738          ENDDO
    695739       ENDIF
    … …  
    771815       ENDIF
    772816       CALL histwrite_phy(o_sicf, zx_tmp_fi2d)
    773        CALL histwrite_phy(o_q2m, zq2m_cor)
    774        CALL histwrite_phy(o_ustar, zustar)
    775        CALL histwrite_phy(o_u10m, zu10m_cor)
    776        CALL histwrite_phy(o_v10m, zv10m_cor)
     817       CALL histwrite_phy(o_q2m, zq2m)
     818       IF (vars_defined) zx_tmp_fi2d = zustar
     819       CALL histwrite_phy(o_ustar, zx_tmp_fi2d)
     820       CALL histwrite_phy(o_u10m, zu10m)
     821       CALL histwrite_phy(o_v10m, zv10m)
    777822
    778823       IF (vars_defined) THEN
    … …  
    14121457!       CALL histwrite_phy(o_rh2m_max, zx_tmp_fi2d)
    14131458
    1414        CALL histwrite_phy(o_qsat2m, zqsat2m_cor)
     1459       CALL histwrite_phy(o_qsat2m, qsat2m)
    14151460       CALL histwrite_phy(o_tpot, tpot)
    14161461       CALL histwrite_phy(o_tpote, tpote)

  • LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd/physiq_mod.F90

    r3666 r3823  
    119119       cdragm, cdragh,                   &
    120120       zustar, zu10m, zv10m, rh2m, qsat2m, &
    121        zq2m, zt2m, weak_inversion, &
    122        zq2m_cor,zt2m_cor,zu10m_cor,zv10m_cor, & ! pour corriger d'un bug
    123        zrh2m_cor,zqsat2m_cor, &
     121       zq2m, zt2m, zn2mout, weak_inversion, &
    124122       zt2m_min_mon, zt2m_max_mon,   &         ! pour calcul_divers.h
    125123       t2m_min_mon, t2m_max_mon,  &            ! pour calcul_divers.h
    … …  
    603601    !$OMP THREADPRIVATE(iflag_alp_wk_cond)
    604602
    605     INTEGER,  SAVE               :: iflag_bug_t2m_ipslcm61=1 !
    606     !$OMP THREADPRIVATE(iflag_bug_t2m_ipslcm61)
    607     INTEGER,  SAVE               :: iflag_bug_t2m_stab_ipslcm61=-1 !
    608     !$OMP THREADPRIVATE(iflag_bug_t2m_stab_ipslcm61)
    609 
    610603    REAL t_w(klon,klev),q_w(klon,klev) ! temperature and moisture profiles in the wake region
    611604    REAL t_x(klon,klev),q_x(klon,klev) ! temperature and moisture profiles in the off-wake region
    … …  
    13141307       tau_gl=86400.*tau_gl
    13151308       WRITE(lunout,*) 'debut physiq_mod tau_gl=',tau_gl
    1316 
    1317        iflag_bug_t2m_ipslcm61 = 1
    1318        CALL getin_p('iflag_bug_t2m_ipslcm61', iflag_bug_t2m_ipslcm61)
    1319        iflag_bug_t2m_stab_ipslcm61 = -1
    1320        CALL getin_p('iflag_bug_t2m_stab_ipslcm61', iflag_bug_t2m_stab_ipslcm61)
    13211309
    13221310       CALL getin_p('iflag_alp_wk_cond', iflag_alp_wk_cond)
    … …  
    25292517                                !albedo SB <<<
    25302518            albsol3_lic,runoff,   snowhgt,   qsnow, to_ice, sissnow, &
    2531             zxtsol,    zxfluxlat, zt2m,    qsat2m, &
     2519            zxtsol,    zxfluxlat, zt2m,    qsat2m, zn2mout, &
    25322520            d_t_vdf,   d_q_vdf,   d_u_vdf, d_v_vdf, d_t_diss, &
    25332521                                !nrlmd<
    … …  
    25762564!>jyg
    25772565       ENDIF
    2578 
    2579 !add limitation for t,q at and wind at 10m
    2580         if ( iflag_bug_t2m_ipslcm61 == 0 ) THEN
    2581           CALL borne_var_surf( klon,klev,nbsrf,                 &
    2582             iflag_bug_t2m_stab_ipslcm61,                        &
    2583             t_seri(:,1),q_seri(:,1),u_seri(:,1),v_seri(:,1),    &
    2584             ftsol,zxqsurf,pctsrf,paprs,                         &
    2585             t2m, q2m, u10m, v10m,                               &
    2586             zt2m_cor, zq2m_cor, zu10m_cor, zv10m_cor,           &
    2587             zrh2m_cor, zqsat2m_cor)
    2588         ELSE
    2589           zt2m_cor(:)=zt2m(:)
    2590           zq2m_cor(:)=zq2m(:)
    2591           zu10m_cor(:)=zu10m(:)
    2592           zv10m_cor(:)=zv10m(:)
    2593           zqsat2m_cor=999.999
    2594         ENDIF
    25952566
    25962567       !---------------------------------------------------------------------
Note: See TracChangeset for help on using the changeset viewer.