- Timestamp:
- 2017-12-26T17:32:56+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icethd_sal.F90
r9019 r9169 12 12 !! 'key_lim3' ESIM sea-ice model 13 13 !!---------------------------------------------------------------------- 14 !! ice_thd_sal 15 !! ice_thd_sal_init 14 !! ice_thd_sal : salinity variations in the ice 15 !! ice_thd_sal_init : initialization 16 16 !!---------------------------------------------------------------------- 17 17 USE dom_oce ! ocean space and time domain … … 32 32 33 33 ! ** namelist (namthd_sal) ** 34 REAL(wp) :: rn_sal_gd 35 REAL(wp) :: rn_time_gd 36 REAL(wp) :: rn_sal_fl 37 REAL(wp) :: rn_time_fl 34 REAL(wp) :: rn_sal_gd ! restoring salinity for gravity drainage [PSU] 35 REAL(wp) :: rn_time_gd ! restoring time constant for gravity drainage (= 20 days) [s] 36 REAL(wp) :: rn_sal_fl ! restoring salinity for flushing [PSU] 37 REAL(wp) :: rn_time_fl ! restoring time constant for gravity drainage (= 10 days) [s] 38 38 39 39 !!---------------------------------------------------------------------- … … 56 56 !!--------------------------------------------------------------------- 57 57 LOGICAL, INTENT(in) :: ld_sal ! gravity drainage and flushing or not 58 ! 58 59 INTEGER :: ji, jk ! dummy loop indices 59 60 REAL(wp) :: iflush, igravdr ! local scalars … … 66 67 ! !---------------------------------------------! 67 68 CASE( 2 ) ! time varying salinity with linear profile ! 68 !!---------------------------------------------!69 ! !---------------------------------------------! 69 70 z1_time_gd = 1._wp / rn_time_gd * rdt_ice 70 71 z1_time_fl = 1._wp / rn_time_fl * rdt_ice 71 72 ! 72 73 DO ji = 1, npti 73 74 ! 74 75 !--------------------------------------------------------- 75 76 ! Update ice salinity from snow-ice and bottom growth … … 82 83 s_i_1d(ji) = s_i_1d(ji) + zs_i_bg + zs_i_si 83 84 ENDIF 84 85 ! 85 86 IF( ld_sal ) THEN 86 87 !--------------------------------------------------------- … … 100 101 ENDIF 101 102 END DO 102 103 ! 103 104 ! Salinity profile 104 105 CALL ice_var_salprof1d 105 106 ! 106 ! !---------------------------------------------!107 CASE( 3 ) ! constant salinity with a fixed profile ! (Schwarzacher (1959) multiyear salinity profile(mean = 2.30)108 ! !---------------------------------------------!107 ! !----------------------------------------! 108 CASE( 3 ) ! constant salinity with a fixed profile ! (Schwarzacher (1959) multiyear salinity profile (mean = 2.30) 109 ! !----------------------------------------! 109 110 CALL ice_var_salprof1d 111 ! 112 END SELECT 110 113 ! 111 END SELECT112 !113 114 END SUBROUTINE ice_thd_sal 114 115 … … 125 126 !! ** input : Namelist namthd_sal 126 127 !!------------------------------------------------------------------- 127 INTEGER :: ios ! Local integer output status for namelist read128 INTEGER :: ios ! Local integer 128 129 !! 129 130 NAMELIST/namthd_sal/ nn_icesal, rn_icesal, rn_sal_gd, rn_time_gd, & … … 133 134 REWIND( numnam_ice_ref ) ! Namelist namthd_sal in reference namelist : Ice salinity 134 135 READ ( numnam_ice_ref, namthd_sal, IOSTAT = ios, ERR = 901) 135 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namthd_sal in reference namelist', lwp ) 136 ! 136 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namthd_sal in reference namelist', lwp ) 137 137 REWIND( numnam_ice_cfg ) ! Namelist namthd_sal in configuration namelist : Ice salinity 138 138 READ ( numnam_ice_cfg, namthd_sal, IOSTAT = ios, ERR = 902 ) 139 902 IF( ios /= 0 )CALL ctl_nam ( ios , 'namthd_sal in configuration namelist', lwp )139 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'namthd_sal in configuration namelist', lwp ) 140 140 IF(lwm) WRITE ( numoni, namthd_sal ) 141 141 !
Note: See TracChangeset
for help on using the changeset viewer.