Changeset 12062 for NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfload.F90

Timestamp:

2019-12-05T10:57:46+01:00 (4 years ago)

Author:

mathiot

Message:

changes required by N.J. review

File:

• NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfload.F90 (modified) (5 diffs)

NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfload.F90

@@ -11 +11 @@
    !!----------------------------------------------------------------------
 
-   USE isf_oce, ONLY: cn_isfload                    ! ice shelf variables
+   USE isf_oce, ONLY: cn_isfload, rn_isfload_T, rn_isfload_S ! ice shelf variables
 
-   USE dom_oce, ONLY: e3w_n, gdept_n, risfdep, mikt ! vertical scale factor
-   USE eosbn2 , ONLY: eos                           ! eos routine
+   USE dom_oce, ONLY: e3w_n, gdept_n, risfdep, mikt          ! vertical scale factor
+   USE eosbn2 , ONLY: eos                                    ! eos routine
 
-   USE lib_mpp, ONLY: ctl_stop                      ! ctl_stop routine
-   USE in_out_manager                               ! 
+   USE lib_mpp, ONLY: ctl_stop                               ! ctl_stop routine
+   USE in_out_manager                                        ! 
 
    IMPLICIT NONE
@@ -44 +44 @@
       ! ice shelf cavity
       SELECT CASE ( cn_isfload )
-      CASE ( 'isomip' )
-         CALL isf_load_isomip ( pisfload )
+      CASE ( 'uniform' )
+         CALL isf_load_uniform ( pisfload )
       CASE DEFAULT
          CALL ctl_stop('STOP','method cn_isfload to compute ice shelf load does not exist (isomip), check your namelist')
@@ -52 +52 @@
    END SUBROUTINE isf_load
 
-   SUBROUTINE isf_load_isomip( pisfload )
+   SUBROUTINE isf_load_uniform( pisfload )
       !!--------------------------------------------------------------------
       !!                  ***  SUBROUTINE isf_load  ***
@@ -77 +77 @@
       znad = 1._wp                     !- To use density and not density anomaly
       !
-      !                                !- assume water displaced by the ice shelf is at T=-1.9 and S=34.4 (rude)
-      zts_top(:,:,jp_tem) = -1.9_wp   ;   zts_top(:,:,jp_sal) = 34.4_wp
+      !                                !- assume water displaced by the ice shelf is at T=rn_isfload_T and S=rn_isfload_S (rude)
+      zts_top(:,:,jp_tem) = rn_isfload_T   ;   zts_top(:,:,jp_sal) = rn_isfload_S
       !
       DO jk = 1, jpk                   !- compute density of the water displaced by the ice shelf 
@@ -111 +111 @@
       END DO
       !
-   END SUBROUTINE isf_load_isomip
+   END SUBROUTINE isf_load_uniform
 
 END MODULE isfload