Changeset 5282 for LMDZ6/trunk/libf/phylmd/ocean_cpl_mod.f90
- Timestamp:
- Oct 28, 2024, 1:11:48 PM (6 hours ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/ocean_cpl_mod.f90
r5274 r5282 8 8 ! 9 9 10 IMPLICIT NONE 10 USE clesphys_mod_h 11 IMPLICIT NONE 11 12 PRIVATE 12 13 … … 40 41 ! Initialize module cpl_init 41 42 CALL cpl_init(dtime, rlon, rlat) 42 43 43 44 END SUBROUTINE ocean_cpl_init 44 45 ! … … 60 61 ! 61 62 ! This subroutine treats the "open ocean", all grid points that are not entierly covered 62 ! by ice. The subroutine first receives fields from coupler, then some calculations at 63 ! by ice. The subroutine first receives fields from coupler, then some calculations at 63 64 ! surface is done and finally it sends some fields to the coupler. 64 65 ! … … 82 83 cpl_send_ocean_fields 83 84 use config_ocean_skin_m, only: activate_ocean_skin 84 85 86 INCLUDE "clesphys.h" 87 ! 88 ! Input arguments 85 USE clesphys_mod_h 86 ! Input arguments 89 87 !**************************************************************************************** 90 88 INTEGER, INTENT(IN) :: itime, knon … … 132 130 REAL, DIMENSION(klon), INTENT(INOUT) :: snow 133 131 REAL, DIMENSION(klon), INTENT(INOUT) :: agesno 134 132 135 133 ! Output arguments 136 134 !**************************************************************************************** … … 139 137 REAL, DIMENSION(klon), INTENT(OUT) :: flux_u1, flux_v1 140 138 REAL, DIMENSION(klon), INTENT(OUT) :: tsurf_new 141 REAL, DIMENSION(klon), INTENT(OUT) :: dflux_s, dflux_l 139 REAL, DIMENSION(klon), INTENT(OUT) :: dflux_s, dflux_l 142 140 REAL, intent(out):: sens_prec_liq(:) ! (knon) 143 141 144 142 REAL, INTENT(OUT):: sss(:) ! (klon) 145 143 ! bulk salinity of the surface layer of the ocean, in ppt 146 144 147 145 148 146 ! Local variables … … 156 154 REAL, DIMENSION(klon) :: u1_lay, v1_lay 157 155 LOGICAL :: check=.FALSE. 158 REAL sens_prec_sol(knon) 159 REAL, DIMENSION(klon) :: lat_prec_liq, lat_prec_sol 156 REAL sens_prec_sol(knon) 157 REAL, DIMENSION(klon) :: lat_prec_liq, lat_prec_sol 160 158 161 159 ! End definitions … … 180 178 agesno(:) = 0. 181 179 lat_prec_liq = 0.; lat_prec_sol = 0. 182 180 183 181 184 182 DO i = 1, knon … … 203 201 204 202 ! assertion: tsurf_new == tsurf_cpl 205 203 206 204 do j = 1, knon 207 205 i = knindex(j) … … 213 211 214 212 215 213 216 214 ! - Flux calculation at first modele level for U and V 217 215 CALL calcul_flux_wind(knon, dtime, & … … 219 217 AcoefU, AcoefV, BcoefU, BcoefV, & 220 218 p1lay, temp_air, & 221 flux_u1, flux_v1) 219 flux_u1, flux_v1) 222 220 223 221 !**************************************************************************************** … … 226 224 !**************************************************************************************** 227 225 fder_new(:) = fder_old(:) + dflux_s(:) + dflux_l(:) 228 226 229 227 iloc = MAXLOC(fder_new(1:klon)) 230 228 IF (check .AND. fder_new(iloc(1))> 0.) THEN … … 262 260 tsurf_new, dflux_s, dflux_l, rhoa) 263 261 ! 264 ! This subroutine treats the ocean where there is ice. The subroutine first receives 265 ! fields from coupler, then some calculations at surface is done and finally sends 262 ! This subroutine treats the ocean where there is ice. The subroutine first receives 263 ! fields from coupler, then some calculations at surface is done and finally sends 266 264 ! some fields to the coupler. 267 ! 265 ! 268 266 USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO & 269 267 , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA & … … 285 283 286 284 287 INCLUDE "clesphys.h"288 285 289 286 ! Input arguments
Note: See TracChangeset
for help on using the changeset viewer.