Changeset 7646 for trunk/NEMOGCM/NEMO/LIM_SRC_3/limupdate1.F90
- Timestamp:
- 2017-02-06T10:25:03+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/LIM_SRC_3/limupdate1.F90
r6403 r7646 15 15 USE sbc_oce ! Surface boundary condition: ocean fields 16 16 USE sbc_ice ! Surface boundary condition: ice fields 17 USE dom_ice18 17 USE dom_oce 19 18 USE phycst ! physical constants … … 22 21 USE limitd_th 23 22 USE limvar 24 USE prtctl ! Print control25 23 USE wrk_nemo ! work arrays 26 24 USE timing ! Timing 27 25 USE limcons ! conservation tests 26 USE limctl ! control prints 28 27 USE lib_mpp ! MPP library 29 28 USE lib_fortran ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) … … 59 58 IF( nn_timing == 1 ) CALL timing_start('limupdate1') 60 59 61 IF( ln_limdyn ) THEN62 63 60 IF( kt == nit000 .AND. lwp ) THEN 64 WRITE(numout,*) ' lim_update1 ' 65 WRITE(numout,*) ' ~~~~~~~~~~~ ' 61 WRITE(numout,*)'' 62 WRITE(numout,*)' lim_update1 ' 63 WRITE(numout,*)' ~~~~~~~~~~~ ' 66 64 ENDIF 67 65 68 66 ! conservation test 69 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) 70 68 71 69 !---------------------------------------------------- … … 137 135 138 136 ! conservation test 139 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) 140 138 141 ! -------------------------------------------------142 139 ! control prints 143 ! ------------------------------------------------- 144 IF(ln_ctl) THEN ! Control print 145 CALL prt_ctl_info(' ') 146 CALL prt_ctl_info(' - Cell values : ') 147 CALL prt_ctl_info(' ~~~~~~~~~~~~~ ') 148 CALL prt_ctl(tab2d_1=e1e2t , clinfo1=' lim_update1 : cell area :') 149 CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_update1 : at_i :') 150 CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_update1 : vt_i :') 151 CALL prt_ctl(tab2d_1=vt_s , clinfo1=' lim_update1 : vt_s :') 152 CALL prt_ctl(tab2d_1=strength , clinfo1=' lim_update1 : strength :') 153 CALL prt_ctl(tab2d_1=u_ice , clinfo1=' lim_update1 : u_ice :', tab2d_2=v_ice , clinfo2=' v_ice :') 154 CALL prt_ctl(tab2d_1=u_ice_b , clinfo1=' lim_update1 : u_ice_b :', tab2d_2=v_ice_b , clinfo2=' v_ice_b :') 140 IF( ln_ctl ) CALL lim_prt3D( 'limupdate1' ) 141 142 IF( nn_timing == 1 ) CALL timing_stop('limupdate1') 155 143 156 DO jl = 1, jpl 157 CALL prt_ctl_info(' ') 158 CALL prt_ctl_info(' - Category : ', ivar1=jl) 159 CALL prt_ctl_info(' ~~~~~~~~~~') 160 CALL prt_ctl(tab2d_1=ht_i (:,:,jl) , clinfo1= ' lim_update1 : ht_i : ') 161 CALL prt_ctl(tab2d_1=ht_s (:,:,jl) , clinfo1= ' lim_update1 : ht_s : ') 162 CALL prt_ctl(tab2d_1=t_su (:,:,jl) , clinfo1= ' lim_update1 : t_su : ') 163 CALL prt_ctl(tab2d_1=t_s (:,:,1,jl) , clinfo1= ' lim_update1 : t_snow : ') 164 CALL prt_ctl(tab2d_1=sm_i (:,:,jl) , clinfo1= ' lim_update1 : sm_i : ') 165 CALL prt_ctl(tab2d_1=o_i (:,:,jl) , clinfo1= ' lim_update1 : o_i : ') 166 CALL prt_ctl(tab2d_1=a_i (:,:,jl) , clinfo1= ' lim_update1 : a_i : ') 167 CALL prt_ctl(tab2d_1=a_i_b (:,:,jl) , clinfo1= ' lim_update1 : a_i_b : ') 168 CALL prt_ctl(tab2d_1=v_i (:,:,jl) , clinfo1= ' lim_update1 : v_i : ') 169 CALL prt_ctl(tab2d_1=v_i_b (:,:,jl) , clinfo1= ' lim_update1 : v_i_b : ') 170 CALL prt_ctl(tab2d_1=v_s (:,:,jl) , clinfo1= ' lim_update1 : v_s : ') 171 CALL prt_ctl(tab2d_1=v_s_b (:,:,jl) , clinfo1= ' lim_update1 : v_s_b : ') 172 CALL prt_ctl(tab2d_1=e_i (:,:,1,jl) , clinfo1= ' lim_update1 : e_i1 : ') 173 CALL prt_ctl(tab2d_1=e_i_b (:,:,1,jl) , clinfo1= ' lim_update1 : e_i1_b : ') 174 CALL prt_ctl(tab2d_1=e_i (:,:,2,jl) , clinfo1= ' lim_update1 : e_i2 : ') 175 CALL prt_ctl(tab2d_1=e_i_b (:,:,2,jl) , clinfo1= ' lim_update1 : e_i2_b : ') 176 CALL prt_ctl(tab2d_1=e_s (:,:,1,jl) , clinfo1= ' lim_update1 : e_snow : ') 177 CALL prt_ctl(tab2d_1=e_s_b (:,:,1,jl) , clinfo1= ' lim_update1 : e_snow_b : ') 178 CALL prt_ctl(tab2d_1=smv_i (:,:,jl) , clinfo1= ' lim_update1 : smv_i : ') 179 CALL prt_ctl(tab2d_1=smv_i_b (:,:,jl) , clinfo1= ' lim_update1 : smv_i_b : ') 180 CALL prt_ctl(tab2d_1=oa_i (:,:,jl) , clinfo1= ' lim_update1 : oa_i : ') 181 CALL prt_ctl(tab2d_1=oa_i_b (:,:,jl) , clinfo1= ' lim_update1 : oa_i_b : ') 144 END SUBROUTINE lim_update1 182 145 183 DO jk = 1, nlay_i184 CALL prt_ctl_info(' - Layer : ', ivar1=jk)185 CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_update1 : t_i : ')186 END DO187 END DO188 189 CALL prt_ctl_info(' ')190 CALL prt_ctl_info(' - Heat / FW fluxes : ')191 CALL prt_ctl_info(' ~~~~~~~~~~~~~~~~~~ ')192 CALL prt_ctl(tab2d_1=sst_m , clinfo1= ' lim_update1 : sst : ', tab2d_2=sss_m , clinfo2= ' sss : ')193 194 CALL prt_ctl_info(' ')195 CALL prt_ctl_info(' - Stresses : ')196 CALL prt_ctl_info(' ~~~~~~~~~~ ')197 CALL prt_ctl(tab2d_1=utau , clinfo1= ' lim_update1 : utau : ', tab2d_2=vtau , clinfo2= ' vtau : ')198 CALL prt_ctl(tab2d_1=utau_ice , clinfo1= ' lim_update1 : utau_ice : ', tab2d_2=vtau_ice , clinfo2= ' vtau_ice : ')199 CALL prt_ctl(tab2d_1=u_oce , clinfo1= ' lim_update1 : u_oce : ', tab2d_2=v_oce , clinfo2= ' v_oce : ')200 ENDIF201 202 ENDIF ! ln_limdyn203 204 IF( nn_timing == 1 ) CALL timing_stop('limupdate1')205 END SUBROUTINE lim_update1206 146 #else 207 147 !!----------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.