Ignore:
Timestamp:
Apr 13, 2015, 10:21:09 AM (9 years ago)
Author:
Laurent Fairhead
Message:

Merged trunk changes 2216:2237 into testing branch

Location:
LMDZ5/branches/testing
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/branches/testing

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

    r2220 r2258  
    181181!!!
    182182       pplay,     paprs,     pctsrf,                  &
    183        ts,        alb1, alb2,ustar, u10m, v10m,wstar, &
     183!albedo SB >>>
     184!       ts,        alb1, alb2,ustar, u10m, v10m,wstar, &
     185       ts,SFRWL,   alb_dir, alb_dif,ustar, u10m, v10m,wstar, &
     186!albedo SB <<<
    184187       cdragh,    cdragm,   zu1,    zv1,              &
    185        alb1_m,    alb2_m,    zxsens,   zxevap,        &
     188!albedo SB >>>
     189!       alb1_m,    alb2_m,    zxsens,   zxevap,        &
     190       alb_dir_m,    alb_dif_m,  zxsens,   zxevap,    &
     191!albedo SB <<<
    186192       alb3_lic,  runoff,    snowhgt,   qsnow,     to_ice,    sissnow,  &
    187193       zxtsol,    zxfluxlat, zt2m,     qsat2m,        &
    … …  
    349355    REAL, DIMENSION(klon, nbsrf), INTENT(INOUT)     :: delta_tsurf !surface temperature difference between
    350356                                                                   !wake and off-wake regions
    351     REAL, DIMENSION(klon, nbsrf), INTENT(INOUT)     :: alb1    ! albedo in visible SW interval
    352     REAL, DIMENSION(klon, nbsrf), INTENT(INOUT)     :: alb2    ! albedo in near infra-red SW interval
     357!albedo SB >>>
     358!    REAL, DIMENSION(klon, nbsrf), INTENT(INOUT)     :: alb1    ! albedo in visible SW interval
     359!    REAL, DIMENSION(klon, nbsrf), INTENT(INOUT)     :: alb2    ! albedo in near infra-red SW interval
     360    REAL, DIMENSIOn(6),intent(in) :: SFRWL
     361    REAL, DIMENSION(klon, nsw, nbsrf), INTENT(INOUT)     :: alb_dir,alb_dif
     362!albedo SB <<<
    353363!jyg Pourquoi ustar et wstar sont-elles INOUT ?
    354364    REAL, DIMENSION(klon, nbsrf), INTENT(INOUT)     :: ustar   ! u* (m/s)
    … …  
    371381    REAL, DIMENSION(klon),        INTENT(OUT)       :: zu1        ! u wind speed in first layer
    372382    REAL, DIMENSION(klon),        INTENT(OUT)       :: zv1        ! v wind speed in first layer
    373     REAL, DIMENSION(klon),        INTENT(OUT)       :: alb1_m     ! mean albedo in visible SW interval
    374     REAL, DIMENSION(klon),        INTENT(OUT)       :: alb2_m     ! mean albedo in near IR SW interval
     383!albedo SB >>>
     384!    REAL, DIMENSION(klon),        INTENT(OUT)       :: alb1_m     ! mean albedo
     385!    in visible SW interval
     386!    REAL, DIMENSION(klon),        INTENT(OUT)       :: alb2_m     ! mean albedo
     387!    in near IR SW interval
     388    REAL, DIMENSION(klon, nsw),        INTENT(OUT)       :: alb_dir_m,alb_dif_m
     389!albedo SB <<<
    375390    ! Martin
    376391    REAL, DIMENSION(klon),        INTENT(OUT)       :: alb3_lic
    … …  
    505520    REAL, DIMENSION(klon)              :: r_co2_ppm     ! taux CO2 atmosphere
    506521    REAL, DIMENSION(klon)              :: yts, yrugos, ypct, yz0_new
    507     REAL, DIMENSION(klon)              :: yalb, yalb1, yalb2
     522!albedo SB >>>
     523!   REAL, DIMENSION(klon)              :: yalb, yalb1, yalb2
     524    REAL, DIMENSION(klon)              :: yalb,yalb_vis
     525!albedo SB <<<
    508526    REAL, DIMENSION(klon)              :: yu1, yv1
    509527    REAL, DIMENSION(klon)              :: ysnow, yqsurf, yagesno, yqsol
    … …  
    535553    REAL, DIMENSION(klon)              :: tair1, qair1, tairsol
    536554    REAL, DIMENSION(klon)              :: psfce, patm
    537     REAL, DIMENSION(klon)              :: qairsol, zgeo1
     555    REAL, DIMENSION(klon)              :: qairsol, zgeo1, speed, zri1, pref !speed, zri1, pref, added by Fuxing WANG, 04/03/2015
    538556    REAL, DIMENSION(klon)              :: rugo1
    539557    REAL, DIMENSION(klon)              :: yfluxsens
    … …  
    542560    REAL, DIMENSION(klon)              :: ypsref
    543561    REAL, DIMENSION(klon)              :: yevap, ytsurf_new, yalb1_new, yalb2_new, yalb3_new
     562!albedo SB >>>
     563    REAL, DIMENSION(klon,nsw)          :: yalb_dir_new, yalb_dif_new
     564!albedo SB <<<
    544565    REAL, DIMENSION(klon)              :: ztsol
    545566    REAL, DIMENSION(klon)              :: alb_m  ! mean albedo for whole SW interval
    … …  
    693714    REAL, DIMENSION(klon)       :: ytrmb3_w
    694715!
    695     REAL, DIMENSION(klon)              :: uzon_x, vmer_x
     716    REAL, DIMENSION(klon)              :: uzon_x, vmer_x, speed_x, zri1_x, pref_x !speed_x, zri1_x, pref_x, added by Fuxing WANG, 04/03/2015
    696717    REAL, DIMENSION(klon)              :: zgeo1_x, tair1_x, qair1_x, tairsol_x
    697718!
    698     REAL, DIMENSION(klon)              :: uzon_w, vmer_w
     719    REAL, DIMENSION(klon)              :: uzon_w, vmer_w, speed_w, zri1_w, pref_w !speed_w, zri1_w, pref_w, added by Fuxing WANG, 04/03/2015
    699720    REAL, DIMENSION(klon)              :: zgeo1_w, tair1_w, qair1_w, tairsol_w
    700721
    … …  
    855876 cdragh(:)=0. ; cdragm(:)=0.
    856877 zu1(:)=0. ; zv1(:)=0.
    857  alb1_m(:)=0. ; alb2_m(:)=0. ; alb3_lic(:)=0.
     878!albedo SB >>>
     879! alb1_m(:)=0. ; alb2_m(:)=0. ; alb3_lic(:)=0.
     880  alb_dir_m=0. ; alb_dif_m=0. ; alb3_lic(:)=0.
     881!albedo SB <<<
    858882 zxsens(:)=0. ; zxevap(:)=0. ; zxtsol(:)=0.
    859883 d_t_w(:,:)=0. ; d_q_w(:,:)=0. ; d_t_x(:,:)=0. ; d_q_x(:,:)=0.
    … …  
    920944    ypct = 0.0    ; yts = 0.0        ; ysnow = 0.0
    921945!!    zv1 = 0.0     ; yqsurf = 0.0     ; yalb1 = 0.0     ; yalb2 = 0.0   
    922     yqsurf = 0.0  ; yalb1 = 0.0      ; yalb2 = 0.0   
     946!albedo SB >>>
     947!    yqsurf = 0.0  ; yalb1 = 0.0      ; yalb2 = 0.0   
     948    yqsurf = 0.0  ; yalb = 0.0 ; yalb_vis = 0.0
     949!albedo SB <<<
    923950    yrain_f = 0.0 ; ysnow_f = 0.0    ; yfder = 0.0     ; ysolsw = 0.0   
    924951    ysollw = 0.0  ; yrugos = 0.0     ; yu1 = 0.0   
    … …  
    10701097! * alb_m  : mean albedo at whole SW interval
    10711098
    1072     alb1_m(:) = 0.0
    1073     alb2_m(:) = 0.0
    1074     DO nsrf = 1, nbsrf
     1099!albedo SB >>>
     1100!    alb1_m(:) = 0.0
     1101!    alb2_m(:) = 0.0
     1102!    DO nsrf = 1, nbsrf
     1103!       DO i = 1, klon
     1104!          alb1_m(i) = alb1_m(i) + alb1(i,nsrf) * pctsrf(i,nsrf)
     1105!          alb2_m(i) = alb2_m(i) + alb2(i,nsrf) * pctsrf(i,nsrf)
     1106!       ENDDO
     1107!    ENDDO
     1108
     1109    alb_dir_m(:,:) = 0.0
     1110    alb_dif_m(:,:) = 0.0
     1111    DO k = 1, nsw
     1112     DO nsrf = 1, nbsrf
    10751113       DO i = 1, klon
    1076           alb1_m(i) = alb1_m(i) + alb1(i,nsrf) * pctsrf(i,nsrf)
    1077           alb2_m(i) = alb2_m(i) + alb2(i,nsrf) * pctsrf(i,nsrf)
     1114          alb_dir_m(i,k) = alb_dir_m(i,k) + alb_dir(i,k,nsrf) * pctsrf(i,nsrf)
     1115          alb_dif_m(i,k) = alb_dif_m(i,k) + alb_dif(i,k,nsrf) * pctsrf(i,nsrf)
    10781116       ENDDO
     1117     ENDDO
    10791118    ENDDO
    10801119
    10811120! We here suppose the fraction f1 of incoming radiance of visible radiance
    10821121! as a fraction of all shortwave radiance
    1083     f1 = 0.5 
     1122    f1 = 0.5
    10841123!    f1 = 1    ! put f1=1 to recreate old calculations
    10851124
    1086     DO nsrf = 1, nbsrf
    1087        DO i = 1, klon
    1088           alb(i,nsrf) = f1*alb1(i,nsrf) + (1-f1)*alb2(i,nsrf)
    1089        ENDDO
     1125!    DO nsrf = 1, nbsrf
     1126!       DO i = 1, klon
     1127!          alb(i,nsrf) = f1*alb1(i,nsrf) + (1-f1)*alb2(i,nsrf)
     1128!       ENDDO
     1129!    ENDDO
     1130!
     1131!    DO i = 1, klon
     1132!       alb_m(i) = f1*alb1_m(i) + (1-f1)*alb2_m(i)
     1133!    END DO
     1134
     1135
     1136!f1 is already included with SFRWL values in each surf files
     1137    alb=0.0
     1138    DO k=1,nsw
     1139      DO nsrf = 1, nbsrf
     1140        DO i = 1, klon
     1141            alb(i,nsrf) = alb(i,nsrf) + alb_dir(i,k,nsrf)*SFRWL(k)
     1142        ENDDO
     1143      ENDDO
    10901144    ENDDO
    10911145
    1092     DO i = 1, klon
    1093        alb_m(i) = f1*alb1_m(i) + (1-f1)*alb2_m(i)
    1094     END DO
     1146    alb_m=0.0
     1147    DO k = 1,nsw
     1148      DO i = 1, klon
     1149        alb_m(i) = alb_m(i) + alb_dir_m(i,k)*SFRWL(k)
     1150      END DO
     1151    ENDDO
     1152!albedo SB <<<
     1153
     1154
    10951155
    10961156! Calculation of mean temperature at surface grid points
    … …  
    11701230          yqsurf(j)  = qsurf(i,nsrf)
    11711231          yalb(j)    = alb(i,nsrf)
    1172           yalb1(j)   = alb1(i,nsrf)
    1173           yalb2(j)   = alb2(i,nsrf)
     1232!albedo SB >>>
     1233!         yalb1(j)   = alb1(i,nsrf)
     1234!         yalb2(j)   = alb2(i,nsrf)
     1235          yalb_vis(j) = alb_dir(i,1,nsrf)
     1236          if(nsw==6)then
     1237            yalb_vis(j)=(alb_dir(i,1,nsrf)*SFRWL(1)+alb_dir(i,2,nsrf)*SFRWL(2) &
     1238              +alb_dir(i,3,nsrf)*SFRWL(3))/(SFRWL(1)+SFRWL(2)+SFRWL(3))
     1239          endif
     1240!albedo SB <<<
    11741241          yrain_f(j) = rain_f(i)
    11751242          ysnow_f(j) = snow_f(i)
    … …  
    12951362!!!
    12961363!!! nrlmd & jyg les 02/05/2011, 13/06/2011, 05/02/2012
    1297         CALL clcdrag( knon, nsrf, ypaprs, ypplay, &
    1298             yu(:,1), yv(:,1), yt(:,1), yq(:,1), &
     1364! Faire disparaitre les lignes commentees fin 2015 (le temps des tests)
     1365!       CALL clcdrag( knon, nsrf, ypaprs, ypplay, &
     1366!           yu(:,1), yv(:,1), yt(:,1), yq(:,1), &
     1367!           yts, yqsurf, yrugos, &
     1368!           ycdragm, ycdragh )
     1369! Fuxing WANG, 04/03/2015, replace the clcdrag by the merged version: cdrag
     1370        DO i = 1, knon
     1371!          print*,'PBL ',i,RD
     1372!          print*,'PBL ',yt(i,1),ypaprs(i,1),ypplay(i,1)
     1373           zgeo1(i) = RD * yt(i,1) / (0.5*(ypaprs(i,1)+ypplay(i,1))) &
     1374                * (ypaprs(i,1)-ypplay(i,1))
     1375           speed(i) = SQRT(yu(i,1)**2+yv(i,1)**2)
     1376        END DO
     1377        CALL cdrag(knon, nsrf, &
     1378            speed, yt(:,1), yq(:,1), zgeo1, ypaprs(:,1),&
    12991379            yts, yqsurf, yrugos, &
    1300             ycdragm, ycdragh )
     1380            ycdragm, ycdragh, zri1, pref )
     1381
    13011382! --- special Dice: on force cdragm ( a defaut de forcer ustar) MPL 05082013
    13021383     IF (ok_prescr_ust) then
    … …  
    13131394        IF (prt_level >=10) print *,'clcdrag -> ycdragh ', ycdragh
    13141395       ELSE  !(iflag_split .eq.0)
    1315         CALL clcdrag( knon, nsrf, ypaprs, ypplay, &
    1316             yu_x(:,1), yv_x(:,1), yt_x(:,1), yq_x(:,1), &
     1396
     1397! Faire disparaitre les lignes commentees fin 2015 (le temps des tests)
     1398!       CALL clcdrag( knon, nsrf, ypaprs, ypplay, &
     1399!           yu_x(:,1), yv_x(:,1), yt_x(:,1), yq_x(:,1), &
     1400!           yts_x, yqsurf, yrugos, &
     1401!           ycdragm_x, ycdragh_x )
     1402! Fuxing WANG, 04/03/2015, replace the clcdrag by the merged version: cdrag
     1403        DO i = 1, knon
     1404           zgeo1_x(i) = RD * yt_x(i,1) / (0.5*(ypaprs(i,1)+ypplay(i,1))) &
     1405                * (ypaprs(i,1)-ypplay(i,1))
     1406           speed_x(i) = SQRT(yu_x(i,1)**2+yv_x(i,1)**2)
     1407        END DO
     1408        CALL cdrag(knon, nsrf, &
     1409            speed_x, yt_x(:,1), yq_x(:,1), zgeo1_x, ypaprs(:,1),&
    13171410            yts_x, yqsurf, yrugos, &
    1318             ycdragm_x, ycdragh_x )
     1411            ycdragm_x, ycdragh_x, zri1_x, pref_x )
     1412
    13191413! --- special Dice. JYG+MPL 25112013
    13201414        IF (ok_prescr_ust) then
    … …  
    17101804               ylwdown, yq2m, yt2m, &
    17111805               ysnow, yqsol, yagesno, ytsoil, &
    1712                yz0_new, yalb1_new, yalb2_new, yevap, yfluxsens, yfluxlat, &
     1806!albedo SB >>>
     1807!              yz0_new, yalb1_new, yalb2_new, yevap, yfluxsens, yfluxlat, &
     1808               yz0_new, SFRWL, yalb_dir_new, yalb_dif_new, yevap, yfluxsens,yfluxlat,&
     1809!albedo SB <<<
    17131810               yqsurf, ytsurf_new, y_dflux_t, y_dflux_q, &
    17141811               y_flux_u1, y_flux_v1 )
    … …  
    17461843               ypsref, yu1, yv1, yrugoro, pctsrf, &
    17471844               ysnow, yqsurf, yqsol, yagesno, &
    1748                ytsoil, yz0_new, yalb1_new, yalb2_new, yevap, yfluxsens, yfluxlat, &
     1845!albedo SB >>>
     1846!              ytsoil, yz0_new, yalb1_new, yalb2_new, yevap, yfluxsens, yfluxlat, &
     1847               ytsoil, yz0_new, SFRWL, yalb_dir_new, yalb_dif_new, yevap,yfluxsens,yfluxlat, &
     1848!albedo SB <<<
    17491849               ytsurf_new, y_dflux_t, y_dflux_q, &
    17501850               yzsig, ycldt, &
    … …  
    17781878         
    17791879       CASE(is_oce)
    1780           CALL surf_ocean(rlon, rlat, ysolsw, ysollw, yalb1, &
     1880!albedo SB >>>
     1881!          CALL surf_ocean(rlon, rlat, ysolsw, ysollw, yalb1, &
     1882           CALL surf_ocean(rlon, rlat, ysolsw, ysollw, yalb_vis, &
     1883!albedo SB <<<
    17811884               yrugos, ywindsp, rmu0, yfder, yts, &
    17821885               itap, dtime, jour, knon, ni, &
    … …  
    17861889               ypsref, yu1, yv1, yrugoro, pctsrf, &
    17871890               ysnow, yqsurf, yagesno, &
    1788                yz0_new, yalb1_new, yalb2_new, yevap, yfluxsens, yfluxlat, &
     1891!albedo SB >>>
     1892!              yz0_new, yalb1_new, yalb2_new, yevap, yfluxsens, yfluxlat, &
     1893               yz0_new, SFRWL,yalb_dir_new, yalb_dif_new, yevap, yfluxsens,yfluxlat,&
     1894!albedo SB <<<
    17891895               ytsurf_new, y_dflux_t, y_dflux_q, slab_wfbils, &
    17901896               y_flux_u1, y_flux_v1)
    … …  
    18071913       CASE(is_sic)
    18081914          CALL surf_seaice( &
    1809                rlon, rlat, ysolsw, ysollw, yalb1, yfder, &
     1915!albedo SB >>>
     1916!               rlon, rlat, ysolsw, ysollw, yalb1, yfder, &
     1917               rlon, rlat, ysolsw, ysollw, yalb_vis, yfder, &
     1918!albedo SB <<<
    18101919               itap, dtime, jour, knon, ni, &
    18111920               lafin, &
    … …  
    18151924               ypsref, yu1, yv1, yrugoro, pctsrf, &
    18161925               ysnow, yqsurf, yqsol, yagesno, ytsoil, &
    1817                yz0_new, yalb1_new, yalb2_new, yevap, yfluxsens, yfluxlat, &
     1926!albedo SB >>>
     1927!               yz0_new, yalb1_new, yalb2_new, yevap, yfluxsens, yfluxlat, &
     1928               yz0_new, SFRWL, yalb_dir_new, yalb_dif_new, yevap, yfluxsens,yfluxlat,&
     1929!albedo SB <<<
    18181930               ytsurf_new, y_dflux_t, y_dflux_q, &
    18191931               y_flux_u1, y_flux_v1)
    … …  
    21852297          evap(i,nsrf) = - flux_q(i,1,nsrf)                  !jyg
    21862298          d_ts(i,nsrf) = y_d_ts(j)
    2187           alb1(i,nsrf) = yalb1_new(j) 
    2188           alb2(i,nsrf) = yalb2_new(j)
     2299!albedo SB >>>
     2300!          alb1(i,nsrf) = yalb1_new(j) 
     2301!          alb2(i,nsrf) = yalb2_new(j)
     2302          do k=1,nsw
     2303          alb_dir(i,k,nsrf) = yalb_dir_new(j,k)
     2304          alb_dif(i,k,nsrf) = yalb_dif_new(j,k)
     2305          enddo
     2306!albedo SB <<<
    21892307          snow(i,nsrf) = ysnow(j) 
    21902308          qsurf(i,nsrf) = yqsurf(j)
    … …  
    29303048!****************************************************************************************
    29313049!
    2932   SUBROUTINE pbl_surface_newfrac(itime, pctsrf_new, pctsrf_old, tsurf, alb1, alb2, ustar, u10m, v10m, tke)
    2933 
     3050
     3051!albedo SB >>>
     3052!  SUBROUTINE pbl_surface_newfrac(itime, pctsrf_new, pctsrf_old, tsurf, alb1, alb2, ustar, u10m, v10m, tke)
     3053SUBROUTINE pbl_surface_newfrac(itime, pctsrf_new, pctsrf_old, tsurf,alb_dir,alb_dif, ustar, u10m, v10m, tke) 
     3054!albedo SB <<<
    29343055    ! Give default values where new fraction has appread
    29353056
    … …  
    29483069!****************************************************************************************
    29493070    REAL, DIMENSION(klon,nbsrf), INTENT(INOUT)        :: tsurf
    2950     REAL, DIMENSION(klon,nbsrf), INTENT(INOUT)        :: alb1, alb2
     3071!albedo SB >>>
     3072!   REAL, DIMENSION(klon,nbsrf), INTENT(INOUT)        :: alb1, alb2
     3073    REAL, DIMENSION(klon,nsw,nbsrf), INTENT(INOUT)       :: alb_dir, alb_dif
     3074    INTEGER :: k
     3075!albedo SB <<<
    29513076    REAL, DIMENSION(klon,nbsrf), INTENT(INOUT)        :: ustar,u10m, v10m
    29523077    REAL, DIMENSION(klon,klev+1,nbsrf+1), INTENT(INOUT) :: tke
    … …  
    29933118                rugos(i,nsrf) = rugos(i,nsrf_comp1)
    29943119                tsurf(i,nsrf) = tsurf(i,nsrf_comp1)
    2995                 alb1(i,nsrf)  = alb1(i,nsrf_comp1)
    2996                 alb2(i,nsrf)  = alb2(i,nsrf_comp1)
     3120!albedo SB >>>
     3121!                alb1(i,nsrf)  = alb1(i,nsrf_comp1)
     3122!                alb2(i,nsrf)  = alb2(i,nsrf_comp1)
     3123                DO k=1,nsw
     3124                 alb_dir(i,k,nsrf)=alb_dir(i,k,nsrf_comp1)
     3125                 alb_dif(i,k,nsrf)=alb_dif(i,k,nsrf_comp1)
     3126                ENDDO
     3127!albedo SB <<<
    29973128                ustar(i,nsrf)  = ustar(i,nsrf_comp1)
    29983129                u10m(i,nsrf)  = u10m(i,nsrf_comp1)
    … …  
    30083139                rugos(i,nsrf) = rugos(i,nsrf_comp2)*pctsrf_old(i,nsrf_comp2) + rugos(i,nsrf_comp3)*pctsrf_old(i,nsrf_comp3)
    30093140                tsurf(i,nsrf) = tsurf(i,nsrf_comp2)*pctsrf_old(i,nsrf_comp2) + tsurf(i,nsrf_comp3)*pctsrf_old(i,nsrf_comp3)
    3010                 alb1(i,nsrf)  = alb1(i,nsrf_comp2) *pctsrf_old(i,nsrf_comp2) + alb1(i,nsrf_comp3) *pctsrf_old(i,nsrf_comp3)
    3011                 alb2(i,nsrf)  = alb2(i,nsrf_comp2) *pctsrf_old(i,nsrf_comp2) + alb2(i,nsrf_comp3) *pctsrf_old(i,nsrf_comp3)
     3141!albedo SB >>>
     3142!                alb1(i,nsrf)  = alb1(i,nsrf_comp2) *pctsrf_old(i,nsrf_comp2) + alb1(i,nsrf_comp3) *pctsrf_old(i,nsrf_comp3)
     3143!                alb2(i,nsrf)  = alb2(i,nsrf_comp2) *pctsrf_old(i,nsrf_comp2) + alb2(i,nsrf_comp3) *pctsrf_old(i,nsrf_comp3)
     3144                DO k=1,nsw
     3145                 alb_dir(i,k,nsrf)=alb_dir(i,k,nsrf_comp2)*pctsrf_old(i,nsrf_comp2)+&
     3146                                        alb_dir(i,k,nsrf_comp3)*pctsrf_old(i,nsrf_comp3)
     3147                 alb_dif(i,k,nsrf)=alb_dif(i,k,nsrf_comp2)*pctsrf_old(i,nsrf_comp2)+&
     3148                                        alb_dif(i,k,nsrf_comp3)*pctsrf_old(i,nsrf_comp3)
     3149                ENDDO
     3150!albedo SB <<<
    30123151                ustar(i,nsrf)  = ustar(i,nsrf_comp2) *pctsrf_old(i,nsrf_comp2) + ustar(i,nsrf_comp3) *pctsrf_old(i,nsrf_comp3)
    30133152                u10m(i,nsrf)  = u10m(i,nsrf_comp2) *pctsrf_old(i,nsrf_comp2) + u10m(i,nsrf_comp3) *pctsrf_old(i,nsrf_comp3)
Note: See TracChangeset for help on using the changeset viewer.