Changeset 10534 for NEMO/trunk/src/ICE/icethd_zdf_bl99.F90
- Timestamp:
- 2019-01-16T17:49:45+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/ICE/icethd_zdf_bl99.F90
r10531 r10534 36 36 CONTAINS 37 37 38 SUBROUTINE ice_thd_zdf_BL99( k_ jules)38 SUBROUTINE ice_thd_zdf_BL99( k_cnd ) 39 39 !!------------------------------------------------------------------- 40 40 !! *** ROUTINE ice_thd_zdf_BL99 *** … … 73 73 !! total ice/snow thickness : h_i_1d, h_s_1d 74 74 !!------------------------------------------------------------------- 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) 76 76 ! 77 77 INTEGER :: ji, jk ! spatial loop index … … 164 164 ! 165 165 ! Store initial temperatures and non solar heat fluxes 166 IF( k_ jules == np_jules_OFF .OR. k_jules == np_jules_EMULE) THEN166 IF( k_cnd == np_cnd_OFF .OR. k_cnd == np_cnd_EMU ) THEN 167 167 ! 168 168 ztsub (1:npti) = t_su_1d(1:npti) ! surface temperature at iteration n-1 … … 332 332 !----------------------------------------! 333 333 ! ! 334 ! JULES COUPLING IS OFF OR EMULATED!334 ! Conduction flux is off or emulated ! 335 335 ! ! 336 336 !----------------------------------------! 337 337 ! 338 IF( k_ jules == np_jules_OFF .OR. k_jules == np_jules_EMULE) THEN338 IF( k_cnd == np_cnd_OFF .OR. k_cnd == np_cnd_EMU ) THEN 339 339 ! 340 340 ! ==> The original BL99 temperature computation is used … … 581 581 !----------------------------------------! 582 582 ! ! 583 ! JULES COUPLING IS ACTIVE!583 ! Conduction flux is on ! 584 584 ! ! 585 585 !----------------------------------------! 586 586 ! 587 ELSEIF( k_ jules == np_jules_ACTIVE) THEN587 ELSEIF( k_cnd == np_cnd_ON ) THEN 588 588 ! 589 589 ! ==> we use a modified BL99 solver with conduction flux (qcn_ice) as forcing term … … 754 754 END DO 755 755 756 ENDIF ! k_ jules756 ENDIF ! k_cnd 757 757 758 758 END DO ! End of the do while iterative procedure … … 781 781 ! --- Diagnose the heat loss due to changing non-solar / conduction flux --- ! 782 782 ! 783 IF( k_ jules == np_jules_OFF .OR. k_jules == np_jules_EMULE) THEN783 IF( k_cnd == np_cnd_OFF .OR. k_cnd == np_cnd_EMU ) THEN 784 784 ! 785 785 DO ji = 1, npti … … 787 787 END DO 788 788 ! 789 ELSEIF( k_ jules == np_jules_ACTIVE) THEN789 ELSEIF( k_cnd == np_cnd_ON ) THEN 790 790 ! 791 791 DO ji = 1, npti … … 798 798 ! --- Diagnose the heat loss due to non-fully converged temperature solution (should not be above 10-4 W-m2) --- ! 799 799 ! 800 IF( k_ jules == np_jules_OFF .OR. k_jules == np_jules_ACTIVE) THEN800 IF( k_cnd == np_cnd_OFF .OR. k_cnd == np_cnd_ON ) THEN 801 801 802 802 CALL ice_var_enthalpy … … 807 807 & SUM( e_s_1d(ji,1:nlay_s) ) * h_s_1d(ji) * r1_nlay_s ) 808 808 809 IF( k_ jules == np_jules_OFF ) THEN809 IF( k_cnd == np_cnd_OFF ) THEN 810 810 811 811 IF( t_su_1d(ji) < rt0 ) THEN ! case T_su < 0degC … … 817 817 ENDIF 818 818 819 ELSEIF( k_ jules == np_jules_ACTIVE) THEN819 ELSEIF( k_cnd == np_cnd_ON ) THEN 820 820 821 821 zhfx_err = ( qcn_ice_top_1d(ji) + qtr_ice_top_1d(ji) - zradtr_i(ji,nlay_i) - qcn_ice_bot_1d(ji) & … … 834 834 ENDIF 835 835 ! 836 !-------------------------------------------------------------------- -------------------837 ! 11) Jules coupling:reset inner snow and ice temperatures, update conduction fluxes838 !-------------------------------------------------------------------- -------------------836 !-------------------------------------------------------------------- 837 ! 11) reset inner snow and ice temperatures, update conduction fluxes 838 !-------------------------------------------------------------------- 839 839 ! effective conductivity and 1st layer temperature (needed by Met Office) 840 840 DO ji = 1, npti … … 851 851 END DO 852 852 ! 853 IF( k_ jules == np_jules_EMULE) THEN853 IF( k_cnd == np_cnd_EMU ) THEN 854 854 ! Restore temperatures to their initial values 855 855 t_s_1d (1:npti,:) = ztsold (1:npti,:)
Note: See TracChangeset
for help on using the changeset viewer.