Changeset 8173


Ignore:
Timestamp:
2017-06-14T15:50:33+02:00 (3 years ago)
Author:
dford
Message:

Add 25 hour mean zooplankton diagnostic.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/CO6_KD490_amm7_oper_fabm/NEMOGCM/NEMO/OPA_SRC/DIA/dia25h.F90

    r8049 r8173  
    195195      fabm3d_25h(:,:,:,10)= fabm_get_bulk_diagnostic_data(model, jp_fabm_xeps)  
    196196      fabm3d_25h(:,:,:,11)= fabm_get_bulk_diagnostic_data(model, jp_fabmdia_phytot)  
     197      fabm3d_25h(:,:,:,12)= trn(:,:,:,jp_fabm_m1+jp_fabm_z4c) + & 
     198         &                  trn(:,:,:,jp_fabm_m1+jp_fabm_z5c) + & 
     199         &                  trn(:,:,:,jp_fabm_m1+jp_fabm_z6c) 
    197200#endif 
    198201#if defined key_lim3 || defined key_lim2 
     
    307310         fabm3d_25h(:,:,:,10)= fabm3d_25h(:,:,:,10) + fabm_get_bulk_diagnostic_data(model, jp_fabm_xeps) ! light attenuation 
    308311         fabm3d_25h(:,:,:,11)= fabm3d_25h(:,:,:,11) + fabm_get_bulk_diagnostic_data(model, jp_fabmdia_phytot) ! total phytoplankton 
     312         fabm3d_25h(:,:,:,12)= fabm3d_25h(:,:,:,12) + trn(:,:,:,jp_fabm_m1+jp_fabm_z4c) + & 
     313            &                                         trn(:,:,:,jp_fabm_m1+jp_fabm_z5c) + & 
     314            &                                         trn(:,:,:,jp_fabm_m1+jp_fabm_z6c)      ! total zooplankton 
    309315#endif 
    310316         cnt_25h = cnt_25h + 1 
     
    386392            zw3d(:,:,:) = fabm3d_25h(:,:,:,11)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
    387393            CALL iom_put( "PhytoC25h", zw3d  )   ! total phytoplankton 
     394            zw3d(:,:,:) = fabm3d_25h(:,:,:,11)*tmask(:,:,:)*0.0832590939745 + zmdi*(1.0-tmask(:,:,:)) 
     395            CALL iom_put( "PhyCMol25h", zw3d  )   ! total phytoplankton in mmol 
     396            zw3d(:,:,:) = fabm3d_25h(:,:,:,12)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
     397            CALL iom_put( "ZooC25h", zw3d  )   ! total zooplankton 
     398            zw3d(:,:,:) = fabm3d_25h(:,:,:,12)*tmask(:,:,:)*0.0832590939745 + zmdi*(1.0-tmask(:,:,:)) 
     399            CALL iom_put( "ZooCMol25h", zw3d  )   ! total zooplankton in mmol 
    388400#endif 
    389401 
     
    442454            fabm3d_25h(:,:,:,10)= fabm_get_bulk_diagnostic_data(model, jp_fabm_xeps) 
    443455            fabm3d_25h(:,:,:,11)= fabm_get_bulk_diagnostic_data(model, jp_fabmdia_phytot) 
     456            fabm3d_25h(:,:,:,12)= trn(:,:,:,jp_fabm_m1+jp_fabm_z4c) + & 
     457               &                  trn(:,:,:,jp_fabm_m1+jp_fabm_z5c) + & 
     458               &                  trn(:,:,:,jp_fabm_m1+jp_fabm_z6c) 
    444459#endif 
    445460            cnt_25h = 1 
Note: See TracChangeset for help on using the changeset viewer.