Changeset 9260
- Timestamp:
- 2018-01-18T14:47:37+01:00 (7 years ago)
- Location:
- branches/NERC/dev_r5518_GO6_conserv_check_up/NEMOGCM/NEMO/TOP_SRC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_GO6_conserv_check_up/NEMOGCM/NEMO/TOP_SRC/trcini.F90
r9208 r9260 218 218 IF ( lwp ) WRITE(numout,*) 219 219 IF ( lwp ) WRITE(numout,*) ' Elemental cycle totals: ' 220 IF ( lwp ) CALL flush(numout)221 220 ! nitrogen 222 221 ztot3d(:,:,:) = trn(:,:,:,jpphn) + trn(:,:,:,jpphd) + trn(:,:,:,jpzmi) + & … … 227 226 cycletot(1) = zsum3d + zsum2d 228 227 IF ( lwp ) WRITE(numout,9010) 'nitrogen', zsum3d, zsum2d, cycletot(1) 229 IF ( lwp ) CALL flush(numout)230 228 ! silicon 231 229 ztot3d(:,:,:) = trn(:,:,:,jppds) + trn(:,:,:,jpsil) … … 235 233 cycletot(2) = zsum3d + zsum2d 236 234 IF ( lwp ) WRITE(numout,9010) 'silicon', zsum3d, zsum2d, cycletot(2) 237 IF ( lwp ) CALL flush(numout)238 235 ! iron 239 236 ztot3d(:,:,:) = ((trn(:,:,:,jpphn) + trn(:,:,:,jpphd) + trn(:,:,:,jpzmi) + & … … 244 241 cycletot(3) = zsum3d + zsum2d 245 242 IF ( lwp ) WRITE(numout,9010) 'iron', zsum3d, zsum2d, cycletot(3) 246 IF ( lwp ) CALL flush(numout)247 243 ! carbon (uses fixed C:N ratios on plankton tracers) 248 244 ztot3d(:,:,:) = (trn(:,:,:,jpphn) * xthetapn) + (trn(:,:,:,jpphd) * xthetapd) + & … … 254 250 cycletot(4) = zsum3d + zsum2d 255 251 IF ( lwp ) WRITE(numout,9010) 'carbon', zsum3d, zsum2d, cycletot(4) 256 IF ( lwp ) CALL flush(numout)257 252 ! alkalinity (note benthic correction) 258 253 ztot3d(:,:,:) = trn(:,:,:,jpalk) … … 262 257 cycletot(5) = zsum3d + zsum2d 263 258 IF ( lwp ) WRITE(numout,9010) 'alkalinity', zsum3d, zsum2d, cycletot(5) 264 IF ( lwp ) CALL flush(numout)265 259 ! oxygen (note no benthic) 266 260 ztot3d(:,:,:) = trn(:,:,:,jpoxy) … … 270 264 cycletot(6) = zsum3d + zsum2d 271 265 IF ( lwp ) WRITE(numout,9010) 'oxygen', zsum3d, zsum2d, cycletot(6) 272 IF ( lwp ) CALL flush(numout)273 266 ! Check 274 267 zsum3d = glob_sum( cvol(:,:,:) ) 275 268 zsum2d = glob_sum( carea(:,:) ) 276 IF ( lwp ) WRITE(numout,*) 277 IF ( lwp ) WRITE(numout,*) ' check : cvol : ', zsum3d 278 IF ( lwp ) WRITE(numout,*) ' check : carea : ', zsum2d 279 IF ( lwp ) WRITE(numout,*) 280 IF ( lwp ) CALL flush(numout) 281 ! 282 # endif 283 284 IF(lwp) WRITE(numout,*) 285 IF(lwp) WRITE(numout,*) 'trc_init : passive tracer set up completed' 286 IF(lwp) WRITE(numout,*) '~~~~~~~' 287 IF(lwp) CALL flush(numout) 269 IF ( lwp ) THEN 270 WRITE(numout,*) 271 WRITE(numout,*) ' check : cvol : ', zsum3d 272 WRITE(numout,*) ' check : carea : ', zsum2d 273 WRITE(numout,*) 274 ENDIF 275 ! 276 # endif 277 278 IF(lwp) THEN 279 WRITE(numout,*) 280 WRITE(numout,*) 'trc_init : passive tracer set up completed' 281 WRITE(numout,*) '~~~~~~~' 282 ENDIF 288 283 # if defined key_debug_medusa 289 284 CALL trc_rst_stat -
branches/NERC/dev_r5518_GO6_conserv_check_up/NEMOGCM/NEMO/TOP_SRC/trcrst.F90
r9208 r9260 844 844 zsum3d = glob_sum( zvol(:,:,:) ) 845 845 zsum2d = glob_sum( zarea(:,:) ) 846 IF ( lwp ) WRITE(numout,*) 847 IF ( lwp ) WRITE(numout,*) ' check : cvol : ', zsum3d 848 IF ( lwp ) WRITE(numout,*) ' check : carea : ', zsum2d 849 IF ( lwp ) WRITE(numout,*) 850 IF ( lwp ) CALL flush(numout) 846 IF ( lwp ) THEN 847 WRITE(numout,*) 848 WRITE(numout,*) ' check : cvol : ', zsum3d 849 WRITE(numout,*) ' check : carea : ', zsum2d 850 WRITE(numout,*) 851 ENDIF 851 852 ! 852 853 9010 FORMAT(' element:',a10, &
Note: See TracChangeset
for help on using the changeset viewer.