Changeset 5022 for LMDZ6/trunk/libf/phylmd/surf_ocean_mod.F90
- Timestamp:
- Jul 5, 2024, 4:38:48 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/surf_ocean_mod.F90
r4526 r5022 21 21 tsurf_new, dflux_s, dflux_l, lmt_bils, & 22 22 flux_u1, flux_v1, delta_sst, delta_sal, ds_ns, dt_ns, dter, dser, & 23 dt_ds, tkt, tks, taur, sss) 23 dt_ds, tkt, tks, taur, sss & 24 #ifdef ISO 25 & ,xtprecip_rain, xtprecip_snow,xtspechum,Roce, & 26 & xtsnow,xtevap,h1 & 27 #endif 28 & ) 24 29 25 30 use albedo, only: alboc, alboc_cd … … 31 36 USE ocean_cpl_mod, ONLY : ocean_cpl_noice 32 37 USE indice_sol_mod, ONLY : nbsrf, is_oce 38 #ifdef ISO 39 USE infotrac_phy, ONLY : ntraciso=>ntiso,niso 40 #ifdef ISOVERIF 41 USE isotopes_mod, ONLY: iso_eau,ridicule 42 USE isotopes_verif_mod 43 #endif 44 #endif 33 45 USE limit_read_mod 34 use config_ocean_skin_m, only: activate_ocean_skin46 USE config_ocean_skin_m, ONLY: activate_ocean_skin 35 47 ! 36 48 ! This subroutine will make a call to ocean_XXX_noice according to the ocean mode (force, … … 68 80 REAL, DIMENSION(klon), INTENT(IN) :: rugoro 69 81 REAL, DIMENSION(klon,nbsrf), INTENT(IN) :: pctsrf 82 #ifdef ISO 83 REAL, DIMENSION(ntraciso,klon), INTENT(IN) :: xtprecip_rain, xtprecip_snow 84 REAL, DIMENSION(ntraciso,klon), INTENT(IN) :: xtspechum 85 #endif 70 86 71 87 ! In/Output variables … … 75 91 REAL, DIMENSION(klon), INTENT(INOUT) :: agesno 76 92 REAL, DIMENSION(klon), INTENT(inOUT) :: z0h 93 #ifdef ISO 94 REAL, DIMENSION(niso,klon), INTENT(IN) :: xtsnow 95 REAL, DIMENSION(niso,klon), INTENT(INOUT):: Roce 96 #endif 77 97 78 98 REAL, intent(inout):: delta_sst(:) ! (knon) … … 136 156 ! size klon because of the coupling machinery.) 137 157 158 #ifdef ISO 159 REAL, DIMENSION(ntraciso,klon), INTENT(out) :: xtevap ! isotopes in surface evaporation flux 160 REAL, DIMENSION(klon), INTENT(out) :: h1 ! just a diagnostic, not useful for the simulation 161 #endif 162 138 163 ! Local variables 139 164 !************************************************************************* … … 146 171 REAL, DIMENSION(klon) :: precip_totsnow 147 172 CHARACTER(len=20),PARAMETER :: modname="surf_ocean" 148 realrhoa(knon) ! density of moist air (kg / m3)173 REAL rhoa(knon) ! density of moist air (kg / m3) 149 174 REAL sens_prec_liq(knon) 150 175 151 176 REAL t_int(knon) ! ocean-air interface temperature, in K 152 reals_int(knon) ! ocean-air interface salinity, in ppt177 REAL s_int(knon) ! ocean-air interface salinity, in ppt 153 178 154 179 !************************************************************************** 155 180 181 #ifdef ISO 182 #ifdef ISOVERIF 183 DO i = 1, knon 184 IF (iso_eau > 0) THEN 185 CALL iso_verif_egalite_choix(xtspechum(iso_eau,i), & 186 & spechum(i),'surf_ocean_mod 117', & 187 & errmax,errmaxrel) 188 CALL iso_verif_egalite_choix(xtsnow(iso_eau,i), & 189 & snow(i),'surf_ocean_mod 127', & 190 & errmax,errmaxrel) 191 ENDIF !IF (iso_eau > 0) then 192 ENDDO !DO i=1,klon 193 #endif 194 #endif 156 195 157 196 !****************************************************************************** … … 230 269 radsol, snow, agesno, & 231 270 qsurf, evap, fluxsens, fluxlat, flux_u1, flux_v1, & 232 tsurf_new, dflux_s, dflux_l, sens_prec_liq, rhoa) 271 tsurf_new, dflux_s, dflux_l, sens_prec_liq, rhoa & 272 #ifdef ISO 273 ,xtprecip_rain, xtprecip_snow, xtspechum,Roce,rlat, & 274 xtsnow,xtevap,h1 & 275 #endif 276 ) 233 277 END SELECT 234 278
Note: See TracChangeset
for help on using the changeset viewer.