Ignore:
Timestamp:
2019-10-29T12:00:04+01:00 (12 months ago)
Author:
mathiot
Message:

rm useless USE statement, option compatibility test + minor changes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/DOM/domvvl.F90

    r11553 r11823  
    1919   !!---------------------------------------------------------------------- 
    2020   USE oce             ! ocean dynamics and tracers 
    21    USE isfcpl 
    22    USE isf 
    2321   USE phycst          ! physical constant 
    2422   USE dom_oce         ! ocean space and time domain 
     
    3836 
    3937   PUBLIC  dom_vvl_init       ! called by domain.F90 
     38   PUBLIC  dom_vvl_zgr        ! called by isfcpl.F90 
    4039   PUBLIC  dom_vvl_sf_nxt     ! called by step.F90 
    4140   PUBLIC  dom_vvl_sf_swp     ! called by step.F90 
     
    117116      !! Reference  : Leclair, M., and G. Madec, 2011, Ocean Modelling. 
    118117      !!---------------------------------------------------------------------- 
     118      ! 
     119      IF(lwp) WRITE(numout,*) 
     120      IF(lwp) WRITE(numout,*) 'dom_vvl_init : Variable volume activated' 
     121      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
     122      ! 
     123      CALL dom_vvl_ctl     ! choose vertical coordinate (z_star, z_tilde or layer) 
     124      ! 
     125      !                    ! Allocate module arrays 
     126      IF( dom_vvl_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'dom_vvl_init : unable to allocate arrays' ) 
     127      ! 
     128      !                    ! Read or initialize e3t_(b/n), tilde_e3t_(b/n) and hdiv_lf 
     129      CALL dom_vvl_rst( nit000, 'READ' ) 
     130      e3t_a(:,:,jpk) = e3t_0(:,:,jpk)  ! last level always inside the sea floor set one for all 
     131      ! 
     132      CALL dom_vvl_zgr() ! interpolation scale factor, depth and water column 
     133      ! 
     134   END SUBROUTINE dom_vvl_init 
     135   ! 
     136   SUBROUTINE dom_vvl_zgr() 
     137      !!---------------------------------------------------------------------- 
     138      !!                ***  ROUTINE dom_vvl_init  *** 
     139      !!                    
     140      !! ** Purpose :  Interpolation of all scale factors,  
     141      !!               depths and water column heights 
     142      !! 
     143      !! ** Method  :  - interpolate scale factors 
     144      !! 
     145      !! ** Action  : - e3t_(n/b) and tilde_e3t_(n/b) 
     146      !!              - Regrid: e3(u/v)_n 
     147      !!                        e3(u/v)_b        
     148      !!                        e3w_n            
     149      !!                        e3(u/v)w_b       
     150      !!                        e3(u/v)w_n       
     151      !!                        gdept_n, gdepw_n and gde3w_n 
     152      !!              - h(t/u/v)_0 
     153      !!              - frq_rst_e3t and frq_rst_hdv 
     154      !! 
     155      !! Reference  : Leclair, M., and G. Madec, 2011, Ocean Modelling. 
     156      !!---------------------------------------------------------------------- 
    119157      INTEGER ::   ji, jj, jk 
    120158      INTEGER ::   ii0, ii1, ij0, ij1 
    121159      REAL(wp)::   zcoef 
    122160      !!---------------------------------------------------------------------- 
    123       ! 
    124       IF(lwp) WRITE(numout,*) 
    125       IF(lwp) WRITE(numout,*) 'dom_vvl_init : Variable volume activated' 
    126       IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
    127       ! 
    128       CALL dom_vvl_ctl     ! choose vertical coordinate (z_star, z_tilde or layer) 
    129       ! 
    130       !                    ! Allocate module arrays 
    131       IF( dom_vvl_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'dom_vvl_init : unable to allocate arrays' ) 
    132       ! 
    133       !                    ! Read or initialize e3t_(b/n), tilde_e3t_(b/n) and hdiv_lf 
    134       CALL dom_vvl_rst( nit000, 'READ' ) 
    135       e3t_a(:,:,jpk) = e3t_0(:,:,jpk)  ! last level always inside the sea floor set one for all 
    136161      ! 
    137162      !                    !== Set of all other vertical scale factors  ==!  (now and before) 
     
    265290      ENDIF 
    266291      ! 
    267    END SUBROUTINE dom_vvl_init 
    268  
     292   END SUBROUTINE dom_vvl_zgr 
    269293 
    270294   SUBROUTINE dom_vvl_sf_nxt( kt, kcall )  
     
    816840            id5 = iom_varid( numror, 'hdiv_lf', ldstop = .FALSE. ) 
    817841            ! 
    818             ! coupling with an ice sheet model (grounding/calving line migration) 
    819             IF ( ln_isfcpl ) THEN 
    820                ! 
    821                ! extrapolate ssh over new water column if needed 
    822                CALL isfcpl_ssh() 
    823                ! 
    824                ! e3t need to be computed where ice shelf draft changed (new ht_0) 
    825                id1 = 0 ; id2 = 0 ; id3 = 0 ; id4 = 0 ; id5 = 0 
    826                ! 
    827             END IF 
    828842            !                             ! --------- ! 
    829843            !                             ! all cases ! 
    830844            !                             ! --------- ! 
     845            ! 
    831846            IF( MIN( id1, id2 ) > 0 ) THEN       ! all required arrays exist 
    832847               CALL iom_get( numror, jpdom_autoglo, 'e3t_b', e3t_b(:,:,:), ldxios = lrxios ) 
     
    10461061      ! 
    10471062      IF( ioptio /= 1 )   CALL ctl_stop( 'Choose ONE vertical coordinate in namelist nam_vvl' ) 
    1048       IF( .NOT. ln_vvl_zstar .AND. ln_isf ) CALL ctl_stop( 'Only vvl_zstar has been tested with ice shelf cavity' ) 
    10491063      ! 
    10501064      IF(lwp) THEN                   ! Print the choice 
Note: See TracChangeset for help on using the changeset viewer.