Changeset 106 for trunk/NEMO
- Timestamp:
- 2004-06-28T10:08:06+02:00 (20 years ago)
- Location:
- trunk/NEMO
- Files:
-
- 51 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/LIM_SRC/icestp.F90
r88 r106 184 184 185 185 186 IF( l_ctl) THEN ! print mean trends (used for debugging)186 IF(l_ctl) THEN ! print mean trends (used for debugging) 187 187 WRITE(numout,*) 'Ice Forcings ' 188 WRITE(numout,*) ' qsr_oce : ', SUM( qsr_oce ( :,:) ), ' qsr_ice : ', SUM( qsr_ice (:,:) )189 WRITE(numout,*) ' qnsr_oce : ', SUM( qnsr_oce( :,:) ), ' qnsr_ice : ', SUM( qnsr_ice(:,:) )190 WRITE(numout,*) ' evap : ', SUM( evap (:,:) )191 WRITE(numout,*) ' precip : ', SUM( tprecip( :,:) ), ' Snow : ', SUM( sprecip (:,:) )192 WRITE(numout,*) ' u-stress : ', SUM( gtaux ( :,:) ), ' v-stress : ', SUM( gtauy (:,:) )193 WRITE(numout,*) ' sst : ', SUM( sst_io ( :,:) ), ' sss : ', SUM( sss_io (:,:) )194 WRITE(numout,*) ' u_io : ', SUM( u_io ( :,:) ), ' v_io : ', SUM( v_io (:,:) )195 WRITE(numout,*) ' hsnif 1 : ', SUM( hsnif ( :,:) ), ' hicnif : ', SUM( hicif (:,:) )196 WRITE(numout,*) ' frld 1 : ', SUM( frld ( :,:) ), ' sist : ', SUM( sist (:,:) )188 WRITE(numout,*) ' qsr_oce : ', SUM( qsr_oce (2:nictl,2:njctl) ), ' qsr_ice : ', SUM( qsr_ice (2:nictl,2:njctl) ) 189 WRITE(numout,*) ' qnsr_oce : ', SUM( qnsr_oce(2:nictl,2:njctl) ), ' qnsr_ice : ', SUM( qnsr_ice(2:nictl,2:njctl) ) 190 WRITE(numout,*) ' evap : ', SUM( evap (2:nictl,2:njctl) ) 191 WRITE(numout,*) ' precip : ', SUM( tprecip(2:nictl,2:njctl) ), ' Snow : ', SUM( sprecip (2:nictl,2:njctl) ) 192 WRITE(numout,*) ' u-stress : ', SUM( gtaux (2:nictl,2:njctl) ), ' v-stress : ', SUM( gtauy (2:nictl,2:njctl) ) 193 WRITE(numout,*) ' sst : ', SUM( sst_io (2:nictl,2:njctl) ), ' sss : ', SUM( sss_io (2:nictl,2:njctl) ) 194 WRITE(numout,*) ' u_io : ', SUM( u_io (2:nictl,2:njctl) ), ' v_io : ', SUM( v_io (2:nictl,2:njctl) ) 195 WRITE(numout,*) ' hsnif 1 : ', SUM( hsnif (2:nictl,2:njctl) ), ' hicnif : ', SUM( hicif (2:nictl,2:njctl) ) 196 WRITE(numout,*) ' frld 1 : ', SUM( frld (2:nictl,2:njctl) ), ' sist : ', SUM( sist (2:nictl,2:njctl) ) 197 197 ENDIF 198 198 … … 203 203 CALL lim_dyn ! Ice dynamics ! ( rheology/dynamics ) 204 204 ! !--------------! 205 IF( l_ctl) THEN206 WRITE(numout,*) ' hsnif 2 : ', SUM( hsnif ( :,:) ), ' hicnif : ', SUM( hicif (:,:) )207 WRITE(numout,*) ' frld 2 : ', SUM( frld ( :,:) ), ' sist : ', SUM( sist (:,:) )205 IF(l_ctl) THEN 206 WRITE(numout,*) ' hsnif 2 : ', SUM( hsnif (2:nictl,2:njctl) ), ' hicnif : ', SUM( hicif (2:nictl,2:njctl) ) 207 WRITE(numout,*) ' frld 2 : ', SUM( frld (2:nictl,2:njctl) ), ' sist : ', SUM( sist (2:nictl,2:njctl) ) 208 208 ENDIF 209 209 … … 212 212 CALL lim_trp ! Ice transport ! ( Advection/diffusion ) 213 213 ! !---------------! 214 IF( l_ctl) THEN215 WRITE(numout,*) ' hsnif 3 : ', SUM( hsnif ( :,:) ), ' hicnif : ', SUM( hicif (:,:) )216 WRITE(numout,*) ' frld 3 : ', SUM( frld ( :,:) ), ' sist : ', SUM( sist (:,:) )214 IF(l_ctl) THEN 215 WRITE(numout,*) ' hsnif 3 : ', SUM( hsnif (2:nictl,2:njctl) ), ' hicnif : ', SUM( hicif (2:nictl,2:njctl) ) 216 WRITE(numout,*) ' frld 3 : ', SUM( frld (2:nictl,2:njctl) ), ' sist : ', SUM( sist (2:nictl,2:njctl) ) 217 217 ENDIF 218 218 … … 221 221 CALL lim_thd ! Ice thermodynamics ! 222 222 ! !--------------------! 223 IF( l_ctl) THEN224 WRITE(numout,*) ' hsnif 4 : ', SUM( hsnif ( :,:) ), ' hicnif : ', SUM( hicif (:,:) )225 WRITE(numout,*) ' frld 4 : ', SUM( frld ( :,:) ), ' sist : ', SUM( sist (:,:) )223 IF(l_ctl) THEN 224 WRITE(numout,*) ' hsnif 4 : ', SUM( hsnif (2:nictl,2:njctl) ), ' hicnif : ', SUM( hicif (2:nictl,2:njctl) ) 225 WRITE(numout,*) ' frld 4 : ', SUM( frld (2:nictl,2:njctl) ), ' sist : ', SUM( sist (2:nictl,2:njctl) ) 226 226 ENDIF 227 227 -
trunk/NEMO/LIM_SRC/limadv.F90
r88 r106 223 223 224 224 IF(l_ctl) THEN 225 WRITE(numout,*) ' lim_adv_x: psm ', SUM( psm ), ' ps0 ', SUM( ps0)226 WRITE(numout,*) ' lim_adv_x: psx ', SUM( psx ), ' psxx ', SUM( psxx)227 WRITE(numout,*) ' lim_adv_x: psy ', SUM( psy ), ' psyy ', SUM( psyy)228 WRITE(numout,*) ' lim_adv_x: psxy ', SUM( psxy )225 WRITE(numout,*) ' lim_adv_x: psm ', SUM( psm (2:nictl,2:njctl) ), ' ps0 ', SUM( ps0 (2:nictl,2:njctl) ) 226 WRITE(numout,*) ' lim_adv_x: psx ', SUM( psx (2:nictl,2:njctl) ), ' psxx ', SUM( psxx(2:nictl,2:njctl) ) 227 WRITE(numout,*) ' lim_adv_x: psy ', SUM( psy (2:nictl,2:njctl) ), ' psyy ', SUM( psyy(2:nictl,2:njctl) ) 228 WRITE(numout,*) ' lim_adv_x: psxy ', SUM( psxy(2:nictl,2:njctl) ) 229 229 ENDIF 230 230 … … 423 423 424 424 IF(l_ctl) THEN 425 WRITE(numout,*) ' lim_adv_y: psm ', SUM( psm ), ' ps0 ', SUM( ps0)426 WRITE(numout,*) ' lim_adv_y: psx ', SUM( psx ), ' psxx ', SUM( psxx)427 WRITE(numout,*) ' lim_adv_y: psy ', SUM( psy ), ' psyy ', SUM( psyy)428 WRITE(numout,*) ' lim_adv_y: psxy ', SUM( psxy )425 WRITE(numout,*) ' lim_adv_y: psm ', SUM( psm (2:nictl,2:njctl) ), ' ps0 ', SUM( ps0 (2:nictl,2:njctl) ) 426 WRITE(numout,*) ' lim_adv_y: psx ', SUM( psx (2:nictl,2:njctl) ), ' psxx ', SUM( psxx(2:nictl,2:njctl) ) 427 WRITE(numout,*) ' lim_adv_y: psy ', SUM( psy (2:nictl,2:njctl) ), ' psyy ', SUM( psyy(2:nictl,2:njctl) ) 428 WRITE(numout,*) ' lim_adv_y: psxy ', SUM( psxy(2:nictl,2:njctl) ) 429 429 ENDIF 430 430 -
trunk/NEMO/LIM_SRC/limdyn.F90
r76 r106 90 90 i_j1 = 1 91 91 i_jpj = jpj 92 IF(l_ctl) WRITE(numout,*) 'lim_ rhg: i_j1 = ', i_j1, ' ij_jpj = ', i_jpj92 IF(l_ctl) WRITE(numout,*) 'lim_dyn : i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 93 93 CALL lim_rhg( i_j1, i_jpj ) 94 94 … … 98 98 zind(jj) = SUM( frld (:,jj ) ) ! = FLOAT(jpj) if ocean everywhere on a j-line 99 99 zmsk(jj) = SUM( tmask(:,jj,1) ) ! = 0 if land everywhere on a j-line 100 !!i write(numout,*) narea, 'lim rhg' , jj, zind(jj), zmsk(jj)100 !!i write(numout,*) narea, 'limdyn' , jj, zind(jj), zmsk(jj) 101 101 END DO 102 102 … … 111 111 END DO 112 112 i_j1 = MAX( 1, i_j1-1 ) 113 IF(l_ctl) WRITE(numout,*) 'lim_ rhg: NH i_j1 = ', i_j1, ' ij_jpj = ', i_jpj113 IF(l_ctl) WRITE(numout,*) 'lim_dyn : NH i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 114 114 115 115 CALL lim_rhg( i_j1, i_jpj ) … … 122 122 END DO 123 123 i_jpj = MIN( jpj, i_jpj+2 ) 124 IF(l_ctl) WRITE(numout,*) 'lim_ rhg: SH i_j1 = ', i_j1, ' ij_jpj = ', i_jpj124 IF(l_ctl) WRITE(numout,*) 'lim_dyn : SH i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 125 125 126 126 CALL lim_rhg( i_j1, i_jpj ) … … 141 141 i_jpj = MIN( jpj, i_jpj+2) 142 142 143 IF(l_ctl) WRITE(numout,*) 'lim_ rhg: one hemisphere: i_j1 = ', i_j1, ' ij_jpj = ', i_jpj143 IF(l_ctl) WRITE(numout,*) 'lim_dyn : one hemisphere: i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 144 144 145 145 CALL lim_rhg( i_j1, i_jpj ) … … 153 153 154 154 IF(l_ctl) THEN 155 WRITE(numout,*) ' lim_dyn : u_oce ', SUM( u_oce ), ' v_oce ', SUM( v_oce)156 WRITE(numout,*) ' lim_dyn : u_ice ', SUM( u_ice ), ' v_ice ', SUM( v_ice)155 WRITE(numout,*) ' lim_dyn : u_oce ', SUM( u_oce(2:nictl,2:njctl) ), ' v_oce ', SUM( v_oce(2:nictl,2:njctl) ) 156 WRITE(numout,*) ' lim_dyn : u_ice ', SUM( u_ice(2:nictl,2:njctl) ), ' v_ice ', SUM( v_ice(2:nictl,2:njctl) ) 157 157 ENDIF 158 158 … … 243 243 244 244 IF(l_ctl) THEN 245 WRITE(numout,*) ' lim_dyn : tio_u ', SUM( tio_u ), ' tio_v ', SUM( tio_v)246 WRITE(numout,*) ' lim_dyn : ust2s ', SUM( ust2s )245 WRITE(numout,*) ' lim_dyn : tio_u ', SUM( tio_u(2:nictl,2:njctl) ), ' tio_v ', SUM( tio_v(2:nictl,2:njctl) ) 246 WRITE(numout,*) ' lim_dyn : ust2s ', SUM( ust2s(2:nictl,2:njctl) ) 247 247 ENDIF 248 248 -
trunk/NEMO/LIM_SRC/limflx.F90
r88 r106 28 28 29 29 !! * Module variables 30 REAL(wp) :: 30 REAL(wp) :: & ! constant values 31 31 epsi16 = 1.e-16 , & 32 32 rzero = 0.e0 , & … … 194 194 IF(l_ctl) THEN 195 195 WRITE(numout,*) ' lim_flx ' 196 WRITE(numout,*) ' fsolar ', SUM( fsolar), ' fnsolar', SUM( fnsolar)197 WRITE(numout,*) ' fmass ', SUM( fmass ), ' fsalt ', SUM(fsalt)198 WRITE(numout,*) ' ftaux ', SUM( ftaux ), ' ftauy ', SUM(ftauy)199 WRITE(numout,*) ' freeze ', SUM( freeze), ' tn_ice ', SUM(tn_ice)196 WRITE(numout,*) ' fsolar ', SUM( fsolar(2:nictl,2:njctl) ), ' fnsolar', SUM( fnsolar(2:nictl,2:njctl) ) 197 WRITE(numout,*) ' fmass ', SUM( fmass (2:nictl,2:njctl) ), ' fsalt ', SUM( fsalt (2:nictl,2:njctl) ) 198 WRITE(numout,*) ' ftaux ', SUM( ftaux (2:nictl,2:njctl) ), ' ftauy ', SUM( ftauy (2:nictl,2:njctl) ) 199 WRITE(numout,*) ' freeze ', SUM( freeze(2:nictl,2:njctl) ), ' tn_ice ', SUM( tn_ice (2:nictl,2:njctl) ) 200 200 ENDIF 201 201 -
trunk/NEMO/LIM_SRC/limhdf.F90
r88 r106 167 167 ptab(:,:) = ptab(:,:) 168 168 169 IF(l_ctl) WRITE(numout,*) ' lim_hdf : ', SUM( ptab-ptab0 ), ' zconv= ', zconv, ' iter= ', iter 169 IF(l_ctl) WRITE(numout,*) ' lim_hdf : ', SUM( ptab(2:nictl,2:njctl)-ptab0(2:nictl,2:njctl) ), & 170 & ' zconv= ', zconv, ' iter= ', iter 170 171 171 172 END SUBROUTINE lim_hdf -
trunk/NEMO/LIM_SRC/limrhg.F90
r77 r106 27 27 28 28 !! * Module variables 29 REAL(wp) :: 29 REAL(wp) :: & ! constant values 30 30 rzero = 0.e0 , & 31 31 rone = 1.e0 … … 631 631 632 632 IF(l_ctl) THEN 633 WRITE(numout,*) ' lim_rhg : res= ', zresm, 'iter= ', jter,' u_ice ', SUM( u_ice ) , ' v_ice ', SUM( v_ice ) 633 WRITE(numout,*) ' lim_rhg : res= ', zresm, 'iter= ', jter,' u_ice ', SUM( u_ice(2:nictl,2:njctl) ) , & 634 & ' v_ice ', SUM( v_ice(2:nictl,2:njctl) ) 634 635 ENDIF 635 636 -
trunk/NEMO/LIM_SRC/limthd.F90
r70 r106 116 116 END DO 117 117 END DO 118 IF(l_ctl) WRITE(numout,*) 'lim_thd : ', SUM( hsnif( :,:) ) , ' hsnif'118 IF(l_ctl) WRITE(numout,*) 'lim_thd : ', SUM( hsnif(2:nictl,2:njctl) ) , ' hsnif' 119 119 120 120 … … 149 149 END DO 150 150 IF(l_ctl) THEN 151 WRITE(numout,*) 'lim_thd: hicif : ', SUM( hicif ), ' hsnif ', SUM( hsnif ) 152 WRITE(numout,*) 'lim_thd: dmgwi : ', SUM( dmgwi ), ' qstoif ', SUM( qstoif ) 153 WRITE(numout,*) 'lim_thd: frld : ', SUM( frld ) 151 WRITE(numout,*) 'lim_thd: hicif : ', SUM( hicif (2:nictl,2:njctl) ) 152 WRITE(numout,*) 'lim_thd: hsnif : ', SUM( hsnif (2:nictl,2:njctl) ) 153 WRITE(numout,*) 'lim_thd: dmgwi : ', SUM( dmgwi (2:nictl,2:njctl) ) 154 WRITE(numout,*) 'lim_thd: qstoif: ', SUM( qstoif(2:nictl,2:njctl) ) 155 WRITE(numout,*) 'lim_thd: frld : ', SUM( frld (2:nictl,2:njctl) ) 154 156 ENDIF 155 157 … … 218 220 END DO 219 221 IF(l_ctl) THEN 220 WRITE(numout,*) 'lim_thd: pfrld ' , SUM( pfrld ), ' thcm ', SUM( thcm)221 WRITE(numout,*) 'lim_thd: fdtcn ' , SUM( fdtcn ), ' qdtcn ', SUM( qdtcn)222 WRITE(numout,*) 'lim_thd: qldif ' , SUM( qldif ), ' zqlbsbq ', SUM( zqlbsbq)223 WRITE(numout,*) 'lim_thd: qcmif ' , SUM( qcmif ), ' fbif ', SUM( fbif)224 WRITE(numout,*) 'lim_thd: qcmif ', SUM( qcmif*tms)225 WRITE(numout,*) 'lim_thd: zhicifp', SUM( zhicifp )222 WRITE(numout,*) 'lim_thd: pfrld ' , SUM( pfrld (2:nictl,2:njctl) ), ' thcm ', SUM( thcm (2:nictl,2:njctl) ) 223 WRITE(numout,*) 'lim_thd: fdtcn ' , SUM( fdtcn (2:nictl,2:njctl) ), ' qdtcn ', SUM( qdtcn (2:nictl,2:njctl) ) 224 WRITE(numout,*) 'lim_thd: qldif ' , SUM( qldif (2:nictl,2:njctl) ), ' zqlbsbq ', SUM( zqlbsbq(2:nictl,2:njctl) ) 225 WRITE(numout,*) 'lim_thd: qcmif ' , SUM( qcmif (2:nictl,2:njctl) ), ' fbif ', SUM( fbif (2:nictl,2:njctl) ) 226 WRITE(numout,*) 'lim_thd: qcmif ', SUM( qcmif (2:nictl,2:njctl)*tms(2:nictl,2:njctl) ) 227 WRITE(numout,*) 'lim_thd: zhicifp', SUM( zhicifp(2:nictl,2:njctl) ) 226 228 WRITE(numout,*) 'limthd : nbpb = ', nbpb 227 229 ENDIF … … 320 322 321 323 IF(l_ctl) THEN 322 WRITE(numout,*) 'lim_thd : phicif ', SUM( phicif ), ' hicif ', SUM( hicif)324 WRITE(numout,*) 'lim_thd : phicif ', SUM( phicif(2:nictl,2:njctl) ), ' hicif ', SUM( hicif(2:nictl,2:njctl) ) 323 325 WRITE(numout,*) 'lim_thd : nbpac = ', nbpac 324 326 ENDIF … … 375 377 IF(l_ctl) THEN 376 378 WRITE(numout,*) ' lim_thd end ' 377 WRITE(numout,*) ' hicif ', SUM( hicif ), ' hsnif ', SUM( hsnif)378 WRITE(numout,*) ' frld ', SUM( frld ), ' hicifp', SUM( hicifp)379 WRITE(numout,*) ' phicif', SUM( phicif ), ' pfrld ', SUM( pfrld)380 WRITE(numout,*) ' sist ', SUM( sist ), ' tbif 1', SUM( tbif (:,:,1) )381 WRITE(numout,*) ' tbif 2', SUM( tbif (:,:,2) ), ' tbif 3', SUM( tbif (:,:,3) )382 WRITE(numout,*) ' fdtcn ', SUM( fdtcn ), ' qdtcn ', SUM( qdtcn)383 WRITE(numout,*) ' qstoif', SUM( qstoif ), ' fsbbq ', SUM( fsbbq)379 WRITE(numout,*) ' hicif ', SUM( hicif (2:nictl,2:njctl) ), ' hsnif ', SUM( hsnif (2:nictl,2:njctl) ) 380 WRITE(numout,*) ' frld ', SUM( frld (2:nictl,2:njctl) ), ' hicifp', SUM( hicifp(2:nictl,2:njctl) ) 381 WRITE(numout,*) ' phicif', SUM( phicif(2:nictl,2:njctl) ), ' pfrld ', SUM( pfrld (2:nictl,2:njctl) ) 382 WRITE(numout,*) ' sist ', SUM( sist (2:nictl,2:njctl) ), ' tbif 1', SUM( tbif (2:nictl,2:njctl,1) ) 383 WRITE(numout,*) ' tbif 2', SUM( tbif (2:nictl,2:njctl,2) ), ' tbif 3', SUM( tbif (2:nictl,2:njctl,3) ) 384 WRITE(numout,*) ' fdtcn ', SUM( fdtcn (2:nictl,2:njctl) ), ' qdtcn ', SUM( qdtcn (2:nictl,2:njctl) ) 385 WRITE(numout,*) ' qstoif', SUM( qstoif(2:nictl,2:njctl) ), ' fsbbq ', SUM( fsbbq (2:nictl,2:njctl) ) 384 386 ENDIF 385 387 -
trunk/NEMO/LIM_SRC/limthd_lac.F90
r88 r106 23 23 24 24 !! * Module variables 25 REAL(wp) :: 25 REAL(wp) :: & ! constant values 26 26 epsi20 = 1.e-20 , & 27 27 epsi13 = 1.e-13 , & -
trunk/NEMO/OPA_SRC/DOM/domain.F90
r93 r106 141 141 USE ioipsl 142 142 #endif 143 NAMELIST/namrun/ no , cexper , ln_rstart , nrstdt , nit000, & 144 & nitend, ndate0 , nleapy , ninist , nstock, & 145 & nprint, nwrite , nrunoff , ln_ctl , nbench 143 NAMELIST/namrun/ no , cexper , ln_rstart , nrstdt , nit000, & 144 & nitend, ndate0 , nleapy , ninist , nstock, & 145 & nprint, nwrite , nrunoff , ln_ctl , nictl , njctl, & 146 & nbench 146 147 NAMELIST/namdom/ ntopo , e3zps_min, e3zps_rat, ngrid , nmsh , & 147 148 & nacc , atfp , rdt , rdtmin , rdtmax, & … … 176 177 WRITE(numout,*) ' runoff option nrunoff = ', nrunoff 177 178 WRITE(numout,*) ' run control (for debugging) ln_ctl = ', ln_ctl 179 WRITE(numout,*) ' Max i indice for SUM control nictl = ', nictl 180 WRITE(numout,*) ' Max j indice for SUM control njctl = ', njctl 178 181 WRITE(numout,*) ' benchmark parameter (0/1) nbench = ', nbench 179 182 ENDIF … … 232 235 WRITE(numout,*) ' nb of seconds per year raass = ', raass, ' s' 233 236 WRITE(numout,*) ' nb of seconds per month rmoss = ', rmoss, ' s' 237 ENDIF 238 239 ! ... Control the Max i and j indices used for the SUM control (i.e. when ln_ctl=.true.) 240 IF ( nictl > jpim1 ) THEN 241 IF(lwp) THEN 242 WRITE(numout,cform_war) 243 WRITE(numout,*) ' nictl = ', nictl, ' must be <= to jpim1 ' 244 WRITE(numout,*) ' nictl forced to be equal to jpim1 ' 245 ENDIF 246 nwarn = nwarn + 1 247 nictl = jpim1 248 ENDIF 249 250 IF ( njctl > jpjm1 ) THEN 251 IF(lwp) THEN 252 WRITE(numout,cform_war) 253 WRITE(numout,*) ' njctl = ', njctl, ' must be <= to jpjm1 ' 254 WRITE(numout,*) ' njctl forced to be equal to jpjm1 ' 255 ENDIF 256 nwarn = nwarn + 1 257 njctl = jpjm1 234 258 ENDIF 235 259 … … 257 281 WRITE(numout,*) ' depth variation tracer step rdth = ', rdth 258 282 IF( lk_ice_lim ) THEN 259 WRITE(numout,*) ' ice model coupling frequency nfice = ', nfice283 IF(lwp)WRITE(numout,*) ' ice model coupling frequency nfice = ', nfice 260 284 nfbulk = nfice 261 285 ENDIF 262 WRITE(numout,*) ' bulk computation frequency nfbulk = ', nfbulk, ' = nfice if ice model used'263 WRITE(numout,*) ' flag closed sea or not nclosea = ', nclosea286 IF(lwp)WRITE(numout,*) ' bulk computation frequency nfbulk = ', nfbulk, ' = nfice if ice model used' 287 IF(lwp)WRITE(numout,*) ' flag closed sea or not nclosea = ', nclosea 264 288 ENDIF 265 289 -
trunk/NEMO/OPA_SRC/DYN/dynhpg.F90
r84 r106 163 163 164 164 IF(l_ctl) THEN ! print sum trends (used for debugging) 165 zuap = SUM( ua(2: jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) )166 zvap = SUM( va(2: jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) )165 zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 166 zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 167 167 WRITE(numout,*) ' hpg - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 168 168 u_ctl = zuap ; v_ctl = zvap … … 318 318 319 319 IF(l_ctl) THEN ! print sum trends (used for debugging) 320 zuap = SUM( ua(2: jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) )321 zvap = SUM( va(2: jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) )320 zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 321 zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 322 322 WRITE(numout,*) ' hpg - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 323 323 u_ctl = zuap ; v_ctl = zvap … … 429 429 430 430 IF(l_ctl) THEN ! print sum trends (used for debugging) 431 zuap = SUM( ua(2: jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) )432 zvap = SUM( va(2: jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) )431 zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 432 zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 433 433 WRITE(numout,*) ' hpg - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 434 434 u_ctl = zuap ; v_ctl = zvap -
trunk/NEMO/OPA_SRC/DYN/dynhpg_atsk.F90
r84 r106 151 151 152 152 IF(l_ctl) THEN ! print sum trends (used for debugging) 153 zuap = SUM( ua(2: jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) )154 zvap = SUM( va(2: jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) )153 zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 154 zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 155 155 WRITE(numout,*) ' hpg - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 156 156 u_ctl = zuap ; v_ctl = zvap … … 297 297 298 298 IF(l_ctl) THEN ! print sum trends (used for debugging) 299 zuap = SUM( ua(2: jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) )300 zvap = SUM( va(2: jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) )299 zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 300 zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 301 301 WRITE(numout,*) ' hpg - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 302 302 u_ctl = zuap ; v_ctl = zvap … … 411 411 412 412 IF(l_ctl) THEN ! print sum trends (used for debugging) 413 WRITE(numout,*) ' hpg - Ua: ', SUM(ua*umask), ' Va: ', SUM(va*vmask) 414 zuap = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 415 zvap = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 413 WRITE(numout,*) ' hpg - Ua: ', SUM(ua(2:nictl,2:njctl,1:jpkm1)*umask(2:nictl,2:njctl,1:jpkm1)) 414 WRITE(numout,*) ' - Va: ', SUM(va(2:nictl,2:njctl,1:jpkm1)*vmask(2:nictl,2:njctl,1:jpkm1)) 415 zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 416 zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 416 417 WRITE(numout,*) ' hpg - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 417 418 u_ctl = zuap ; v_ctl = zvap -
trunk/NEMO/OPA_SRC/DYN/dynkeg.F90
r84 r106 120 120 121 121 IF(l_ctl) THEN ! print sum trends (used for debugging) 122 zua = SUM( ua(2: jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) )123 zva = SUM( va(2: jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) )122 zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 123 zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 124 124 WRITE(numout,*) ' keg - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 125 125 u_ctl = zua ; v_ctl = zva -
trunk/NEMO/OPA_SRC/DYN/dynldf_bilap.F90
r84 r106 221 221 222 222 IF(l_ctl) THEN ! print sum trends (used for debugging) 223 zua = SUM( ua(2: jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) )224 zva = SUM( va(2: jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) )223 zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 224 zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 225 225 WRITE(numout,*) ' ldf - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 226 226 u_ctl = zua ; v_ctl = zva -
trunk/NEMO/OPA_SRC/DYN/dynldf_bilapg.F90
r84 r106 129 129 130 130 IF(l_ctl) THEN ! print sum trends (used for debugging) 131 zua = SUM( ua(2: jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) )132 zva = SUM( va(2: jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) )131 zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 132 zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 133 133 WRITE(numout,*) ' ldf - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 134 134 u_ctl = zua ; v_ctl = zva -
trunk/NEMO/OPA_SRC/DYN/dynldf_lap.F90
r84 r106 119 119 120 120 IF(l_ctl) THEN ! print sum trends (used for debugging) 121 zua = SUM( ua(2: jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) )122 zva = SUM( va(2: jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) )121 zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 122 zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 123 123 WRITE(numout,*) ' ldf - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 124 124 u_ctl = zua ; v_ctl = zva -
trunk/NEMO/OPA_SRC/DYN/dynnxt.F90
r84 r106 133 133 ! ! =============== 134 134 135 IF(l_ctl) WRITE(numout,*) ' nxt - Un: ', SUM(un*umask), ' Vn: ', SUM(vn*vmask) 135 IF(l_ctl) WRITE(numout,*) ' nxt - Un: ', SUM(un(2:nictl,2:njctl,1:jpkm1)*umask(2:nictl,2:njctl,1:jpkm1)), & 136 & ' Vn: ', SUM(vn(2:nictl,2:njctl,1:jpkm1)*vmask(2:nictl,2:njctl,1:jpkm1)) 136 137 137 138 END SUBROUTINE dyn_nxt -
trunk/NEMO/OPA_SRC/DYN/dynspg_fsc.F90
r84 r106 266 266 ! --------------------- 267 267 rnorme =0. 268 DO jj = 1, jpj 269 DO ji = 1, jpi 270 zgwgt = gcdmat(ji,jj) * gcb(ji,jj) 271 rnorme = rnorme + gcb(ji,jj) * zgwgt 272 END DO 273 END DO 268 rnorme = SUM( gcb(:,:) * gcdmat(:,:) * gcb(:,:) ) 274 269 IF( lk_mpp ) CALL mpp_sum( rnorme ) ! sum over the global domain 275 270 … … 348 343 349 344 IF(l_ctl) THEN ! print sum trends (used for debugging) 350 WRITE(numout,*) ' spg - Ua: ', SUM( ua(2: jpim1,2:jpjm1,1:jpkm1)*umask(2:jpim1,2:jpjm1,1:jpkm1) ), &351 & ' Va: ', SUM( va(2: jpim1,2:jpjm1,1:jpkm1)*vmask(2:jpim1,2:jpjm1,1:jpkm1) )345 WRITE(numout,*) ' spg - Ua: ', SUM( ua(2:nictl,2:njctl,1:jpkm1)*umask(2:nictl,2:njctl,1:jpkm1) ), & 346 & ' Va: ', SUM( va(2:nictl,2:njctl,1:jpkm1)*vmask(2:nictl,2:njctl,1:jpkm1) ) 352 347 ENDIF 353 348 … … 382 377 CALL lbc_lnk( sshn, 'T', 1. ) 383 378 379 IF(l_ctl) THEN ! print sum trends (used for debugging) 380 WRITE(numout,*) ' spg - ssh:', SUM( sshn(2:nictl,2:njctl)*tmask(2:nictl,2:njctl,1) ) 381 ENDIF 382 383 384 384 END SUBROUTINE dyn_spg_fsc 385 385 -
trunk/NEMO/OPA_SRC/DYN/dynspg_fsc_atsk.F90
r84 r106 363 363 364 364 IF(l_ctl) THEN ! print sum trends (used for debugging) 365 WRITE(numout,*) ' spg - Ua: ', SUM( ua(2: jpim1,2:jpjm1,1:jpkm1)*umask(2:jpim1,2:jpjm1,1:jpkm1) ), &366 & ' Va: ', SUM( va(2: jpim1,2:jpjm1,1:jpkm1)*vmask(2:jpim1,2:jpjm1,1:jpkm1) )365 WRITE(numout,*) ' spg - Ua: ', SUM( ua(2:nictl,2:njctl,1:jpkm1)*umask(2:nictl,2:njctl,1:jpkm1) ), & 366 & ' Va: ', SUM( va(2:nictl,2:njctl,1:jpkm1)*vmask(2:nictl,2:njctl,1:jpkm1) ) 367 367 ENDIF 368 368 … … 395 395 CALL lbc_lnk( sshn, 'T', 1. ) 396 396 397 IF(l_ctl) THEN ! print sum trends (used for debugging) 398 WRITE(numout,*) ' spg - ssh:', SUM( sshn(2:nictl,2:njctl)*tmask(2:nictl,2:njctl,1) ) 399 ENDIF 400 401 397 402 END SUBROUTINE dyn_spg_fsc_atsk 398 403 -
trunk/NEMO/OPA_SRC/DYN/dynvor.F90
r52 r106 146 146 147 147 IF(l_ctl) THEN ! print sum trends (used for debugging) 148 zua = SUM( ua(2: jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) )149 zva = SUM( va(2: jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) )148 zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 149 zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 150 150 WRITE(numout,*) ' vor - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 151 151 u_ctl = zua ; v_ctl = zva … … 261 261 262 262 IF(l_ctl) THEN ! print sum trends (used for debugging) 263 zua = SUM( ua(2: jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) )264 zva = SUM( va(2: jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) )263 zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 264 zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 265 265 WRITE(numout,*) ' vor - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 266 266 u_ctl = zua ; v_ctl = zva … … 395 395 396 396 IF(l_ctl) THEN ! print sum trends (used for debugging) 397 zua = SUM( ua(2: jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) )398 zva = SUM( va(2: jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) )397 zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 398 zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 399 399 WRITE(numout,*) ' vor - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 400 400 u_ctl = zua ; v_ctl = zva -
trunk/NEMO/OPA_SRC/DYN/dynzad.F90
r84 r106 127 127 128 128 IF(l_ctl) THEN ! print sum trends (used for debugging) 129 zua = SUM( ua(2: jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) )130 zva = SUM( va(2: jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) )129 zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 130 zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 131 131 WRITE(numout,*) ' zad - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 132 132 u_ctl = zua ; v_ctl = zva … … 228 228 229 229 IF(l_ctl) THEN ! print sum trends (used for debugging) 230 zua = SUM( ua(2: jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) )231 zva = SUM( va(2: jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) )230 zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 231 zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 232 232 WRITE(numout,*) ' zad - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 233 233 u_ctl = zua ; v_ctl = zva -
trunk/NEMO/OPA_SRC/DYN/dynzdf_imp.F90
r84 r106 343 343 344 344 IF(l_ctl) THEN ! print sum trends (used for debugging) 345 zua = SUM( ua(2: jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) )346 zva = SUM( va(2: jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) )345 zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 346 zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 347 347 WRITE(numout,*) ' zdf - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 348 348 u_ctl = zua ; v_ctl = zva -
trunk/NEMO/OPA_SRC/DYN/dynzdf_imp_atsk.F90
r84 r106 314 314 315 315 IF(l_ctl) THEN ! print sum trends (used for debugging) 316 zua = SUM( ua(2: jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) )317 zva = SUM( va(2: jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) )316 zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 317 zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 318 318 WRITE(numout,*) ' zdf - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 319 319 u_ctl = zua ; v_ctl = zva -
trunk/NEMO/OPA_SRC/DYN/wzvmod.F90
r3 r106 79 79 80 80 IF(l_ctl) THEN ! print mean trends (used for debugging) 81 WRITE(numout,*) ' w**2 - : ', SUM( wn( :,:,:) * wn (:,:,:) )81 WRITE(numout,*) ' w**2 - : ', SUM( wn(2:nictl,2:njctl,1:jpkm1) * wn (2:nictl,2:njctl,1:jpkm1) ) 82 82 ENDIF 83 83 … … 130 130 131 131 IF(l_ctl) THEN ! print mean trends (used for debugging) 132 WRITE(numout,*) ' w**2 - : ', SUM( wn( :,:,:) * wn (:,:,:) )132 WRITE(numout,*) ' w**2 - : ', SUM( wn(2:nictl,2:njctl,1:jpkm1) * wn (2:nictl,2:njctl,1:jpkm1) ) 133 133 ENDIF 134 134 -
trunk/NEMO/OPA_SRC/LDF/ldfeiv.F90
r80 r106 158 158 CALL lbc_lnk( aeiv, 'V', 1. ) 159 159 160 IF(l_ctl) WRITE(numout,*) ' eiv - u: ', SUM( aeiu( :,:) ), ' v: ', SUM( aeiv(:,:) )160 IF(l_ctl) WRITE(numout,*) ' eiv - u: ', SUM( aeiu(1:nictl+1,1:njctl+1) ), ' v: ', SUM( aeiv(1:nictl+1,1:njctl+1) ) 161 161 162 162 … … 174 174 END DO 175 175 END DO 176 IF(l_ctl) WRITE(numout,*) ' aht - u: ', SUM( ahtu(:,:) ), ' v: ', SUM( ahtv(:,:) ), ' w: ', SUM( ahtw(:,:) ) 176 IF(l_ctl) THEN 177 WRITE(numout,*)' aht - u: ', SUM( ahtu(1:nictl+1,1:njctl+1) ) 178 WRITE(numout,*)' - v: ', SUM( ahtv(1:nictl+1,1:njctl+1) ) 179 WRITE(numout,*)' - w: ', SUM( ahtw(1:nictl+1,1:njctl+1) ) 180 ENDIF 177 181 ENDIF 178 182 … … 319 323 CALL lbc_lnk( aeiv, 'V', 1. ) 320 324 321 IF(l_ctl) WRITE(numout,*) ' eiv - u: ', SUM( aeiu( :,:) ), ' v: ', SUM( aeiv(:,:) )325 IF(l_ctl) WRITE(numout,*) ' eiv - u: ', SUM( aeiu(1:nictl+1,1:njctl+1) ), ' v: ', SUM( aeiv(1:nictl+1,1:njctl+1) ) 322 326 323 327 ! ORCA R05: add a space variation on aht (=aeiv except at the equator and river mouth) … … 334 338 END DO 335 339 END DO 336 IF(l_ctl) WRITE(numout,*) ' aht - u: ', SUM( ahtu(:,:) ), ' v: ', SUM( ahtv(:,:) ), ' w: ', SUM( ahtw(:,:) ) 340 IF(l_ctl) THEN 341 WRITE(numout,*) ' aht - u: ', SUM( ahtu(1:nictl+1,1:njctl+1) ) 342 WRITE(numout,*) ' - v: ', SUM( ahtv(1:nictl+1,1:njctl+1) ) 343 WRITE(numout,*) ' - w: ', SUM( ahtw(1:nictl+1,1:njctl+1) ) 344 ENDIF 337 345 ENDIF 338 346 -
trunk/NEMO/OPA_SRC/LDF/ldfslp.F90
r49 r106 479 479 480 480 IF(l_ctl) THEN 481 WRITE(numout,*) ' slp - u : ', SUM( uslp ( :,:,:) ), ' v : ', SUM( vslp (:,:,:) )482 WRITE(numout,*) ' - wi: ', SUM( wslpi( :,:,:) ), ' wj: ', SUM( wslpj(:,:,:) )481 WRITE(numout,*) ' slp - u : ', SUM( uslp (1:nictl,1:njctl,:) ), ' v : ', SUM( vslp (1:nictl,1:njctl,:) ) 482 WRITE(numout,*) ' - wi: ', SUM( wslpi(1:nictl,1:njctl,:) ), ' wj: ', SUM( wslpj(1:nictl,1:njctl,:) ) 483 483 ENDIF 484 484 -
trunk/NEMO/OPA_SRC/SBC/bulk.F90
r84 r106 97 97 IF(l_ctl) THEN ! print mean trends (used for debugging) 98 98 WRITE(numout,*) ' Forcings ' 99 WRITE(numout,*) ' qsr_oce : ', SUM( qsr_oce ( :,:) * tmask(:,:,1) )100 WRITE(numout,*) ' qsr_ice : ', SUM( qsr_ice ( :,:) * tmask(:,:,1) )101 WRITE(numout,*) ' qnsr_oce : ', SUM( qnsr_oce( :,:) * tmask(:,:,1) )102 WRITE(numout,*) ' qnsr_ice : ', SUM( qnsr_ice( :,:) * tmask(:,:,1) )103 WRITE(numout,*) ' evap : ', SUM( evap ( :,:) * tmask(:,:,1) )104 WRITE(numout,*) ' precip : ', SUM( tprecip ( :,:) * tmask(:,:,1) ) / rday105 WRITE(numout,*) ' Snow : ', SUM( sprecip ( :,:) * tmask(:,:,1) ) / rday106 WRITE(numout,*) ' u-stress : ', SUM( taux ( :,:) * umask(:,:,1) )107 WRITE(numout,*) ' v-stress : ', SUM( tauy ( :,:) * vmask(:,:,1) )108 WRITE(numout,*) ' sst : ', SUM( zsst ( :,:) * tmask(:,:,1) )109 WRITE(numout,*) ' sss : ', SUM( zsss ( :,:) * tmask(:,:,1) )99 WRITE(numout,*) ' qsr_oce : ', SUM( qsr_oce (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 100 WRITE(numout,*) ' qsr_ice : ', SUM( qsr_ice (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 101 WRITE(numout,*) ' qnsr_oce : ', SUM( qnsr_oce(1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 102 WRITE(numout,*) ' qnsr_ice : ', SUM( qnsr_ice(1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 103 WRITE(numout,*) ' evap : ', SUM( evap (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 104 WRITE(numout,*) ' precip : ', SUM( tprecip (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) / rday 105 WRITE(numout,*) ' Snow : ', SUM( sprecip (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) / rday 106 WRITE(numout,*) ' u-stress : ', SUM( taux (1:nictl+1,1:njctl+1) * umask(1:nictl+1,1:njctl+1,1) ) 107 WRITE(numout,*) ' v-stress : ', SUM( tauy (1:nictl+1,1:njctl+1) * vmask(1:nictl+1,1:njctl+1,1) ) 108 WRITE(numout,*) ' sst : ', SUM( zsst (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 109 WRITE(numout,*) ' sss : ', SUM( zsss (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 110 110 WRITE(numout,*) 111 111 ENDIF -
trunk/NEMO/OPA_SRC/TRA/traadv_cen2.F90
r74 r106 301 301 302 302 IF(l_ctl) THEN 303 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )304 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )303 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 304 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 305 305 WRITE(numout,*) ' had - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 306 306 t_ctl = zta ; s_ctl = zsa … … 396 396 397 397 IF(l_ctl) THEN 398 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )399 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )398 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 399 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 400 400 WRITE(numout,*) ' zad - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl, ' centered2' 401 401 t_ctl = zta ; s_ctl = zsa -
trunk/NEMO/OPA_SRC/TRA/traadv_cen2_atsk.h90
r74 r106 249 249 ! ! =============== 250 250 251 IF( l_ctl .AND. lwp) THEN252 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )253 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )251 IF(l_ctl) THEN 252 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 253 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 254 254 WRITE(numout,*) ' had - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 255 255 t_ctl = zta ; s_ctl = zsa … … 342 342 ! ! =============== 343 343 344 IF( l_ctl .AND. lwp) THEN345 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )346 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )344 IF(l_ctl) THEN 345 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 346 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 347 347 WRITE(numout,*) ' zad - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl, ' centered2 autotsk' 348 348 t_ctl = zta ; s_ctl = zsa -
trunk/NEMO/OPA_SRC/TRA/traadv_muscl.F90
r97 r106 268 268 269 269 IF(l_ctl) THEN 270 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )271 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )270 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 271 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 272 272 WRITE(numout,*) ' had - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl, ' muscl' 273 273 t_ctl = zta ; s_ctl = zsa … … 394 394 END DO 395 395 396 IF( l_ctl .AND. lwp) THEN397 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )398 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )396 IF(l_ctl) THEN 397 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 398 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 399 399 WRITE(numout,*) ' zad - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl, ' muscl' 400 400 t_ctl = zta ; s_ctl = zsa -
trunk/NEMO/OPA_SRC/TRA/traadv_muscl2.F90
r97 r106 329 329 330 330 IF(l_ctl) THEN 331 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )332 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )331 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 332 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 333 333 WRITE(numout,*) ' had - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl, ' muscl' 334 334 t_ctl = zta ; s_ctl = zsa … … 468 468 END DO 469 469 470 IF( l_ctl .AND. lwp) THEN471 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )472 zsa = SUM( sa(2:jpim1,2: jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )470 IF(l_ctl) THEN 471 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 472 zsa = SUM( sa(2:jpim1,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 473 473 WRITE(numout,*) ' zad - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl, ' muscl2' 474 474 t_ctl = zta ; s_ctl = zsa -
trunk/NEMO/OPA_SRC/TRA/traadv_tvd.F90
r74 r106 260 260 END DO 261 261 262 IF( l_ctl .AND. lwp) THEN263 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )264 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )262 IF(l_ctl) THEN 263 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 264 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 265 265 WRITE(numout,*) ' zad - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl, ' tvd' 266 266 t_ctl = zta ; s_ctl = zsa -
trunk/NEMO/OPA_SRC/TRA/trabbc.F90
r32 r106 104 104 END DO 105 105 #endif 106 IF( l_ctl .AND. lwp ) THEN 107 zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 106 107 IF(l_ctl) THEN 108 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 108 109 WRITE(numout,*) ' bbc - Ta: ', zta-t_ctl 109 110 t_ctl = zta -
trunk/NEMO/OPA_SRC/TRA/trabbl.F90
r87 r106 340 340 341 341 IF(l_ctl) THEN 342 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )343 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )342 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 343 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 344 344 WRITE(numout,*) ' bbl - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 345 345 t_ctl = zta ; s_ctl = zsa -
trunk/NEMO/OPA_SRC/TRA/trabbl_adv.h90
r21 r106 407 407 END DO 408 408 409 IF( l_ctl .AND. lwp) THEN410 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )411 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )409 IF(l_ctl) THEN 410 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 411 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 412 412 WRITE(numout,*) ' bbl - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 413 413 t_ctl = zta ; s_ctl = zsa -
trunk/NEMO/OPA_SRC/TRA/tradmp.F90
r61 r106 178 178 179 179 IF(l_ctl) THEN ! print mean trends (used for debugging) 180 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )181 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )180 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 181 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 182 182 WRITE(numout,*) ' dmp - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 183 183 t_ctl = zta ; s_ctl = zsa -
trunk/NEMO/OPA_SRC/TRA/traldf_iso.F90
r74 r106 232 232 ! ! =============== 233 233 234 IF( l_ctl) THEN ! print mean trends (used for debugging)235 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )236 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )234 IF(l_ctl) THEN ! print mean trends (used for debugging) 235 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 236 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 237 237 WRITE(numout,*) ' ldf - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 238 238 t_ctl = zta ; s_ctl = zsa -
trunk/NEMO/OPA_SRC/TRA/traldf_iso_zps.F90
r63 r106 268 268 ! ! =============== 269 269 270 IF( l_ctl) THEN ! print mean trends (used for debugging)271 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )272 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )270 IF(l_ctl) THEN ! print mean trends (used for debugging) 271 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 272 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 273 273 WRITE(numout,*) ' ldf - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 274 274 t_ctl = zta ; s_ctl = zsa -
trunk/NEMO/OPA_SRC/TRA/traldf_lap.F90
r74 r106 142 142 ! ! ============= 143 143 144 IF( l_ctl) THEN ! print mean trends (used for debugging)145 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )146 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )144 IF(l_ctl) THEN ! print mean trends (used for debugging) 145 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 146 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 147 147 WRITE(numout,*) ' ldf - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 148 148 t_ctl = zta ; s_ctl = zsa -
trunk/NEMO/OPA_SRC/TRA/tranxt.F90
r87 r106 164 164 165 165 IF(l_ctl) THEN ! print mean field (used for debugging) 166 WRITE(numout,*) ' nxt - Tn: ', SUM(tn*tmask), ' Sn: ', SUM(sn*tmask) 166 WRITE(numout,*) ' nxt - Tn: ', SUM(tn(2:nictl,2:njctl,1:jpkm1)*tmask(2:nictl,2:njctl,1:jpkm1)), & 167 & ' Sn: ', SUM(sn(2:nictl,2:njctl,1:jpkm1)*tmask(2:nictl,2:njctl,1:jpkm1)) 167 168 ENDIF 168 169 -
trunk/NEMO/OPA_SRC/TRA/traqsr.F90
r87 r106 173 173 174 174 IF(l_ctl) THEN ! print mean trends (used for debugging) 175 zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 176 zta = SUM( ta * tmask ) 175 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 177 176 WRITE(numout,*) ' qsr - Ta: ', zta-t_ctl 178 177 t_ctl = zta -
trunk/NEMO/OPA_SRC/TRA/trasbc.F90
r87 r106 122 122 123 123 IF(l_ctl) THEN ! print mean trends (used for debugging) 124 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )125 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )124 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 125 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 126 126 WRITE(numout,*) ' sbc - Ta: ', zta, ' Sa: ', zsa 127 127 t_ctl = zta ; s_ctl = zsa -
trunk/NEMO/OPA_SRC/TRA/trazdf_iso.F90
r87 r106 371 371 372 372 IF(l_ctl) THEN ! print mean trends (used for debugging) 373 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )374 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )373 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 374 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 375 375 WRITE(numout,*) ' zdf 1- Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 376 376 t_ctl = zta ; s_ctl = zsa … … 534 534 535 535 IF(l_ctl) THEN ! print mean trends (used for debugging) 536 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )537 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )536 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 537 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 538 538 WRITE(numout,*) ' zdf 2- Ta: ', zta, ' Sa: ', zsa 539 539 ENDIF -
trunk/NEMO/OPA_SRC/TRA/trazdf_iso_vopt.F90
r67 r106 81 81 82 82 IF(l_ctl) THEN ! print mean trends (used for debugging) 83 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )84 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )83 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 84 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 85 85 WRITE(numout,*) ' zdf 1- Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 86 86 t_ctl = zta ; s_ctl = zsa … … 93 93 94 94 IF(l_ctl) THEN ! print mean trends (used for debugging) 95 zta = SUM( ta(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )96 zsa = SUM( sa(2: jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) )97 WRITE(numout,*) ' zdf 1- Ta: ', zta, ' Sa: ', zsa95 zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 96 zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 97 WRITE(numout,*) ' zdf 2- Ta: ', zta, ' Sa: ', zsa 98 98 ENDIF 99 99 -
trunk/NEMO/OPA_SRC/ZDF/zdfbfr.F90
r49 r106 168 168 CALL lbc_lnk( avmv, 'V', 1. ) 169 169 170 IF( l_ctl ) WRITE(numout,*) ' bfr u : ', SUM( avmu(:,:,:) ), ' v : ', SUM( avmv(:,:,:) )170 IF(l_ctl) WRITE(numout,*) ' bfr u : ', SUM( avmu(1:nictl+1,1:njctl+1,:) ), ' v : ', SUM( avmv(1:nictl+1,1:njctl+1,:) ) 171 171 172 172 END SUBROUTINE zdf_bfr -
trunk/NEMO/OPA_SRC/ZDF/zdfddm.F90
r49 r106 199 199 CALL lbc_lnk( avmv, 'V', 1. ) 200 200 201 IF( l_ctl) THEN202 WRITE(numout,*) ' ddm t : ', SUM( avt ( :,:,:) ), ' s : ', SUM( avs (:,:,:) )203 WRITE(numout,*) ' u : ', SUM( avmu( :,:,:) ), ' v : ', SUM( avmv(:,:,:) )201 IF(l_ctl) THEN 202 WRITE(numout,*) ' ddm t : ', SUM( avt (1:nictl+1,1:njctl+1,:) ), ' s : ', SUM( avs (1:nictl+1,1:njctl+1,:) ) 203 WRITE(numout,*) ' u : ', SUM( avmu(1:nictl+1,1:njctl+1,:) ), ' v : ', SUM( avmv(1:nictl+1,1:njctl+1,:) ) 204 204 ENDIF 205 205 -
trunk/NEMO/OPA_SRC/ZDF/zdfmxl.F90
r49 r106 132 132 ! ! =============== 133 133 134 IF( l_ctl ) WRITE(numout,*) ' nmln : ', SUM( nmln(:,:) ), ' hmld: ', SUM( hmld(:,:) )134 IF(l_ctl) WRITE(numout,*) ' nmln : ', SUM( nmln(1:nictl+1,1:njctl+1) ), ' hmld: ', SUM( hmld(1:nictl+1,1:njctl+1) ) 135 135 136 136 END SUBROUTINE zdf_mxl … … 175 175 IF(lwp) WRITE(numout,*) 'zdf_mxl : mixed layer depth' 176 176 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 177 IF(lwp) WRITE(numout,*)178 177 ENDIF 179 178 … … 234 233 END DO 235 234 236 IF( l_ctl ) WRITE(numout,*) ' nmln : ', SUM( nmln(:,:) ), ' hmld: ', SUM( hmld(:,:) )235 IF(l_ctl) WRITE(numout,*) ' nmln : ', SUM( nmln(1:nictl+1,1:njctl+1) ), ' hmld: ', SUM( hmld(1:nictl+1,1:njctl+1) ) 237 236 238 237 END SUBROUTINE zdf_mxl -
trunk/NEMO/OPA_SRC/ZDF/zdftke.F90
r49 r106 578 578 CALL lbc_lnk( avt, 'W', 1. ) 579 579 580 IF( l_ctl) THEN581 WRITE(numout,*) ' tke e : ', SUM( en ( :,:,:) ), ' t : ', SUM( avt (:,:,:) )582 WRITE(numout,*) ' u : ', SUM( avmu( :,:,:) ), ' v : ', SUM( avmv(:,:,:) )580 IF(l_ctl) THEN 581 WRITE(numout,*) ' tke e : ', SUM( en (1:nictl+1,1:njctl+1,:) ), ' t : ', SUM( avt (1:nictl+1,1:njctl+1,:) ) 582 WRITE(numout,*) ' u : ', SUM( avmu(1:nictl+1,1:njctl+1,:) ), ' v : ', SUM( avmv(1:nictl+1,1:njctl+1,:) ) 583 583 ENDIF 584 584 -
trunk/NEMO/OPA_SRC/ZDF/zdftke_atsk.h90
r3 r106 476 476 CALL lbc_lnk( avt, 'W', 1. ) 477 477 478 IF(l_ctl) THEN 479 WRITE(numout,*) ' tke e : ', SUM( en (1:nictl+1,1:njctl+1,:) ), ' t : ', SUM( avt (1:nictl+1,1:njctl+1,:) ) 480 WRITE(numout,*) ' u : ', SUM( avmu(1:nictl+1,1:njctl+1,:) ), ' v : ', SUM( avmv(1:nictl+1,1:njctl+1,:) ) 481 ENDIF 482 483 478 484 END SUBROUTINE zdf_tke -
trunk/NEMO/OPA_SRC/eosbn2.F90
r99 r106 222 222 END SELECT 223 223 224 IF(l_ctl) WRITE(numout,*) ' eos : ', SUM( prd( :,:,:) )224 IF(l_ctl) WRITE(numout,*) ' eos : ', SUM( prd(1:nictl+1,1:njctl+1,:) ) 225 225 226 226 END SUBROUTINE eos_insitu … … 408 408 END SELECT 409 409 410 IF(l_ctl) WRITE(numout,*) ' eos-p: ', SUM( prd( :,:,:) ), ' pot : ', SUM( prhop(:,:,:) )410 IF(l_ctl) WRITE(numout,*) ' eos-p: ', SUM( prd(1:nictl+1,1:njctl+1,:) ), ' pot : ', SUM( prhop(1:nictl+1,1:njctl+1,:) ) 411 411 412 412 END SUBROUTINE eos_insitu_pot … … 584 584 END SELECT 585 585 586 IF(l_ctl) WRITE(numout,*) ' eos2d: ', SUM( prd(1: jpim1,1:jpjm1) )586 IF(l_ctl) WRITE(numout,*) ' eos2d: ', SUM( prd(1:nictl+1,1:njctl+1) ) 587 587 588 588 END SUBROUTINE eos_insitu_2d … … 759 759 END SELECT 760 760 761 IF(l_ctl) WRITE(numout,*) ' bn2 : ', SUM( pn2 ( :,:,:) )761 IF(l_ctl) WRITE(numout,*) ' bn2 : ', SUM( pn2 (1:nictl+1,1:njctl+1,:) ) 762 762 #if defined key_zdfddm 763 IF(l_ctl) WRITE(numout,*) ' rrau : ', SUM( rrau( :,:,:) )763 IF(l_ctl) WRITE(numout,*) ' rrau : ', SUM( rrau(1:nictl+1,1:njctl+1,:) ) 764 764 #endif 765 765 -
trunk/NEMO/OPA_SRC/in_out_manager.F90
r100 r106 3 3 USE lib_print ! formated print library 4 4 USE par_kind 5 USE par_oce 5 6 6 7 PUBLIC … … 49 50 50 51 INTEGER :: & !: 51 nstock = 10 , & !: restart file frequency 52 nprint = 0 , & !: level of print (0 no print) 53 nwrite = 10 !: restart file frequency 52 nstock = 10 , & !: restart file frequency 53 nprint = 0 , & !: level of print (0 no print) 54 nwrite = 10 , & !: restart file frequency 55 nictl = jpim1, & !: max i indice for the SUM control 56 njctl = jpjm1 !: max j indice for the SUM control 54 57 !!---------------------------------------------------------------------- 55 58 !! logical units -
trunk/NEMO/OPA_SRC/step.F90
r75 r106 214 214 215 215 IF(l_ctl) THEN ! print mean trends (used for debugging) 216 WRITE(numout,*) ' emp - : ', SUM( emp * tmask(:,:,1) )217 WRITE(numout,*) ' emps - : ', SUM( emps * tmask(:,:,1) )218 WRITE(numout,*) ' qt - : ', SUM( qt * tmask(:,:,1) )219 WRITE(numout,*) ' qsr - : ', SUM( qsr * tmask(:,:,1) )220 WRITE(numout,*) ' runoff : ', SUM( runoff * tmask(:,:,1) )221 WRITE(numout,*) ' tmask : ', SUM( tmask 222 WRITE(numout,*) ' sst - : ', SUM( tn (:,:,1) * tmask(:,:,1) )223 WRITE(numout,*) ' sss - : ', SUM( sn (:,:,1) * tmask(:,:,1) )224 WRITE(numout,*) ' tau - x : ', SUM( taux ), ' - y : ', SUM( tauy)216 WRITE(numout,*) ' emp - : ', SUM( emp (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 217 WRITE(numout,*) ' emps - : ', SUM( emps (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 218 WRITE(numout,*) ' qt - : ', SUM( qt (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 219 WRITE(numout,*) ' qsr - : ', SUM( qsr (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 220 WRITE(numout,*) ' runoff : ', SUM( runoff(1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 221 WRITE(numout,*) ' tmask : ', SUM( tmask (1:nictl+1,1:njctl+1,:) ) 222 WRITE(numout,*) ' sst - : ', SUM( tn (1:nictl+1,1:njctl+1,1) * tmask(1:nictl+1,1:njctl+1,1) ) 223 WRITE(numout,*) ' sss - : ', SUM( sn (1:nictl+1,1:njctl+1,1) * tmask(1:nictl+1,1:njctl+1,1) ) 224 WRITE(numout,*) ' tau - x : ', SUM( taux (1:nictl+1,1:njctl+1) ), ' - y : ', SUM( tauy(1:nictl+1,1:njctl+1) ) 225 225 ENDIF 226 226
Note: See TracChangeset
for help on using the changeset viewer.