Changeset 14053 for NEMO/trunk/src/OCE/oce.F90

Timestamp:

2020-12-03T14:48:38+01:00 (3 years ago)

Author:

techene

Message:

#2385 added to the trunk

File:

• NEMO/trunk/src/OCE/oce.F90 (modified) (3 diffs)

NEMO/trunk/src/OCE/oce.F90

@@ -16 +16 @@
    PRIVATE
 
-   PUBLIC oce_alloc ! routine called by nemo_init in nemogcm.F90
+   PUBLIC oce_alloc       ! routine called by nemo_init in     nemogcm.F90
+   PUBLIC oce_SWE_alloc   ! routine called by nemo_init in SWE/nemogcm.F90 (Shallow Water Eq. case)
 
    !! dynamics and tracer fields
@@ -68 +69 @@
    INTEGER, PUBLIC, DIMENSION(2) :: noce_array                             !: unused array but seems to be needed to prevent agrif from creating an empty module
 
+   !! Shallow Water Eq. case (SWE)
+   LOGICAL, PUBLIC ::   lk_SWE = .FALSE.                                   !: shallow water flag =T in SWE configurations only
+
+   !! Stand-Alone Surface module (SAS)
+   LOGICAL, PUBLIC ::   l_SAS = .FALSE.                                    !: SAS flag =T in SAS configurations only
+   
+   
    !!----------------------------------------------------------------------
    !! NEMO/OCE 4.0 , NEMO Consortium (2018)
@@ -115 +123 @@
    END FUNCTION oce_alloc
 
+
+   INTEGER FUNCTION oce_SWE_alloc()
+      !!----------------------------------------------------------------------
+      !!                   ***  FUNCTION oce_SWE_alloc  ***
+      !!----------------------------------------------------------------------
+      INTEGER :: ierr(2)
+      !!----------------------------------------------------------------------
+      !
+      lk_SWE  = .TRUE.                   ! =T SWE case 
+      !
+      ierr(:) = 0 
+      ALLOCATE( uu(jpi,jpj,jpk,jpt) , vv  (jpi,jpj,jpk,jpt) ,     &          
+         &      ww(jpi,jpj,jpk)     , hdiv(jpi,jpj,jpk)     , ssh(jpi,jpj,jpt) , STAT=ierr(1) )
+         !
+      ALLOCATE(   ts(jpi,jpj,jpk,jpts,jpt) , fraqsr_1lev(jpi,jpj) ,  &
+         &      uu_b(jpi,jpj,jpt) , vv_b(jpi,jpj,jpt)       , rn2(jpi,jpj,jpk) , STAT=ierr(2) )
+         !
+      oce_SWE_alloc = MAXVAL( ierr )
+      IF( oce_SWE_alloc /= 0 )   CALL ctl_stop( 'STOP', 'oce_SWE_alloc: failed to allocate arrays' )
+      !
+   END FUNCTION oce_SWE_alloc
+
    !!======================================================================
 END MODULE oce