- Timestamp:
- Dec 6, 2022, 12:01:16 AM (22 months ago)
- Location:
- LMDZ6/branches/Ocean_skin
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Ocean_skin
- Property svn:mergeinfo changed
-
LMDZ6/branches/Ocean_skin/libf/phylmd/pbl_surface_mod.F90
r4020 r4368 14 14 USE mod_grid_phy_lmdz, ONLY : klon_glo 15 15 USE ioipsl 16 USE surface_data, ONLY : type_ocean, ok_veget 16 USE surface_data, ONLY : type_ocean, ok_veget, landice_opt 17 17 USE surf_land_mod, ONLY : surf_land 18 18 USE surf_landice_mod, ONLY : surf_landice … … 2060 2060 CASE(is_lic) 2061 2061 ! Martin 2062 CALL surf_landice(itap, dtime, knon, ni, & 2063 rlon, rlat, debut, lafin, & 2064 yrmu0, ylwdown, yalb, zgeo1, & 2065 ysolsw, ysollw, yts, ypplay(:,1), & 2066 !!jyg ycdragh, ycdragm, yrain_f, ysnow_f, yt(:,1), yq(:,1),& 2067 ycdragh, ycdragm, yrain_f, ysnow_f, yt1, yq1,& 2068 AcoefH, AcoefQ, BcoefH, BcoefQ, & 2069 AcoefU, AcoefV, BcoefU, BcoefV, & 2070 ypsref, yu1, yv1, ygustiness, yrugoro, pctsrf, & 2071 ysnow, yqsurf, yqsol, yagesno, & 2072 ytsoil, yz0m, yz0h, SFRWL, yalb_dir_new, yalb_dif_new, yevap,yfluxsens,yfluxlat, & 2073 ytsurf_new, y_dflux_t, y_dflux_q, & 2074 yzmea, yzsig, ycldt, & 2075 ysnowhgt, yqsnow, ytoice, ysissnow, & 2076 yalb3_new, yrunoff, & 2077 y_flux_u1, y_flux_v1) 2078 2079 !jyg< 2080 !! alb3_lic(:)=0. 2081 !>jyg 2082 DO j = 1, knon 2083 i = ni(j) 2084 alb3_lic(i) = yalb3_new(j) 2085 snowhgt(i) = ysnowhgt(j) 2086 qsnow(i) = yqsnow(j) 2087 to_ice(i) = ytoice(j) 2088 sissnow(i) = ysissnow(j) 2089 runoff(i) = yrunoff(j) 2090 ENDDO 2091 ! Martin 2092 ! Special DICE MPL 05082013 puis BOMEX MPL 20150410 2093 IF (ok_prescr_ust) THEN 2094 DO j=1,knon 2095 y_flux_u1(j)=ycdragm(j)*(1.+sqrt(yu(j,1)*yu(j,1)+yv(j,1)*yv(j,1)))*yu(j,1)*ypplay(j,1)/RD/yt(j,1) 2096 y_flux_v1(j)=ycdragm(j)*(1.+sqrt(yu(j,1)*yu(j,1)+yv(j,1)*yv(j,1)))*yv(j,1)*ypplay(j,1)/RD/yt(j,1) 2097 ENDDO 2098 ENDIF 2062 2063 IF (landice_opt .LT. 2) THEN 2064 ! Land ice is treated by LMDZ and not by ORCHIDEE 2065 CALL surf_landice(itap, dtime, knon, ni, & 2066 rlon, rlat, debut, lafin, & 2067 yrmu0, ylwdown, yalb, zgeo1, & 2068 ysolsw, ysollw, yts, ypplay(:,1), & 2069 !!jyg ycdragh, ycdragm, yrain_f, ysnow_f, yt(:,1), yq(:,1),& 2070 ycdragh, ycdragm, yrain_f, ysnow_f, yt1, yq1,& 2071 AcoefH, AcoefQ, BcoefH, BcoefQ, & 2072 AcoefU, AcoefV, BcoefU, BcoefV, & 2073 ypsref, yu1, yv1, ygustiness, yrugoro, pctsrf, & 2074 ysnow, yqsurf, yqsol, yagesno, & 2075 ytsoil, yz0m, yz0h, SFRWL, yalb_dir_new, yalb_dif_new, yevap,yfluxsens,yfluxlat, & 2076 ytsurf_new, y_dflux_t, y_dflux_q, & 2077 yzmea, yzsig, ycldt, & 2078 ysnowhgt, yqsnow, ytoice, ysissnow, & 2079 yalb3_new, yrunoff, & 2080 y_flux_u1, y_flux_v1) 2081 2082 !jyg< 2083 !! alb3_lic(:)=0. 2084 !>jyg 2085 DO j = 1, knon 2086 i = ni(j) 2087 alb3_lic(i) = yalb3_new(j) 2088 snowhgt(i) = ysnowhgt(j) 2089 qsnow(i) = yqsnow(j) 2090 to_ice(i) = ytoice(j) 2091 sissnow(i) = ysissnow(j) 2092 runoff(i) = yrunoff(j) 2093 ENDDO 2094 ! Martin 2095 ! Special DICE MPL 05082013 puis BOMEX MPL 20150410 2096 IF (ok_prescr_ust) THEN 2097 DO j=1,knon 2098 y_flux_u1(j)=ycdragm(j)*(1.+sqrt(yu(j,1)*yu(j,1)+yv(j,1)*yv(j,1)))*yu(j,1)*ypplay(j,1)/RD/yt(j,1) 2099 y_flux_v1(j)=ycdragm(j)*(1.+sqrt(yu(j,1)*yu(j,1)+yv(j,1)*yv(j,1)))*yv(j,1)*ypplay(j,1)/RD/yt(j,1) 2100 ENDDO 2101 ENDIF 2102 2103 END IF 2099 2104 2100 2105 CASE(is_oce)
Note: See TracChangeset
for help on using the changeset viewer.