Ignore:
Timestamp:
2019-11-04T19:00:27+01:00 (13 months ago)
Author:
mathiot
Message:

ENHANCE-02_ISF_nemo: fix WED025 restartability, finish removing useless USE, remove useless lbc_lnk

File:
1 edited

Legend:

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

    r11823 r11852  
    1111   !!---------------------------------------------------------------------- 
    1212 
    13    USE isf             ! ice shelf variables 
    14    USE dom_oce         ! vertical scale factor 
     13   USE isf, ONLY: cn_isfload   ! ice shelf variables 
     14 
     15   USE dom_oce, ONLY: e3w_n, gdept_n, risfdep, mikt ! vertical scale factor 
     16   USE eosbn2 , ONLY: eos                           ! eos routine 
     17 
     18   USE lib_mpp, ONLY: ctl_stop ! ctl_stop routine 
    1519   USE in_out_manager  !  
    16    USE eosbn2          ! eos routine 
    17    USE lib_mpp         ! ctl_stop routine 
    1820 
    1921   IMPLICIT NONE 
     
    8688      ! 
    8789      !                                !- Surface value + ice shelf gradient 
    88       risfload(:,:) = 0._wp                       ! compute pressure due to ice shelf load  
     90      pisfload(:,:) = 0._wp                       ! compute pressure due to ice shelf load  
    8991      DO jj = 1, jpj                         ! (used to compute hpgi/j for all the level from 1 to miku/v) 
    9092         DO ji = 1, jpi                      ! divided by 2 later 
     
    9496               ! 
    9597               ! top layer of the ice shelf 
    96                risfload(ji,jj) = risfload(ji,jj) + (znad + zrhd(ji,jj,1) ) * e3w_n(ji,jj,1) 
     98               pisfload(ji,jj) = pisfload(ji,jj) + (znad + zrhd(ji,jj,1) ) * e3w_n(ji,jj,1) 
    9799               ! 
    98100               ! core layers of the ice shelf 
    99101               DO jk = 2, ikt-1 
    100                   risfload(ji,jj) = risfload(ji,jj) + (2._wp * znad + zrhd(ji,jj,jk-1) + zrhd(ji,jj,jk)) * e3w_n(ji,jj,jk) 
     102                  pisfload(ji,jj) = pisfload(ji,jj) + (2._wp * znad + zrhd(ji,jj,jk-1) + zrhd(ji,jj,jk)) * e3w_n(ji,jj,jk) 
    101103               END DO 
    102104               ! 
    103105               ! deepest part of the ice shelf (between deepest T point and ice/ocean interface 
    104                risfload(ji,jj) = risfload(ji,jj) + (2._wp * znad + zrhdtop_isf(ji,jj) + zrhd(ji,jj,ikt-1)) & 
     106               pisfload(ji,jj) = pisfload(ji,jj) + (2._wp * znad + zrhdtop_isf(ji,jj) + zrhd(ji,jj,ikt-1)) & 
    105107                  &                                              * ( risfdep(ji,jj) - gdept_n(ji,jj,ikt-1) ) 
    106108               ! 
Note: See TracChangeset for help on using the changeset viewer.