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

Timestamp:

Oct 10, 2016, 10:57:24 AM (8 years ago)

Author:

Ehouarn Millour

Message:

Making the slab work:

• added a slab_heat_transp_mod module for horizontal diffusion and Ekman transport
• added storage and output of relevent variables in phyredem, phyetat0, phy_output_ctrlout_mod, phys_output_write_mod
• moved nslay (number of slab layers) out of dimphy into ocean_slab_mod.

FC

File:

• LMDZ5/trunk/libf/phylmd/phyredem.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/phyredem.F90

@@ -30 +30 @@
   USE indice_sol_mod, ONLY: nbsrf, is_oce, is_sic, is_ter, is_lic, epsfra
   USE surface_data, ONLY: type_ocean, version_ocean
-  USE ocean_slab_mod, ONLY : tslab, seaice, tice, fsic
+  USE ocean_slab_mod, ONLY : nslay, tslab, seaice, tice, fsic
   USE time_phylmdz_mod, ONLY: annee_ref, day_end, itau_phy, pdtphys
 
@@ -58 +58 @@
 
   INTEGER isoil, nsrf,isw
-  CHARACTER (len=7) :: str7
+  CHARACTER (len=2) :: str2
   CHARACTER (len=256) :: nam, lnam
   INTEGER           :: it, iiq
@@ -303 +303 @@
   ! Restart variables for Slab ocean
   IF (type_ocean == 'slab') THEN
-      CALL put_field("tslab", "Slab ocean temperature", tslab)
+      IF (nslay.EQ.1) THEN
+        CALL put_field("tslab", "Slab ocean temperature", tslab)
+      ELSE
+        DO it=1,nslay
+          WRITE(str2,'(i2.2)') it
+          CALL put_field("tslab"//str2, "Slab ocean temperature", tslab(:,it))
+        END DO
+      END IF
       IF (version_ocean == 'sicINT') THEN
           CALL put_field("seaice", "Slab seaice (kg/m2)", seaice)