[781] | 1 | ! |
---|
[2227] | 2 | ! $Id: surf_seaice_mod.F90 3411 2018-11-05 14:24:59Z dcugnet $ |
---|
| 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 | |
---|
| 40 | ! Input arguments |
---|
| 41 | !**************************************************************************************** |
---|
| 42 | INTEGER, INTENT(IN) :: itime, jour, knon |
---|
| 43 | INTEGER, DIMENSION(klon), INTENT(IN) :: knindex |
---|
[1067] | 44 | LOGICAL, INTENT(IN) :: lafin |
---|
[781] | 45 | REAL, INTENT(IN) :: dtime |
---|
| 46 | REAL, DIMENSION(klon), INTENT(IN) :: rlon, rlat |
---|
[888] | 47 | REAL, DIMENSION(klon), INTENT(IN) :: swnet ! net shortwave radiation at surface |
---|
| 48 | REAL, DIMENSION(klon), INTENT(IN) :: lwnet ! net longwave radiation at surface |
---|
| 49 | REAL, DIMENSION(klon), INTENT(IN) :: alb1 ! albedo in visible SW interval |
---|
[781] | 50 | REAL, DIMENSION(klon), INTENT(IN) :: fder |
---|
| 51 | REAL, DIMENSION(klon), INTENT(IN) :: tsurf |
---|
| 52 | REAL, DIMENSION(klon), INTENT(IN) :: p1lay |
---|
[1072] | 53 | REAL, DIMENSION(klon), INTENT(IN) :: cdragh, cdragm |
---|
[781] | 54 | REAL, DIMENSION(klon), INTENT(IN) :: precip_rain, precip_snow |
---|
| 55 | REAL, DIMENSION(klon), INTENT(IN) :: temp_air, spechum |
---|
[1067] | 56 | REAL, DIMENSION(klon), INTENT(IN) :: AcoefH, AcoefQ, BcoefH, BcoefQ |
---|
| 57 | REAL, DIMENSION(klon), INTENT(IN) :: AcoefU, AcoefV, BcoefU, BcoefV |
---|
[781] | 58 | REAL, DIMENSION(klon), INTENT(IN) :: ps |
---|
[2240] | 59 | REAL, DIMENSION(klon), INTENT(IN) :: u1, v1, gustiness |
---|
[781] | 60 | REAL, DIMENSION(klon,nbsrf), INTENT(IN) :: pctsrf |
---|
| 61 | |
---|
| 62 | ! In/Output arguments |
---|
| 63 | !**************************************************************************************** |
---|
| 64 | REAL, DIMENSION(klon), INTENT(INOUT) :: snow, qsurf, qsol |
---|
| 65 | REAL, DIMENSION(klon), INTENT(INOUT) :: agesno |
---|
| 66 | REAL, DIMENSION(klon, nsoilmx), INTENT(INOUT) :: tsoil |
---|
| 67 | |
---|
| 68 | ! Output arguments |
---|
| 69 | !**************************************************************************************** |
---|
[2243] | 70 | REAL, DIMENSION(klon), INTENT(OUT) :: z0m, z0h |
---|
[2227] | 71 | !albedo SB >>> |
---|
| 72 | ! REAL, DIMENSION(klon), INTENT(OUT) :: alb1_new ! new albedo in visible SW interval |
---|
| 73 | ! REAL, DIMENSION(klon), INTENT(OUT) :: alb2_new ! new albedo in near IR interval |
---|
| 74 | REAL, DIMENSION(6), INTENT(IN) :: SFRWL |
---|
| 75 | REAL, DIMENSION(klon,nsw), INTENT(OUT) :: alb_dir_new,alb_dif_new |
---|
| 76 | !albedo SB <<< |
---|
[781] | 77 | REAL, DIMENSION(klon), INTENT(OUT) :: evap, fluxsens, fluxlat |
---|
[888] | 78 | REAL, DIMENSION(klon), INTENT(OUT) :: tsurf_new |
---|
[781] | 79 | REAL, DIMENSION(klon), INTENT(OUT) :: dflux_s, dflux_l |
---|
[1067] | 80 | REAL, DIMENSION(klon), INTENT(OUT) :: flux_u1, flux_v1 |
---|
[781] | 81 | |
---|
[888] | 82 | ! Local arguments |
---|
| 83 | !**************************************************************************************** |
---|
| 84 | REAL, DIMENSION(klon) :: radsol |
---|
[1067] | 85 | |
---|
[2227] | 86 | !albedo SB >>> |
---|
| 87 | REAL, DIMENSION(klon) :: alb1_new,alb2_new |
---|
| 88 | !albedo SB <<< |
---|
[888] | 89 | ! |
---|
[781] | 90 | ! End definitions |
---|
| 91 | !**************************************************************************************** |
---|
| 92 | |
---|
[888] | 93 | |
---|
[781] | 94 | !**************************************************************************************** |
---|
[888] | 95 | ! Calculate total net radiance at surface |
---|
| 96 | ! |
---|
| 97 | !**************************************************************************************** |
---|
| 98 | radsol(:) = 0.0 |
---|
| 99 | radsol(1:knon) = swnet(1:knon) + lwnet(1:knon) |
---|
| 100 | |
---|
| 101 | !**************************************************************************************** |
---|
[781] | 102 | ! Switch according to type of ocean (couple, slab or forced) |
---|
| 103 | ! |
---|
| 104 | !**************************************************************************************** |
---|
[996] | 105 | IF (type_ocean == 'couple') THEN |
---|
| 106 | |
---|
[781] | 107 | CALL ocean_cpl_ice( & |
---|
[996] | 108 | rlon, rlat, swnet, lwnet, alb1, & |
---|
| 109 | fder, & |
---|
| 110 | itime, dtime, knon, knindex, & |
---|
| 111 | lafin,& |
---|
[1067] | 112 | p1lay, cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum,& |
---|
| 113 | AcoefH, AcoefQ, BcoefH, BcoefQ, & |
---|
| 114 | AcoefU, AcoefV, BcoefU, BcoefV, & |
---|
[2240] | 115 | ps, u1, v1, gustiness, pctsrf, & |
---|
[996] | 116 | radsol, snow, qsurf, & |
---|
[1067] | 117 | alb1_new, alb2_new, evap, fluxsens, fluxlat, flux_u1, flux_v1, & |
---|
[996] | 118 | tsurf_new, dflux_s, dflux_l) |
---|
| 119 | |
---|
[2209] | 120 | ELSE IF (type_ocean == 'slab'.AND.version_ocean=='sicINT') THEN |
---|
| 121 | CALL ocean_slab_ice( & |
---|
| 122 | itime, dtime, jour, knon, knindex, & |
---|
| 123 | tsurf, p1lay, cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum,& |
---|
| 124 | AcoefH, AcoefQ, BcoefH, BcoefQ, & |
---|
| 125 | AcoefU, AcoefV, BcoefU, BcoefV, & |
---|
[2240] | 126 | ps, u1, v1, gustiness, & |
---|
[2209] | 127 | radsol, snow, qsurf, qsol, agesno, & |
---|
| 128 | alb1_new, alb2_new, evap, fluxsens, fluxlat, flux_u1, flux_v1, & |
---|
| 129 | tsurf_new, dflux_s, dflux_l, swnet) |
---|
| 130 | |
---|
[2057] | 131 | ELSE ! type_ocean=force or slab +sicOBS or sicNO |
---|
[1067] | 132 | CALL ocean_forced_ice( & |
---|
| 133 | itime, dtime, jour, knon, knindex, & |
---|
| 134 | tsurf, p1lay, cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum,& |
---|
| 135 | AcoefH, AcoefQ, BcoefH, BcoefQ, & |
---|
| 136 | AcoefU, AcoefV, BcoefU, BcoefV, & |
---|
[2240] | 137 | ps, u1, v1, gustiness, & |
---|
[996] | 138 | radsol, snow, qsol, agesno, tsoil, & |
---|
[1067] | 139 | qsurf, alb1_new, alb2_new, evap, fluxsens, fluxlat, flux_u1, flux_v1, & |
---|
[996] | 140 | tsurf_new, dflux_s, dflux_l) |
---|
[781] | 141 | |
---|
[996] | 142 | END IF |
---|
| 143 | |
---|
[781] | 144 | !**************************************************************************************** |
---|
| 145 | ! Calculate rugosity |
---|
| 146 | ! |
---|
| 147 | !**************************************************************************************** |
---|
| 148 | |
---|
[2243] | 149 | z0m=z0m_seaice |
---|
| 150 | z0h = z0h_seaice |
---|
[2227] | 151 | |
---|
| 152 | !albedo SB >>> |
---|
| 153 | select case(NSW) |
---|
| 154 | case(2) |
---|
| 155 | alb_dir_new(1:knon,1)=alb1_new(1:knon) |
---|
| 156 | alb_dir_new(1:knon,2)=alb2_new(1:knon) |
---|
| 157 | case(4) |
---|
| 158 | alb_dir_new(1:knon,1)=alb1_new(1:knon) |
---|
| 159 | alb_dir_new(1:knon,2)=alb2_new(1:knon) |
---|
| 160 | alb_dir_new(1:knon,3)=alb2_new(1:knon) |
---|
| 161 | alb_dir_new(1:knon,4)=alb2_new(1:knon) |
---|
| 162 | case(6) |
---|
| 163 | alb_dir_new(1:knon,1)=alb1_new(1:knon) |
---|
| 164 | alb_dir_new(1:knon,2)=alb1_new(1:knon) |
---|
| 165 | alb_dir_new(1:knon,3)=alb1_new(1:knon) |
---|
| 166 | alb_dir_new(1:knon,4)=alb2_new(1:knon) |
---|
| 167 | alb_dir_new(1:knon,5)=alb2_new(1:knon) |
---|
| 168 | alb_dir_new(1:knon,6)=alb2_new(1:knon) |
---|
| 169 | end select |
---|
| 170 | alb_dif_new=alb_dir_new |
---|
| 171 | !albedo SB <<< |
---|
| 172 | |
---|
| 173 | |
---|
| 174 | |
---|
| 175 | |
---|
[781] | 176 | END SUBROUTINE surf_seaice |
---|
| 177 | ! |
---|
| 178 | !**************************************************************************************** |
---|
| 179 | ! |
---|
| 180 | END MODULE surf_seaice_mod |
---|
| 181 | |
---|