New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 8648 for branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc_v2/NEMOGCM/NEMO/OPA_SRC/ASM/asmlogchlbal_hadocc.F90 – NEMO

Ignore:
Timestamp:
2017-10-20T15:31:59+02:00 (7 years ago)
Author:
dford
Message:

Change source of MLD for direct initialisation of logchl and pCO2 to work.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc_v2/NEMOGCM/NEMO/OPA_SRC/ASM/asmlogchlbal_hadocc.F90

    r8495 r8648  
    6969 
    7070   SUBROUTINE asm_logchl_bal_hadocc( logchl_bkginc, aincper, mld_choice_bgc, & 
    71       &                              k_maxchlinc, ld_logchlbal,              & 
     71      &                              k_maxchlinc, ld_logchlbal, ld_asmdin,   & 
    7272      &                              pgrow_avg_bkg, ploss_avg_bkg,           & 
    7373      &                              phyt_avg_bkg, mld_max_bkg,              & 
    … …  
    9393      REAL(wp), INTENT(in   )                               :: k_maxchlinc    ! Max chl increment 
    9494      LOGICAL,  INTENT(in   )                               :: ld_logchlbal   ! Balancing y/n 
     95      LOGICAL,  INTENT(in   )                               :: ld_asmdin      ! Direct initialisation y/n 
    9596      REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj)           :: pgrow_avg_bkg  ! Avg phyto growth 
    9697      REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj)           :: ploss_avg_bkg  ! Avg phyto loss 
    … …  
    136137       
    137138      ! Select mixed layer 
    138       SELECT CASE( mld_choice_bgc ) 
    139       CASE ( 1 )                   ! Turbocline/mixing depth [W points] 
    140          zmld(:,:) = hmld(:,:) 
    141       CASE ( 2 )                   ! Density criterion (0.01 kg/m^3 change from 10m) [W points] 
    142          zmld(:,:) = hmlp(:,:) 
    143       CASE ( 3 )                   ! Kara MLD [Interpolated] 
     139      IF ( ld_asmdin ) THEN 
     140         CALL ctl_warn( ' Doing direct initialisation of HadOCC with chlorophyll assimilation', & 
     141            &           ' Mixed layer depth taken to be background maximum mld_max_bkg' ) 
     142         zmld(:,:) = mld_max_bkg(:,:) 
     143      ELSE 
     144         SELECT CASE( mld_choice_bgc ) 
     145         CASE ( 1 )                   ! Turbocline/mixing depth [W points] 
     146            zmld(:,:) = hmld(:,:) 
     147         CASE ( 2 )                   ! Density criterion (0.01 kg/m^3 change from 10m) [W points] 
     148            zmld(:,:) = hmlp(:,:) 
     149         CASE ( 3 )                   ! Kara MLD [Interpolated] 
    144150#if defined key_karaml 
    145          IF ( ln_kara ) THEN 
    146             zmld(:,:) = hmld_kara(:,:) 
    147          ELSE 
     151            IF ( ln_kara ) THEN 
     152               zmld(:,:) = hmld_kara(:,:) 
     153            ELSE 
     154               CALL ctl_stop( ' Kara mixed layer requested for LogChl assimilation,', & 
     155                  &           ' but ln_kara=.false.' ) 
     156            ENDIF 
     157#else 
    148158            CALL ctl_stop( ' Kara mixed layer requested for LogChl assimilation,', & 
    149                &           ' but ln_kara=.false.' ) 
    150          ENDIF 
    151 #else 
    152          CALL ctl_stop( ' Kara mixed layer requested for LogChl assimilation,', & 
    153             &           ' but is not defined' ) 
     159               &           ' but is not defined' ) 
    154160#endif 
    155       CASE ( 4 )                   ! Temperature criterion (0.2 K change from surface) [T points] 
    156          !zmld(:,:) = hmld_tref(:,:) 
    157          CALL ctl_stop( ' hmld_tref mixed layer requested for LogChl assimilation,', & 
    158             &           ' but is not available in this version' ) 
    159       CASE ( 5 )                   ! Density criterion (0.01 kg/m^3 change from 10m) [T points] 
    160          zmld(:,:) = hmlpt(:,:) 
    161       END SELECT 
     161         CASE ( 4 )                   ! Temperature criterion (0.2 K change from surface) [T points] 
     162            !zmld(:,:) = hmld_tref(:,:) 
     163            CALL ctl_stop( ' hmld_tref mixed layer requested for LogChl assimilation,', & 
     164               &           ' but is not available in this version' ) 
     165         CASE ( 5 )                   ! Density criterion (0.01 kg/m^3 change from 10m) [T points] 
     166            zmld(:,:) = hmlpt(:,:) 
     167         END SELECT 
     168      ENDIF 
    162169       
    163170      IF ( ld_logchlbal ) THEN   ! Nitrogen balancing 
Note: See TracChangeset for help on using the changeset viewer.