Ignore:
Timestamp:
2019-12-05T10:57:46+01:00 (12 months ago)
Author:
mathiot
Message:

changes required by N.J. review

Location:
NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/ISOMIP+
Files:
1 added
1 deleted
3 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/ISOMIP+/EXPREF/namelist_cfg

    r11896 r12062  
    165165   ! ---------------- ice shelf load ------------------------------- 
    166166   ! 
    167    cn_isfload = 'isomip'      ! scheme to compute ice shelf load (ln_isfcav = .true. in domain_cfg.nc) 
     167   cn_isfload = 'uniform'      ! scheme to compute ice shelf load (ln_isfcav = .true. in domain_cfg.nc) 
     168      rn_isfload_T = -1.0 
     169      rn_isfload_S =  34.2 
    168170   ! 
    169171   ! ---------------- ice shelf melt formulation ------------------------------- 
     
    182184         !                       ! oasis = fwfisf is given by oasis and pattern by file sn_isfcav_fwf 
    183185         !              !  cn_isfcav_mlt = 2eq or 3eq cases: 
    184          cn_gammablk = 'ad15'    ! scheme to compute gammat/s (spe,ad15,hj99) 
    185          !                       ! ad15 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
    186          !                       ! hj99 = velocity and stability dependent Gamma    (Holland et al. 1999) 
     186         cn_gammablk = 'vel'     ! scheme to compute gammat/s (spe,ad15,hj99) 
     187         !                       ! spe      = constant transfert velocity (rn_gammat0, rn_gammas0) 
     188         !                       ! vel      = velocity dependent transfert velocity (u* * gammat/s) (Asay-Davis et al. 2016 for a short description) 
     189         !                       ! vel_stab = velocity and stability dependent transfert coeficient (Holland et al. 1999 for a complete description) 
    187190         rn_gammat0  = 0.0215    ! gammat coefficient used in blk formula 
    188191         rn_gammas0  = 0.614e-3  ! gammas coefficient used in blk formula 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/ISOMIP+/MY_SRC/isfcavgam.F90

    r11987 r12062  
    6565         ! gamma is constant (specified in namelist) 
    6666         ! nothing to do 
    67       CASE ('ad15', 'hj99') 
     67      CASE ('vel', 'vel_stab') 
    6868         ! compute velocity in tbl 
    6969         CALL isf_tbl(un(:,:,:) ,zutbl(:,:),'U', miku, rhisf_tbl_cav) 
     
    8989         pgt(:,:) = rn_gammat0 
    9090         pgs(:,:) = rn_gammas0 
    91       CASE ( 'ad15' ) ! gamma is proportional to u* 
    92          CALL gammats_AD15 (              zutbl, zvtbl, rCd0_top, rn_vtide**2,               pgt, pgs ) 
    93       CASE ( 'hj99' ) ! gamma depends of stability of boundary layer and u* 
    94          CALL gammats_HJ99 (pttbl, pstbl, zutbl, zvtbl, rCd0_top, r_ke0_top, pqoce, pqfwf, pgt, pgs ) 
     91      CASE ( 'vel' ) ! gamma is proportional to u* 
     92         CALL gammats_vel      (              zutbl, zvtbl, rCd0_top, rn_vtide**2,             pgt, pgs ) 
     93      CASE ( 'vel_stab' ) ! gamma depends of stability of boundary layer and u* 
     94         CALL gammats_vel_stab (pttbl, pstbl, zutbl, zvtbl, rCd0_top, r_ke0_top, pqoce, pqfwf, pgt, pgs ) 
    9595      CASE DEFAULT 
    9696         CALL ctl_stop('STOP','method to compute gamma (cn_gammablk) is unknown (should not see this)') 
     
    115115   !!----------------------------------------------------------------------------------------------------- 
    116116   ! 
    117    SUBROUTINE gammats_AD15(putbl, pvtbl, pCd, pke2, &   ! <<== in 
     117   SUBROUTINE gammats_vel( putbl, pvtbl, pCd, pke2, &   ! <<== in 
    118118      &                                  pgt, pgs   )   ! ==>> out gammats [m/s] 
    119119      !!---------------------------------------------------------------------- 
     
    145145      CALL iom_put('isfustar',zustar(:,:)) 
    146146      ! 
    147    END SUBROUTINE gammats_AD15 
    148  
    149    SUBROUTINE gammats_HJ99( pttbl, pstbl, putbl, pvtbl, pCd, pke2, pqoce, pqfwf, &  ! <<== in 
    150       &                                                            pgt  , pgs    )  ! ==>> out gammats [m/s] 
     147   END SUBROUTINE gammats_vel 
     148 
     149   SUBROUTINE gammats_vel_stab( pttbl, pstbl, putbl, pvtbl, pCd, pke2, pqoce, pqfwf, &  ! <<== in 
     150      &                                                                pgt  , pgs    )  ! ==>> out gammats [m/s] 
    151151      !!---------------------------------------------------------------------- 
    152152      !! ** Purpose    : compute the coefficient echange coefficient  
     
    252252      END DO 
    253253 
    254    END SUBROUTINE gammats_HJ99 
     254   END SUBROUTINE gammats_vel_stab 
    255255 
    256256END MODULE isfcavgam 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/ISOMIP+/MY_SRC/isfstp.F90

    r11987 r12062  
    220220         END IF 
    221221 
    222          IF (ln_isfcav) WRITE(numout,*) '      Ice shelf load method                   cn_isfload        = ', TRIM(cn_isfload) 
     222         IF (ln_isfcav) THEN 
     223            WRITE(numout,*) '      Ice shelf load method                   cn_isfload        = ', TRIM(cn_isfload) 
     224            WRITE(numout,*) '         Temperature used to compute the ice shelf load            = ', rn_isfload_T 
     225            WRITE(numout,*) '         Salinity    used to compute the ice shelf load            = ', rn_isfload_S 
     226         END IF 
    223227         WRITE(numout,*) '' 
    224228         FLUSH(numout) 
     
    276280      INTEGER               :: ios                  ! Local integer output status for namelist read 
    277281      !!---------------------------------------------------------------------- 
    278       NAMELIST/namisf/ ln_isf       ,                                                                               &  
    279          &             ln_isfcav_mlt, cn_isfcav_mlt, cn_gammablk, rn_gammat0, rn_gammas0, rn_htbl, sn_isfcav_fwf,   & 
    280          &             ln_isfpar_mlt, cn_isfpar_mlt, sn_isfpar_fwf, sn_isfpar_zmin, sn_isfpar_zmax, sn_isfpar_Leff, & 
    281          &             ln_isfcpl    , nn_drown, ln_isfcpl_cons, ln_isfdebug, rn_vtide,                              & 
    282          &             cn_isfload   , cn_isfdir 
     282      NAMELIST/namisf/ ln_isf        ,                                                           &  
     283         &             cn_gammablk   , rn_gammat0    , rn_gammas0    , rn_htbl, sn_isfcav_fwf,   & 
     284         &             ln_isfcav_mlt , cn_isfcav_mlt , sn_isfcav_fwf ,                           & 
     285         &             ln_isfpar_mlt , cn_isfpar_mlt , sn_isfpar_fwf ,                           & 
     286         &             sn_isfpar_zmin, sn_isfpar_zmax, sn_isfpar_Leff,                           & 
     287         &             ln_isfcpl     , nn_drown      , ln_isfcpl_cons, ln_isfdebug, rn_vtide,    & 
     288         &             cn_isfload    , rn_isfload_T  , rn_isfload_S  , cn_isfdir 
    283289      !!---------------------------------------------------------------------- 
    284290      ! 
Note: See TracChangeset for help on using the changeset viewer.