Changeset 8167


Ignore:
Timestamp:
2017-06-13T17:16:50+02:00 (3 years ago)
Author:
dford
Message:

Alterations to 25 hour mean diagnostics for ERSEM.

File:
1 edited

Legend:

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

    r8049 r8167  
    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 
     
    378384            zw3d(:,:,:) = fabm3d_25h(:,:,:,7)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
    379385            CALL iom_put( "pCO2_25h", zw3d  )   ! pCO2 
     386            zw2d(:,:)   = fabm3d_25h(:,:,1,7)*tmask(:,:,1)*0.101325 + zmdi*(1.0-tmask(:,:,1)) 
     387            CALL iom_put( "spCO2_25h", zw2d  )  ! surface pCO2 in Pa 
    380388            zw3d(:,:,:) = fabm3d_25h(:,:,:,8)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
    381389            CALL iom_put( "CHL25h", zw3d  )   ! total Chl 
     
    386394            zw3d(:,:,:) = fabm3d_25h(:,:,:,11)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
    387395            CALL iom_put( "PhytoC25h", zw3d  )   ! total phytoplankton 
     396            zw3d(:,:,:) = fabm3d_25h(:,:,:,11)*tmask(:,:,:)*0.0832590939745 + zmdi*(1.0-tmask(:,:,:)) 
     397            CALL iom_put( "PhyCMol25h", zw3d  )   ! total phytoplankton in mmol 
     398            zw3d(:,:,:) = fabm3d_25h(:,:,:,12)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
     399            CALL iom_put( "ZooC25h", zw3d  )   ! total zooplankton 
     400            zw3d(:,:,:) = fabm3d_25h(:,:,:,12)*tmask(:,:,:)*0.0832590939745 + zmdi*(1.0-tmask(:,:,:)) 
     401            CALL iom_put( "ZooCMol25h", zw3d  )   ! total zooplankton in mmol 
    388402#endif 
    389403 
     
    442456            fabm3d_25h(:,:,:,10)= fabm_get_bulk_diagnostic_data(model, jp_fabm_xeps) 
    443457            fabm3d_25h(:,:,:,11)= fabm_get_bulk_diagnostic_data(model, jp_fabmdia_phytot) 
     458            fabm3d_25h(:,:,:,12)= trn(:,:,:,jp_fabm_m1+jp_fabm_z4c) + & 
     459               &                  trn(:,:,:,jp_fabm_m1+jp_fabm_z5c) + & 
     460               &                  trn(:,:,:,jp_fabm_m1+jp_fabm_z6c) 
    444461#endif 
    445462            cnt_25h = 1 
Note: See TracChangeset for help on using the changeset viewer.