Ignore:
Timestamp:
2019-08-08T16:02:49+02:00 (15 months ago)
Author:
mathiot
Message:

ENHANCE-02_ISF_nemo : add UKESM ice sheet coupling method (ticket #2142)

File:
1 edited

Legend:

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

    r11403 r11423  
    9090      ! 
    9191      !                                !- Surface value + ice shelf gradient 
    92       risfload = 0._wp                       ! compute pressure due to ice shelf load  
     92      risfload(:,:) = 0._wp                       ! compute pressure due to ice shelf load  
    9393      DO jj = 1, jpj                         ! (used to compute hpgi/j for all the level from 1 to miku/v) 
    9494         DO ji = 1, jpi                      ! divided by 2 later 
    9595            ikt = mikt(ji,jj) 
    9696            ! 
    97             ! top layer of the ice shelf 
    98             risfload(ji,jj) = risfload(ji,jj) + (znad + zrhd(ji,jj,1) ) * e3w_n(ji,jj,1) 
    99             ! 
    100             ! core layers of the ice shelf 
    101             DO jk = 2, ikt-1 
    102                risfload(ji,jj) = risfload(ji,jj) + (2._wp * znad + zrhd(ji,jj,jk-1) + zrhd(ji,jj,jk)) * e3w_n(ji,jj,jk) 
    103             END DO 
    104             ! 
    105             ! deepest part of the ice shelf (between deepest T point and ice/ocean interface 
    106             risfload(ji,jj) = risfload(ji,jj) + (2._wp * znad + zrhdtop_isf(ji,jj) + zrhd(ji,jj,ikt-1)) & 
    107                &                                              * ( risfdep(ji,jj) - gdept_n(ji,jj,ikt-1) ) 
    108             ! 
     97            IF ( ikt > 1 ) THEN 
     98               ! 
     99               ! top layer of the ice shelf 
     100               risfload(ji,jj) = risfload(ji,jj) + (znad + zrhd(ji,jj,1) ) * e3w_n(ji,jj,1) 
     101               ! 
     102               ! core layers of the ice shelf 
     103               DO jk = 2, ikt-1 
     104                  risfload(ji,jj) = risfload(ji,jj) + (2._wp * znad + zrhd(ji,jj,jk-1) + zrhd(ji,jj,jk)) * e3w_n(ji,jj,jk) 
     105               END DO 
     106               ! 
     107               ! deepest part of the ice shelf (between deepest T point and ice/ocean interface 
     108               risfload(ji,jj) = risfload(ji,jj) + (2._wp * znad + zrhdtop_isf(ji,jj) + zrhd(ji,jj,ikt-1)) & 
     109                  &                                              * ( risfdep(ji,jj) - gdept_n(ji,jj,ikt-1) ) 
     110               ! 
     111            END IF 
    109112         END DO 
    110113      END DO 
    111       ! 
    112       ! mask it  
    113       risfload(:,:) = risfload(:,:) * mskisf_cav(:,:)  ! need to be saved for diaar5 
    114114      ! 
    115115   END SUBROUTINE isf_load_isomip 
Note: See TracChangeset for help on using the changeset viewer.