Changeset 2209 for LMDZ5/trunk/libf/phylmd/surf_seaice_mod.F90

Timestamp:

Feb 16, 2015, 8:34:28 AM (10 years ago)

Author:

Ehouarn Millour

Message:

Update of the slab ocean by Francis Codron. There are now 3 possibilities for the "version_ocean" slab type:
sicOBS = prescribed ice fraction. Water temperature nearby is set to -1.8°C and cannot become lower.
sicNO = ignore sea ice. One can prescribe a fraction, but the nearby ocean evolves freely, depending on surface fluxes: temperature can go below freezing point or above...
sicINT = interactive sea ice.
EM

File:

• LMDZ5/trunk/libf/phylmd/surf_seaice_mod.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/surf_seaice_mod.F90

@@ -25 +25 @@
   USE ocean_forced_mod, ONLY : ocean_forced_ice
   USE ocean_cpl_mod, ONLY    : ocean_cpl_ice
+  USE ocean_slab_mod, ONLY   : ocean_slab_ice
   USE indice_sol_mod
 
@@ -108 +109 @@
             tsurf_new, dflux_s, dflux_l)
        
-!    ELSE IF (type_ocean == 'slab'.AND.version_ocean=='sicINT') THEN
-!       CALL ocean_slab_ice( & 
-!          itime, dtime, jour, knon, knindex, &
-!          debut, &
-!          tsurf, p1lay, cdragh, precip_rain, precip_snow, temp_air, spechum,&
-!          AcoefH, AcoefQ, BcoefH, BcoefQ, &
-!          ps, u1, v1, &
-!          radsol, snow, qsurf, qsol, agesno, tsoil, &
-!          alb1_new, alb2_new, evap, fluxsens, fluxlat, &
-!          tsurf_new, dflux_s, dflux_l)
-!
+    ELSE IF (type_ocean == 'slab'.AND.version_ocean=='sicINT') THEN
+       CALL ocean_slab_ice( & 
+          itime, dtime, jour, knon, knindex, &
+          tsurf, p1lay, cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum,&
+          AcoefH, AcoefQ, BcoefH, BcoefQ, &
+            AcoefU, AcoefV, BcoefU, BcoefV, &
+          ps, u1, v1, &
+          radsol, snow, qsurf, qsol, agesno, &
+          alb1_new, alb2_new, evap, fluxsens, fluxlat, flux_u1, flux_v1, &
+          tsurf_new, dflux_s, dflux_l, swnet)
+
       ELSE ! type_ocean=force or slab +sicOBS or sicNO
        CALL ocean_forced_ice( &