Changeset 3458 for LMDZ6/branches/Ocean_skin/libf/phylmd/surf_ocean_mod.F90

Timestamp:

Jan 30, 2019, 5:55:48 PM (5 years ago)

Author:

lguez

Message:

Introduce variable activate_ocean_skin in module config_ocean_skin_m.

Bug fix in phys_state_var_end: we need to deallocate variables for
lmdz1d (although it is useless for a 3D run).

File:

• LMDZ6/branches/Ocean_skin/libf/phylmd/surf_ocean_mod.F90 (modified) (3 diffs)

LMDZ6/branches/Ocean_skin/libf/phylmd/surf_ocean_mod.F90

@@ -33 +33 @@
     USE limit_read_mod
     use sens_heat_rain_m, only: sens_heat_rain
+    use config_ocean_skin_m, only: activate_ocean_skin
     !
     ! This subroutine will make a call to ocean_XXX_noice according to the ocean mode (force, 
@@ -103 +104 @@
     REAL, intent(out):: t_int(:) ! (knon) interface temperature, in K
     real, intent(out):: s_int(:) ! (knon) interface salinity, in ppt
-    
+
     REAL, intent(out):: dter(:) ! (knon)
     ! Temperature variation in the diffusive microlayer, that is
@@ -311 +312 @@
     ENDIF
 
-    rhoa = PS(:KNON) / (Rd * temp_air(:knon) * (1. + retv * spechum(:knon)))
-    xlv = rlvtt
-    precip_tot = precip_rain(:knon) + precip_snow(:knon)
-    rf =  sens_heat_rain(precip_tot, temp_air(:knon), spechum(:knon), rhoa, &
-         xlv, tsurf_new(:knon), ps(:knon))
-    s1 = 35.
-    call bulk_flux(tkt, tks, taur, dter, dser, t_int, s_int, ds_ns, dt_ns, &
-         u = windsp(:knon), t_ocean_1 = tsurf_new(:knon), s1 = s1, &
-         rain = precip_tot, hf = - fluxsens(:knon), hlb = - fluxlat(:knon), &
-         rnl = - lwnet(:knon), &
-         tau = sqrt(flux_u1(:knon)**2 + flux_v1(:knon)**2), rhoa = rhoa, &
-         xlv = xlv, rf = rf, dtime = dtime, rns = swnet(:knon))
+    if (activate_ocean_skin >= 1) then
+       rhoa = PS(:KNON) / (Rd * temp_air(:knon) * (1. + retv * spechum(:knon)))
+       xlv = rlvtt
+       precip_tot = precip_rain(:knon) + precip_snow(:knon)
+       rf =  sens_heat_rain(precip_tot, temp_air(:knon), spechum(:knon), rhoa, &
+            xlv, tsurf_new(:knon), ps(:knon))
+       s1 = 35.
+       call bulk_flux(tkt, tks, taur, dter, dser, t_int, s_int, ds_ns, dt_ns, &
+            u = windsp(:knon), t_ocean_1 = tsurf_new(:knon), s1 = s1, &
+            rain = precip_tot, hf = - fluxsens(:knon), hlb = - fluxlat(:knon), &
+            rnl = - lwnet(:knon), &
+            tau = sqrt(flux_u1(:knon)**2 + flux_v1(:knon)**2), rhoa = rhoa, &
+            xlv = xlv, rf = rf, dtime = dtime, rns = swnet(:knon))
+    end if
+    
+    if (activate_ocean_skin == 2) tsurf_new(:knon) = t_int
 
   END SUBROUTINE surf_ocean
-  !******************************************************************************
+  !****************************************************************************
   !
 END MODULE surf_ocean_mod