- Timestamp:
- 2013-02-08T12:13:47+01:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_MERGE_2012/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsms.F90
r3557 r3780 373 373 374 374 IF( ln_check_mass .AND. kt == nitend ) THEN ! Compute the budget of NO3, ALK, Si, Fer 375 no3budget = glob_sum( ( trn(:,:,:,jpno3) + trn(:,:,:,jpnh4) & 376 & + trn(:,:,:,jpphy) + trn(:,:,:,jpdia) & 377 & + trn(:,:,:,jpzoo) + trn(:,:,:,jpmes) & 375 no3budget = glob_sum( ( trn(:,:,:,jpno3) + trn(:,:,:,jpnh4) & 376 & + trn(:,:,:,jpphy) + trn(:,:,:,jpdia) & 377 & + trn(:,:,:,jpzoo) + trn(:,:,:,jpmes) & 378 & + trn(:,:,:,jppoc) & 378 379 #if ! defined key_kriest 379 & + trn(:,:,:,jppoc) + trn(:,:,:,jpgoc) & 380 #else 381 & + trn(:,:,:,jppoc) & 380 & + trn(:,:,:,jpgoc) & 382 381 #endif 383 & 382 & + trn(:,:,:,jpdoc) ) * cvol(:,:,:) ) 384 383 ! 385 384 silbudget = glob_sum( ( trn(:,:,:,jpsil) + trn(:,:,:,jpgsi) & 386 & 385 & + trn(:,:,:,jpdsi) ) * cvol(:,:,:) ) 387 386 ! 388 387 alkbudget = glob_sum( ( trn(:,:,:,jpno3) * rno3 & 389 & 390 & 388 & + trn(:,:,:,jptal) & 389 & + trn(:,:,:,jpcal) * 2. ) * cvol(:,:,:) ) 391 390 ! 392 391 ferbudget = glob_sum( ( trn(:,:,:,jpfer) + trn(:,:,:,jpnfe) & 393 & + trn(:,:,:,jpdfe) + trn(:,:,:,jpbfe) & 394 & + trn(:,:,:,jpsfe) & 395 & + trn(:,:,:,jpzoo) & 396 & + trn(:,:,:,jpmes) * ferat3 ) * cvol(:,:,:) ) 392 & + trn(:,:,:,jpdfe) & 393 #if ! defined key_kriest 394 & + trn(:,:,:,jpbfe) & 395 #endif 396 & + trn(:,:,:,jpsfe) & 397 & + trn(:,:,:,jpzoo) & 398 & + trn(:,:,:,jpmes) * ferat3 ) * cvol(:,:,:) ) 397 399 398 400 !
Note: See TracChangeset
for help on using the changeset viewer.