Changeset 9023 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/SBC/tideini.F90

Timestamp:

2017-12-13T18:08:50+01:00 (6 years ago)

Author:

timgraham

Message:

Merged METO_MERCATOR branch and resolved all conflicts in OPA_SRC

File:

• branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/SBC/tideini.F90 (modified) (6 diffs)

branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/SBC/tideini.F90

@@ -8 +8 @@
    USE oce            ! ocean dynamics and tracers variables
    USE dom_oce        ! ocean space and time domain
-   USE phycst         ! physical constant
-   USE daymod         ! calandar
+   USE phycst         ! physical constants
+   USE daymod         ! calendar
    USE tide_mod       ! 
    !
@@ -27 +27 @@
    LOGICAL , PUBLIC ::   ln_tide         !:
    LOGICAL , PUBLIC ::   ln_tide_pot     !:
+   LOGICAL , PUBLIC ::   ln_read_load    !:
+   LOGICAL , PUBLIC ::   ln_scal_load    !:
    LOGICAL , PUBLIC ::   ln_tide_ramp    !:
    INTEGER , PUBLIC ::   nb_harmo        !:
    INTEGER , PUBLIC ::   kt_tide         !:
    REAL(wp), PUBLIC ::   rdttideramp     !:
-   
+   REAL(wp), PUBLIC ::   rn_scal_load    !:
+   CHARACTER(lc), PUBLIC ::   cn_tide_load   !: 
+
    INTEGER , PUBLIC, ALLOCATABLE, DIMENSION(:) ::   ntide   !:
 
@@ -49 +53 @@
       INTEGER  ::   ios                 ! Local integer output status for namelist read
       !
-      NAMELIST/nam_tide/ln_tide, ln_tide_pot, ln_tide_ramp, rdttideramp, clname
+      NAMELIST/nam_tide/ln_tide, ln_tide_pot, ln_scal_load, ln_read_load, cn_tide_load, &
+                  &     ln_tide_ramp, rn_scal_load, rdttideramp, clname
       !!----------------------------------------------------------------------
       !
@@ -69 +74 @@
             WRITE(numout,*) '   Namelist nam_tide'
             WRITE(numout,*) '              Use tidal components : ln_tide      = ', ln_tide
-            WRITE(numout,*) '      Apply astronomical potential : ln_tide_pot  = ', ln_tide_pot
-            WRITE(numout,*) '                                     nb_harmo     = ', nb_harmo
-            WRITE(numout,*) '                                     ln_tide_ramp = ', ln_tide_ramp
-            WRITE(numout,*) '                                     rdttideramp  = ', rdttideramp
+            WRITE(numout,*) '      Apply astronomical potential          : ln_tide_pot  = ', ln_tide_pot
+            WRITE(numout,*) '      Use scalar approx. for load potential : ln_scal_load = ', ln_scal_load
+            WRITE(numout,*) '      Read load potential from file         : ln_read_load = ', ln_read_load
+            WRITE(numout,*) '      Apply ramp on tides at startup        : ln_tide_ramp = ', ln_tide_ramp
+            WRITE(numout,*) '      Fraction of SSH used in scal. approx. : rn_scal_load = ', rn_scal_load
+            WRITE(numout,*) '      Duration (days) of ramp               : rdttideramp  = ', rdttideramp
          ENDIF
       ELSE
@@ -93 +100 @@
       IF( nb_harmo == 0 )   CALL ctl_stop( 'tide_init : No tidal components set in nam_tide' )
       !
+      IF( ln_read_load.AND.(.NOT.ln_tide_pot) ) &
+          &   CALL ctl_stop('ln_read_load requires ln_tide_pot')
+      IF( ln_scal_load.AND.(.NOT.ln_tide_pot) ) &
+          &   CALL ctl_stop('ln_scal_load requires ln_tide_pot')
+      IF( ln_scal_load.AND.ln_read_load ) &
+          &   CALL ctl_stop('Choose between ln_scal_load and ln_read_load')
       IF( ln_tide_ramp.AND.((nitend-nit000+1)*rdt/rday < rdttideramp) )   &
          &   CALL ctl_stop('rdttideramp must be lower than run duration')
@@ -112 +125 @@
       kt_tide = nit000
       !
+      IF (.NOT.ln_scal_load ) rn_scal_load = 0._wp
+      !
    END SUBROUTINE tide_init