Ignore:
Timestamp:
Feb 22, 2021, 12:44:07 PM (4 years ago)
Author:
dcugnet
Message:

Update the branch to the current trunk.

Location:
LMDZ6/branches/LMDZ-tracers
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/LMDZ-tracers

  • LMDZ6/branches/LMDZ-tracers/libf/phylmd/pbl_surface_mod.F90

    r3435 r3851  
    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, &
    2730!!                                  wx_pbl_fuse, wx_pbl_split
    2831                                  wx_pbl0_fuse, wx_pbl0_split
     32  use config_ocean_skin_m, only: activate_ocean_skin
    2933
    3034  IMPLICIT NONE
    … …  
    4246  INTEGER, SAVE :: iflag_pbl_surface_t2m_bug
    4347  !$OMP THREADPRIVATE(iflag_pbl_surface_t2m_bug)
     48  integer, save :: iflag_new_t2mq2m
     49  !$OMP THREADPRIVATE(iflag_new_t2mq2m)
     50
    4451!FC
    4552!  integer, save :: iflag_frein
    … …  
    164171       rlon,      rlat,      rugoro,   rmu0,          &
    165172       zsig,      lwdown_m,  pphi,     cldt,          &
    166        rain_f,    snow_f,    solsw_m,  sollw_m,       &
     173       rain_f,    snow_f,    solsw_m,  solswfdiff_m, sollw_m,       &
    167174       gustiness,                                     &
    168175       t,         q,         u,        v,             &
    … …  
    177184       alb_dir_m,    alb_dif_m,  zxsens,   zxevap,    &
    178185       alb3_lic,  runoff,    snowhgt,   qsnow,     to_ice,    sissnow,  &
    179        zxtsol,    zxfluxlat, zt2m,     qsat2m,       &
     186       zxtsol,    zxfluxlat, zt2m,     qsat2m, zn2mout, &
    180187       d_t,       d_q,       d_u,      d_v, d_t_diss, &
    181188!!! nrlmd+jyg le 02/05/2011 et le 20/02/2012
    … …  
    206213!jyg<
    207214!!       zxfluxt,   zxfluxq,   q2m,      flux_q, tke,   &
    208        zxfluxt,   zxfluxq,   q2m,      flux_q, tke_x,   &
     215       zxfluxt,   zxfluxq,   q2m,      flux_q, tke_x,  &
    209216!>jyg
    210217!!! nrlmd+jyg le 02/05/2011 et le 20/02/2012
    … …  
    286293    USE carbon_cycle_mod,   ONLY : carbon_cycle_cpl, carbon_cycle_tr, level_coupling_esm
    287294    USE carbon_cycle_mod,   ONLY : co2_send, nbcf_out, fields_out, yfields_out, cfname_out
     295    use hbtm_mod, only: hbtm
    288296    USE indice_sol_mod
    289297    USE time_phylmdz_mod,   ONLY : day_ini,annee_ref,itau_phy
    … …  
    291299    USE print_control_mod,  ONLY : prt_level,lunout
    292300    USE ioipsl_getin_p_mod, ONLY : getin_p
     301    use phys_state_var_mod, only: ds_ns, dt_ns, delta_sst, delta_sal
     302    use phys_output_var_mod, only: dter, dser, tkt, tks, taur, sss
     303#ifdef CPP_XIOS
     304    USE wxios, ONLY: missing_val
     305#else
     306    use netcdf, only: missing_val => nf90_fill_real
     307#endif
    293308
    294309    IMPLICIT NONE
    … …  
    318333    REAL, DIMENSION(klon),        INTENT(IN)        :: snow_f  ! snow fall
    319334    REAL, DIMENSION(klon),        INTENT(IN)        :: solsw_m ! net shortwave radiation at mean surface
     335    REAL, DIMENSION(klon),        INTENT(IN)        :: solswfdiff_m ! diffuse fraction fordownward shortwave radiation at mean surface
    320336    REAL, DIMENSION(klon),        INTENT(IN)        :: sollw_m ! net longwave radiation at mean surface
    321337    REAL, DIMENSION(klon,klev),   INTENT(IN)        :: t       ! temperature (K)
    … …  
    394410    REAL, DIMENSION(klon),        INTENT(OUT)       :: zxfluxlat  ! latent flux, mean for each grid point
    395411    REAL, DIMENSION(klon),        INTENT(OUT)       :: zt2m       ! temperature at 2m, mean for each grid point
     412    INTEGER, DIMENSION(klon, 6),  INTENT(OUT)       :: zn2mout    ! number of times the 2m temperature is out of the [tsol,temp]
    396413    REAL, DIMENSION(klon),        INTENT(OUT)       :: qsat2m
    397414    REAL, DIMENSION(klon, klev),  INTENT(OUT)       :: d_t        ! change in temperature
    … …  
    453470    REAL, DIMENSION(klon, klev),  INTENT(OUT)       :: zxfluxu    ! u wind tension, mean for each grid point
    454471    REAL, DIMENSION(klon, klev),  INTENT(OUT)       :: zxfluxv    ! v wind tension, mean for each grid point
    455     REAL, DIMENSION(klon, nbsrf+1), INTENT(INOUT)     :: z0m,z0h      ! rugosity length (m)
    456     REAL, DIMENSION(klon, nbsrf), INTENT(INOUT)       :: agesno   ! age of snow at surface
     472    REAL, DIMENSION(klon, nbsrf+1), INTENT(INOUT)   :: z0m,z0h      ! rugosity length (m)
     473    REAL, DIMENSION(klon, nbsrf), INTENT(INOUT)     :: agesno   ! age of snow at surface
    457474    REAL, DIMENSION(klon, nbsrf), INTENT(OUT)       :: solsw      ! net shortwave radiation at surface
    458475    REAL, DIMENSION(klon, nbsrf), INTENT(OUT)       :: sollw      ! net longwave radiation at surface
    459476    REAL, DIMENSION(klon, nbsrf), INTENT(OUT)       :: d_ts       ! change in temperature at surface
    460     REAL, DIMENSION(klon, nbsrf), INTENT(INOUT)       :: evap     ! evaporation at surface
     477    REAL, DIMENSION(klon, nbsrf), INTENT(INOUT)     :: evap     ! evaporation at surface
    461478    REAL, DIMENSION(klon, nbsrf), INTENT(OUT)       :: fluxlat    ! latent flux
    462479    REAL, DIMENSION(klon, nbsrf), INTENT(OUT)       :: t2m        ! temperature at 2 meter height
    … …  
    537554    REAL, DIMENSION(klon)              :: y_flux_u1, y_flux_v1
    538555    REAL, DIMENSION(klon)              :: yt2m, yq2m, yu10m
     556    INTEGER, DIMENSION(klon, nbsrf, 6) :: yn2mout, yn2mout_x, yn2mout_w
     557    INTEGER, DIMENSION(klon, nbsrf, 6) :: n2mout, n2mout_x, n2mout_w
    539558    REAL, DIMENSION(klon)              :: yustar
    540559    REAL, DIMENSION(klon)              :: ywstar
    … …  
    822841    ! Martin
    823842
     843    real, DIMENSION(klon):: ydelta_sst, ydelta_sal, yds_ns, ydt_ns, ydter, ydser, &
     844         ytkt, ytks, ytaur, ysss
     845    ! compression of delta_sst, delta_sal, ds_ns, dt_ns, dter, dser, tkt, tks,
     846    ! taur, sss on ocean points
     847
    824848!****************************************************************************************
    825849! End of declarations
    … …  
    838862
    839863    IF (first_call) THEN
     864
     865       iflag_new_t2mq2m=1
     866       CALL getin_p('iflag_new_t2mq2m',iflag_new_t2mq2m)
     867       print*,'pbl_iflag_new_t2mq2m=',iflag_new_t2mq2m
     868
    840869       print*,'PBL SURFACE AVEC GUSTINESS'
    841870       first_call=.FALSE.
    … …  
    843872       ! Initialize ok_flux_surf (for 1D model)
    844873       if (klon_glo>1) ok_flux_surf=.FALSE.
     874       if (klon_glo>1) ok_forc_tsurf=.FALSE.
    845875
    846876       ! intialize beta_land
    … …  
    903933 zxfluxlat(:)=0.
    904934 zt2m(:)=0. ; zq2m(:)=0. ; qsat2m(:)=0. ; rh2m(:)=0.
     935 zn2mout(:,:)=0 ;
    905936 d_t(:,:)=0. ; d_t_diss(:,:)=0. ; d_q(:,:)=0. ; d_u(:,:)=0. ; d_v(:,:)=0.
    906937 zcoefh(:,:,:)=0. ; zcoefm(:,:,:)=0.
    … …  
    945976!!    tke(:,:,is_ave)=0.
    946977    tke_x(:,:,is_ave)=0.
     978
    947979    wake_dltke(:,:,is_ave)=0.
    948980!>jyg
    … …  
    955987!!! jyg le 10/02/2012
    956988    rh2m_x(:) = 0. ; qsat2m_x(:) = 0. ; rh2m_w(:) = 0. ; qsat2m_w(:) = 0.
    957 !!!
    958989
    959990! 2b) Initialization of all local variables that will be compressed later
    … …  
    965996    yqsurf = 0.0  ; yalb = 0.0 ; yalb_vis = 0.0
    966997!albedo SB <<<
    967     yrain_f = 0.0 ; ysnow_f = 0.0    ; yfder = 0.0     ; ysolsw = 0.0   
     998    yrain_f = 0.0 ; ysnow_f = 0.0    ; yfder = 0.0     ; ysolsw = 0.0
    968999    ysollw = 0.0  ; yz0m = 0.0 ; yz0h = 0.0    ; yu1 = 0.0   
    9691000    yv1 = 0.0     ; ypaprs = 0.0     ; ypplay = 0.0
    … …  
    9761007!!    d_t_diss= 0.0 ;d_u = 0.0     ; d_v = 0.0
    9771008    yqsol = 0.0   
    978     ytherm = 0.0  ; ytke=0.
     1009
     1010    ytke=0.
    9791011!FC
    9801012    y_treedrg=0.
    … …  
    11831215!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    11841216
     1217!--OB this line is not satisfactory because alb is the direct albedo not total albedo
    11851218          solsw(i,nsrf) = solsw_m(i) * (1.-alb(i,nsrf)) / (1.-alb_m(i))
    11861219       ENDDO
    … …  
    12051238!>al1
    12061239
     1240!--OB add diffuse fraction of SW down
     1241   DO n=1,nbcf_out
     1242       IF (cfname_out(n) == "swdownfdiff" ) fields_out(:,n) = solswfdiff_m(:)
     1243   ENDDO
    12071244! >> PC
    12081245   IF (carbon_cycle_cpl .AND. carbon_cycle_tr .AND. nbcf_out.GT.0 ) THEN
    … …  
    13811418             ytke_w(j,k)      = tke_x(i,k,nsrf)+wake_dltke(i,k,nsrf)
    13821419             ywake_dltke(j,k) = wake_dltke(i,k,nsrf)
     1420           
    13831421!>jyg
    13841422          ENDDO
    … …  
    14151453          ENDDO
    14161454       ENDIF
     1455
     1456       if (nsrf == is_oce .and. activate_ocean_skin >= 1) then
     1457          if (activate_ocean_skin == 2 .and. type_ocean == "couple") then
     1458             ydelta_sal(:knon) = delta_sal(ni(:knon))
     1459             ydelta_sst(:knon) = delta_sst(ni(:knon))
     1460          end if
     1461         
     1462          yds_ns(:knon) = ds_ns(ni(:knon))
     1463          ydt_ns(:knon) = dt_ns(ni(:knon))
     1464       end if
    14171465       
    14181466!****************************************************************************************
    … …  
    14551503      ENDDO
    14561504     ENDIF
     1505
    14571506        IF (prt_level >=10) print *,'clcdrag -> ycdragh ', ycdragh
    14581507       ELSE  !(iflag_split .eq.0)
    … …  
    15381587      print *,' args coef_diff_turb: ycdragh ', ycdragh
    15391588      print *,' args coef_diff_turb: ytke ', ytke
     1589
    15401590       ENDIF
    15411591        CALL coef_diff_turb(dtime, nsrf, knon, ni,  &
    … …  
    15671617      print *,' args coef_diff_turb: ycdragh_x ', ycdragh_x
    15681618      print *,' args coef_diff_turb: ytke_x ', ytke_x
     1619
    15691620       ENDIF
    15701621        CALL coef_diff_turb(dtime, nsrf, knon, ni,  &
    … …  
    18211872
    18221873          ! Calculate the temperature et relative humidity at 2m and the wind at 10m
     1874          IF (iflag_new_t2mq2m==1) THEN
     1875           CALL stdlevvarn(klon, knon, is_ter, zxli, &
     1876               yu(:,1), yv(:,1), yt(:,1), yq(:,1), zgeo1, &
     1877               yts, yqsurf, yz0m, yz0h, ypaprs(:,1), ypplay(:,1), &
     1878               yt2m, yq2m, yt10m, yq10m, yu10m, yustar, &
     1879               yn2mout(:, nsrf, :))
     1880          ELSE
    18231881          CALL stdlevvar(klon, knon, is_ter, zxli, &
    18241882               yu(:,1), yv(:,1), yt(:,1), yq(:,1), zgeo1, &
    18251883               yts, yqsurf, yz0m, yz0h, ypaprs(:,1), ypplay(:,1), &
    18261884               yt2m, yq2m, yt10m, yq10m, yu10m, yustar)
     1885          ENDIF
    18271886         
    18281887       ENDIF
    … …  
    19351994               yz0m, yz0h, SFRWL,yalb_dir_new, yalb_dif_new, yevap, yfluxsens,yfluxlat,&
    19361995               ytsurf_new, y_dflux_t, y_dflux_q, slab_wfbils, &
    1937                y_flux_u1, y_flux_v1)
     1996               y_flux_u1, y_flux_v1, ydelta_sst(:knon), ydelta_sal(:knon), &
     1997               yds_ns(:knon), ydt_ns(:knon), ydter(:knon), ydser(:knon), &
     1998               ytkt(:knon), ytks(:knon), ytaur(:knon), ysss)
    19381999      IF (prt_level >=10) THEN
    19392000          print *,'arg de surf_ocean: ycdragh ',ycdragh
    … …  
    20132074!
    20142075!****************************************************************************************
    2015 
    2016 !!!
    2017 !!! jyg le 10/04/2013
     2076!!
     2077!!!
     2078!!! jyg le 10/04/2013 et EV 10/2020
     2079
     2080        IF (ok_forc_tsurf) THEN
     2081            DO j=1,knon
     2082                ytsurf_new(j)=tg
     2083                y_d_ts(j) = ytsurf_new(j) - yts(j)
     2084            ENDDO
     2085        ENDIF ! ok_forc_tsurf
     2086
    20182087!!!
    20192088        IF (ok_flux_surf) THEN
    … …  
    24442513              tke_x(i,k,nsrf)    = ytke(j,k)
    24452514              tke_x(i,k,is_ave) = tke_x(i,k,is_ave) + ytke(j,k)*ypct(j)
     2515
    24462516!>jyg
    24472517           ENDDO
    … …  
    24572527!!            tke(i,k,is_ave) = tke(i,k,is_ave) + tke(i,k,nsrf)*ypct(j)
    24582528            tke_x(i,k,nsrf)   = ytke_x(j,k)
    2459             tke_x(i,k,is_ave)   = tke_x(i,k,is_ave) + tke_x(i,k,nsrf)*ypct(j)
     2529            tke_x(i,k,is_ave)   = tke_x(i,k,is_ave) + tke_x(i,k,nsrf)*ypct(j)       
    24602530            wake_dltke(i,k,is_ave)   = wake_dltke(i,k,is_ave) + wake_dltke(i,k,nsrf)*ypct(j)
     2531           
    24612532
    24622533!>jyg
    … …  
    25382609          d_t_w(:,1), d_t_x(:,1), d_t(:,1)
    25392610       ENDIF
     2611
     2612       if (nsrf == is_oce .and. activate_ocean_skin >= 1) then
     2613          delta_sal = missing_val
     2614          ds_ns = missing_val
     2615          dt_ns = missing_val
     2616          delta_sst = missing_val
     2617          dter = missing_val
     2618          dser = missing_val
     2619          tkt = missing_val
     2620          tks = missing_val
     2621          taur = missing_val
     2622          sss = missing_val
     2623         
     2624          delta_sal(ni(:knon)) = ydelta_sal(:knon)
     2625          ds_ns(ni(:knon)) = yds_ns(:knon)
     2626          dt_ns(ni(:knon)) = ydt_ns(:knon)
     2627          delta_sst(ni(:knon)) = ydelta_sst(:knon)
     2628          dter(ni(:knon)) = ydter(:knon)
     2629          dser(ni(:knon)) = ydser(:knon)
     2630          tkt(ni(:knon)) = ytkt(:knon)
     2631          tks(ni(:knon)) = ytks(:knon)
     2632          taur(ni(:knon)) = ytaur(:knon)
     2633          sss(ni(:knon)) = ysss(:knon)
     2634       end if
    25402635
    25412636!****************************************************************************************
    … …  
    26142709!!! jyg le 07/02/2012
    26152710       IF (iflag_split .eq.0) THEN
     2711        IF (iflag_new_t2mq2m==1) THEN
     2712         CALL stdlevvarn(klon, knon, nsrf, zxli, &
     2713            uzon, vmer, tair1, qair1, zgeo1, &
     2714            tairsol, qairsol, yz0m, yz0h_oupas, psfce, patm, &
     2715            yt2m, yq2m, yt10m, yq10m, yu10m, yustar, &
     2716            yn2mout(:, nsrf, :))
     2717        ELSE
    26162718        CALL stdlevvar(klon, knon, nsrf, zxli, &
    26172719            uzon, vmer, tair1, qair1, zgeo1, &
    26182720            tairsol, qairsol, yz0m, yz0h_oupas, psfce, patm, &
    26192721            yt2m, yq2m, yt10m, yq10m, yu10m, yustar)
     2722        ENDIF
    26202723       ELSE  !(iflag_split .eq.0)
     2724        IF (iflag_new_t2mq2m==1) THEN
     2725         CALL stdlevvarn(klon, knon, nsrf, zxli, &
     2726            uzon_x, vmer_x, tair1_x, qair1_x, zgeo1_x, &
     2727            tairsol_x, qairsol, yz0m, yz0h_oupas, psfce, patm, &
     2728            yt2m_x, yq2m_x, yt10m_x, yq10m_x, yu10m_x, yustar_x, &
     2729            yn2mout_x(:, nsrf, :))
     2730         CALL stdlevvarn(klon, knon, nsrf, zxli, &
     2731            uzon_w, vmer_w, tair1_w, qair1_w, zgeo1_w, &
     2732            tairsol_w, qairsol, yz0m, yz0h_oupas, psfce, patm, &
     2733            yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yustar_w, &
     2734            yn2mout_w(:, nsrf, :))
     2735        ELSE
    26212736        CALL stdlevvar(klon, knon, nsrf, zxli, &
    26222737            uzon_x, vmer_x, tair1_x, qair1_x, zgeo1_x, &
    … …  
    26272742            tairsol_w, qairsol, yz0m, yz0h_oupas, psfce, patm, &
    26282743            yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yustar_w)
     2744        ENDIF
    26292745!!!
    26302746       ENDIF  ! (iflag_split .eq.0)
    … …  
    26402756          u10m(i,nsrf)=(yu10m(j) * uzon(j))/SQRT(uzon(j)**2+vmer(j)**2)
    26412757          v10m(i,nsrf)=(yu10m(j) * vmer(j))/SQRT(uzon(j)**2+vmer(j)**2)
     2758!
     2759          DO k = 1, 6
     2760           n2mout(i,nsrf,k) = yn2mout(j,nsrf,k)
     2761          END DO 
     2762!
    26422763        ENDDO
    26432764       ELSE  !(iflag_split .eq.0)
    … …  
    26502771          u10m_x(i,nsrf)=(yu10m_x(j) * uzon_x(j))/SQRT(uzon_x(j)**2+vmer_x(j)**2)
    26512772          v10m_x(i,nsrf)=(yu10m_x(j) * vmer_x(j))/SQRT(uzon_x(j)**2+vmer_x(j)**2)
     2773!
     2774          DO k = 1, 6
     2775           n2mout_x(i,nsrf,k) = yn2mout_x(j,nsrf,k)
     2776          END DO 
     2777!
    26522778        ENDDO
    26532779        DO j=1, knon
    … …  
    26632789          u10m(i,nsrf) = u10m_x(i,nsrf) + wake_s(i)*(u10m_w(i,nsrf)-u10m_x(i,nsrf))
    26642790          v10m(i,nsrf) = v10m_x(i,nsrf) + wake_s(i)*(v10m_w(i,nsrf)-v10m_x(i,nsrf))
     2791!
     2792          DO k = 1, 6
     2793           n2mout_w(i,nsrf,k) = yn2mout_w(j,nsrf,k)
     2794          END DO 
     2795!
    26652796        ENDDO
    26662797!!!
    … …  
    29193050!
    29203051    zxtsol(:) = 0.0  ; zxfluxlat(:) = 0.0
    2921     zt2m(:) = 0.0    ; zq2m(:) = 0.0
     3052    zt2m(:) = 0.0    ; zq2m(:) = 0.0 ; zn2mout(:,:) = 0
    29223053    zustar(:)=0.0 ; zu10m(:) = 0.0   ; zv10m(:) = 0.0
    29233054    s_pblh(:) = 0.0  ; s_plcl(:) = 0.0
    … …  
    29723103          zt2m(i)  = zt2m(i)  + t2m(i,nsrf)  * pctsrf(i,nsrf)
    29733104          zq2m(i)  = zq2m(i)  + q2m(i,nsrf)  * pctsrf(i,nsrf)
     3105!
     3106          DO k = 1, 6
     3107           zn2mout(i,k)  = zn2mout(i,k)  + n2mout(i,nsrf,k)  * pctsrf(i,nsrf)
     3108          ENDDO 
     3109!
    29743110          zustar(i) = zustar(i) + ustar(i,nsrf) * pctsrf(i,nsrf)
    29753111          wstar(i,is_ave)=wstar(i,is_ave)+wstar(i,nsrf)*pctsrf(i,nsrf)
    … …  
    30033139          zt2m(i)  = zt2m(i)  + (t2m_x(i,nsrf)+wake_s(i)*(t2m_w(i,nsrf)-t2m_x(i,nsrf))) * pctsrf(i,nsrf)
    30043140          zq2m(i)  = zq2m(i)  + q2m_x(i,nsrf)  * pctsrf(i,nsrf)
     3141!
     3142          DO k = 1, 6
     3143           zn2mout(i,k)  = zn2mout(i,k)  + n2mout_x(i,nsrf,k)  * pctsrf(i,nsrf)
     3144          ENDDO
     3145!
    30053146          zustar(i) = zustar(i) + ustar_x(i,nsrf) * pctsrf(i,nsrf)
    30063147          wstar(i,is_ave)=wstar(i,is_ave)+wstar_x(i,nsrf)*pctsrf(i,nsrf)
    … …  
    30813222    DO nsrf = 1, nbsrf
    30823223       DO i = 1, klon
    3083           zxqsurf(i) = zxqsurf(i) + qsurf(i,nsrf) * pctsrf(i,nsrf)
     3224          zxqsurf(i) = zxqsurf(i) + MAX(qsurf(i,nsrf),0.0) * pctsrf(i,nsrf)
    30843225          zxsnow(i)  = zxsnow(i)  + snow(i,nsrf)  * pctsrf(i,nsrf)
    30853226       ENDDO
    … …  
    31423283
    31433284!albedo SB >>>
    3144 SUBROUTINE pbl_surface_newfrac(itime, pctsrf_new, pctsrf_old, &
    3145      evap, z0m, z0h, agesno,                                  &
    3146      tsurf,alb_dir,alb_dif, ustar, u10m, v10m, tke) 
    3147 !albedo SB <<<
     3285  SUBROUTINE pbl_surface_newfrac(itime, pctsrf_new, pctsrf_old, &
     3286       evap, z0m, z0h, agesno,                                  &
     3287       tsurf,alb_dir,alb_dif, ustar, u10m, v10m, tke) 
     3288    !albedo SB <<<
    31483289    ! Give default values where new fraction has appread
    31493290
    31503291    USE indice_sol_mod
     3292    use phys_state_var_mod, only: delta_sal, ds_ns, dt_ns, delta_sst
     3293    use config_ocean_skin_m, only: activate_ocean_skin
    31513294
    31523295    INCLUDE "dimsoil.h"
    … …  
    32353378                      alb_dif(i,k,nsrf) = 0.06
    32363379                   ENDDO
     3380                   if (activate_ocean_skin >= 1) then
     3381                      if (activate_ocean_skin == 2 &
     3382                           .and. type_ocean == "couple") then
     3383                         delta_sal(i) = 0.
     3384                         delta_sst(i) = 0.
     3385                      end if
     3386                     
     3387                      ds_ns(i) = 0.
     3388                      dt_ns(i) = 0.
     3389                   end if
    32373390                ELSE IF (nsrf.EQ.is_sic) THEN
    32383391                   tsurf(i,nsrf) = 271.35
Note: See TracChangeset for help on using the changeset viewer.