- Timestamp:
- Jul 24, 2024, 2:54:37 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/Ocean_skin/bulk_flux_m.F90
r5113 r5116 8 8 u, t_ocean_1, s1, rain, hf, hlb, rnl, tau, rhoa, xlv, rf, dtime, rns) 9 9 10 use config_ocean_skin_m, only: jwarm, jcool, rain_effect11 use Microlayer_m, only: Microlayer12 use mom_flux_rain_m, only: mom_flux_rain13 use Near_Surface_m, only: Near_Surface, depth14 use therm_expans_m, only: therm_expans10 use config_ocean_skin_m, ONLY: jwarm, jcool, rain_effect 11 use Microlayer_m, ONLY: Microlayer 12 use mom_flux_rain_m, ONLY: mom_flux_rain 13 use Near_Surface_m, ONLY: Near_Surface, depth 14 use therm_expans_m, ONLY: therm_expans 15 15 16 16 real, intent(out):: tkt(:) … … 97 97 !------------------------------------------------------------------- 98 98 99 if (rain_effect) then99 if (rain_effect) THEN 100 100 taur = mom_flux_rain(u, rain) 101 101 else … … 106 106 if (jwarm .or. jcool) tau_with_min = tau + tau_0 * (1. - exp(- tau_0 / tau)) 107 107 108 if (Jwarm) then109 if (rain_effect) then108 if (Jwarm) THEN 109 if (rain_effect) THEN 110 110 CALL Near_Surface(al, t_subskin, s_subskin, ds_ns, dt_ns, & 111 111 tau_with_min, taur, hlb, rhoa, xlv, dtime, t_ocean_1, s1, rain, & … … 122 122 end if 123 123 124 if (Jcool) then124 if (Jcool) THEN 125 125 ! First guess: 126 126 tkt = 0.001 … … 131 131 dels = rns * (0.065 + 11. * tkt - 6.6e-5 / tkt & 132 132 * (1. - exp(- tkt / 8e-4))) ! equation 16 Ohlmann 133 if (rain_effect) then133 if (rain_effect) THEN 134 134 CALL Microlayer(dter, dser, tkt, tks, hlb, tau_with_min, & 135 135 s_subskin, al, xlv, taur, rf, rain, &
Note: See TracChangeset
for help on using the changeset viewer.