Ignore:
Timestamp:
2017-09-04T18:50:10+02:00 (3 years ago)
Author:
dford
Message:

Merge in changes from dev_r5518_GO6_package_asm_surf_bgc, and adapt to the updated MEDUSA structure.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc_v2/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.F90

    r8442 r8495  
    4545                                   f3_co3, f3_h2co3, f3_hco3,           & 
    4646                                   f3_omarg, f3_omcal, f3_pH,           & 
     47# if defined key_foam_medusa 
     48                                   mld_max, pgrow_avg,                  & 
     49                                   ploss_avg, phyt_avg,                 & 
     50# endif 
    4751                                   za_sed_c, za_sed_ca, za_sed_fe,      & 
    4852                                   za_sed_n, za_sed_si,                 & 
     
    5357      USE trc,               ONLY: med_diag, nittrc000  
    5458      USE trcnam_trp,        ONLY: ln_trcadv_cen2, ln_trcadv_tvd 
     59# if defined key_foam_medusa 
     60      USE zdfmxl,            ONLY: hmld 
     61# endif 
    5562  
    5663      !! time (integer timestep) 
     
    218225      endif 
    219226 
     227# if defined key_foam_medusa 
     228!!---------------------------------------------------------------------- 
     229!! Diagnostics required for ocean colour assimilation: 
     230!! Mixed layer average phytoplankton growth, loss and concentration 
     231!! Maximum mixed layer depth 
     232!!---------------------------------------------------------------------- 
     233!! 
     234      DO jj = 2,jpjm1 
     235         DO ji = 2,jpim1 
     236            pgrow_avg(ji,jj) = pgrow_avg(ji,jj) / hmld(ji,jj) 
     237            ploss_avg(ji,jj) = ploss_avg(ji,jj) / hmld(ji,jj) 
     238            phyt_avg(ji,jj)  = phyt_avg(ji,jj)  / hmld(ji,jj) 
     239            IF ( hmld(ji,jj) .GT. mld_max(ji,jj) ) THEN 
     240               mld_max(ji,jj) = hmld(ji,jj) 
     241            ENDIF 
     242         END DO 
     243      END DO 
     244# endif 
     245 
    220246#  if defined key_debug_medusa 
    221247         !! AXY (12/07/17) 
Note: See TracChangeset for help on using the changeset viewer.