- Timestamp:
- 2019-10-29T12:00:04+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfstp.F90
r11553 r11823 15 15 USE oce ! ocean dynamics and tracers 16 16 USE dom_oce ! ocean space and time domain 17 USE domvvl, ONLY : ln_vvl_zstar 17 18 USE phycst ! physical constants 18 19 USE eosbn2 ! equation of state … … 120 121 !! ** Purpose : Initialisation of the ice shelf public variables 121 122 !! 122 !! ** Method : Read the namsbc namelist and set derived parameters 123 !! Call init routines for all other SBC modules that have one 124 !! 125 !! ** Action : - read namsbc parameters 123 !! ** Method : Read the namisf namelist, check option compatibility and set derived parameters 124 !! 125 !! ** Action : - read namisf parameters 126 126 !! - allocate memory 127 !! - call cav/param init routine 127 !! - output print 128 !! - ckeck option compatibility 129 !! - call cav/param/isfcpl init routine 128 130 !!---------------------------------------------------------------------- 129 131 ! … … 144 146 ! terminate routine now if no ice shelf melt formulation specify 145 147 IF ( ln_isf ) THEN 146 !147 !---------------------------------------------------------------------------------------------------------------------148 ! initialisation melt in the cavity149 IF ( ln_isfcav_mlt ) THEN150 !151 ! initialisation of cav variable152 CALL isf_cav_init()153 !154 ! read cav variable from restart155 IF ( ln_rstart ) CALL isfrst_read('cav', risf_cav_tsc, fwfisf_cav, risf_cav_tsc_b, fwfisf_cav_b)156 !157 END IF158 !159 !---------------------------------------------------------------------------------------------------------------------160 ! initialisation parametrised melt161 IF ( ln_isfpar_mlt ) THEN162 !163 ! initialisation of par variable164 CALL isf_par_init()165 !166 ! read par variable from restart167 IF ( ln_rstart ) CALL isfrst_read('par', risf_par_tsc, fwfisf_par, risf_par_tsc_b, fwfisf_par_b)168 !169 END IF170 148 ! 171 149 !--------------------------------------------------------------------------------------------------------------------- … … 173 151 IF( ln_isfcpl ) THEN 174 152 ! 153 ! start on an euler time step 154 neuler = 0 155 ! 156 ! ice sheet coupling: extrapolation of restart to fill new wet cell and compute divergence correction 157 CALL isfcpl_init() 158 ! 159 ! all before fields set to now values 160 tsb (:,:,:,:) = tsn (:,:,:,:) 161 ub (:,:,:) = un (:,:,:) 162 vb (:,:,:) = vn (:,:,:) 163 sshb (:,:) = sshn (:,:) 164 e3t_b(:,:,:) = e3t_n(:,:,:) 165 175 166 ! prepare writing restart 176 167 IF( lwxios ) CALL iom_set_rstw_var_active('ssmask') … … 182 173 END IF 183 174 ! 175 !--------------------------------------------------------------------------------------------------------------------- 176 ! initialisation melt in the cavity 177 IF ( ln_isfcav_mlt ) THEN 178 ! 179 ! initialisation of cav variable 180 CALL isf_cav_init() 181 ! 182 ! read cav variable from restart 183 IF ( ln_rstart ) CALL isfrst_read('cav', risf_cav_tsc, fwfisf_cav, risf_cav_tsc_b, fwfisf_cav_b) 184 ! 185 END IF 186 ! 187 !--------------------------------------------------------------------------------------------------------------------- 188 ! initialisation parametrised melt 189 IF ( ln_isfpar_mlt ) THEN 190 ! 191 ! initialisation of par variable 192 CALL isf_par_init() 193 ! 194 ! read par variable from restart 195 IF ( ln_rstart ) CALL isfrst_read('par', risf_par_tsc, fwfisf_par, risf_par_tsc_b, fwfisf_par_b) 196 ! 197 END IF 198 ! 184 199 END IF 185 200 186 201 END SUBROUTINE isf_init 187 202 188 SUBROUTINE isf_ctl()! 203 SUBROUTINE isf_ctl() 204 !!--------------------------------------------------------------------- 205 !! *** ROUTINE isf_ctl *** 206 !! 207 !! ** Purpose : output print and option compatibility check 208 !! 209 !!---------------------------------------------------------------------- 189 210 IF (lwp) THEN 190 211 WRITE(numout,*) … … 272 293 END IF 273 294 ! 295 ! compatibility ice shelf and vvl 296 IF( .NOT. ln_vvl_zstar .AND. ln_isf ) CALL ctl_stop( 'Only vvl_zstar has been tested with ice shelf cavity' ) 297 ! 274 298 END IF 275 299 END SUBROUTINE isf_ctl
Note: See TracChangeset
for help on using the changeset viewer.