[782] | 1 | ! |
---|
[2538] | 2 | ! $Id: ocean_cpl_mod.F90 3463 2019-02-08 13:28:03Z lguez $ |
---|
| 3 | ! |
---|
[782] | 4 | MODULE ocean_cpl_mod |
---|
| 5 | ! |
---|
| 6 | ! This module is used both for the sub-surface ocean and sea-ice for the case of a |
---|
| 7 | ! coupled model configuration, ocean=couple. |
---|
| 8 | ! |
---|
| 9 | |
---|
| 10 | IMPLICIT NONE |
---|
| 11 | PRIVATE |
---|
| 12 | |
---|
[996] | 13 | PUBLIC :: ocean_cpl_init, ocean_cpl_noice, ocean_cpl_ice |
---|
[782] | 14 | |
---|
[2538] | 15 | |
---|
[782] | 16 | !**************************************************************************************** |
---|
| 17 | ! |
---|
| 18 | CONTAINS |
---|
| 19 | ! |
---|
| 20 | !**************************************************************************************** |
---|
| 21 | ! |
---|
| 22 | SUBROUTINE ocean_cpl_init(dtime, rlon, rlat) |
---|
| 23 | ! |
---|
| 24 | ! Allocate fields for this module and initailize the module mod_cpl |
---|
| 25 | ! |
---|
[1067] | 26 | USE dimphy, ONLY : klon |
---|
| 27 | USE cpl_mod |
---|
| 28 | |
---|
[782] | 29 | ! Input arguments |
---|
| 30 | !************************************************************************************* |
---|
| 31 | REAL, INTENT(IN) :: dtime |
---|
| 32 | REAL, DIMENSION(klon), INTENT(IN) :: rlon, rlat |
---|
| 33 | |
---|
| 34 | ! Local variables |
---|
| 35 | !************************************************************************************* |
---|
| 36 | INTEGER :: error |
---|
| 37 | CHARACTER (len = 80) :: abort_message |
---|
| 38 | CHARACTER (len = 20) :: modname = 'ocean_cpl_init' |
---|
| 39 | |
---|
| 40 | ! Initialize module cpl_init |
---|
| 41 | CALL cpl_init(dtime, rlon, rlat) |
---|
| 42 | |
---|
| 43 | END SUBROUTINE ocean_cpl_init |
---|
| 44 | ! |
---|
| 45 | !**************************************************************************************** |
---|
| 46 | ! |
---|
| 47 | SUBROUTINE ocean_cpl_noice( & |
---|
[888] | 48 | swnet, lwnet, alb1, & |
---|
[1067] | 49 | windsp, fder_old, & |
---|
[782] | 50 | itime, dtime, knon, knindex, & |
---|
[2254] | 51 | p1lay, cdragh, cdragq, cdragm, precip_rain, precip_snow, temp_air, spechum, & |
---|
[1067] | 52 | AcoefH, AcoefQ, BcoefH, BcoefQ, & |
---|
| 53 | AcoefU, AcoefV, BcoefU, BcoefV, & |
---|
[3463] | 54 | ps, u1, v1, gustiness, tsurf_in, & |
---|
[888] | 55 | radsol, snow, agesno, & |
---|
[1067] | 56 | qsurf, evap, fluxsens, fluxlat, flux_u1, flux_v1, & |
---|
[996] | 57 | tsurf_new, dflux_s, dflux_l) |
---|
[1067] | 58 | |
---|
[782] | 59 | ! |
---|
| 60 | ! This subroutine treats the "open ocean", all grid points that are not entierly covered |
---|
| 61 | ! by ice. The subroutine first receives fields from coupler, then some calculations at |
---|
| 62 | ! surface is done and finally it sends some fields to the coupler. |
---|
| 63 | ! |
---|
[1067] | 64 | USE dimphy, ONLY : klon |
---|
| 65 | USE cpl_mod |
---|
| 66 | USE calcul_fluxs_mod |
---|
[1785] | 67 | USE indice_sol_mod |
---|
[2538] | 68 | USE phys_output_var_mod, ONLY : sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o |
---|
| 69 | USE cpl_mod, ONLY : gath2cpl |
---|
[3463] | 70 | use config_ocean_skin_m, only: activate_ocean_skin |
---|
[1067] | 71 | |
---|
[793] | 72 | INCLUDE "YOMCST.h" |
---|
[2254] | 73 | INCLUDE "clesphys.h" |
---|
[782] | 74 | ! |
---|
| 75 | ! Input arguments |
---|
| 76 | !**************************************************************************************** |
---|
| 77 | INTEGER, INTENT(IN) :: itime, knon |
---|
| 78 | INTEGER, DIMENSION(klon), INTENT(IN) :: knindex |
---|
| 79 | REAL, INTENT(IN) :: dtime |
---|
[888] | 80 | REAL, DIMENSION(klon), INTENT(IN) :: swnet |
---|
| 81 | REAL, DIMENSION(klon), INTENT(IN) :: lwnet |
---|
| 82 | REAL, DIMENSION(klon), INTENT(IN) :: alb1 ! albedo in visible SW interval |
---|
[782] | 83 | REAL, DIMENSION(klon), INTENT(IN) :: windsp |
---|
| 84 | REAL, DIMENSION(klon), INTENT(IN) :: fder_old |
---|
| 85 | REAL, DIMENSION(klon), INTENT(IN) :: p1lay |
---|
[2254] | 86 | REAL, DIMENSION(klon), INTENT(IN) :: cdragh, cdragq, cdragm |
---|
[782] | 87 | REAL, DIMENSION(klon), INTENT(IN) :: precip_rain, precip_snow |
---|
| 88 | REAL, DIMENSION(klon), INTENT(IN) :: temp_air, spechum |
---|
[1067] | 89 | REAL, DIMENSION(klon), INTENT(IN) :: AcoefH, AcoefQ, BcoefH, BcoefQ |
---|
| 90 | REAL, DIMENSION(klon), INTENT(IN) :: AcoefU, AcoefV, BcoefU, BcoefV |
---|
[782] | 91 | REAL, DIMENSION(klon), INTENT(IN) :: ps |
---|
[2240] | 92 | REAL, DIMENSION(klon), INTENT(IN) :: u1, v1, gustiness |
---|
[3463] | 93 | REAL, DIMENSION(klon), INTENT(IN) :: tsurf_in |
---|
[782] | 94 | |
---|
| 95 | ! In/Output arguments |
---|
| 96 | !**************************************************************************************** |
---|
| 97 | REAL, DIMENSION(klon), INTENT(INOUT) :: radsol |
---|
| 98 | REAL, DIMENSION(klon), INTENT(INOUT) :: snow |
---|
| 99 | REAL, DIMENSION(klon), INTENT(INOUT) :: agesno |
---|
| 100 | |
---|
| 101 | ! Output arguments |
---|
| 102 | !**************************************************************************************** |
---|
| 103 | REAL, DIMENSION(klon), INTENT(OUT) :: qsurf |
---|
| 104 | REAL, DIMENSION(klon), INTENT(OUT) :: evap, fluxsens, fluxlat |
---|
[1067] | 105 | REAL, DIMENSION(klon), INTENT(OUT) :: flux_u1, flux_v1 |
---|
[782] | 106 | REAL, DIMENSION(klon), INTENT(OUT) :: tsurf_new |
---|
| 107 | REAL, DIMENSION(klon), INTENT(OUT) :: dflux_s, dflux_l |
---|
[2538] | 108 | |
---|
[782] | 109 | |
---|
| 110 | ! Local variables |
---|
| 111 | !**************************************************************************************** |
---|
[2538] | 112 | INTEGER :: i, j |
---|
[782] | 113 | INTEGER, DIMENSION(1) :: iloc |
---|
| 114 | REAL, DIMENSION(klon) :: cal, beta, dif_grnd |
---|
| 115 | REAL, DIMENSION(klon) :: fder_new |
---|
| 116 | REAL, DIMENSION(klon) :: tsurf_cpl |
---|
[1067] | 117 | REAL, DIMENSION(klon) :: u0_cpl, v0_cpl |
---|
| 118 | REAL, DIMENSION(klon) :: u1_lay, v1_lay |
---|
[782] | 119 | LOGICAL :: check=.FALSE. |
---|
[2538] | 120 | REAL, DIMENSION(klon) :: sens_prec_liq, sens_prec_sol |
---|
| 121 | REAL, DIMENSION(klon) :: lat_prec_liq, lat_prec_sol |
---|
[782] | 122 | |
---|
| 123 | ! End definitions |
---|
| 124 | !**************************************************************************************** |
---|
| 125 | |
---|
| 126 | IF (check) WRITE(*,*)' Entering ocean_cpl_noice' |
---|
| 127 | |
---|
| 128 | !**************************************************************************************** |
---|
[996] | 129 | ! Receive sea-surface temperature(tsurf_cpl) from coupler |
---|
[782] | 130 | ! |
---|
| 131 | !**************************************************************************************** |
---|
[1067] | 132 | CALL cpl_receive_ocean_fields(knon, knindex, tsurf_cpl, u0_cpl, v0_cpl) |
---|
[782] | 133 | |
---|
| 134 | !**************************************************************************************** |
---|
| 135 | ! Calculate fluxes at surface |
---|
| 136 | ! |
---|
| 137 | !**************************************************************************************** |
---|
| 138 | cal = 0. |
---|
| 139 | beta = 1. |
---|
| 140 | dif_grnd = 0. |
---|
| 141 | agesno(:) = 0. |
---|
[2538] | 142 | sens_prec_liq = 0.; sens_prec_sol = 0.; lat_prec_liq = 0.; lat_prec_sol = 0. |
---|
| 143 | |
---|
[1067] | 144 | |
---|
| 145 | DO i = 1, knon |
---|
| 146 | u1_lay(i) = u1(i) - u0_cpl(i) |
---|
| 147 | v1_lay(i) = v1(i) - v0_cpl(i) |
---|
| 148 | END DO |
---|
| 149 | |
---|
[782] | 150 | CALL calcul_fluxs(knon, is_oce, dtime, & |
---|
[3463] | 151 | merge(tsurf_in, tsurf_cpl, activate_ocean_skin == 2), p1lay, cal, & |
---|
| 152 | beta, cdragh, cdragq, ps, & |
---|
[782] | 153 | precip_rain, precip_snow, snow, qsurf, & |
---|
[2240] | 154 | radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, gustiness, & |
---|
[2254] | 155 | f_qsat_oce,AcoefH, AcoefQ, BcoefH, BcoefQ, & |
---|
[2538] | 156 | tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l, & |
---|
| 157 | sens_prec_liq, sens_prec_sol, lat_prec_liq, lat_prec_sol) |
---|
[3463] | 158 | if (activate_ocean_skin == 2) tsurf_new = tsurf_cpl |
---|
| 159 | |
---|
[2538] | 160 | do j = 1, knon |
---|
| 161 | i = knindex(j) |
---|
| 162 | sens_prec_liq_o(i,1) = sens_prec_liq(j) |
---|
| 163 | sens_prec_sol_o(i,1) = sens_prec_sol(j) |
---|
| 164 | lat_prec_liq_o(i,1) = lat_prec_liq(j) |
---|
| 165 | lat_prec_sol_o(i,1) = lat_prec_sol(j) |
---|
| 166 | enddo |
---|
| 167 | |
---|
| 168 | |
---|
[782] | 169 | |
---|
[1067] | 170 | ! - Flux calculation at first modele level for U and V |
---|
| 171 | CALL calcul_flux_wind(knon, dtime, & |
---|
[2240] | 172 | u0_cpl, v0_cpl, u1, v1, gustiness, cdragm, & |
---|
[1067] | 173 | AcoefU, AcoefV, BcoefU, BcoefV, & |
---|
| 174 | p1lay, temp_air, & |
---|
| 175 | flux_u1, flux_v1) |
---|
[782] | 176 | |
---|
| 177 | !**************************************************************************************** |
---|
| 178 | ! Calculate fder : flux derivative (sensible and latente) |
---|
| 179 | ! |
---|
| 180 | !**************************************************************************************** |
---|
| 181 | fder_new(:) = fder_old(:) + dflux_s(:) + dflux_l(:) |
---|
| 182 | |
---|
| 183 | iloc = MAXLOC(fder_new(1:klon)) |
---|
| 184 | IF (check .AND. fder_new(iloc(1))> 0.) THEN |
---|
| 185 | WRITE(*,*)'**** Debug fder****' |
---|
| 186 | WRITE(*,*)'max fder(',iloc(1),') = ',fder_new(iloc(1)) |
---|
| 187 | WRITE(*,*)'fder_old, dflux_s, dflux_l',fder_old(iloc(1)), & |
---|
| 188 | dflux_s(iloc(1)), dflux_l(iloc(1)) |
---|
| 189 | ENDIF |
---|
| 190 | |
---|
| 191 | !**************************************************************************************** |
---|
| 192 | ! Send and cumulate fields to the coupler |
---|
| 193 | ! |
---|
| 194 | !**************************************************************************************** |
---|
| 195 | |
---|
| 196 | CALL cpl_send_ocean_fields(itime, knon, knindex, & |
---|
[888] | 197 | swnet, lwnet, fluxlat, fluxsens, & |
---|
[2872] | 198 | precip_rain, precip_snow, evap, tsurf_new, fder_new, alb1, flux_u1, flux_v1, windsp,& |
---|
| 199 | sens_prec_liq, sens_prec_sol, lat_prec_liq, lat_prec_sol) |
---|
[782] | 200 | |
---|
| 201 | |
---|
| 202 | END SUBROUTINE ocean_cpl_noice |
---|
| 203 | ! |
---|
| 204 | !**************************************************************************************** |
---|
| 205 | ! |
---|
| 206 | SUBROUTINE ocean_cpl_ice( & |
---|
[888] | 207 | rlon, rlat, swnet, lwnet, alb1, & |
---|
[782] | 208 | fder_old, & |
---|
| 209 | itime, dtime, knon, knindex, & |
---|
| 210 | lafin, & |
---|
[1067] | 211 | p1lay, cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum, & |
---|
| 212 | AcoefH, AcoefQ, BcoefH, BcoefQ, & |
---|
| 213 | AcoefU, AcoefV, BcoefU, BcoefV, & |
---|
[2240] | 214 | ps, u1, v1, gustiness, pctsrf, & |
---|
[782] | 215 | radsol, snow, qsurf, & |
---|
[1067] | 216 | alb1_new, alb2_new, evap, fluxsens, fluxlat, flux_u1, flux_v1, & |
---|
[996] | 217 | tsurf_new, dflux_s, dflux_l) |
---|
[782] | 218 | ! |
---|
| 219 | ! This subroutine treats the ocean where there is ice. The subroutine first receives |
---|
| 220 | ! fields from coupler, then some calculations at surface is done and finally sends |
---|
| 221 | ! some fields to the coupler. |
---|
| 222 | ! |
---|
[1067] | 223 | USE dimphy, ONLY : klon |
---|
| 224 | USE cpl_mod |
---|
| 225 | USE calcul_fluxs_mod |
---|
[1785] | 226 | USE indice_sol_mod |
---|
[2538] | 227 | USE phys_output_var_mod, ONLY : sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o |
---|
[1067] | 228 | |
---|
[793] | 229 | INCLUDE "YOMCST.h" |
---|
[2254] | 230 | INCLUDE "clesphys.h" |
---|
[782] | 231 | |
---|
| 232 | ! Input arguments |
---|
| 233 | !**************************************************************************************** |
---|
| 234 | INTEGER, INTENT(IN) :: itime, knon |
---|
| 235 | INTEGER, DIMENSION(klon), INTENT(IN) :: knindex |
---|
| 236 | LOGICAL, INTENT(IN) :: lafin |
---|
| 237 | REAL, INTENT(IN) :: dtime |
---|
| 238 | REAL, DIMENSION(klon), INTENT(IN) :: rlon, rlat |
---|
[888] | 239 | REAL, DIMENSION(klon), INTENT(IN) :: swnet |
---|
| 240 | REAL, DIMENSION(klon), INTENT(IN) :: lwnet |
---|
| 241 | REAL, DIMENSION(klon), INTENT(IN) :: alb1 ! albedo in visible SW interval |
---|
[782] | 242 | REAL, DIMENSION(klon), INTENT(IN) :: fder_old |
---|
| 243 | REAL, DIMENSION(klon), INTENT(IN) :: p1lay |
---|
[1067] | 244 | REAL, DIMENSION(klon), INTENT(IN) :: cdragh, cdragm |
---|
[782] | 245 | REAL, DIMENSION(klon), INTENT(IN) :: precip_rain, precip_snow |
---|
| 246 | REAL, DIMENSION(klon), INTENT(IN) :: temp_air, spechum |
---|
[1067] | 247 | REAL, DIMENSION(klon), INTENT(IN) :: AcoefH, AcoefQ, BcoefH, BcoefQ |
---|
| 248 | REAL, DIMENSION(klon), INTENT(IN) :: AcoefU, AcoefV, BcoefU, BcoefV |
---|
[782] | 249 | REAL, DIMENSION(klon), INTENT(IN) :: ps |
---|
[2240] | 250 | REAL, DIMENSION(klon), INTENT(IN) :: u1, v1, gustiness |
---|
[996] | 251 | REAL, DIMENSION(klon,nbsrf), INTENT(IN) :: pctsrf |
---|
[782] | 252 | |
---|
| 253 | ! In/output arguments |
---|
| 254 | !**************************************************************************************** |
---|
| 255 | REAL, DIMENSION(klon), INTENT(INOUT) :: radsol |
---|
| 256 | REAL, DIMENSION(klon), INTENT(INOUT) :: snow |
---|
| 257 | |
---|
| 258 | ! Output arguments |
---|
| 259 | !**************************************************************************************** |
---|
| 260 | REAL, DIMENSION(klon), INTENT(OUT) :: qsurf |
---|
[888] | 261 | REAL, DIMENSION(klon), INTENT(OUT) :: alb1_new, alb2_new |
---|
[782] | 262 | REAL, DIMENSION(klon), INTENT(OUT) :: evap, fluxsens, fluxlat |
---|
[1067] | 263 | REAL, DIMENSION(klon), INTENT(OUT) :: flux_u1, flux_v1 |
---|
[888] | 264 | REAL, DIMENSION(klon), INTENT(OUT) :: tsurf_new |
---|
[782] | 265 | REAL, DIMENSION(klon), INTENT(OUT) :: dflux_s, dflux_l |
---|
[2538] | 266 | |
---|
[782] | 267 | |
---|
| 268 | ! Local variables |
---|
| 269 | !**************************************************************************************** |
---|
[2538] | 270 | INTEGER :: i, j |
---|
[782] | 271 | INTEGER, DIMENSION(1) :: iloc |
---|
| 272 | LOGICAL :: check=.FALSE. |
---|
| 273 | REAL, PARAMETER :: t_grnd=271.35 |
---|
| 274 | REAL, DIMENSION(klon) :: cal, beta, dif_grnd |
---|
| 275 | REAL, DIMENSION(klon) :: tsurf_cpl, fder_new |
---|
[888] | 276 | REAL, DIMENSION(klon) :: alb_cpl |
---|
[1067] | 277 | REAL, DIMENSION(klon) :: u0, v0 |
---|
| 278 | REAL, DIMENSION(klon) :: u1_lay, v1_lay |
---|
[2538] | 279 | REAL, DIMENSION(klon) :: sens_prec_liq, sens_prec_sol |
---|
| 280 | REAL, DIMENSION(klon) :: lat_prec_liq, lat_prec_sol |
---|
[782] | 281 | |
---|
| 282 | ! End definitions |
---|
| 283 | !**************************************************************************************** |
---|
| 284 | |
---|
| 285 | IF (check) WRITE(*,*)'Entering surface_seaice, knon=',knon |
---|
| 286 | |
---|
[2538] | 287 | sens_prec_liq = 0.; sens_prec_sol = 0.; lat_prec_liq = 0.; lat_prec_sol = 0. |
---|
| 288 | |
---|
[782] | 289 | !**************************************************************************************** |
---|
[996] | 290 | ! Receive ocean temperature(tsurf_cpl) and albedo(alb_new) from coupler |
---|
[782] | 291 | ! |
---|
| 292 | !**************************************************************************************** |
---|
| 293 | |
---|
| 294 | CALL cpl_receive_seaice_fields(knon, knindex, & |
---|
[1146] | 295 | tsurf_cpl, alb_cpl, u0, v0) |
---|
[888] | 296 | |
---|
| 297 | alb1_new(1:knon) = alb_cpl(1:knon) |
---|
| 298 | alb2_new(1:knon) = alb_cpl(1:knon) |
---|
| 299 | |
---|
[782] | 300 | |
---|
| 301 | !**************************************************************************************** |
---|
| 302 | ! Calculate fluxes at surface |
---|
| 303 | ! |
---|
| 304 | !**************************************************************************************** |
---|
| 305 | cal = 0. |
---|
| 306 | dif_grnd = 0. |
---|
| 307 | beta = 1.0 |
---|
| 308 | |
---|
[1146] | 309 | DO i = 1, knon |
---|
| 310 | u1_lay(i) = u1(i) - u0(i) |
---|
| 311 | v1_lay(i) = v1(i) - v0(i) |
---|
| 312 | END DO |
---|
[782] | 313 | |
---|
| 314 | CALL calcul_fluxs(knon, is_sic, dtime, & |
---|
[2254] | 315 | tsurf_cpl, p1lay, cal, beta, cdragh, cdragh, ps, & |
---|
[782] | 316 | precip_rain, precip_snow, snow, qsurf, & |
---|
[2240] | 317 | radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, gustiness, & |
---|
[2254] | 318 | f_qsat_oce,AcoefH, AcoefQ, BcoefH, BcoefQ, & |
---|
[2538] | 319 | tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l, & |
---|
| 320 | sens_prec_liq, sens_prec_sol, lat_prec_liq, lat_prec_sol) |
---|
| 321 | do j = 1, knon |
---|
| 322 | i = knindex(j) |
---|
| 323 | sens_prec_liq_o(i,2) = sens_prec_liq(j) |
---|
| 324 | sens_prec_sol_o(i,2) = sens_prec_sol(j) |
---|
| 325 | lat_prec_liq_o(i,2) = lat_prec_liq(j) |
---|
| 326 | lat_prec_sol_o(i,2) = lat_prec_sol(j) |
---|
| 327 | enddo |
---|
[782] | 328 | |
---|
[1067] | 329 | |
---|
| 330 | ! - Flux calculation at first modele level for U and V |
---|
| 331 | CALL calcul_flux_wind(knon, dtime, & |
---|
[2240] | 332 | u0, v0, u1, v1, gustiness, cdragm, & |
---|
[1067] | 333 | AcoefU, AcoefV, BcoefU, BcoefV, & |
---|
| 334 | p1lay, temp_air, & |
---|
| 335 | flux_u1, flux_v1) |
---|
| 336 | |
---|
[782] | 337 | !**************************************************************************************** |
---|
| 338 | ! Calculate fder : flux derivative (sensible and latente) |
---|
| 339 | ! |
---|
| 340 | !**************************************************************************************** |
---|
| 341 | fder_new(:) = fder_old(:) + dflux_s(:) + dflux_l(:) |
---|
| 342 | |
---|
| 343 | iloc = MAXLOC(fder_new(1:klon)) |
---|
| 344 | IF (check .AND. fder_new(iloc(1))> 0.) THEN |
---|
| 345 | WRITE(*,*)'**** Debug fder ****' |
---|
| 346 | WRITE(*,*)'max fder(',iloc(1),') = ',fder_new(iloc(1)) |
---|
| 347 | WRITE(*,*)'fder_old, dflux_s, dflux_l',fder_old(iloc(1)), & |
---|
| 348 | dflux_s(iloc(1)), dflux_l(iloc(1)) |
---|
| 349 | ENDIF |
---|
| 350 | |
---|
| 351 | !**************************************************************************************** |
---|
| 352 | ! Send and cumulate fields to the coupler |
---|
| 353 | ! |
---|
| 354 | !**************************************************************************************** |
---|
| 355 | |
---|
| 356 | CALL cpl_send_seaice_fields(itime, dtime, knon, knindex, & |
---|
[996] | 357 | pctsrf, lafin, rlon, rlat, & |
---|
[888] | 358 | swnet, lwnet, fluxlat, fluxsens, & |
---|
[2872] | 359 | precip_rain, precip_snow, evap, tsurf_new, fder_new, alb1, flux_u1, flux_v1,& |
---|
| 360 | sens_prec_liq, sens_prec_sol, lat_prec_liq, lat_prec_sol) |
---|
[2538] | 361 | |
---|
[782] | 362 | |
---|
| 363 | |
---|
| 364 | END SUBROUTINE ocean_cpl_ice |
---|
| 365 | ! |
---|
| 366 | !**************************************************************************************** |
---|
| 367 | ! |
---|
| 368 | END MODULE ocean_cpl_mod |
---|