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

Timestamp:

Feb 16, 2015, 8:34:28 AM (9 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/phyredem.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/phyredem.F90

@@ -16 +16 @@
   USE indice_sol_mod
   USE surface_data
-  USE ocean_slab_mod, ONLY : tslab
+  USE ocean_slab_mod, ONLY : tslab, seaice, tice, fsic
 
   IMPLICIT none
@@ -106 +106 @@
 
   ! BP ajout des fraction de chaque sous-surface
+
+  ! Get last fractions from slab ocean
+  IF (type_ocean == 'slab' .AND. version_ocean == "sicINT") THEN
+      WHERE (1.-zmasq(:).GT.EPSFRA)
+          pctsrf(:,is_oce)=(1.-fsic(:))*(1.-zmasq(:))
+          pctsrf(:,is_sic)=fsic(:)*(1.-zmasq(:))
+      END WHERE
+  END IF
 
   ! 1. fraction de terre 
@@ -352 +360 @@
   IF (type_ocean == 'slab') THEN
       CALL put_field("tslab", "Slab ocean temperature", tslab)
+      IF (version_ocean == 'sicINT') THEN
+          CALL put_field("seaice", "Slab seaice (kg/m2)", seaice)
+          CALL put_field("slab_tice", "Slab sea ice temperature", tice)
+      END IF
   END IF