New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 11823 for NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfstp.F90 – NEMO

Ignore:
Timestamp:
2019-10-29T12:00:04+01:00 (4 years 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/ISF/isfstp.F90

    r11553 r11823  
    1515   USE oce            ! ocean dynamics and tracers 
    1616   USE dom_oce        ! ocean space and time domain 
     17   USE domvvl, ONLY : ln_vvl_zstar 
    1718   USE phycst         ! physical constants 
    1819   USE eosbn2         ! equation of state 
     
    120121      !! ** Purpose :   Initialisation of the ice shelf public variables 
    121122      !! 
    122       !! ** Method  :   Read the namsbc namelist and set derived parameters 
    123       !!                Call init routines for all other SBC modules that have one 
    124       !! 
    125       !! ** Action  : - read namsbc parameters 
     123      !! ** Method  :   Read the namisf namelist, check option compatibility and set derived parameters 
     124      !! 
     125      !! ** Action  : - read namisf parameters 
    126126      !!              - allocate memory 
    127       !!              - call cav/param init routine 
     127      !!              - output print 
     128      !!              - ckeck option compatibility 
     129      !!              - call cav/param/isfcpl init routine 
    128130      !!---------------------------------------------------------------------- 
    129131      ! 
     
    144146      ! terminate routine now if no ice shelf melt formulation specify 
    145147      IF ( ln_isf ) THEN 
    146          ! 
    147          !--------------------------------------------------------------------------------------------------------------------- 
    148          ! initialisation melt in the cavity 
    149          IF ( ln_isfcav_mlt ) THEN 
    150             ! 
    151             ! initialisation  of cav variable 
    152             CALL isf_cav_init() 
    153             ! 
    154             ! read cav variable from restart 
    155             IF ( ln_rstart ) CALL isfrst_read('cav', risf_cav_tsc, fwfisf_cav, risf_cav_tsc_b, fwfisf_cav_b) 
    156             ! 
    157          END IF 
    158          ! 
    159          !--------------------------------------------------------------------------------------------------------------------- 
    160          ! initialisation parametrised melt 
    161          IF ( ln_isfpar_mlt ) THEN 
    162             ! 
    163             ! initialisation  of par variable 
    164             CALL isf_par_init() 
    165             ! 
    166             ! read par variable from restart 
    167             IF ( ln_rstart ) CALL isfrst_read('par', risf_par_tsc, fwfisf_par, risf_par_tsc_b, fwfisf_par_b) 
    168             ! 
    169          END IF 
    170148         ! 
    171149         !--------------------------------------------------------------------------------------------------------------------- 
     
    173151         IF( ln_isfcpl ) THEN  
    174152            ! 
     153            ! start on an euler time step 
     154            neuler = 0 
     155            ! 
     156            ! ice sheet coupling: extrapolation of restart to fill new wet cell and compute divergence correction 
     157            CALL isfcpl_init() 
     158            ! 
     159            ! all before fields set to now values 
     160            tsb  (:,:,:,:) = tsn  (:,:,:,:) 
     161            ub   (:,:,:)   = un   (:,:,:) 
     162            vb   (:,:,:)   = vn   (:,:,:) 
     163            sshb (:,:)     = sshn (:,:) 
     164            e3t_b(:,:,:)   = e3t_n(:,:,:) 
     165  
    175166            ! prepare writing restart 
    176167            IF( lwxios ) CALL iom_set_rstw_var_active('ssmask') 
     
    182173         END IF 
    183174         ! 
     175         !--------------------------------------------------------------------------------------------------------------------- 
     176         ! initialisation melt in the cavity 
     177         IF ( ln_isfcav_mlt ) THEN 
     178            ! 
     179            ! initialisation  of cav variable 
     180            CALL isf_cav_init() 
     181            ! 
     182            ! read cav variable from restart 
     183            IF ( ln_rstart ) CALL isfrst_read('cav', risf_cav_tsc, fwfisf_cav, risf_cav_tsc_b, fwfisf_cav_b) 
     184            ! 
     185         END IF 
     186         ! 
     187         !--------------------------------------------------------------------------------------------------------------------- 
     188         ! initialisation parametrised melt 
     189         IF ( ln_isfpar_mlt ) THEN 
     190            ! 
     191            ! initialisation  of par variable 
     192            CALL isf_par_init() 
     193            ! 
     194            ! read par variable from restart 
     195            IF ( ln_rstart ) CALL isfrst_read('par', risf_par_tsc, fwfisf_par, risf_par_tsc_b, fwfisf_par_b) 
     196            ! 
     197         END IF 
     198         ! 
    184199      END IF 
    185200          
    186201  END SUBROUTINE isf_init 
    187202 
    188   SUBROUTINE isf_ctl()! 
     203  SUBROUTINE isf_ctl() 
     204      !!--------------------------------------------------------------------- 
     205      !!                  ***  ROUTINE isf_ctl  *** 
     206      !! 
     207      !! ** Purpose :   output print and option compatibility check 
     208      !! 
     209      !!---------------------------------------------------------------------- 
    189210      IF (lwp) THEN 
    190211         WRITE(numout,*) 
     
    272293         END IF 
    273294         ! 
     295         ! compatibility ice shelf and vvl 
     296         IF( .NOT. ln_vvl_zstar .AND. ln_isf ) CALL ctl_stop( 'Only vvl_zstar has been tested with ice shelf cavity' ) 
     297         ! 
    274298      END IF 
    275299   END SUBROUTINE isf_ctl 
Note: See TracChangeset for help on using the changeset viewer.