Changeset 9336


Ignore:
Timestamp:
2018-02-19T18:37:36+01:00 (2 years ago)
Author:
dford
Message:

Add growth and loss diagnostics.

Location:
branches/UKMO/CO6_KD490_amm7_oper_fabm_hem08/NEMOGCM/NEMO/TOP_SRC/FABM
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/CO6_KD490_amm7_oper_fabm_hem08/NEMOGCM/NEMO/TOP_SRC/FABM/par_fabm.F90

    r9330 r9336  
    3636 
    3737   INTEGER, PUBLIC :: jp_fabmdia_chltot, jp_fabmdia_netpp, jp_fabmdia_phytot 
     38   INTEGER, PUBLIC :: jp_fabmdia_pgrow,  jp_fabmdia_ploss 
    3839 
    3940#if defined key_fabm 
  • branches/UKMO/CO6_KD490_amm7_oper_fabm_hem08/NEMOGCM/NEMO/TOP_SRC/FABM/trcini_fabm.F90

    r9330 r9336  
    121121      jp_fabm_xeps = fabm_diag_index( 'light_xEPS' ) 
    122122      jp_fabmdia_phytot = fabm_diag_index( 'total_phytoplankton_result' ) 
     123      jp_fabmdia_pgrow = fabm_diag_index( 'p_grow_sum_result' ) 
     124      jp_fabmdia_ploss = fabm_diag_index( 'p_loss_sum_result' ) 
     125       
    123126 
    124127      IF(lwp) WRITE(numout,*) 'DAF: jp_fabm_n5s = ', jp_fabm_n5s 
  • branches/UKMO/CO6_KD490_amm7_oper_fabm_hem08/NEMOGCM/NEMO/TOP_SRC/FABM/trcsms_fabm.F90

    r9330 r9336  
    136136   SUBROUTINE asmdiags_fabm 
    137137      INTEGER :: ji,jj,jk,jkmax 
     138      REAL(wp), DIMENSION(jpi,jpj,jpk) :: pgrow_3d, ploss_3d 
    138139    
    139       PGROW_AVG(:,:) = 1.0 
    140       PLOSS_AVG(:,:) = 1.0 
     140      PGROW_AVG(:,:) = 0.0 
     141      PLOSS_AVG(:,:) = 0.0 
    141142      PHYT_AVG(:,:)  = 0.0 
     143       
     144      pgrow_3d(:,:,:) = fabm_get_bulk_diagnostic_data(model, jp_fabmdia_pgrow) 
     145      ploss_3d(:,:,:) = fabm_get_bulk_diagnostic_data(model, jp_fabmdia_ploss) 
    142146 
    143147      DO jj = 1, jpj 
     
    159163                  &              trn(ji,jj,jk,jp_fabm_p3n) + & 
    160164                  &              trn(ji,jj,jk,jp_fabm_p4n) 
     165               IF ( pgrow_3d(ji,jj,jk) .GT. 0.0 ) THEN 
     166                  PGROW_AVG(ji,jj) = PGROW_AVG(ji,jj) + & 
     167                     &               pgrow_3d(ji,jj,jk) 
     168               ENDIF 
     169               IF ( ploss_3d(ji,jj,jk) .GT. 0.0 ) THEN 
     170                  PLOSS_AVG(ji,jj) = PLOSS_AVG(ji,jj) + & 
     171                     &               ploss_3d(ji,jj,jk) 
     172               ENDIF 
    161173            END DO 
    162174             
    163             PHYT_AVG(ji,jj) = PHYT_AVG(ji,jj) / REAL(jkmax) 
     175            PHYT_AVG(ji,jj)  = PHYT_AVG(ji,jj)  / REAL(jkmax) 
     176            PGROW_AVG(ji,jj) = PGROW_AVG(ji,jj) / REAL(jkmax) 
     177            PLOSS_AVG(ji,jj) = PLOSS_AVG(ji,jj) / REAL(jkmax) 
    164178 
    165179            IF ( hmld_tref(ji,jj) .GT. MLD_MAX(ji,jj) ) THEN 
Note: See TracChangeset for help on using the changeset viewer.