New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 8889 for branches/2017/dev_CNRS_2017/NEMOGCM/NEMO – NEMO

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

change namelist nn_iceflx into the more explicit nn_flxdist

Location:
branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/LIM_SRC_3
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • 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.