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 10534 for NEMO/trunk/src/ICE/icethd_zdf_bl99.F90 – NEMO

Ignore:
Timestamp:
2019-01-16T17:49:45+01:00 (5 years ago)
Author:
clem
Message:

change some more names in the ice

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/ICE/icethd_zdf_bl99.F90

    r10531 r10534  
    3636CONTAINS 
    3737 
    38    SUBROUTINE ice_thd_zdf_BL99( k_jules ) 
     38   SUBROUTINE ice_thd_zdf_BL99( k_cnd ) 
    3939      !!------------------------------------------------------------------- 
    4040      !!                ***  ROUTINE ice_thd_zdf_BL99  *** 
     
    7373      !!           total ice/snow thickness         : h_i_1d, h_s_1d 
    7474      !!------------------------------------------------------------------- 
    75       INTEGER, INTENT(in) ::   k_jules     ! Jules coupling (0=OFF, 1=EMULATED, 2=ACTIVE) 
     75      INTEGER, INTENT(in) ::   k_cnd     ! conduction flux (off, on, emulated) 
    7676      ! 
    7777      INTEGER ::   ji, jk         ! spatial loop index 
     
    164164      ! 
    165165      ! Store initial temperatures and non solar heat fluxes 
    166       IF( k_jules == np_jules_OFF .OR. k_jules == np_jules_EMULE ) THEN 
     166      IF( k_cnd == np_cnd_OFF .OR. k_cnd == np_cnd_EMU ) THEN 
    167167         ! 
    168168         ztsub      (1:npti) = t_su_1d(1:npti)                          ! surface temperature at iteration n-1 
     
    332332         !----------------------------------------! 
    333333         !                                        ! 
    334          !   JULES COUPLING IS OFF OR EMULATED    ! 
     334         !   Conduction flux is off or emulated   ! 
    335335         !                                        ! 
    336336         !----------------------------------------! 
    337337         ! 
    338          IF( k_jules == np_jules_OFF .OR. k_jules == np_jules_EMULE ) THEN 
     338         IF( k_cnd == np_cnd_OFF .OR. k_cnd == np_cnd_EMU ) THEN 
    339339            ! 
    340340            ! ==> The original BL99 temperature computation is used 
     
    581581         !----------------------------------------! 
    582582         !                                        ! 
    583          !      JULES COUPLING IS ACTIVE          ! 
     583         !      Conduction flux is on             ! 
    584584         !                                        ! 
    585585         !----------------------------------------! 
    586586         ! 
    587          ELSEIF( k_jules == np_jules_ACTIVE ) THEN 
     587         ELSEIF( k_cnd == np_cnd_ON ) THEN 
    588588            ! 
    589589            ! ==> we use a modified BL99 solver with conduction flux (qcn_ice) as forcing term 
     
    754754            END DO 
    755755 
    756          ENDIF ! k_jules 
     756         ENDIF ! k_cnd 
    757757          
    758758      END DO  ! End of the do while iterative procedure 
     
    781781      ! --- Diagnose the heat loss due to changing non-solar / conduction flux --- ! 
    782782      ! 
    783       IF( k_jules == np_jules_OFF .OR. k_jules == np_jules_EMULE ) THEN 
     783      IF( k_cnd == np_cnd_OFF .OR. k_cnd == np_cnd_EMU ) THEN 
    784784         ! 
    785785         DO ji = 1, npti 
     
    787787         END DO 
    788788         ! 
    789       ELSEIF( k_jules == np_jules_ACTIVE ) THEN 
     789      ELSEIF( k_cnd == np_cnd_ON ) THEN 
    790790         ! 
    791791         DO ji = 1, npti 
     
    798798      ! --- Diagnose the heat loss due to non-fully converged temperature solution (should not be above 10-4 W-m2) --- ! 
    799799      ! 
    800       IF( k_jules == np_jules_OFF .OR. k_jules == np_jules_ACTIVE ) THEN 
     800      IF( k_cnd == np_cnd_OFF .OR. k_cnd == np_cnd_ON ) THEN 
    801801          
    802802         CALL ice_var_enthalpy        
     
    807807               &                   SUM( e_s_1d(ji,1:nlay_s) ) * h_s_1d(ji) * r1_nlay_s ) 
    808808             
    809             IF( k_jules == np_jules_OFF ) THEN 
     809            IF( k_cnd == np_cnd_OFF ) THEN 
    810810                
    811811               IF( t_su_1d(ji) < rt0 ) THEN  ! case T_su < 0degC 
     
    817817               ENDIF 
    818818                
    819             ELSEIF( k_jules == np_jules_ACTIVE ) THEN 
     819            ELSEIF( k_cnd == np_cnd_ON ) THEN 
    820820             
    821821               zhfx_err    = ( qcn_ice_top_1d(ji) + qtr_ice_top_1d(ji) - zradtr_i(ji,nlay_i) - qcn_ice_bot_1d(ji)  & 
     
    834834      ENDIF 
    835835      ! 
    836       !--------------------------------------------------------------------------------------- 
    837       ! 11) Jules coupling: reset inner snow and ice temperatures, update conduction fluxes 
    838       !--------------------------------------------------------------------------------------- 
     836      !-------------------------------------------------------------------- 
     837      ! 11) reset inner snow and ice temperatures, update conduction fluxes 
     838      !-------------------------------------------------------------------- 
    839839      ! effective conductivity and 1st layer temperature (needed by Met Office) 
    840840      DO ji = 1, npti 
     
    851851      END DO 
    852852      ! 
    853       IF( k_jules == np_jules_EMULE ) THEN 
     853      IF( k_cnd == np_cnd_EMU ) THEN 
    854854         ! Restore temperatures to their initial values 
    855855         t_s_1d    (1:npti,:) = ztsold        (1:npti,:) 
Note: See TracChangeset for help on using the changeset viewer.