Ignore:
Timestamp:
Feb 5, 2021, 1:23:18 PM (3 years ago)
Author:
musat
Message:

Nouveaux calculs 2m/10m activables par iflag_new_t2mq2m=1
Valeur par defaut iflag_new_t2mq2m=0 (<=> anciens calculs CMIP6)

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/IPSLCM6.0.15/libf/phylmd/pbl_surface_mod.F90

    r3200 r3826  
    2323  USE climb_wind_mod,      ONLY : climb_wind_down, climb_wind_up
    2424  USE coef_diff_turb_mod,  ONLY : coef_diff_turb
    25 
     25  USE ioipsl_getin_p_mod,  ONLY : getin_p
     26  USE cdrag_mod
     27  USE stdlevvar_mod
    2628
    2729  IMPLICIT NONE
    … …  
    3941  integer, save :: iflag_pbl_surface_t2m_bug
    4042  !$OMP THREADPRIVATE(iflag_pbl_surface_t2m_bug)
     43  integer, save :: iflag_new_t2mq2m
     44  !$OMP THREADPRIVATE(iflag_new_t2mq2m)
    4145!FC
    4246!  integer, save :: iflag_frein
    … …  
    168172       alb_dir_m,    alb_dif_m,  zxsens,   zxevap,    &
    169173       alb3_lic,  runoff,    snowhgt,   qsnow,     to_ice,    sissnow,  &
    170        zxtsol,    zxfluxlat, zt2m,     qsat2m,       &
     174       zxtsol,    zxfluxlat, zt2m,     qsat2m, zn2mout, &
    171175       d_t,       d_q,       d_u,      d_v, d_t_diss, &
    172176!!! nrlmd+jyg le 02/05/2011 et le 20/02/2012
    … …  
    382386    REAL, DIMENSION(klon),        INTENT(OUT)       :: zxfluxlat  ! latent flux, mean for each grid point
    383387    REAL, DIMENSION(klon),        INTENT(OUT)       :: zt2m       ! temperature at 2m, mean for each grid point
     388    INTEGER, DIMENSION(klon, 6),  INTENT(OUT)       :: zn2mout ! nb of times t2m is out of the [tsol,temp]
    384389    REAL, DIMENSION(klon),        INTENT(OUT)       :: qsat2m
    385390    REAL, DIMENSION(klon, klev),  INTENT(OUT)       :: d_t        ! change in temperature
    … …  
    519524    REAL, DIMENSION(klon)              :: y_flux_u1, y_flux_v1
    520525    REAL, DIMENSION(klon)              :: yt2m, yq2m, yu10m
     526    INTEGER, DIMENSION(klon, nbsrf, 6) :: yn2mout, yn2mout_x, yn2mout_w
     527    INTEGER, DIMENSION(klon, nbsrf, 6) :: n2mout, n2mout_x, n2mout_w
    521528    REAL, DIMENSION(klon)              :: yustar
    522529    REAL, DIMENSION(klon)              :: ywstar
    … …  
    819826
    820827    IF (first_call) THEN
     828
     829       iflag_new_t2mq2m=0
     830       CALL getin_p('iflag_new_t2mq2m',iflag_new_t2mq2m)
     831       print*,'pbl_iflag_new_t2mq2m=',iflag_new_t2mq2m
     832
    821833       print*,'PBL SURFACE AVEC GUSTINESS'
    822834       first_call=.FALSE.
    … …  
    880892 zxfluxlat(:)=0.
    881893 zt2m(:)=0. ; zq2m(:)=0. ; qsat2m(:)=0. ; rh2m(:)=0.
     894 zn2mout(:,:)=0 ;
    882895 d_t(:,:)=0. ; d_t_diss(:,:)=0. ; d_q(:,:)=0. ; d_u(:,:)=0. ; d_v(:,:)=0.
    883896 zxsens_x(:)=0. ; zxsens_w(:)=0. ; zxfluxlat_x(:)=0. ; zxfluxlat_w(:)=0.
    … …  
    17941807
    17951808          ! Calculate the temperature et relative humidity at 2m and the wind at 10m
     1809          IF (iflag_new_t2mq2m==1) THEN
     1810           CALL stdlevvarn(klon, knon, is_ter, zxli, &
     1811               yu(:,1), yv(:,1), yt(:,1), yq(:,1), zgeo1, &
     1812               yts, yqsurf, yz0m, yz0h, ypaprs(:,1), ypplay(:,1), &
     1813               yt2m, yq2m, yt10m, yq10m, yu10m, yn2mout(:, nsrf, :))
     1814          ELSE
    17961815          CALL stdlevvar(klon, knon, is_ter, zxli, &
    17971816               yu(:,1), yv(:,1), yt(:,1), yq(:,1), zgeo1, &
    17981817               yts, yqsurf, yz0m, yz0h, ypaprs(:,1), ypplay(:,1), &
    17991818               yt2m, yq2m, yt10m, yq10m, yu10m, yustar)
    1800          
     1819          ENDIF
    18011820       END IF
    18021821
    … …  
    25812600!!! jyg le 07/02/2012
    25822601       IF (iflag_split .eq.0) THEN
     2602        IF (iflag_new_t2mq2m==1) THEN
     2603         CALL stdlevvarn(klon, knon, nsrf, zxli, &
     2604            uzon, vmer, tair1, qair1, zgeo1, &
     2605            tairsol, qairsol, yz0m, yz0h_oupas, psfce, patm, &
     2606            yt2m, yq2m, yt10m, yq10m, yu10m, yn2mout(:, nsrf, :))
     2607        ELSE
    25832608        CALL stdlevvar(klon, knon, nsrf, zxli, &
    25842609            uzon, vmer, tair1, qair1, zgeo1, &
    25852610            tairsol, qairsol, yz0m, yz0h_oupas, psfce, patm, &
    25862611            yt2m, yq2m, yt10m, yq10m, yu10m, yustar)
     2612        ENDIF
    25872613       ELSE  !(iflag_split .eq.0)
     2614        IF (iflag_new_t2mq2m==1) THEN
     2615         CALL stdlevvarn(klon, knon, nsrf, zxli, &
     2616            uzon_x, vmer_x, tair1_x, qair1_x, zgeo1_x, &
     2617            tairsol_x, qairsol, yz0m, yz0h_oupas, psfce, patm, &
     2618            yt2m_x, yq2m_x, yt10m_x, yq10m_x, yu10m_x, yn2mout_x(:, nsrf, :))
     2619         CALL stdlevvarn(klon, knon, nsrf, zxli, &
     2620            uzon_w, vmer_w, tair1_w, qair1_w, zgeo1_w, &
     2621            tairsol_w, qairsol, yz0m, yz0h_oupas, psfce, patm, &
     2622            yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yn2mout_w(:, nsrf, :))
     2623        ELSE
    25882624        CALL stdlevvar(klon, knon, nsrf, zxli, &
    25892625            uzon_x, vmer_x, tair1_x, qair1_x, zgeo1_x, &
    … …  
    25942630            tairsol_w, qairsol, yz0m, yz0h_oupas, psfce, patm, &
    25952631            yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yustar_w)
     2632        ENDIF
    25962633!!!
    25972634       ENDIF  ! (iflag_split .eq.0)
    … …  
    26072644          u10m(i,nsrf)=(yu10m(j) * uzon(j))/SQRT(uzon(j)**2+vmer(j)**2)
    26082645          v10m(i,nsrf)=(yu10m(j) * vmer(j))/SQRT(uzon(j)**2+vmer(j)**2)
     2646!
     2647          DO k = 1, 6
     2648           n2mout(i,nsrf,k) = yn2mout(j,nsrf,k)
     2649          END DO
     2650!
    26092651        END DO
    26102652       ELSE  !(iflag_split .eq.0)
    … …  
    26172659          u10m_x(i,nsrf)=(yu10m_x(j) * uzon_x(j))/SQRT(uzon_x(j)**2+vmer_x(j)**2)
    26182660          v10m_x(i,nsrf)=(yu10m_x(j) * vmer_x(j))/SQRT(uzon_x(j)**2+vmer_x(j)**2)
     2661!
     2662          DO k = 1, 6
     2663           n2mout_x(i,nsrf,k) = yn2mout_x(j,nsrf,k)
     2664          END DO
     2665!
    26192666        END DO
    26202667        DO j=1, knon
    … …  
    26302677          u10m(i,nsrf) = u10m_x(i,nsrf) + wake_s(i)*(u10m_w(i,nsrf)-u10m_x(i,nsrf))
    26312678          v10m(i,nsrf) = v10m_x(i,nsrf) + wake_s(i)*(v10m_w(i,nsrf)-v10m_x(i,nsrf))
     2679!
     2680          DO k = 1, 6
     2681           n2mout_w(i,nsrf,k) = yn2mout_w(j,nsrf,k)
     2682          END DO
     2683!
    26322684        END DO
    26332685!!!
    … …  
    28862938!
    28872939    zxtsol(:) = 0.0  ; zxfluxlat(:) = 0.0
    2888     zt2m(:) = 0.0    ; zq2m(:) = 0.0
     2940    zt2m(:) = 0.0    ; zq2m(:) = 0.0 ; zn2mout(:,:) = 0
    28892941    zustar(:)=0.0 ; zu10m(:) = 0.0   ; zv10m(:) = 0.0
    28902942    s_pblh(:) = 0.0  ; s_plcl(:) = 0.0
    … …  
    29392991          zt2m(i)  = zt2m(i)  + t2m(i,nsrf)  * pctsrf(i,nsrf)
    29402992          zq2m(i)  = zq2m(i)  + q2m(i,nsrf)  * pctsrf(i,nsrf)
     2993!
     2994          DO k = 1, 6
     2995           zn2mout(i,k)  = zn2mout(i,k)  + n2mout(i,nsrf,k)  * pctsrf(i,nsrf)
     2996          ENDDO
     2997!
    29412998          zustar(i) = zustar(i) + ustar(i,nsrf) * pctsrf(i,nsrf)
    29422999          wstar(i,is_ave)=wstar(i,is_ave)+wstar(i,nsrf)*pctsrf(i,nsrf)
    … …  
    29703027          zt2m(i)  = zt2m(i)  + (t2m_x(i,nsrf)+wake_s(i)*(t2m_w(i,nsrf)-t2m_x(i,nsrf))) * pctsrf(i,nsrf)
    29713028          zq2m(i)  = zq2m(i)  + q2m_x(i,nsrf)  * pctsrf(i,nsrf)
     3029!
     3030          DO k = 1, 6
     3031           zn2mout(i,k)  = zn2mout(i,k)  + n2mout_x(i,nsrf,k)  * pctsrf(i,nsrf)
     3032          ENDDO
     3033!
    29723034          zustar(i) = zustar(i) + ustar_x(i,nsrf) * pctsrf(i,nsrf)
    29733035          wstar(i,is_ave)=wstar(i,is_ave)+wstar_x(i,nsrf)*pctsrf(i,nsrf)
    … …  
    30483110    DO nsrf = 1, nbsrf
    30493111       DO i = 1, klon
    3050           zxqsurf(i) = zxqsurf(i) + qsurf(i,nsrf) * pctsrf(i,nsrf)
     3112          zxqsurf(i) = zxqsurf(i) + MAX(qsurf(i,nsrf),0.0) * pctsrf(i,nsrf)
    30513113          zxsnow(i)  = zxsnow(i)  + snow(i,nsrf)  * pctsrf(i,nsrf)
    30523114       END DO
Note: See TracChangeset for help on using the changeset viewer.