[781] | 1 | ! |
---|
[2227] | 2 | ! $Id: surf_seaice_mod.F90 3815 2021-02-01 14:30:57Z jyg $ |
---|
| 3 | ! |
---|
[781] | 4 | MODULE surf_seaice_mod |
---|
| 5 | |
---|
| 6 | IMPLICIT NONE |
---|
| 7 | |
---|
| 8 | CONTAINS |
---|
| 9 | ! |
---|
| 10 | !**************************************************************************************** |
---|
| 11 | ! |
---|
| 12 | SUBROUTINE surf_seaice( & |
---|
[888] | 13 | rlon, rlat, swnet, lwnet, alb1, fder, & |
---|
[781] | 14 | itime, dtime, jour, knon, knindex, & |
---|
[1067] | 15 | lafin, & |
---|
| 16 | tsurf, p1lay, cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum, & |
---|
| 17 | AcoefH, AcoefQ, BcoefH, BcoefQ, & |
---|
| 18 | AcoefU, AcoefV, BcoefU, BcoefV, & |
---|
[2244] | 19 | ps, u1, v1, gustiness, pctsrf, & |
---|
[888] | 20 | snow, qsurf, qsol, agesno, tsoil, & |
---|
[2243] | 21 | z0m, z0h, SFRWL, alb_dir_new, alb_dif_new, evap, fluxsens, fluxlat, & |
---|
[1067] | 22 | tsurf_new, dflux_s, dflux_l, & |
---|
| 23 | flux_u1, flux_v1) |
---|
| 24 | |
---|
| 25 | USE dimphy |
---|
| 26 | USE surface_data |
---|
| 27 | USE ocean_forced_mod, ONLY : ocean_forced_ice |
---|
| 28 | USE ocean_cpl_mod, ONLY : ocean_cpl_ice |
---|
[2209] | 29 | USE ocean_slab_mod, ONLY : ocean_slab_ice |
---|
[1785] | 30 | USE indice_sol_mod |
---|
[1067] | 31 | |
---|
[781] | 32 | ! |
---|
| 33 | ! This subroutine will make a call to ocean_XXX_ice according to the ocean mode (force, |
---|
| 34 | ! slab or couple). The calculation of rugosity for the sea-ice surface is also done |
---|
| 35 | ! in here because it is the same calculation for the different modes of ocean. |
---|
| 36 | ! |
---|
| 37 | INCLUDE "dimsoil.h" |
---|
[2227] | 38 | INCLUDE "clesphys.h" |
---|
[781] | 39 | |
---|
[3815] | 40 | INCLUDE "YOMCST.h" |
---|
| 41 | ! for rd and retv |
---|
| 42 | |
---|
[781] | 43 | ! Input arguments |
---|
| 44 | !**************************************************************************************** |
---|
| 45 | INTEGER, INTENT(IN) :: itime, jour, knon |
---|
| 46 | INTEGER, DIMENSION(klon), INTENT(IN) :: knindex |
---|
[1067] | 47 | LOGICAL, INTENT(IN) :: lafin |
---|
[781] | 48 | REAL, INTENT(IN) :: dtime |
---|
| 49 | REAL, DIMENSION(klon), INTENT(IN) :: rlon, rlat |
---|
[888] | 50 | REAL, DIMENSION(klon), INTENT(IN) :: swnet ! net shortwave radiation at surface |
---|
| 51 | REAL, DIMENSION(klon), INTENT(IN) :: lwnet ! net longwave radiation at surface |
---|
| 52 | REAL, DIMENSION(klon), INTENT(IN) :: alb1 ! albedo in visible SW interval |
---|
[781] | 53 | REAL, DIMENSION(klon), INTENT(IN) :: fder |
---|
| 54 | REAL, DIMENSION(klon), INTENT(IN) :: tsurf |
---|
| 55 | REAL, DIMENSION(klon), INTENT(IN) :: p1lay |
---|
[1072] | 56 | REAL, DIMENSION(klon), INTENT(IN) :: cdragh, cdragm |
---|
[781] | 57 | REAL, DIMENSION(klon), INTENT(IN) :: precip_rain, precip_snow |
---|
| 58 | REAL, DIMENSION(klon), INTENT(IN) :: temp_air, spechum |
---|
[1067] | 59 | REAL, DIMENSION(klon), INTENT(IN) :: AcoefH, AcoefQ, BcoefH, BcoefQ |
---|
| 60 | REAL, DIMENSION(klon), INTENT(IN) :: AcoefU, AcoefV, BcoefU, BcoefV |
---|
[781] | 61 | REAL, DIMENSION(klon), INTENT(IN) :: ps |
---|
[2240] | 62 | REAL, DIMENSION(klon), INTENT(IN) :: u1, v1, gustiness |
---|
[781] | 63 | REAL, DIMENSION(klon,nbsrf), INTENT(IN) :: pctsrf |
---|
| 64 | |
---|
| 65 | ! In/Output arguments |
---|
| 66 | !**************************************************************************************** |
---|
| 67 | REAL, DIMENSION(klon), INTENT(INOUT) :: snow, qsurf, qsol |
---|
| 68 | REAL, DIMENSION(klon), INTENT(INOUT) :: agesno |
---|
| 69 | REAL, DIMENSION(klon, nsoilmx), INTENT(INOUT) :: tsoil |
---|
| 70 | |
---|
| 71 | ! Output arguments |
---|
| 72 | !**************************************************************************************** |
---|
[2243] | 73 | REAL, DIMENSION(klon), INTENT(OUT) :: z0m, z0h |
---|
[2227] | 74 | !albedo SB >>> |
---|
| 75 | ! REAL, DIMENSION(klon), INTENT(OUT) :: alb1_new ! new albedo in visible SW interval |
---|
| 76 | ! REAL, DIMENSION(klon), INTENT(OUT) :: alb2_new ! new albedo in near IR interval |
---|
| 77 | REAL, DIMENSION(6), INTENT(IN) :: SFRWL |
---|
| 78 | REAL, DIMENSION(klon,nsw), INTENT(OUT) :: alb_dir_new,alb_dif_new |
---|
| 79 | !albedo SB <<< |
---|
[781] | 80 | REAL, DIMENSION(klon), INTENT(OUT) :: evap, fluxsens, fluxlat |
---|
[888] | 81 | REAL, DIMENSION(klon), INTENT(OUT) :: tsurf_new |
---|
[781] | 82 | REAL, DIMENSION(klon), INTENT(OUT) :: dflux_s, dflux_l |
---|
[1067] | 83 | REAL, DIMENSION(klon), INTENT(OUT) :: flux_u1, flux_v1 |
---|
[781] | 84 | |
---|
[888] | 85 | ! Local arguments |
---|
| 86 | !**************************************************************************************** |
---|
| 87 | REAL, DIMENSION(klon) :: radsol |
---|
[1067] | 88 | |
---|
[2227] | 89 | !albedo SB >>> |
---|
| 90 | REAL, DIMENSION(klon) :: alb1_new,alb2_new |
---|
| 91 | !albedo SB <<< |
---|
[3815] | 92 | |
---|
| 93 | real rhoa(knon) ! density of moist air (kg / m3) |
---|
| 94 | |
---|
[781] | 95 | ! End definitions |
---|
| 96 | !**************************************************************************************** |
---|
| 97 | |
---|
[888] | 98 | |
---|
[781] | 99 | !**************************************************************************************** |
---|
[888] | 100 | ! Calculate total net radiance at surface |
---|
| 101 | ! |
---|
| 102 | !**************************************************************************************** |
---|
| 103 | radsol(:) = 0.0 |
---|
| 104 | radsol(1:knon) = swnet(1:knon) + lwnet(1:knon) |
---|
| 105 | |
---|
[3815] | 106 | rhoa = PS(:KNON) / (Rd * temp_air(:knon) * (1. + retv * spechum(:knon))) |
---|
| 107 | |
---|
[888] | 108 | !**************************************************************************************** |
---|
[781] | 109 | ! Switch according to type of ocean (couple, slab or forced) |
---|
| 110 | ! |
---|
| 111 | !**************************************************************************************** |
---|
[996] | 112 | IF (type_ocean == 'couple') THEN |
---|
| 113 | |
---|
[781] | 114 | CALL ocean_cpl_ice( & |
---|
[996] | 115 | rlon, rlat, swnet, lwnet, alb1, & |
---|
| 116 | fder, & |
---|
| 117 | itime, dtime, knon, knindex, & |
---|
| 118 | lafin,& |
---|
[1067] | 119 | p1lay, cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum,& |
---|
| 120 | AcoefH, AcoefQ, BcoefH, BcoefQ, & |
---|
| 121 | AcoefU, AcoefV, BcoefU, BcoefV, & |
---|
[2240] | 122 | ps, u1, v1, gustiness, pctsrf, & |
---|
[996] | 123 | radsol, snow, qsurf, & |
---|
[1067] | 124 | alb1_new, alb2_new, evap, fluxsens, fluxlat, flux_u1, flux_v1, & |
---|
[3815] | 125 | tsurf_new, dflux_s, dflux_l, rhoa) |
---|
[996] | 126 | |
---|
[2209] | 127 | ELSE IF (type_ocean == 'slab'.AND.version_ocean=='sicINT') THEN |
---|
| 128 | CALL ocean_slab_ice( & |
---|
| 129 | itime, dtime, jour, knon, knindex, & |
---|
| 130 | tsurf, p1lay, cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum,& |
---|
| 131 | AcoefH, AcoefQ, BcoefH, BcoefQ, & |
---|
| 132 | AcoefU, AcoefV, BcoefU, BcoefV, & |
---|
[2240] | 133 | ps, u1, v1, gustiness, & |
---|
[2209] | 134 | radsol, snow, qsurf, qsol, agesno, & |
---|
| 135 | alb1_new, alb2_new, evap, fluxsens, fluxlat, flux_u1, flux_v1, & |
---|
| 136 | tsurf_new, dflux_s, dflux_l, swnet) |
---|
| 137 | |
---|
[2057] | 138 | ELSE ! type_ocean=force or slab +sicOBS or sicNO |
---|
[1067] | 139 | CALL ocean_forced_ice( & |
---|
| 140 | itime, dtime, jour, knon, knindex, & |
---|
| 141 | tsurf, p1lay, cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum,& |
---|
| 142 | AcoefH, AcoefQ, BcoefH, BcoefQ, & |
---|
| 143 | AcoefU, AcoefV, BcoefU, BcoefV, & |
---|
[2240] | 144 | ps, u1, v1, gustiness, & |
---|
[996] | 145 | radsol, snow, qsol, agesno, tsoil, & |
---|
[1067] | 146 | qsurf, alb1_new, alb2_new, evap, fluxsens, fluxlat, flux_u1, flux_v1, & |
---|
[3815] | 147 | tsurf_new, dflux_s, dflux_l, rhoa) |
---|
[781] | 148 | |
---|
[996] | 149 | END IF |
---|
| 150 | |
---|
[781] | 151 | !**************************************************************************************** |
---|
| 152 | ! Calculate rugosity |
---|
| 153 | ! |
---|
| 154 | !**************************************************************************************** |
---|
| 155 | |
---|
[2243] | 156 | z0m=z0m_seaice |
---|
| 157 | z0h = z0h_seaice |
---|
[2227] | 158 | |
---|
| 159 | !albedo SB >>> |
---|
| 160 | select case(NSW) |
---|
| 161 | case(2) |
---|
| 162 | alb_dir_new(1:knon,1)=alb1_new(1:knon) |
---|
| 163 | alb_dir_new(1:knon,2)=alb2_new(1:knon) |
---|
| 164 | case(4) |
---|
| 165 | alb_dir_new(1:knon,1)=alb1_new(1:knon) |
---|
| 166 | alb_dir_new(1:knon,2)=alb2_new(1:knon) |
---|
| 167 | alb_dir_new(1:knon,3)=alb2_new(1:knon) |
---|
| 168 | alb_dir_new(1:knon,4)=alb2_new(1:knon) |
---|
| 169 | case(6) |
---|
| 170 | alb_dir_new(1:knon,1)=alb1_new(1:knon) |
---|
| 171 | alb_dir_new(1:knon,2)=alb1_new(1:knon) |
---|
| 172 | alb_dir_new(1:knon,3)=alb1_new(1:knon) |
---|
| 173 | alb_dir_new(1:knon,4)=alb2_new(1:knon) |
---|
| 174 | alb_dir_new(1:knon,5)=alb2_new(1:knon) |
---|
| 175 | alb_dir_new(1:knon,6)=alb2_new(1:knon) |
---|
| 176 | end select |
---|
| 177 | alb_dif_new=alb_dir_new |
---|
| 178 | !albedo SB <<< |
---|
| 179 | |
---|
| 180 | |
---|
| 181 | |
---|
| 182 | |
---|
[781] | 183 | END SUBROUTINE surf_seaice |
---|
| 184 | ! |
---|
| 185 | !**************************************************************************************** |
---|
| 186 | ! |
---|
| 187 | END MODULE surf_seaice_mod |
---|
| 188 | |
---|