- Timestamp:
- 2020-09-14T17:40:34+02:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11351_fldread_with_XIOS
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11351_fldread_with_XIOS
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEADext/AGRIF5 ^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL 8 9 # SETTE 10 ^/utils/CI/sette@13382 sette
-
- Property svn:externals
-
NEMO/branches/2019/dev_r11351_fldread_with_XIOS/src/ICE/icestp.F90
r10994 r13463 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 241 ! 242 #if defined key_agrif 243 CALL Agrif_Declare_Var_ice ! " " " " " Sea ice 244 #endif 240 245 ! 241 246 ! ! Allocate the ice arrays (sbc_ice already allocated in sbc_init) … … 254 259 IF( .NOT. ln_rstart ) THEN ! start from rest: sea-ice deduced from sst 255 260 CALL ice_istate_init 256 CALL ice_istate 261 CALL ice_istate( nit000, Kbb, Kmm, Kaa ) 257 262 ELSE ! start from a restart file 258 CALL ice_rst_read 263 CALL ice_rst_read( Kbb, Kmm, Kaa ) 259 264 ENDIF 260 265 CALL ice_var_glo2eqv … … 301 306 !!------------------------------------------------------------------- 302 307 ! 303 REWIND( numnam_ice_ref ) ! Namelist nampar in reference namelist : Parameters for ice304 308 READ ( numnam_ice_ref, nampar, IOSTAT = ios, ERR = 901) 305 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampar in reference namelist', lwp ) 306 REWIND( numnam_ice_cfg ) ! Namelist nampar in configuration namelist : Parameters for ice 309 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampar in reference namelist' ) 307 310 READ ( numnam_ice_cfg, nampar, IOSTAT = ios, ERR = 902 ) 308 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'nampar in configuration namelist' , lwp)311 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'nampar in configuration namelist' ) 309 312 IF(lwm) WRITE( numoni, nampar ) 310 313 ! … … 339 342 IF( ln_bdy .AND. ln_icediachk ) CALL ctl_warn('par_init: online conservation check does not work with BDY') 340 343 ! 341 r dt_ice = REAL(nn_fsbc) * rdt !--- sea-ice timestep and its inverse342 r1_ rdtice = 1._wp / rdt_ice344 rDt_ice = REAL(nn_fsbc) * rn_Dt !--- sea-ice timestep and its inverse 345 r1_Dt_ice = 1._wp / rDt_ice 343 346 IF(lwp) WRITE(numout,*) 344 IF(lwp) WRITE(numout,*) ' ice timestep r dt_ice = nn_fsbc*rdt = ', rdt_ice347 IF(lwp) WRITE(numout,*) ' ice timestep rDt_ice = nn_fsbc*rn_Dt = ', rDt_ice 345 348 ! 346 349 r1_nlay_i = 1._wp / REAL( nlay_i, wp ) !--- inverse of nlay_i and nlay_s … … 425 428 wfx_err_sub(:,:) = 0._wp 426 429 ! 427 afx_tot(:,:) = 0._wp ;428 !429 430 diag_heat(:,:) = 0._wp ; diag_sice(:,:) = 0._wp 430 431 diag_vice(:,:) = 0._wp ; diag_vsnw(:,:) = 0._wp
Note: See TracChangeset
for help on using the changeset viewer.