Changeset 6994
- Timestamp:
- 2016-10-05T15:07:10+02:00 (8 years ago)
- Location:
- branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM
- Files:
-
- 19 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/CONFIG/SHARED/namelist_ice_lim3_ref
r6989 r6994 36 36 &namicediag ! Diagnostics 37 37 !------------------------------------------------------------------------------ 38 ln_limdia hsb= .false. ! check online the heat, mass & salt budgets (T) or not (F)39 ln_limdia out= .false. ! output the heat, mass & salt budgets (T) or not (F)40 ln_ icectl = .false. ! ice points output for debug (T or F)38 ln_limdiachk = .false. ! check online the heat, mass & salt budgets (T) or not (F) 39 ln_limdiahsb = .false. ! output the heat, mass & salt budgets (T) or not (F) 40 ln_limctl = .false. ! ice points output for debug (T or F) 41 41 iiceprt = 10 ! i-index for debug 42 42 jiceprt = 10 ! j-index for debug -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/ice.F90
r6989 r6994 173 173 174 174 ! !!** ice-diagnostics namelist (namicediag) ** 175 LOGICAL , PUBLIC :: ln_limdiachk !: flag for ice diag (T) or not (F) 175 176 LOGICAL , PUBLIC :: ln_limdiahsb !: flag for ice diag (T) or not (F) 176 LOGICAL , PUBLIC :: ln_limdiaout !: flag for ice diag (T) or not (F) 177 LOGICAL , PUBLIC :: ln_icectl !: flag for sea-ice points output (T) or not (F) 177 LOGICAL , PUBLIC :: ln_limctl !: flag for sea-ice points output (T) or not (F) 178 178 INTEGER , PUBLIC :: iiceprt !: debug i-point 179 179 INTEGER , PUBLIC :: jiceprt !: debug j-point -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limcons.F90
r6970 r6994 164 164 !! + test if ice concentration and volume are > 0 165 165 !! 166 !! ** Method : This is an online diagnostics which can be activated with ln_limdia hsb=true166 !! ** Method : This is an online diagnostics which can be activated with ln_limdiachk=true 167 167 !! It prints in ocean.output if there is a violation of conservation at each time-step 168 168 !! The thresholds (zv_sill, zs_sill, zh_sill) which determine violations are set to … … 274 274 !! ** Purpose : Test the conservation of heat, salt and mass at the end of each ice time-step 275 275 !! 276 !! ** Method : This is an online diagnostics which can be activated with ln_limdia hsb=true276 !! ** Method : This is an online diagnostics which can be activated with ln_limdiachk=true 277 277 !! It prints in ocean.output if there is a violation of conservation at each time-step 278 278 !! The thresholds (zv_sill, zs_sill, zh_sill) which determine the violation are set to -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limctl.F90
r6746 r6994 5 5 !!====================================================================== 6 6 !! History : 3.5 ! 2015-01 (M. Vancoppenolle) Original code 7 !! 3.7 ! 2016-10 (C. Rousset) Add routine lim_prt3D 7 8 !!---------------------------------------------------------------------- 8 9 #if defined key_lim3 … … 12 13 !! lim_ctl : control prints in case of crash 13 14 !! lim_prt : ice control print at a given grid point 15 !! lim_prt3D : control prints of ice arrays 14 16 !!---------------------------------------------------------------------- 15 17 USE oce ! ocean dynamics and tracers … … 34 36 PUBLIC lim_ctl 35 37 PUBLIC lim_prt 38 PUBLIC lim_prt3D 36 39 37 40 !! * Substitutions … … 444 447 END SUBROUTINE lim_prt 445 448 449 SUBROUTINE lim_prt3D( cd_routine ) 450 !!--------------------------------------------------------------------------------------------------------- 451 !! *** ROUTINE lim_prt3D *** 452 !! 453 !! ** Purpose : CTL prints of ice arrays in case ln_ctl is activated 454 !! 455 !!--------------------------------------------------------------------------------------------------------- 456 CHARACTER(len=*), INTENT(in) :: cd_routine ! name of the routine 457 INTEGER :: jk, jl ! dummy loop indices 458 459 CALL prt_ctl_info(' ========== ') 460 CALL prt_ctl_info( cd_routine ) 461 CALL prt_ctl_info(' ========== ') 462 CALL prt_ctl_info(' - Cell values : ') 463 CALL prt_ctl_info(' ~~~~~~~~~~~~~ ') 464 CALL prt_ctl(tab2d_1=e12t , clinfo1=' cell area :') 465 CALL prt_ctl(tab2d_1=at_i , clinfo1=' at_i :') 466 CALL prt_ctl(tab2d_1=ato_i , clinfo1=' ato_i :') 467 CALL prt_ctl(tab2d_1=vt_i , clinfo1=' vt_i :') 468 CALL prt_ctl(tab2d_1=vt_s , clinfo1=' vt_s :') 469 CALL prt_ctl(tab2d_1=divu_i , clinfo1=' divu_i :') 470 CALL prt_ctl(tab2d_1=delta_i , clinfo1=' delta_i :') 471 CALL prt_ctl(tab2d_1=stress1_i , clinfo1=' stress1_i :') 472 CALL prt_ctl(tab2d_1=stress2_i , clinfo1=' stress2_i :') 473 CALL prt_ctl(tab2d_1=stress12_i , clinfo1=' stress12_i :') 474 CALL prt_ctl(tab2d_1=strength , clinfo1=' strength :') 475 CALL prt_ctl(tab2d_1=delta_i , clinfo1=' delta_i :') 476 CALL prt_ctl(tab2d_1=u_ice , clinfo1=' u_ice :', tab2d_2=v_ice , clinfo2=' v_ice :') 477 478 DO jl = 1, jpl 479 CALL prt_ctl_info(' ') 480 CALL prt_ctl_info(' - Category : ', ivar1=jl) 481 CALL prt_ctl_info(' ~~~~~~~~~~') 482 CALL prt_ctl(tab2d_1=ht_i (:,:,jl) , clinfo1= ' ht_i : ') 483 CALL prt_ctl(tab2d_1=ht_s (:,:,jl) , clinfo1= ' ht_s : ') 484 CALL prt_ctl(tab2d_1=t_su (:,:,jl) , clinfo1= ' t_su : ') 485 CALL prt_ctl(tab2d_1=t_s (:,:,1,jl) , clinfo1= ' t_snow : ') 486 CALL prt_ctl(tab2d_1=sm_i (:,:,jl) , clinfo1= ' sm_i : ') 487 CALL prt_ctl(tab2d_1=o_i (:,:,jl) , clinfo1= ' o_i : ') 488 CALL prt_ctl(tab2d_1=a_i (:,:,jl) , clinfo1= ' a_i : ') 489 CALL prt_ctl(tab2d_1=v_i (:,:,jl) , clinfo1= ' v_i : ') 490 CALL prt_ctl(tab2d_1=v_s (:,:,jl) , clinfo1= ' v_s : ') 491 CALL prt_ctl(tab2d_1=e_i (:,:,1,jl) , clinfo1= ' e_i1 : ') 492 CALL prt_ctl(tab2d_1=e_s (:,:,1,jl) , clinfo1= ' e_snow : ') 493 CALL prt_ctl(tab2d_1=smv_i (:,:,jl) , clinfo1= ' smv_i : ') 494 CALL prt_ctl(tab2d_1=oa_i (:,:,jl) , clinfo1= ' oa_i : ') 495 496 DO jk = 1, nlay_i 497 CALL prt_ctl_info(' - Layer : ', ivar1=jk) 498 CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' t_i : ') 499 END DO 500 END DO 501 502 CALL prt_ctl_info(' ') 503 CALL prt_ctl_info(' - Heat / FW fluxes : ') 504 CALL prt_ctl_info(' ~~~~~~~~~~~~~~~~~~ ') 505 CALL prt_ctl(tab2d_1=sst_m , clinfo1= ' sst : ', tab2d_2=sss_m , clinfo2= ' sss : ') 506 CALL prt_ctl(tab2d_1=qsr , clinfo1= ' qsr : ', tab2d_2=qns , clinfo2= ' qns : ') 507 CALL prt_ctl(tab2d_1=emp , clinfo1= ' emp : ', tab2d_2=sfx , clinfo2= ' sfx : ') 508 509 CALL prt_ctl_info(' ') 510 CALL prt_ctl_info(' - Stresses : ') 511 CALL prt_ctl_info(' ~~~~~~~~~~ ') 512 CALL prt_ctl(tab2d_1=utau , clinfo1= ' utau : ', tab2d_2=vtau , clinfo2= ' vtau : ') 513 CALL prt_ctl(tab2d_1=utau_ice , clinfo1= ' utau_ice : ', tab2d_2=vtau_ice , clinfo2= ' vtau_ice : ') 514 CALL prt_ctl(tab2d_1=u_oce , clinfo1= ' u_oce : ', tab2d_2=v_oce , clinfo2= ' v_oce : ') 515 516 END SUBROUTINE lim_prt3D 517 446 518 #else 447 519 !!-------------------------------------------------------------------------- … … 453 525 SUBROUTINE lim_prt ! Empty routine 454 526 END SUBROUTINE lim_prt 527 SUBROUTINE lim_prt3D ! Empty routine 528 END SUBROUTINE lim_prt3D 455 529 #endif 456 530 !!====================================================================== -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limdyn.F90
r6989 r6994 24 24 USE wrk_nemo ! work arrays 25 25 USE in_out_manager ! I/O manager 26 USE prtctl ! Print control27 26 USE lib_fortran ! glob_sum 28 27 USE timing ! Timing 29 28 USE limcons ! conservation tests 29 USE limctl ! control prints 30 30 USE limvar 31 31 … … 67 67 ! 68 68 ! conservation test 69 IF( ln_limdia hsb) CALL lim_cons_hsm(0, 'limdyn', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b)69 IF( ln_limdiachk ) CALL lim_cons_hsm(0, 'limdyn', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 70 70 71 71 ! ice velocities before rheology … … 85 85 CALL lim_rhg 86 86 ! 87 ! 87 ! conservation test 88 IF( ln_limdiachk ) CALL lim_cons_hsm(1, 'limdyn', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 89 88 90 ! Control prints 89 IF(ln_ctl) THEN 90 CALL prt_ctl_info(' ') 91 CALL prt_ctl_info(' - Cell values : ') 92 CALL prt_ctl_info(' ~~~~~~~~~~~~~ ') 93 CALL prt_ctl(tab2d_1=divu_i , clinfo1=' lim_dyn : divu_i :') 94 CALL prt_ctl(tab2d_1=delta_i , clinfo1=' lim_dyn : delta_i :') 95 CALL prt_ctl(tab2d_1=strength , clinfo1=' lim_dyn : strength :') 96 CALL prt_ctl(tab2d_1=e12t , clinfo1=' lim_dyn : cell area :') 97 CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_dyn : at_i :') 98 CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_dyn : vt_i :') 99 CALL prt_ctl(tab2d_1=vt_s , clinfo1=' lim_dyn : vt_s :') 100 CALL prt_ctl(tab2d_1=stress1_i , clinfo1=' lim_dyn : stress1_i :') 101 CALL prt_ctl(tab2d_1=stress2_i , clinfo1=' lim_dyn : stress2_i :') 102 CALL prt_ctl(tab2d_1=stress12_i, clinfo1=' lim_dyn : stress12_i:') 103 DO jl = 1, jpl 104 CALL prt_ctl_info(' ') 105 CALL prt_ctl_info(' - Category : ', ivar1=jl) 106 CALL prt_ctl_info(' ~~~~~~~~~~') 107 CALL prt_ctl(tab2d_1=a_i (:,:,jl) , clinfo1= ' lim_dyn : a_i : ') 108 CALL prt_ctl(tab2d_1=ht_i (:,:,jl) , clinfo1= ' lim_dyn : ht_i : ') 109 CALL prt_ctl(tab2d_1=ht_s (:,:,jl) , clinfo1= ' lim_dyn : ht_s : ') 110 CALL prt_ctl(tab2d_1=v_i (:,:,jl) , clinfo1= ' lim_dyn : v_i : ') 111 CALL prt_ctl(tab2d_1=v_s (:,:,jl) , clinfo1= ' lim_dyn : v_s : ') 112 CALL prt_ctl(tab2d_1=e_s (:,:,1,jl) , clinfo1= ' lim_dyn : e_s : ') 113 CALL prt_ctl(tab2d_1=t_su (:,:,jl) , clinfo1= ' lim_dyn : t_su : ') 114 CALL prt_ctl(tab2d_1=t_s (:,:,1,jl) , clinfo1= ' lim_dyn : t_snow : ') 115 CALL prt_ctl(tab2d_1=sm_i (:,:,jl) , clinfo1= ' lim_dyn : sm_i : ') 116 CALL prt_ctl(tab2d_1=smv_i (:,:,jl) , clinfo1= ' lim_dyn : smv_i : ') 117 DO jk = 1, nlay_i 118 CALL prt_ctl_info(' ') 119 CALL prt_ctl_info(' - Layer : ', ivar1=jk) 120 CALL prt_ctl_info(' ~~~~~~~') 121 CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_dyn : t_i : ') 122 CALL prt_ctl(tab2d_1=e_i(:,:,jk,jl) , clinfo1= ' lim_dyn : e_i : ') 123 END DO 124 END DO 125 ENDIF 126 ! 127 ! conservation test 128 IF( ln_limdiahsb ) CALL lim_cons_hsm(1, 'limdyn', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 91 IF( ln_ctl ) CALL lim_prt3D( 'limdyn' ) 129 92 ! 130 93 IF( nn_timing == 1 ) CALL timing_stop('limdyn') -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90
r6796 r6994 22 22 USE lib_mpp ! MPP library 23 23 USE wrk_nemo ! work arrays 24 USE prtctl ! Print control25 24 26 25 USE in_out_manager ! I/O manager 27 26 USE iom ! I/O manager 28 27 USE lib_fortran ! glob_sum 29 USE limdiahsb30 28 USE timing ! Timing 31 29 USE limcons ! conservation tests 30 USE limctl ! control prints 32 31 33 32 IMPLICIT NONE … … 125 124 CALL wrk_alloc( jpi,jpj, closing_net, divu_adv, opning, closing_gross ) 126 125 127 IF(ln_ctl) THEN128 CALL prt_ctl(tab2d_1=ato_i , clinfo1=' lim_itd_me: ato_i : ', tab2d_2=at_i , clinfo2=' at_i : ')129 CALL prt_ctl(tab2d_1=divu_i, clinfo1=' lim_itd_me: divu_i : ', tab2d_2=delta_i, clinfo2=' delta_i : ')130 ENDIF131 132 126 ! conservation test 133 IF( ln_limdia hsb) CALL lim_cons_hsm(0, 'limitd_me', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b)127 IF( ln_limdiachk ) CALL lim_cons_hsm(0, 'limitd_me', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 134 128 135 129 !-----------------------------------------------------------------------------! … … 289 283 ! control prints 290 284 !-----------------------------------------------------------------------------! 291 IF(ln_ctl) THEN292 CALL lim_var_glo2eqv293 294 CALL prt_ctl_info(' ')295 CALL prt_ctl_info(' - Cell values : ')296 CALL prt_ctl_info(' ~~~~~~~~~~~~~ ')297 CALL prt_ctl(tab2d_1=e12t , clinfo1=' lim_itd_me : cell area :')298 CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_itd_me : at_i :')299 CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_itd_me : vt_i :')300 CALL prt_ctl(tab2d_1=vt_s , clinfo1=' lim_itd_me : vt_s :')301 DO jl = 1, jpl302 CALL prt_ctl_info(' ')303 CALL prt_ctl_info(' - Category : ', ivar1=jl)304 CALL prt_ctl_info(' ~~~~~~~~~~')305 CALL prt_ctl(tab2d_1=a_i (:,:,jl) , clinfo1= ' lim_itd_me : a_i : ')306 CALL prt_ctl(tab2d_1=ht_i (:,:,jl) , clinfo1= ' lim_itd_me : ht_i : ')307 CALL prt_ctl(tab2d_1=ht_s (:,:,jl) , clinfo1= ' lim_itd_me : ht_s : ')308 CALL prt_ctl(tab2d_1=v_i (:,:,jl) , clinfo1= ' lim_itd_me : v_i : ')309 CALL prt_ctl(tab2d_1=v_s (:,:,jl) , clinfo1= ' lim_itd_me : v_s : ')310 CALL prt_ctl(tab2d_1=e_s (:,:,1,jl) , clinfo1= ' lim_itd_me : e_s : ')311 CALL prt_ctl(tab2d_1=t_su (:,:,jl) , clinfo1= ' lim_itd_me : t_su : ')312 CALL prt_ctl(tab2d_1=t_s (:,:,1,jl) , clinfo1= ' lim_itd_me : t_snow : ')313 CALL prt_ctl(tab2d_1=sm_i (:,:,jl) , clinfo1= ' lim_itd_me : sm_i : ')314 CALL prt_ctl(tab2d_1=smv_i (:,:,jl) , clinfo1= ' lim_itd_me : smv_i : ')315 DO jk = 1, nlay_i316 CALL prt_ctl_info(' ')317 CALL prt_ctl_info(' - Layer : ', ivar1=jk)318 CALL prt_ctl_info(' ~~~~~~~')319 CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_itd_me : t_i : ')320 CALL prt_ctl(tab2d_1=e_i(:,:,jk,jl) , clinfo1= ' lim_itd_me : e_i : ')321 END DO322 END DO323 ENDIF324 325 285 ! conservation test 326 IF( ln_limdiahsb ) CALL lim_cons_hsm(1, 'limitd_me', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 286 IF( ln_limdiachk ) CALL lim_cons_hsm(1, 'limitd_me', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 287 288 ! control prints 289 IF( ln_ctl ) CALL lim_prt3D( 'limitd_me' ) 327 290 328 291 CALL wrk_dealloc( jpi, jpj, closing_net, divu_adv, opning, closing_gross ) -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limrhg.F90
r6989 r6994 38 38 USE bdyice_lim 39 39 #endif 40 !!!clem41 !! USE iom42 40 43 41 IMPLICIT NONE … … 703 701 ENDIF 704 702 ENDIF 705 ! 706 707 !! ! velocity 708 !! IF ( iom_use( "uice_ipa" ) .OR. iom_use( "vice_ipa" ) .OR. iom_use( "icevel" ) ) THEN 709 !! CALL iom_put( "uice_ipa" , u_ice ) ! ice velocity u component 710 !! CALL iom_put( "vice_ipa" , v_ice ) ! ice velocity v component 711 !! ENDIF 712 !! IF ( iom_use( "iceconc_cat" ) ) CALL iom_put( "iceconc_cat" , a_i ) ! area for categories 713 !! IF ( iom_use( "icethic_cat" ) ) CALL iom_put( "icethic_cat" , ht_i ) ! thickness for categories 714 703 ! 715 704 716 705 CALL wrk_dealloc( jpi,jpj, z1_e1t0, z1_e2t0, zp_delt ) -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limrst.F90
r6989 r6994 93 93 ENDIF 94 94 ! 95 IF( ln_ icectl ) CALL lim_prt( kt, iiceprt, jiceprt, 1, ' - Beginning the time step - ' ) ! control print95 IF( ln_limctl ) CALL lim_prt( kt, iiceprt, jiceprt, 1, ' - Beginning the time step - ' ) ! control print 96 96 END SUBROUTINE lim_rst_opn 97 97 -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90
r6989 r6994 36 36 USE wrk_nemo ! work arrays 37 37 USE in_out_manager ! I/O manager 38 USE prtctl ! Print control39 38 USE lib_fortran ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) 40 39 USE traqsr ! add penetration of solar flux in the calculation of heat budget … … 215 214 216 215 ! conservation test 217 IF( ln_limdia hsb )CALL lim_cons_final( 'limsbc' )216 IF( ln_limdiachk ) CALL lim_cons_final( 'limsbc' ) 218 217 219 218 ! control prints 220 IF( ln_icectl ) CALL lim_prt( kt, iiceprt, jiceprt, 3, ' - Final state lim_sbc - ' ) 221 222 IF(ln_ctl) THEN 223 CALL prt_ctl( tab2d_1=qsr , clinfo1=' lim_sbc: qsr : ', tab2d_2=qns , clinfo2=' qns : ' ) 224 CALL prt_ctl( tab2d_1=emp , clinfo1=' lim_sbc: emp : ', tab2d_2=sfx , clinfo2=' sfx : ' ) 225 CALL prt_ctl( tab2d_1=fr_i , clinfo1=' lim_sbc: fr_i : ' ) 226 CALL prt_ctl( tab3d_1=tn_ice, clinfo1=' lim_sbc: tn_ice : ', kdim=jpl ) 227 ENDIF 219 IF( ln_limctl ) CALL lim_prt( kt, iiceprt, jiceprt, 3, ' - Final state lim_sbc - ' ) 220 IF( ln_ctl ) CALL lim_prt3D( 'limsbc' ) 228 221 229 222 END SUBROUTINE lim_sbc_flx … … 300 293 CALL lbc_lnk_multi( utau, 'U', -1., vtau, 'V', -1. ) ! lateral boundary condition 301 294 ! 302 IF(ln_ctl) CALL prt_ctl( tab2d_1=utau, clinfo1=' lim_sbc: utau : ', mask1=umask, &303 & tab2d_2=vtau, clinfo2=' vtau : ' , mask2=vmask )304 295 ! 305 296 END SUBROUTINE lim_sbc_tau -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90
r6746 r6994 39 39 USE wrk_nemo ! work arrays 40 40 USE in_out_manager ! I/O manager 41 USE prtctl ! Print control42 41 USE lib_fortran ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) 43 42 USE timing ! Timing … … 104 103 105 104 ! conservation test 106 IF( ln_limdia hsb) CALL lim_cons_hsm(0, 'limthd', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b)105 IF( ln_limdiachk ) CALL lim_cons_hsm(0, 'limthd', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 107 106 108 107 CALL lim_var_glo2eqv … … 244 243 ! debug point to follow 245 244 jiindex_1d = 0 246 IF( ln_ icectl ) THEN245 IF( ln_limctl ) THEN 247 246 DO ji = mi0(iiceprt), mi1(iiceprt) 248 247 DO jj = mj0(jiceprt), mj1(jiceprt) … … 316 315 317 316 ! control checks 318 IF( ln_ icectl )CALL lim_prt( kt, iiceprt, jiceprt, 1, ' - ice thermodyn. - ' ) ! control print319 ! 320 IF( ln_limdia hsb) CALL lim_cons_hsm(1, 'limthd', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b)317 IF( ln_limctl ) CALL lim_prt( kt, iiceprt, jiceprt, 1, ' - ice thermodyn. - ' ) ! control print 318 ! 319 IF( ln_limdiachk ) CALL lim_cons_hsm(1, 'limthd', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 321 320 322 321 !------------------------------------------------! … … 324 323 !------------------------------------------------! 325 324 ! Given thermodynamic growth rates, transport ice between thickness categories. 326 IF( ln_limdia hsb) CALL lim_cons_hsm(0, 'limitd_th_rem', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b)325 IF( ln_limdiachk ) CALL lim_cons_hsm(0, 'limitd_th_rem', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 327 326 328 327 IF( jpl > 1 ) CALL lim_itd_th_rem( 1, jpl, kt ) 329 328 330 IF( ln_limdia hsb) CALL lim_cons_hsm(1, 'limitd_th_rem', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b)329 IF( ln_limdiachk ) CALL lim_cons_hsm(1, 'limitd_th_rem', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 331 330 332 331 !------------------------------------------------! 333 332 ! Add frazil ice growing in leads 334 333 !------------------------------------------------! 335 IF( ln_limdia hsb) CALL lim_cons_hsm(0, 'limthd_lac', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b)334 IF( ln_limdiachk ) CALL lim_cons_hsm(0, 'limthd_lac', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 336 335 337 336 IF( ln_limdO ) CALL lim_thd_lac 338 337 339 IF( ln_limdiahsb ) CALL lim_cons_hsm(1, 'limthd_lac', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 338 ! conservation test 339 IF( ln_limdiachk ) CALL lim_cons_hsm(1, 'limthd_lac', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 340 340 341 341 ! Control print 342 IF(ln_ctl) THEN 343 CALL lim_var_glo2eqv 344 345 CALL prt_ctl_info(' ') 346 CALL prt_ctl_info(' - Cell values : ') 347 CALL prt_ctl_info(' ~~~~~~~~~~~~~ ') 348 CALL prt_ctl(tab2d_1=e12t , clinfo1=' lim_itd_th : cell area :') 349 CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_itd_th : at_i :') 350 CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_itd_th : vt_i :') 351 CALL prt_ctl(tab2d_1=vt_s , clinfo1=' lim_itd_th : vt_s :') 352 DO jl = 1, jpl 353 CALL prt_ctl_info(' ') 354 CALL prt_ctl_info(' - Category : ', ivar1=jl) 355 CALL prt_ctl_info(' ~~~~~~~~~~') 356 CALL prt_ctl(tab2d_1=a_i (:,:,jl) , clinfo1= ' lim_itd_th : a_i : ') 357 CALL prt_ctl(tab2d_1=ht_i (:,:,jl) , clinfo1= ' lim_itd_th : ht_i : ') 358 CALL prt_ctl(tab2d_1=ht_s (:,:,jl) , clinfo1= ' lim_itd_th : ht_s : ') 359 CALL prt_ctl(tab2d_1=v_i (:,:,jl) , clinfo1= ' lim_itd_th : v_i : ') 360 CALL prt_ctl(tab2d_1=v_s (:,:,jl) , clinfo1= ' lim_itd_th : v_s : ') 361 CALL prt_ctl(tab2d_1=e_s (:,:,1,jl) , clinfo1= ' lim_itd_th : e_s : ') 362 CALL prt_ctl(tab2d_1=t_su (:,:,jl) , clinfo1= ' lim_itd_th : t_su : ') 363 CALL prt_ctl(tab2d_1=t_s (:,:,1,jl) , clinfo1= ' lim_itd_th : t_snow : ') 364 CALL prt_ctl(tab2d_1=sm_i (:,:,jl) , clinfo1= ' lim_itd_th : sm_i : ') 365 CALL prt_ctl(tab2d_1=smv_i (:,:,jl) , clinfo1= ' lim_itd_th : smv_i : ') 366 DO jk = 1, nlay_i 367 CALL prt_ctl_info(' ') 368 CALL prt_ctl_info(' - Layer : ', ivar1=jk) 369 CALL prt_ctl_info(' ~~~~~~~') 370 CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_itd_th : t_i : ') 371 CALL prt_ctl(tab2d_1=e_i(:,:,jk,jl) , clinfo1= ' lim_itd_th : e_i : ') 372 END DO 373 END DO 374 ENDIF 342 IF( ln_ctl ) CALL lim_prt3D( 'limthd' ) 375 343 ! 376 344 CALL wrk_dealloc( jpi,jpj, zu_io, zv_io, zfric ) -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90
r6515 r6994 613 613 hfx_out(ii,ij) = hfx_out(ii,ij) + ( zq_rema(ji) * a_i_1d(ji) ) * r1_rdtice 614 614 615 IF( ln_ icectl .AND. zq_rema(ji) < 0. .AND. lwp ) WRITE(numout,*) 'ALERTE zq_rema <0 = ', zq_rema(ji)615 IF( ln_limctl .AND. zq_rema(ji) < 0. .AND. lwp ) WRITE(numout,*) 'ALERTE zq_rema <0 = ', zq_rema(ji) 616 616 END DO 617 617 -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90
r5512 r6994 734 734 END DO ! End of the do while iterative procedure 735 735 736 IF( ln_ icectl .AND. lwp ) THEN736 IF( ln_limctl .AND. lwp ) THEN 737 737 WRITE(numout,*) ' zerritmax : ', zerritmax 738 738 WRITE(numout,*) ' nconv : ', nconv -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd_lac.F90
r6746 r6994 254 254 ! debug point to follow 255 255 jiindex_1d = 0 256 IF( ln_ icectl ) THEN256 IF( ln_limctl ) THEN 257 257 DO ji = mi0(iiceprt), mi1(iiceprt) 258 258 DO jj = mj0(jiceprt), mj1(jiceprt) … … 264 264 ENDIF 265 265 266 IF( ln_ icectl ) WRITE(numout,*) 'lim_thd_lac : nbpac = ', nbpac266 IF( ln_limctl ) WRITE(numout,*) 'lim_thd_lac : nbpac = ', nbpac 267 267 268 268 !------------------------------ -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limtrp.F90
r6989 r6994 109 109 110 110 ! conservation test 111 IF( ln_limdia hsb) CALL lim_cons_hsm(0, 'limtrp', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b)111 IF( ln_limdiachk ) CALL lim_cons_hsm(0, 'limtrp', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 112 112 113 113 ! store old values for diag … … 517 517 518 518 ! conservation test 519 IF( ln_limdia hsb) CALL lim_cons_hsm(1, 'limtrp', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b)519 IF( ln_limdiachk ) CALL lim_cons_hsm(1, 'limtrp', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 520 520 521 521 ! ------------------------------------------------- 522 522 ! control prints 523 523 ! ------------------------------------------------- 524 IF( ln_ icectl ) CALL lim_prt( kt, iiceprt, jiceprt,-1, ' - ice dyn & trp - ' )524 IF( ln_limctl ) CALL lim_prt( kt, iiceprt, jiceprt,-1, ' - ice dyn & trp - ' ) 525 525 ! 526 526 CALL wrk_dealloc( jpi,jpj, zatold, zeiold, zesold, zsmvold ) -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limupdate1.F90
r6746 r6994 21 21 USE limitd_th 22 22 USE limvar 23 USE prtctl ! Print control24 23 USE wrk_nemo ! work arrays 25 24 USE timing ! Timing 26 25 USE limcons ! conservation tests 26 USE limctl ! control prints 27 27 USE lib_mpp ! MPP library 28 28 USE lib_fortran ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) … … 65 65 66 66 ! conservation test 67 IF( ln_limdia hsb) CALL lim_cons_hsm(0, 'limupdate1', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b)67 IF( ln_limdiachk ) CALL lim_cons_hsm(0, 'limupdate1', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 68 68 69 69 !---------------------------------------------------- … … 135 135 136 136 ! conservation test 137 IF( ln_limdia hsb) CALL lim_cons_hsm(1, 'limupdate1', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b)137 IF( ln_limdiachk ) CALL lim_cons_hsm(1, 'limupdate1', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 138 138 139 ! -------------------------------------------------140 139 ! control prints 141 ! ------------------------------------------------- 142 IF(ln_ctl) THEN ! Control print 143 CALL prt_ctl_info(' ') 144 CALL prt_ctl_info(' - Cell values : ') 145 CALL prt_ctl_info(' ~~~~~~~~~~~~~ ') 146 CALL prt_ctl(tab2d_1=e12t , clinfo1=' lim_update1 : cell area :') 147 CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_update1 : at_i :') 148 CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_update1 : vt_i :') 149 CALL prt_ctl(tab2d_1=vt_s , clinfo1=' lim_update1 : vt_s :') 150 CALL prt_ctl(tab2d_1=strength , clinfo1=' lim_update1 : strength :') 151 CALL prt_ctl(tab2d_1=u_ice , clinfo1=' lim_update1 : u_ice :', tab2d_2=v_ice , clinfo2=' v_ice :') 152 CALL prt_ctl(tab2d_1=u_ice_b , clinfo1=' lim_update1 : u_ice_b :', tab2d_2=v_ice_b , clinfo2=' v_ice_b :') 153 154 DO jl = 1, jpl 155 CALL prt_ctl_info(' ') 156 CALL prt_ctl_info(' - Category : ', ivar1=jl) 157 CALL prt_ctl_info(' ~~~~~~~~~~') 158 CALL prt_ctl(tab2d_1=ht_i (:,:,jl) , clinfo1= ' lim_update1 : ht_i : ') 159 CALL prt_ctl(tab2d_1=ht_s (:,:,jl) , clinfo1= ' lim_update1 : ht_s : ') 160 CALL prt_ctl(tab2d_1=t_su (:,:,jl) , clinfo1= ' lim_update1 : t_su : ') 161 CALL prt_ctl(tab2d_1=t_s (:,:,1,jl) , clinfo1= ' lim_update1 : t_snow : ') 162 CALL prt_ctl(tab2d_1=sm_i (:,:,jl) , clinfo1= ' lim_update1 : sm_i : ') 163 CALL prt_ctl(tab2d_1=o_i (:,:,jl) , clinfo1= ' lim_update1 : o_i : ') 164 CALL prt_ctl(tab2d_1=a_i (:,:,jl) , clinfo1= ' lim_update1 : a_i : ') 165 CALL prt_ctl(tab2d_1=a_i_b (:,:,jl) , clinfo1= ' lim_update1 : a_i_b : ') 166 CALL prt_ctl(tab2d_1=v_i (:,:,jl) , clinfo1= ' lim_update1 : v_i : ') 167 CALL prt_ctl(tab2d_1=v_i_b (:,:,jl) , clinfo1= ' lim_update1 : v_i_b : ') 168 CALL prt_ctl(tab2d_1=v_s (:,:,jl) , clinfo1= ' lim_update1 : v_s : ') 169 CALL prt_ctl(tab2d_1=v_s_b (:,:,jl) , clinfo1= ' lim_update1 : v_s_b : ') 170 CALL prt_ctl(tab2d_1=e_i (:,:,1,jl) , clinfo1= ' lim_update1 : e_i1 : ') 171 CALL prt_ctl(tab2d_1=e_i_b (:,:,1,jl) , clinfo1= ' lim_update1 : e_i1_b : ') 172 CALL prt_ctl(tab2d_1=e_i (:,:,2,jl) , clinfo1= ' lim_update1 : e_i2 : ') 173 CALL prt_ctl(tab2d_1=e_i_b (:,:,2,jl) , clinfo1= ' lim_update1 : e_i2_b : ') 174 CALL prt_ctl(tab2d_1=e_s (:,:,1,jl) , clinfo1= ' lim_update1 : e_snow : ') 175 CALL prt_ctl(tab2d_1=e_s_b (:,:,1,jl) , clinfo1= ' lim_update1 : e_snow_b : ') 176 CALL prt_ctl(tab2d_1=smv_i (:,:,jl) , clinfo1= ' lim_update1 : smv_i : ') 177 CALL prt_ctl(tab2d_1=smv_i_b (:,:,jl) , clinfo1= ' lim_update1 : smv_i_b : ') 178 CALL prt_ctl(tab2d_1=oa_i (:,:,jl) , clinfo1= ' lim_update1 : oa_i : ') 179 CALL prt_ctl(tab2d_1=oa_i_b (:,:,jl) , clinfo1= ' lim_update1 : oa_i_b : ') 180 181 DO jk = 1, nlay_i 182 CALL prt_ctl_info(' - Layer : ', ivar1=jk) 183 CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_update1 : t_i : ') 184 END DO 185 END DO 186 187 CALL prt_ctl_info(' ') 188 CALL prt_ctl_info(' - Heat / FW fluxes : ') 189 CALL prt_ctl_info(' ~~~~~~~~~~~~~~~~~~ ') 190 CALL prt_ctl(tab2d_1=sst_m , clinfo1= ' lim_update1 : sst : ', tab2d_2=sss_m , clinfo2= ' sss : ') 191 192 CALL prt_ctl_info(' ') 193 CALL prt_ctl_info(' - Stresses : ') 194 CALL prt_ctl_info(' ~~~~~~~~~~ ') 195 CALL prt_ctl(tab2d_1=utau , clinfo1= ' lim_update1 : utau : ', tab2d_2=vtau , clinfo2= ' vtau : ') 196 CALL prt_ctl(tab2d_1=utau_ice , clinfo1= ' lim_update1 : utau_ice : ', tab2d_2=vtau_ice , clinfo2= ' vtau_ice : ') 197 CALL prt_ctl(tab2d_1=u_oce , clinfo1= ' lim_update1 : u_oce : ', tab2d_2=v_oce , clinfo2= ' v_oce : ') 198 ENDIF 140 IF( ln_ctl ) CALL lim_prt3D( 'limupdate1' ) 199 141 200 142 IF( nn_timing == 1 ) CALL timing_stop('limupdate1') 143 201 144 END SUBROUTINE lim_update1 145 202 146 #else 203 147 !!---------------------------------------------------------------------- … … 207 151 SUBROUTINE lim_update1 ! Empty routine 208 152 END SUBROUTINE lim_update1 209 210 153 #endif 211 154 -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limupdate2.F90
r6746 r6994 21 21 USE limitd_th 22 22 USE limvar 23 USE prtctl ! Print control24 23 USE lbclnk ! lateral boundary condition - MPP exchanges 25 24 USE wrk_nemo ! work arrays … … 67 66 68 67 ! conservation test 69 IF( ln_limdia hsb) CALL lim_cons_hsm(0, 'limupdate2', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b)68 IF( ln_limdiachk ) CALL lim_cons_hsm(0, 'limupdate2', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 70 69 71 70 !---------------------------------------------------------------------- … … 176 175 177 176 ! conservation test 178 IF( ln_limdia hsb) CALL lim_cons_hsm(1, 'limupdate2', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b)177 IF( ln_limdiachk ) CALL lim_cons_hsm(1, 'limupdate2', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 179 178 180 ! -------------------------------------------------181 179 ! control prints 182 ! ------------------------------------------------- 183 IF( ln_icectl ) CALL lim_prt( kt, iiceprt, jiceprt, 2, ' - Final state - ' ) ! control print 184 185 IF(ln_ctl) THEN ! Control print 186 CALL prt_ctl_info(' ') 187 CALL prt_ctl_info(' - Cell values : ') 188 CALL prt_ctl_info(' ~~~~~~~~~~~~~ ') 189 CALL prt_ctl(tab2d_1=e12t , clinfo1=' lim_update2 : cell area :') 190 CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_update2 : at_i :') 191 CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_update2 : vt_i :') 192 CALL prt_ctl(tab2d_1=vt_s , clinfo1=' lim_update2 : vt_s :') 193 CALL prt_ctl(tab2d_1=strength , clinfo1=' lim_update2 : strength :') 194 CALL prt_ctl(tab2d_1=u_ice , clinfo1=' lim_update2 : u_ice :', tab2d_2=v_ice , clinfo2=' v_ice :') 195 CALL prt_ctl(tab2d_1=u_ice_b , clinfo1=' lim_update2 : u_ice_b :', tab2d_2=v_ice_b , clinfo2=' v_ice_b :') 196 197 DO jl = 1, jpl 198 CALL prt_ctl_info(' ') 199 CALL prt_ctl_info(' - Category : ', ivar1=jl) 200 CALL prt_ctl_info(' ~~~~~~~~~~') 201 CALL prt_ctl(tab2d_1=ht_i (:,:,jl) , clinfo1= ' lim_update2 : ht_i : ') 202 CALL prt_ctl(tab2d_1=ht_s (:,:,jl) , clinfo1= ' lim_update2 : ht_s : ') 203 CALL prt_ctl(tab2d_1=t_su (:,:,jl) , clinfo1= ' lim_update2 : t_su : ') 204 CALL prt_ctl(tab2d_1=t_s (:,:,1,jl) , clinfo1= ' lim_update2 : t_snow : ') 205 CALL prt_ctl(tab2d_1=sm_i (:,:,jl) , clinfo1= ' lim_update2 : sm_i : ') 206 CALL prt_ctl(tab2d_1=o_i (:,:,jl) , clinfo1= ' lim_update2 : o_i : ') 207 CALL prt_ctl(tab2d_1=a_i (:,:,jl) , clinfo1= ' lim_update2 : a_i : ') 208 CALL prt_ctl(tab2d_1=a_i_b (:,:,jl) , clinfo1= ' lim_update2 : a_i_b : ') 209 CALL prt_ctl(tab2d_1=v_i (:,:,jl) , clinfo1= ' lim_update2 : v_i : ') 210 CALL prt_ctl(tab2d_1=v_i_b (:,:,jl) , clinfo1= ' lim_update2 : v_i_b : ') 211 CALL prt_ctl(tab2d_1=v_s (:,:,jl) , clinfo1= ' lim_update2 : v_s : ') 212 CALL prt_ctl(tab2d_1=v_s_b (:,:,jl) , clinfo1= ' lim_update2 : v_s_b : ') 213 CALL prt_ctl(tab2d_1=e_i (:,:,1,jl) , clinfo1= ' lim_update2 : e_i1 : ') 214 CALL prt_ctl(tab2d_1=e_i_b (:,:,1,jl) , clinfo1= ' lim_update2 : e_i1_b : ') 215 CALL prt_ctl(tab2d_1=e_i (:,:,2,jl) , clinfo1= ' lim_update2 : e_i2 : ') 216 CALL prt_ctl(tab2d_1=e_i_b (:,:,2,jl) , clinfo1= ' lim_update2 : e_i2_b : ') 217 CALL prt_ctl(tab2d_1=e_s (:,:,1,jl) , clinfo1= ' lim_update2 : e_snow : ') 218 CALL prt_ctl(tab2d_1=e_s_b (:,:,1,jl) , clinfo1= ' lim_update2 : e_snow_b : ') 219 CALL prt_ctl(tab2d_1=smv_i (:,:,jl) , clinfo1= ' lim_update2 : smv_i : ') 220 CALL prt_ctl(tab2d_1=smv_i_b (:,:,jl) , clinfo1= ' lim_update2 : smv_i_b : ') 221 CALL prt_ctl(tab2d_1=oa_i (:,:,jl) , clinfo1= ' lim_update2 : oa_i : ') 222 CALL prt_ctl(tab2d_1=oa_i_b (:,:,jl) , clinfo1= ' lim_update2 : oa_i_b : ') 223 224 DO jk = 1, nlay_i 225 CALL prt_ctl_info(' - Layer : ', ivar1=jk) 226 CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_update2 : t_i : ') 227 END DO 228 END DO 229 230 CALL prt_ctl_info(' ') 231 CALL prt_ctl_info(' - Heat / FW fluxes : ') 232 CALL prt_ctl_info(' ~~~~~~~~~~~~~~~~~~ ') 233 CALL prt_ctl(tab2d_1=sst_m , clinfo1= ' lim_update2 : sst : ', tab2d_2=sss_m , clinfo2= ' sss : ') 234 235 CALL prt_ctl_info(' ') 236 CALL prt_ctl_info(' - Stresses : ') 237 CALL prt_ctl_info(' ~~~~~~~~~~ ') 238 CALL prt_ctl(tab2d_1=utau , clinfo1= ' lim_update2 : utau : ', tab2d_2=vtau , clinfo2= ' vtau : ') 239 CALL prt_ctl(tab2d_1=utau_ice , clinfo1= ' lim_update2 : utau_ice : ', tab2d_2=vtau_ice , clinfo2= ' vtau_ice : ') 240 CALL prt_ctl(tab2d_1=u_oce , clinfo1= ' lim_update2 : u_oce : ', tab2d_2=v_oce , clinfo2= ' v_oce : ') 241 ENDIF 180 IF( ln_limctl ) CALL lim_prt( kt, iiceprt, jiceprt, 2, ' - Final state - ' ) 181 IF( ln_ctl ) CALL lim_prt3D( 'limupdate2' ) 242 182 243 183 IF( nn_timing == 1 ) CALL timing_stop('limupdate2') 244 184 245 185 END SUBROUTINE lim_update2 186 246 187 #else 247 188 !!---------------------------------------------------------------------- … … 251 192 SUBROUTINE lim_update2 ! Empty routine 252 193 END SUBROUTINE lim_update2 253 254 194 #endif 255 195 -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/BDY/bdyice_lim.F90
r6746 r6994 28 28 USE ice ! LIM_3 ice variables 29 29 USE limvar 30 USE limctl 30 31 #endif 31 32 USE par_oce ! ocean parameters … … 79 80 80 81 #if defined key_lim3 81 CALL lim_var_zapsmall 82 CALL lim_var_agg(1) 82 CALL lim_var_zapsmall 83 CALL lim_var_agg(1) 84 IF( ln_limctl ) CALL lim_prt( kt, iiceprt, jiceprt, 1, ' - ice thermo bdy - ' ) 83 85 #endif 84 86 -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/DIA/diahsb.F90
r6970 r6994 240 240 SUBROUTINE dia_hsb_rst( kt, cdrw ) 241 241 !!--------------------------------------------------------------------- 242 !! *** ROUTINE limdia_rst ***242 !! *** ROUTINE dia_hsb_rst *** 243 243 !! 244 244 !! ** Purpose : Read or write DIA file in restart file -
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90
r6989 r6994 195 195 #if defined key_bdy 196 196 IF( ln_limthd ) CALL bdy_ice_lim( kt ) ! -- bdy ice thermo 197 IF( ln_icectl ) CALL lim_prt( kt, iiceprt, jiceprt, 1, ' - ice thermo bdy - ' )198 197 #endif 199 198 … … 260 259 CALL lim_sbc_flx( kt ) ! -- Update surface ocean mass, heat and salt fluxes 261 260 ! 262 IF( ln_limdiaout)CALL lim_diahsb( kt ) ! -- Diagnostics and outputs261 IF( ln_limdiahsb ) CALL lim_diahsb( kt ) ! -- Diagnostics and outputs 263 262 ! 264 263 CALL lim_wri( 1 ) ! -- Ice outputs … … 269 268 IF( lrst_ice ) CALL lim_rst_write( kt ) ! -- Ice restart file 270 269 ! 271 IF( ln_ icectl ) CALL lim_ctl( kt ) ! alerts in case of model crash270 IF( ln_limctl ) CALL lim_ctl( kt ) ! alerts in case of model crash 272 271 ! 273 272 ENDIF ! End sea-ice time step only … … 346 345 CALL lim_sbc_init ! ice surface boundary condition 347 346 ! 348 IF( ln_limdia out) CALL lim_diahsb_init ! initialization for diags347 IF( ln_limdiahsb) CALL lim_diahsb_init ! initialization for diags 349 348 ! 350 349 fr_i(:,:) = at_i(:,:) ! initialisation of sea-ice fraction … … 386 385 NAMELIST/namicerun/ jpl, nlay_i, nlay_s, rn_amax_n, rn_amax_s, cn_icerst_in, cn_icerst_indir, & 387 386 & cn_icerst_out, cn_icerst_outdir, ln_limthd, ln_limdyn, nn_limdyn, rn_uice, rn_vice 388 NAMELIST/namicediag/ ln_limdia hsb, ln_limdiaout, ln_icectl, iiceprt, jiceprt387 NAMELIST/namicediag/ ln_limdiachk, ln_limdiahsb, ln_limctl, iiceprt, jiceprt 389 388 !!------------------------------------------------------------------- 390 389 ! … … 427 426 WRITE(numout,*) '...and ice diagnostics' 428 427 WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~' 429 WRITE(numout,*) ' Diagnose heat/mass/salt budget or not ln_limdiahsb = ', ln_limdiahsb430 WRITE(numout,*) ' Output heat/mass/salt budget or not ln_limdiaout = ', ln_limdiaout431 WRITE(numout,*) ' control prints in ocean.out for (i,j)=(iiceprt,jiceprt) = ', ln_ icectl432 WRITE(numout,*) ' i-index for control prints (ln_ icectl=true) = ', iiceprt433 WRITE(numout,*) ' j-index for control prints (ln_ icectl=true) = ', jiceprt428 WRITE(numout,*) ' Diagnose online heat/mass/salt budget ln_limdiachk = ', ln_limdiachk 429 WRITE(numout,*) ' Output heat/mass/salt budget ln_limdiahsb = ', ln_limdiahsb 430 WRITE(numout,*) ' control prints in ocean.out for (i,j)=(iiceprt,jiceprt) = ', ln_limctl 431 WRITE(numout,*) ' i-index for control prints (ln_limctl=true) = ', iiceprt 432 WRITE(numout,*) ' j-index for control prints (ln_limctl=true) = ', jiceprt 434 433 ENDIF 435 434 ! … … 443 442 ! 444 443 #if defined key_bdy 445 IF( lwp .AND. ln_limdia hsb) CALL ctl_warn('online conservation check activated but it does not work with BDY')444 IF( lwp .AND. ln_limdiachk ) CALL ctl_warn('online conservation check activated but it does not work with BDY') 446 445 #endif 447 446 !
Note: See TracChangeset
for help on using the changeset viewer.