Changeset 14200 for NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/DOM/domain.F90

Timestamp:

2020-12-17T15:36:44+01:00 (3 years ago)

Author:

mcastril

Message:

Merging r14117 through r14199 into dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final

File:

• NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/DOM/domain.F90 (modified) (7 diffs)

NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/DOM/domain.F90

@@ -28 +28 @@
    USE oce            ! ocean variables
    USE dom_oce        ! domain: ocean
+   USE domtile        ! tiling utilities
 #if defined key_qco
-   USE domqco         ! quasi-eulerian
+   USE domqco         ! quasi-eulerian coord.
+#elif defined key_linssh
+   !                  ! fix in time coord.
 #else
-   USE domvvl         ! variable volume
-#endif
-   USE sshwzv  , ONLY : ssh_init_rst   ! set initial ssh
+   USE domvvl         ! variable volume coord.
+#endif
+#if defined key_agrif
+   USE agrif_oce_interp, ONLY : Agrif_istate_ssh ! ssh interpolated from parent
+#endif
    USE sbc_oce        ! surface boundary condition: ocean
    USE trc_oce        ! shared ocean & passive tracers variab
@@ -39 +44 @@
    USE domhgr         ! domain: set the horizontal mesh
    USE domzgr         ! domain: set the vertical mesh
-   USE domtile
    USE dommsk         ! domain: set the mask system
    USE domwri         ! domain: write the meshmask file
@@ -51 +55 @@
    USE lbclnk         ! ocean lateral boundary condition (or mpp link)
    USE lib_mpp        ! distributed memory computing library
-   USE restart        ! only for lrst_oce
+   USE restart        ! only for lrst_oce and rst_read_ssh
 
    IMPLICIT NONE
@@ -169 +173 @@
          DO_2D( 1, 1, 1, 1 )
             IF( ht_0(ji,jj) < 0._wp .AND. ssmask(ji,jj) == 1._wp ) THEN
-               CALL ctl_stop( 'ssh_init_rst : ht_0 must be positive at potentially wet points' )
+               CALL ctl_stop( 'dom_init : ht_0 must be positive at potentially wet points' )
             ENDIF
          END_2D
@@ -177 +181 @@
       !
       !                                 != ssh initialization
-      IF( .NOT.l_offline .AND. .NOT.l_SAS ) THEN
-         CALL ssh_init_rst( Kbb, Kmm, Kaa )
-      ELSE
+      !
+      IF( l_offline .OR. l_SAS ) THEN        !* No ocean dynamics calculation : set to 0
          ssh(:,:,:) = 0._wp
-      ENDIF
-      !
+#if defined key_agrif
+      ELSEIF( .NOT.Agrif_root() .AND.    &
+         &     ln_init_chfrpar ) THEN        !* Interpolate initial ssh from parent
+         CALL Agrif_istate_ssh( Kbb, Kmm, Kaa )
+#endif
+      ELSE                                   !* Read in restart file or set by user
+         CALL rst_read_ssh( Kbb, Kmm, Kaa )
+      ENDIF
+      !     
 #if defined key_qco
       !                                 != Quasi-Euerian coordinate case
       !
       IF( .NOT.l_offline )   CALL dom_qco_init( Kbb, Kmm, Kaa )
+#elif defined key_linssh
+      !                                 != Fix in time : key_linssh case, set through domzgr_substitute.h90
 #else
       !
@@ -345 +357 @@
       IF(lwm) WRITE( numond, namdom )
       !
+#if defined key_linssh
+      ln_linssh = lk_linssh      ! overwrite ln_linssh with the logical associated with key_linssh
+#endif
+      !
 #if defined key_agrif
       IF( .NOT. Agrif_Root() ) THEN    ! AGRIF child, subdivide the Parent timestep
@@ -371 +387 @@
       !
 #if defined key_qco
-      IF( ln_linssh )   CALL ctl_stop( 'STOP','domain: key_qco and ln_linssh = T are incompatible' )
+      IF( ln_linssh )   CALL ctl_stop( 'STOP','domain: key_qco and ln_linssh=T or key_linssh are incompatible' )
 #endif
       !