Changeset 8889


Ignore:
Timestamp:
2017-12-04T12:27:10+01:00 (3 years ago)
Author:
clem
Message:

change namelist nn_iceflx into the more explicit nn_flxdist

Location:
branches/2017/dev_CNRS_2017/NEMOGCM
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_CNRS_2017/NEMOGCM/CONFIG/SHARED/namelist_ice_lim3_ref

    r8881 r8889  
    113113   rn_blow_s        =   0.66          !  mesure of snow blowing into the leads 
    114114                                      !     = 1 => no snow blowing, < 1 => some snow blowing 
    115    nn_iceflx        =  -1             !  Redistribute heat flux over ice categories 
     115   nn_flxdist       =  -1             !  Redistribute heat flux over ice categories 
    116116                                      !     =-1  Do nothing (needs N(cat) fluxes) 
    117                                       !          ==> forced mode only 
    118117                                      !     = 0  Average N(cat) fluxes then apply the average over the N(cat) ice 
    119                                       !          ==> forced and coupled modes 
    120118                                      !     = 1  Average N(cat) fluxes then redistribute over the N(cat) ice using T-ice and albedo sensitivity 
    121                                       !          ==> forced mode only 
    122119                                      !     = 2  Redistribute a single flux over categories 
    123                                       !          ==> coupled mode only 
    124120   nice_jules       = 0               !  Jules coupling (0=OFF, 1=EMULATED, 2=ACTIVE) 
    125121/ 
  • branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/LIM_SRC_3/ice.F90

    r8882 r8889  
    189189                                          ! -- icethd -- ! 
    190190   REAL(wp), PUBLIC ::   rn_cio           !: drag coefficient for oceanic stress 
    191    INTEGER , PUBLIC ::   nn_iceflx        !: Redistribute heat flux over ice categories 
     191   INTEGER , PUBLIC ::   nn_flxdist       !: Redistribute heat flux over ice categories 
    192192   !                                      !   =-1  Do nothing (needs N(cat) fluxes) 
    193193   !                                      !   = 0  Average N(cat) fluxes then apply the average over the N(cat) ice  
  • branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/LIM_SRC_3/iceforcing.F90

    r8884 r8889  
    142142      CASE( jp_blk )                !--- bulk formulation 
    143143                                CALL blk_ice_flx( t_su, h_s, h_i, alb_ice )    !  
    144          IF( ln_mixcpl      )   CALL sbc_cpl_ice_flx( picefr=at_i_b, palbi=alb_ice, psst=sst_m, pist=t_su, phs=h_s, phi=h_i ) 
    145          IF( nn_iceflx /= 2 )   CALL ice_flx_dist( t_su, alb_ice, qns_ice, qsr_ice, dqns_ice, evap_ice, devap_ice, nn_iceflx ) 
     144         IF( ln_mixcpl       )  CALL sbc_cpl_ice_flx( picefr=at_i_b, palbi=alb_ice, psst=sst_m, pist=t_su, phs=h_s, phi=h_i ) 
     145         IF( nn_flxdist /= 2 )  CALL ice_flx_dist( t_su, alb_ice, qns_ice, qsr_ice, dqns_ice, evap_ice, devap_ice, nn_flxdist ) 
    146146         ! 
    147147      CASE ( jp_purecpl )           !--- coupled formulation 
    148148                                CALL sbc_cpl_ice_flx( picefr=at_i_b, palbi=alb_ice, psst=sst_m, pist=t_su, phs=h_s, phi=h_i ) 
    149          IF( nn_iceflx == 2 )   CALL ice_flx_dist( t_su, alb_ice, qns_ice, qsr_ice, dqns_ice, evap_ice, devap_ice, nn_iceflx ) 
     149         IF( nn_flxdist == 2 )  CALL ice_flx_dist( t_su, alb_ice, qns_ice, qsr_ice, dqns_ice, evap_ice, devap_ice, nn_flxdist ) 
    150150         ! 
    151151      END SELECT 
     
    155155         ELSEWHERE                   ;   zalb(:,:) = SUM( alb_ice * a_i_b, dim=3 ) / at_i_b 
    156156         END WHERE 
    157          CALL iom_put( "icealb" , zalb(:,:) )   ! ice albedo 
    158       ENDIF 
    159  
    160       IF( iom_use('albedo') ) THEN  !--- surface albedo 
     157         CALL iom_put( "icealb" , zalb(:,:) ) 
     158      ENDIF 
     159 
     160      IF( iom_use('albedo') ) THEN    !--- output surface albedo 
    161161         zalb(:,:) = SUM( alb_ice * a_i_b, dim=3 ) + rn_alb_oce * ( 1._wp - at_i_b ) 
    162162         CALL iom_put( "albedo" , zalb(:,:) ) 
     
    267267      INTEGER ::   ios, ioptio   ! Local integer output status for namelist read 
    268268      !! 
    269       NAMELIST/namforcing/ rn_cio, rn_blow_s, nn_iceflx, nice_jules 
     269      NAMELIST/namforcing/ rn_cio, rn_blow_s, nn_flxdist, nice_jules 
    270270      !!------------------------------------------------------------------- 
    271271      ! 
     
    286286         WRITE(numout,*) '      drag coefficient for oceanic stress              rn_cio     = ', rn_cio 
    287287         WRITE(numout,*) '      coefficient for ice-lead partition of snowfall   rn_blow_s  = ', rn_blow_s 
    288          WRITE(numout,*) '      Multicategory heat flux formulation              nn_iceflx  = ', nn_iceflx 
     288         WRITE(numout,*) '      Multicategory heat flux formulation              nn_flxdist = ', nn_flxdist 
    289289         WRITE(numout,*) '      Jules coupling (0=no, 1=emulated, 2=active)      nice_jules = ', nice_jules 
    290290      ENDIF 
    291291      ! 
    292292      IF(lwp) WRITE(numout,*) 
    293       SELECT CASE( nn_iceflx )         ! ESIM Multi-category heat flux formulation 
     293      SELECT CASE( nn_flxdist )         ! ESIM Multi-category heat flux formulation 
    294294      CASE( -1  ) 
    295          IF(lwp) WRITE(numout,*) '   ESIM: use per-category fluxes (nn_iceflx = -1) ' 
    296          IF( ln_cpl )   CALL ctl_stop( 'ice_thd_init : the chosen nn_iceflx for ESIM in coupled mode must be 0 or 2' ) 
     295         IF(lwp) WRITE(numout,*) '   ESIM: use per-category fluxes (nn_flxdist = -1) ' 
     296         IF( ln_cpl )   CALL ctl_stop( 'ice_thd_init : the chosen nn_flxdist for ESIM in coupled mode must be 0 or 2' ) 
    297297      CASE(  0  ) 
    298          IF(lwp) WRITE(numout,*) '   ESIM: use average per-category fluxes (nn_iceflx = 0) ' 
     298         IF(lwp) WRITE(numout,*) '   ESIM: use average per-category fluxes (nn_flxdist = 0) ' 
    299299      CASE(  1  ) 
    300          IF(lwp) WRITE(numout,*) '   ESIM: use average then redistribute per-category fluxes (nn_iceflx = 1) ' 
    301          IF( ln_cpl )   CALL ctl_stop( 'ice_thd_init : the chosen nn_iceflx for ESIM in coupled mode must be 0 or 2' ) 
     300         IF(lwp) WRITE(numout,*) '   ESIM: use average then redistribute per-category fluxes (nn_flxdist = 1) ' 
     301         IF( ln_cpl )   CALL ctl_stop( 'ice_thd_init : the chosen nn_flxdist for ESIM in coupled mode must be 0 or 2' ) 
    302302      CASE(  2  ) 
    303          IF(lwp) WRITE(numout,*) '   ESIM: Redistribute a single flux over categories (nn_iceflx = 2) ' 
    304          IF( .NOT. ln_cpl )   CALL ctl_stop( 'ice_thd_init : the chosen nn_iceflx for ESIM in forced mode cannot be 2' ) 
     303         IF(lwp) WRITE(numout,*) '   ESIM: Redistribute a single flux over categories (nn_flxdist = 2) ' 
     304         IF( .NOT. ln_cpl )   CALL ctl_stop( 'ice_thd_init : the chosen nn_flxdist for ESIM in forced mode cannot be 2' ) 
    305305      CASE DEFAULT 
    306          CALL ctl_stop( 'ice_thd_init: ESIM option, nn_iceflx, should be between -1 and 2' ) 
     306         CALL ctl_stop( 'ice_thd_init: ESIM option, nn_flxdist, should be between -1 and 2' ) 
    307307      END SELECT 
    308308      ! 
Note: See TracChangeset for help on using the changeset viewer.