Changeset 12377 for NEMO/trunk/src/ICE/icestp.F90
- Timestamp:
- 2020-02-12T15:39:06+01:00 (4 years ago)
- Location:
- NEMO/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEAD ext/AGRIF5 ^/vendors/AGRIF/dev_r11615_ENHANCE-04_namelists_as_internalfiles_agrif@HEAD ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL
-
- Property svn:externals
-
NEMO/trunk/src/ICE/icestp.F90
r11536 r12377 86 86 PUBLIC ice_init ! called by sbcmod.F90 87 87 88 !! * Substitutions89 # include "vectopt_loop_substitute.h90"90 88 !!---------------------------------------------------------------------- 91 89 !! NEMO/ICE 4.0 , NEMO Consortium (2018) … … 95 93 CONTAINS 96 94 97 SUBROUTINE ice_stp( kt, ksbc )95 SUBROUTINE ice_stp( kt, Kbb, Kmm, ksbc ) 98 96 !!--------------------------------------------------------------------- 99 97 !! *** ROUTINE ice_stp *** … … 115 113 !! utau, vtau, taum, wndm, qns , qsr, emp , sfx 116 114 !!--------------------------------------------------------------------- 117 INTEGER, INTENT(in) :: kt ! ocean time step 118 INTEGER, INTENT(in) :: ksbc ! flux formulation (user defined, bulk, or Pure Coupled) 115 INTEGER, INTENT(in) :: kt ! ocean time step 116 INTEGER, INTENT(in) :: Kbb, Kmm ! ocean time level indices 117 INTEGER, INTENT(in) :: ksbc ! flux formulation (user defined, bulk, or Pure Coupled) 119 118 ! 120 119 INTEGER :: jl ! dummy loop index … … 160 159 ! 161 160 IF( ln_icedyn .AND. .NOT.lk_c1d ) & 162 & CALL ice_dyn( kt )! -- Ice dynamics161 & CALL ice_dyn( kt, Kmm ) ! -- Ice dynamics 163 162 ! 164 163 ! !== lateral boundary conditions ==! … … 209 208 ! --- Ocean time step --- ! 210 209 !-------------------------! 211 IF( ln_icedyn ) CALL ice_update_tau( kt, u b(:,:,1), vb(:,:,1) ) ! -- update surface ocean stresses210 IF( ln_icedyn ) CALL ice_update_tau( kt, uu(:,:,1,Kbb), vv(:,:,1,Kbb) ) ! -- update surface ocean stresses 212 211 !!gm remark, the ocean-ice stress is not saved in ice diag call above ..... find a solution!!! 213 212 ! … … 217 216 218 217 219 SUBROUTINE ice_init 218 SUBROUTINE ice_init( Kbb, Kmm, Kaa ) 220 219 !!---------------------------------------------------------------------- 221 220 !! *** ROUTINE ice_init *** … … 223 222 !! ** purpose : Initialize sea-ice parameters 224 223 !!---------------------------------------------------------------------- 224 INTEGER, INTENT(in) :: Kbb, Kmm, Kaa 225 ! 225 226 INTEGER :: ji, jj, ierr 226 227 !!---------------------------------------------------------------------- … … 232 233 IF(lwp) WRITE(numout,*) '~~~~~~~~' 233 234 ! 234 ! ! Open the reference and configuration namelist files andnamelist output file235 CALL ctl_opn( numnam_ice_ref, 'namelist_ice_ref', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp)236 CALL ctl_opn( numnam_ice_cfg, 'namelist_ice_cfg', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp)237 IF(lwm) CALL ctl_opn( numoni , 'output.namelist.ice', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, 1 )235 ! ! Load the reference and configuration namelist files and open namelist output file 236 CALL load_nml( numnam_ice_ref, 'namelist_ice_ref', numout, lwm ) 237 CALL load_nml( numnam_ice_cfg, 'namelist_ice_cfg', numout, lwm ) 238 IF(lwm) CALL ctl_opn( numoni , 'output.namelist.ice', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, 1 ) 238 239 ! 239 240 CALL par_init ! set some ice run parameters … … 254 255 IF( .NOT. ln_rstart ) THEN ! start from rest: sea-ice deduced from sst 255 256 CALL ice_istate_init 256 CALL ice_istate( nit000 )257 CALL ice_istate( nit000, Kbb, Kmm, Kaa ) 257 258 ELSE ! start from a restart file 258 CALL ice_rst_read 259 CALL ice_rst_read( Kbb, Kmm, Kaa ) 259 260 ENDIF 260 261 CALL ice_var_glo2eqv … … 301 302 !!------------------------------------------------------------------- 302 303 ! 303 REWIND( numnam_ice_ref ) ! Namelist nampar in reference namelist : Parameters for ice304 304 READ ( numnam_ice_ref, nampar, IOSTAT = ios, ERR = 901) 305 305 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampar in reference namelist' ) 306 REWIND( numnam_ice_cfg ) ! Namelist nampar in configuration namelist : Parameters for ice307 306 READ ( numnam_ice_cfg, nampar, IOSTAT = ios, ERR = 902 ) 308 307 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'nampar in configuration namelist' )
Note: See TracChangeset
for help on using the changeset viewer.