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

LMDZ5/trunk/libf/phylmd/limit_slab.F90

@@ -1 +1 @@
 ! $Header$
 
-SUBROUTINE limit_slab(itime, dtime, jour, lmt_bils, lmt_foce, diff_sst)
+SUBROUTINE limit_slab(itime, dtime, jour, lmt_bils, diff_sst)
 
   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, lmt_foce, diff_sst
+  REAL, DIMENSION(klon), INTENT(OUT) :: lmt_bils, diff_sst
 
 ! Locals variables with attribute SAVE
 !****************************************************************************************
-  REAL, DIMENSION(:), ALLOCATABLE, SAVE :: bils_save, foce_save
-!$OMP THREADPRIVATE(bils_save, foce_save)
+  REAL, DIMENSION(:), ALLOCATABLE, SAVE :: bils_save, diff_sst_save
+!$OMP THREADPRIVATE(bils_save, diff_sst_save)
 
 ! Locals variables
@@ -32 +32 @@
   INTEGER                  :: nvarid, nid, ierr, i
   INTEGER, DIMENSION(2)    :: start, epais 
-  REAL, DIMENSION(klon_glo):: bils_glo, foce_glo, sst_l_glo, sst_lp1_glo, diff_sst_glo
+  REAL, DIMENSION(klon_glo):: bils_glo, sst_l_glo, sst_lp1_glo, diff_sst_glo
   CHARACTER (len = 20)     :: modname = 'limit_slab'
+  LOGICAL                  :: read_bils,read_sst
 
 ! End declaration
@@ -41 +42 @@
   lmt_pas = NINT(86400./dtime)
   
+  ! F. Codron 5/14: add defaults for bils, diff_sst (0)
   IF (MOD(itime-1, lmt_pas) == 0) THEN   ! time to read
      !$OMP MASTER  ! Only master thread
      IF (is_mpi_root) THEN ! Only master processus
         print*,'in limit_slab time to read, itime=',itime
+        read_bils=.TRUE.
+        read_sst=.TRUE.
         
         ierr = NF90_OPEN ('limit_slab.nc', NF90_NOWRITE, nid)
-        IF (ierr /= NF90_NOERR) CALL abort_gcm(modname,&
-             'Pb in opening file limit_slab.nc',1)
+        IF (ierr /= NF90_NOERR) THEN
+            read_bils=.FALSE.
+            read_sst=.FALSE.
+        ELSE ! read file
         
         ! La tranche de donnees a lire:
@@ -57 +63 @@
 
 !****************************************************************************************
-! 2) Read bils and ocean fraction
+! 2) Read bils and SST tendency
 !
 !****************************************************************************************
@@ -63 +69 @@
 ! Read bils_glo
         ierr = NF90_INQ_VARID(nid, 'BILS_OCE', nvarid)
-        IF (ierr /= NF90_NOERR) CALL abort_gcm(modname,'The variable <BILS_OCE> is abstent',1)
-
-        ierr = NF90_GET_VAR(nid,nvarid,bils_glo,start,epais)
-        IF (ierr /= NF90_NOERR) CALL abort_gcm(modname,'Reading of <BILS_OCE> failed',1)
-!
-! Read foce_glo
-        ierr = NF90_INQ_VARID(nid, 'FOCE', nvarid)
-        IF (ierr /= NF90_NOERR) CALL abort_gcm(modname,'The variable <FOCE> is abstent',1)
-
-        ierr = NF90_GET_VAR(nid,nvarid,foce_glo,start,epais)
-        IF (ierr /= NF90_NOERR) CALL abort_gcm(modname,'Reading of <FOCE> failed',1)
-!
+        IF (ierr /= NF90_NOERR) THEN
+            read_bils=.FALSE.
+        ELSE
+            ierr = NF90_GET_VAR(nid,nvarid,bils_glo,start,epais)
+            IF (ierr /= NF90_NOERR) read_bils=.FALSE.
+        END IF
 ! Read sst_glo for this day
         ierr = NF90_INQ_VARID(nid, 'SST', nvarid)
-        IF (ierr /= NF90_NOERR) CALL abort_gcm(modname,'The variable <SST> is abstent',1)
-
-        ierr = NF90_GET_VAR(nid,nvarid,sst_l_glo,start,epais)
-        IF (ierr /= NF90_NOERR) CALL abort_gcm(modname,'Reading of <SST> failed',1)
-
+        IF (ierr /= NF90_NOERR)  THEN
+            read_sst=.FALSE.
+        ELSE
+            ierr = NF90_GET_VAR(nid,nvarid,sst_l_glo,start,epais)
+            IF (ierr /= NF90_NOERR) read_sst=.FALSE.
 ! Read sst_glo for one day ahead
-        start(2) = jour + 1
-        IF (start(2) > 360) start(2)=1
-        ierr = NF90_GET_VAR(nid,nvarid,sst_lp1_glo,start,epais)
-        IF (ierr /= NF90_NOERR) CALL abort_gcm(modname,'Reading of <SST> day+1 failed',1)
-
-! 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)
+            start(2) = jour + 1
+            IF (start(2) > 360) start(2)=1
+            ierr = NF90_GET_VAR(nid,nvarid,sst_lp1_glo,start,epais)
+            IF (ierr /= NF90_NOERR) read_sst=.FALSE.
+        END IF
 
 !****************************************************************************************
-! 5) Close file and distribuate variables to all processus
+! 5) Close file and distribute variables to all processus
 !
 !****************************************************************************************
         ierr = NF90_CLOSE(nid)
         IF (ierr /= NF90_NOERR) CALL abort_gcm(modname,'Pb when closing file', 1)
+        END IF ! Read File
+        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
      ENDIF ! is_mpi_root
 
@@ -105 +111 @@
        
      IF (.NOT. ALLOCATED(bils_save)) THEN
-        ALLOCATE(bils_save(klon), foce_save(klon), stat=ierr)
+        ALLOCATE(bils_save(klon), diff_sst_save(klon), stat=ierr)
         IF (ierr /= 0) CALL abort_gcm('limit_slab', 'pb in allocation',1)
      END IF
 
-     CALL Scatter(bils_glo, bils_save)
-     CALL Scatter(foce_glo, foce_save)
-     CALL Scatter(diff_sst_glo, diff_sst)
+! Giveddefault values if needed
+     IF (read_bils) THEN
+         CALL Scatter(bils_glo, bils_save)
+     ELSE
+         bils_save(:)=0.
+     END IF
+     IF (read_sst) THEN
+         CALL Scatter(diff_sst_glo, diff_sst_save)
+     ELSE
+         diff_sst_save(:)=0.
+     END IF
      
-  ELSE ! not time to read
-     diff_sst(:) = 0.
   ENDIF ! time to read
 
   lmt_bils(:) = bils_save(:)
-  lmt_foce(:) = foce_save(:)
+  diff_sst(:) = diff_sst_save(:)
   
 END SUBROUTINE limit_slab