Changeset 7045


Ignore:
Timestamp:
2016-10-19T17:46:53+02:00 (4 years ago)
Author:
dford
Message:

Update default MLD definition for BGC assimilation, and temporarily comment out code that was causing issues.

Location:
branches/UKMO/dev_r5518_v3.4_asm_nemovar_community_bgc_ersem/NEMOGCM/NEMO/OPA_SRC/ASM
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_v3.4_asm_nemovar_community_bgc_ersem/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

    r6983 r7045  
    188188                                      !: 4) hmld_tref - Temperature criterion (0.2 K change from surface) [T points] 
    189189 
    190    INTEGER :: mld_choice_bgc    = 5   !: choice of mld criteria to use for physics assimilation 
     190   INTEGER :: mld_choice_bgc    = 4   !: choice of mld criteria to use for physics assimilation 
    191191                                      !: 1) hmld      - Turbocline/mixing depth                           [W points] 
    192192                                      !: 2) hmlp      - Density criterion (0.01 kg/m^3 change from 10m)   [W points] 
  • branches/UKMO/dev_r5518_v3.4_asm_nemovar_community_bgc_ersem/NEMOGCM/NEMO/OPA_SRC/ASM/asmlogchlbal_ersem.F90

    r6983 r7045  
    206206      ! Maintain existing background ratios 
    207207      ! 
    208       DO jk = 1, jpk 
    209          DO jj = 1, jpj 
    210             DO ji = 1, jpi 
    211                IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) > 0.0 ) THEN 
    212                   logchl_balinc_ersem_p1c(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk)   * & 
    213                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1c) / & 
    214                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 
    215                   logchl_balinc_ersem_p1n(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk)   * & 
    216                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) / & 
    217                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 
    218                   logchl_balinc_ersem_p1p(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk)   * & 
    219                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1p) / & 
    220                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 
    221                   logchl_balinc_ersem_p1s(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk)   * & 
    222                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1s) / & 
    223                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 
    224                ENDIF 
    225                IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) > 0.0 ) THEN 
    226                   logchl_balinc_ersem_p2c(ji,jj,jk) = logchl_balinc_ersem_chl2(ji,jj,jk)   * & 
    227                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p2c) / & 
    228                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) 
    229                   logchl_balinc_ersem_p2n(ji,jj,jk) = logchl_balinc_ersem_chl2(ji,jj,jk)   * & 
    230                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) / & 
    231                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) 
    232                   logchl_balinc_ersem_p2p(ji,jj,jk) = logchl_balinc_ersem_chl2(ji,jj,jk)   * & 
    233                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p2p) / & 
    234                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) 
    235                ENDIF 
    236                IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) > 0.0 ) THEN 
    237                   logchl_balinc_ersem_p3c(ji,jj,jk) = logchl_balinc_ersem_chl3(ji,jj,jk)   * & 
    238                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p3c) / & 
    239                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) 
    240                   logchl_balinc_ersem_p3n(ji,jj,jk) = logchl_balinc_ersem_chl3(ji,jj,jk)   * & 
    241                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) / & 
    242                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) 
    243                   logchl_balinc_ersem_p3p(ji,jj,jk) = logchl_balinc_ersem_chl3(ji,jj,jk)   * & 
    244                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p3p) / & 
    245                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) 
    246                ENDIF 
    247                IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) > 0.0 ) THEN 
    248                   logchl_balinc_ersem_p4c(ji,jj,jk) = logchl_balinc_ersem_chl4(ji,jj,jk)   * & 
    249                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p4c) / & 
    250                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) 
    251                   logchl_balinc_ersem_p4n(ji,jj,jk) = logchl_balinc_ersem_chl4(ji,jj,jk)   * & 
    252                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) / & 
    253                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) 
    254                   logchl_balinc_ersem_p4p(ji,jj,jk) = logchl_balinc_ersem_chl4(ji,jj,jk)   * & 
    255                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p4p) / & 
    256                      &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) 
    257                ENDIF 
    258             END DO 
    259          END DO 
    260       END DO 
     208      !DO jk = 1, jpk 
     209      !   DO jj = 1, jpj 
     210      !      DO ji = 1, jpi 
     211      !         IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) > 0.0 ) THEN 
     212      !            logchl_balinc_ersem_p1c(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk)   * & 
     213      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1c) / & 
     214      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 
     215      !            logchl_balinc_ersem_p1n(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk)   * & 
     216      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) / & 
     217      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 
     218      !            logchl_balinc_ersem_p1p(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk)   * & 
     219      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1p) / & 
     220      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 
     221      !            logchl_balinc_ersem_p1s(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk)   * & 
     222      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1s) / & 
     223      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 
     224      !         ENDIF 
     225      !         IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) > 0.0 ) THEN 
     226      !            logchl_balinc_ersem_p2c(ji,jj,jk) = logchl_balinc_ersem_chl2(ji,jj,jk)   * & 
     227      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p2c) / & 
     228      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) 
     229      !            logchl_balinc_ersem_p2n(ji,jj,jk) = logchl_balinc_ersem_chl2(ji,jj,jk)   * & 
     230      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) / & 
     231      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) 
     232      !            logchl_balinc_ersem_p2p(ji,jj,jk) = logchl_balinc_ersem_chl2(ji,jj,jk)   * & 
     233      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p2p) / & 
     234      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) 
     235      !         ENDIF 
     236      !         IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) > 0.0 ) THEN 
     237      !            logchl_balinc_ersem_p3c(ji,jj,jk) = logchl_balinc_ersem_chl3(ji,jj,jk)   * & 
     238      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p3c) / & 
     239      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) 
     240      !            logchl_balinc_ersem_p3n(ji,jj,jk) = logchl_balinc_ersem_chl3(ji,jj,jk)   * & 
     241      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) / & 
     242      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) 
     243      !            logchl_balinc_ersem_p3p(ji,jj,jk) = logchl_balinc_ersem_chl3(ji,jj,jk)   * & 
     244      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p3p) / & 
     245      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) 
     246      !         ENDIF 
     247      !         IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) > 0.0 ) THEN 
     248      !            logchl_balinc_ersem_p4c(ji,jj,jk) = logchl_balinc_ersem_chl4(ji,jj,jk)   * & 
     249      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p4c) / & 
     250      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) 
     251      !            logchl_balinc_ersem_p4n(ji,jj,jk) = logchl_balinc_ersem_chl4(ji,jj,jk)   * & 
     252      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) / & 
     253      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) 
     254      !            logchl_balinc_ersem_p4p(ji,jj,jk) = logchl_balinc_ersem_chl4(ji,jj,jk)   * & 
     255      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p4p) / & 
     256      !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) 
     257      !         ENDIF 
     258      !      END DO 
     259      !   END DO 
     260      !END DO 
    261261 
    262262   END SUBROUTINE asm_logchl_bal_ersem 
Note: See TracChangeset for help on using the changeset viewer.