Ignore:
Timestamp:
2019-11-04T19:00:27+01:00 (12 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/isfpar.F90

    r11541 r11852  
    1414   !!   isfpar       : compute ice shelf melt using a prametrisation of ice shelf cavities 
    1515   !!---------------------------------------------------------------------- 
    16    USE oce            ! ocean dynamics and tracers 
     16   !USE oce            ! ocean dynamics and tracers 
    1717   USE isf            ! ice shelf 
    18    USE isfutils       ! 
    19    USE isfparmlt      ! ice shelf parametrisation 
    20    USE isftbl         ! ice shelf depth average 
    21    USE isfdiags       ! ice shelf diagnostics 
    22    USE dom_oce        ! ocean space and time domain 
    23    USE phycst         ! physical constants 
    24    USE eosbn2         ! equation of state 
     18   ! 
     19   USE isfrst   , ONLY: isfrst_write, isfrst_read ! ice shelf restart read/write subroutine 
     20   USE isftbl   , ONLY: isf_tbl_ktop, isf_tbl_lvl ! ice shelf top boundary layer properties subroutine 
     21   USE isfutils , ONLY: debug, read_2dcstdta      ! ice shelf debug subroutine 
     22   USE isfparmlt, ONLY: isfpar_mlt     ! ice shelf melt formulation subroutine 
     23   USE isfdiags , ONLY: isf_diags_flx  ! ice shelf diags subroutine 
     24   ! 
     25   USE dom_oce  , ONLY: jpi,jpj        ! ocean space and time domain 
     26   USE phycst   , ONLY: r1_rau0_rcp    ! physical constants 
    2527   ! 
    2628   USE in_out_manager ! I/O manager 
    2729   USE iom            ! I/O library 
    2830   USE fldread        ! read input field at current time step 
    29    USE lbclnk         ! 
     31   USE lbclnk         ! lbc_lnk 
    3032 
    3133   IMPLICIT NONE 
     
    8082      ptsc(:,:,jp_tem) = zqh(:,:) * r1_rau0_rcp 
    8183      ! 
     84      ! write restart variables (qoceisf, qhcisf, fwfisf for now and before) 
     85      IF (lrst_oce) CALL isfrst_write(kt, 'par', ptsc, pqfwf) 
     86      ! 
     87      IF ( ln_isfdebug ) THEN 
     88         CALL debug('isf_par: ptsc T',ptsc(:,:,1)) 
     89         CALL debug('isf_par: ptsc S',ptsc(:,:,2)) 
     90         CALL debug('isf_par: pqfwf fwf',pqfwf(:,:)) 
     91      END IF 
     92      ! 
    8293   END SUBROUTINE isf_par 
    8394 
     
    100111      rhisf_tbl_par(:,:)  = 1e-20     ; rfrac_tbl_par(:,:)    = 0.0_wp 
    101112      ! 
    102       mskisf_par(:,:) = 0 
    103       ! 
    104113      ! define isf tbl tickness, top and bottom indice 
    105114      CALL read_2dcstdta(TRIM(sn_isfpar_zmax%clname), TRIM(sn_isfpar_zmax%clvar), ztblmax) 
     
    110119      ztblmin(:,:) = ztblmin(:,:) * ssmask(:,:) 
    111120      ! 
    112       ! if param used under an ice shelf overwrite ztblmax by the ice shelf draft 
     121      ! if param used under an ice shelf overwrite ztblmin by the ice shelf draft 
    113122      WHERE ( risfdep > 0._wp .AND. ztblmin > 0._wp ) 
    114123         ztblmin(:,:) = risfdep(:,:) 
     
    127136      END WHERE 
    128137      ! 
    129       ! compute misfkb_par, rhisf_tbl 
    130       rhisf_tbl_par(:,:) = rhisf0_tbl_par(:,:) 
    131       CALL isf_tbl_lvl( ht_n * mskisf_par, e3t_n, misfkt_par, misfkb_par, rhisf_tbl_par, rfrac_tbl_par ) 
     138      ! read par variable from restart 
     139      IF ( ln_rstart ) CALL isfrst_read('par', risf_par_tsc, fwfisf_par, risf_par_tsc_b, fwfisf_par_b) 
    132140      ! 
    133141      SELECT CASE ( TRIM(cn_isfpar_mlt) ) 
Note: See TracChangeset for help on using the changeset viewer.