Changeset 2057 for LMDZ5/trunk/libf/phylmd/phyetat0.F90

Timestamp:

Jun 16, 2014, 4:33:38 PM (10 years ago)

Author:

Ehouarn Millour

Message:

Preparatory stuff to fix and improve the slab ocean model.
FC

File:

• LMDZ5/trunk/libf/phylmd/phyetat0.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/phyetat0.F90

@@ -3 +3 @@
 SUBROUTINE phyetat0 (fichnom, clesphy0, tabcntr0)
 
-  USE dimphy, only: klon, zmasq, klev
+  USE dimphy, only: klon, zmasq, klev, nslay
   USE iophy, ONLY : init_iophy_new
   USE ocean_cpl_mod,    ONLY : ocean_cpl_init
@@ -21 +21 @@
   USE carbon_cycle_mod, ONLY : carbon_cycle_tr, carbon_cycle_cpl, co2_send
   USE indice_sol_mod, only: nbsrf, is_ter, epsfra, is_lic, is_oce, is_sic
+  USE ocean_slab_mod, ONLY: tslab, ocean_slab_init
 
   IMPLICIT none
@@ -40 +41 @@
 
   REAL tsoil(klon, nsoilmx, nbsrf)
-  REAL tslab(klon), seaice(klon)
   REAL qsurf(klon, nbsrf)
   REAL qsol(klon)
@@ -1045 +1045 @@
   end if
 
+  ! Initialize Slab variables
+  IF ( type_ocean == 'slab' ) THEN
+      ALLOCATE(tslab(klon,nslay), stat=ierr)
+      IF (ierr /= 0) CALL abort_gcm &
+         ('phyetat0', 'pb allocation tslab', 1)
+      CALL get_field("tslab", tslab, found)
+      IF (.NOT. found) THEN 
+          PRINT*, "phyetat0: Le champ <tslab> est absent"
+          PRINT*, "Initialisation a tsol_oce"
+          DO i=1,nslay
+              tslab(:,i)=ftsol(:,is_oce)
+          END DO
+      END IF 
+      print*, "calling slab_init"
+      CALL ocean_slab_init(dtime, pctsrf)
+  END IF ! Slab        
+
   ! on ferme le fichier
   CALL close_startphy