Changeset 2209 for LMDZ5/trunk/libf/phylmd/limit_slab.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/limit_slab.F90 (modified) (11 diffs)

LMDZ5/trunk/libf/phylmd/limit_slab.F90

@@ -1 +1 @@
 ! $Header$
 
-SUBROUTINE limit_slab(itime, dtime, jour, lmt_bils, diff_sst)
+SUBROUTINE limit_slab(itime, dtime, jour, lmt_bils, diff_sst, diff_siv)
 
   USE dimphy
@@ -20 +20 @@
   INTEGER, INTENT(IN) :: jour    ! jour a lire dans l'annee
   REAL   , INTENT(IN) :: dtime   ! pas de temps de la physique (en s)
-  REAL, DIMENSION(klon), INTENT(OUT) :: lmt_bils, diff_sst
+  REAL, DIMENSION(klon), INTENT(OUT) :: lmt_bils, diff_sst, diff_siv
 
 ! Locals variables with attribute SAVE
 !****************************************************************************************
   REAL, DIMENSION(:), ALLOCATABLE, SAVE :: bils_save, diff_sst_save
-!$OMP THREADPRIVATE(bils_save, diff_sst_save)
+  REAL, DIMENSION(:), ALLOCATABLE, SAVE :: diff_siv_save
+!$OMP THREADPRIVATE(bils_save, diff_sst_save, diff_siv_save)
 
 ! Locals variables
@@ -33 +34 @@
   INTEGER, DIMENSION(2)    :: start, epais 
   REAL, DIMENSION(klon_glo):: bils_glo, sst_l_glo, sst_lp1_glo, diff_sst_glo
+  REAL, DIMENSION(klon_glo):: siv_l_glo, siv_lp1_glo, diff_siv_glo
   CHARACTER (len = 20)     :: modname = 'limit_slab'
-  LOGICAL                  :: read_bils,read_sst
+  LOGICAL                  :: read_bils,read_sst,read_siv
 
 ! End declaration
@@ -49 +51 @@
         read_bils=.TRUE.
         read_sst=.TRUE.
+        read_siv=.TRUE.
         
         ierr = NF90_OPEN ('limit_slab.nc', NF90_NOWRITE, nid)
@@ -54 +57 @@
             read_bils=.FALSE.
             read_sst=.FALSE.
+            read_siv=.FALSE.
         ELSE ! read file
         
@@ -63 +67 @@
 
 !****************************************************************************************
-! 2) Read bils and SST tendency
+! 2) Read bils and SST/ ice volume tendency
 !
 !****************************************************************************************
@@ -89 +93 @@
         END IF
 
+! Read siv_glo for this day
+        ierr = NF90_INQ_VARID(nid, 'SICV', nvarid)
+        IF (ierr /= NF90_NOERR)  THEN
+            read_siv=.FALSE.
+        ELSE
+            start(2) = jour
+            ierr = NF90_GET_VAR(nid,nvarid,siv_l_glo,start,epais)
+            IF (ierr /= NF90_NOERR) read_siv=.FALSE.
+! Read siv_glo for one day ahead
+            start(2) = jour + 1
+            IF (start(2) > 360) start(2)=1
+            ierr = NF90_GET_VAR(nid,nvarid,siv_lp1_glo,start,epais)
+            IF (ierr /= NF90_NOERR) read_siv=.FALSE.
+        END IF
+
 !****************************************************************************************
 ! 5) Close file and distribute variables to all processus
@@ -98 +117 @@
         IF (read_sst) THEN
 ! Calculate difference in temperature between this day and one ahead
-!            DO i=1, klon_glo-1
-!               diff_sst_glo(i) = sst_lp1_glo(i) - sst_l_glo(i)
-!            END DO
-!            diff_sst_glo(klon_glo) = sst_lp1_glo(klon_glo) - sst_l_glo(1)
             DO i=1, klon_glo
                diff_sst_glo(i) = sst_lp1_glo(i) - sst_l_glo(i)
             END DO
         END IF !read_sst
+        IF (read_siv) THEN
+! Calculate difference in temperature between this day and one ahead
+            DO i=1, klon_glo
+               diff_siv_glo(i) = siv_lp1_glo(i) - siv_l_glo(i)
+            END DO
+        END IF !read_siv
      ENDIF ! is_mpi_root
 
@@ -111 +132 @@
        
      IF (.NOT. ALLOCATED(bils_save)) THEN
-        ALLOCATE(bils_save(klon), diff_sst_save(klon), stat=ierr)
+        ALLOCATE(bils_save(klon), diff_sst_save(klon), diff_siv_save(klon), stat=ierr)
         IF (ierr /= 0) CALL abort_gcm('limit_slab', 'pb in allocation',1)
      END IF
 
-! Giveddefault values if needed
+! Give default values if needed
      IF (read_bils) THEN
          CALL Scatter(bils_glo, bils_save)
@@ -126 +147 @@
          diff_sst_save(:)=0.
      END IF
+     IF (read_siv) THEN
+         CALL Scatter(diff_siv_glo, diff_siv_save)
+     ELSE
+         diff_siv_save(:)=0.
+     END IF
      
   ENDIF ! time to read
@@ -131 +157 @@
   lmt_bils(:) = bils_save(:)
   diff_sst(:) = diff_sst_save(:)
+  diff_siv(:) = diff_siv_save(:)
+
   
 END SUBROUTINE limit_slab