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 7798 for branches/UKMO/dev_r5518_v3.6_asm_nemovar_community_maxchlinc/NEMOGCM/NEMO/OPA_SRC/ASM/asmlogchlbal_ersem.F90 – NEMO

Ignore:
Timestamp:
2017-03-16T11:00:52+01:00 (7 years ago)
Author:
dford
Message:

Update variable naming.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_v3.6_asm_nemovar_community_maxchlinc/NEMOGCM/NEMO/OPA_SRC/ASM/asmlogchlbal_ersem.F90

    r7791 r7798  
    2929CONTAINS 
    3030 
    31    SUBROUTINE asm_logchl_bal_ersem( ln_logchlpftinc, npfts, mld_choice_bgc, & 
    32       &                             rn_maxchlinc, logchl_bkginc, logchl_balinc ) 
     31   SUBROUTINE asm_logchl_bal_ersem( ld_logchlpftinc, npfts, mld_choice_bgc, & 
     32      &                             k_maxchlinc, logchl_bkginc, logchl_balinc ) 
    3333      !!--------------------------------------------------------------------------- 
    3434      !!                    ***  ROUTINE asm_logchl_bal_ersem  *** 
     
    4646      !!--------------------------------------------------------------------------- 
    4747      !! 
    48       LOGICAL,  INTENT(in   )                               :: ln_logchlpftinc 
     48      LOGICAL,  INTENT(in   )                               :: ld_logchlpftinc 
    4949      INTEGER,  INTENT(in   )                               :: npfts 
    5050      INTEGER,  INTENT(in   )                               :: mld_choice_bgc 
    51       REAL(wp), INTENT(in   )                               :: rn_maxchlinc 
     51      REAL(wp), INTENT(in   )                               :: k_maxchlinc 
    5252      REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj,npfts)     :: logchl_bkginc 
    5353      REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk,jptra) :: logchl_balinc 
     
    6666      ! 2) Take 10^log10(analysis) to get analysis 
    6767      ! 3) Subtract background from analysis to get chl incs 
    68       ! If rn_maxchlinc > 0 then cap total absolute chlorophyll increment at that value 
     68      ! If k_maxchlinc > 0 then cap total absolute chlorophyll increment at that value 
    6969      ! 
    7070      ! Only apply increments if all of Chl1-4 background values are > 0 
     
    7272      ! that something's going wrong which the assimilation might make worse 
    7373      ! 
    74       IF ( ln_logchlpftinc ) THEN 
     74      IF ( ld_logchlpftinc ) THEN 
    7575         ! 
    7676         ! Assimilating separate PFTs, so separately transform each from LogChl to Chl 
     
    105105                        &                                             trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) 
    106106                  ENDIF 
    107                   IF (rn_maxchlinc > 0.0) THEN 
     107                  IF (k_maxchlinc > 0.0) THEN 
    108108                     chl_inc = logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) + & 
    109109                               logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) + & 
    110110                               logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) + & 
    111111                               logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) 
    112                      IF ( ABS(chl_inc) > rn_maxchlinc ) THEN 
    113                         chl_tot = ABS(chl_inc) / rn_maxchlinc 
     112                     IF ( ABS(chl_inc) > k_maxchlinc ) THEN 
     113                        chl_tot = ABS(chl_inc) / k_maxchlinc 
    114114                        logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) = logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) / chl_tot 
    115115                        logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) = logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) / chl_tot 
     
    139139                     &      trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) + trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) 
    140140                  chl_inc = 10**( LOG10( chl_tot ) + logchl_bkginc(ji,jj,1) ) - chl_tot 
    141                   IF (rn_maxchlinc > 0.0) THEN 
    142                      chl_inc = MAX( -1.0 * rn_maxchlinc, MIN( chl_inc, rn_maxchlinc ) ) 
     141                  IF (k_maxchlinc > 0.0) THEN 
     142                     chl_inc = MAX( -1.0 * k_maxchlinc, MIN( chl_inc, k_maxchlinc ) ) 
    143143                  ENDIF 
    144144                  logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) = chl_inc * trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) / chl_tot 
     
    210210   !!---------------------------------------------------------------------- 
    211211CONTAINS 
    212    SUBROUTINE asm_logchl_bal_ersem( ln_logchlpftinc, npfts, mld_choice_bgc, & 
    213       &                             logchl_bkginc, logchl_balinc ) 
    214       LOGICAL :: ln_logchlpftinc 
     212   SUBROUTINE asm_logchl_bal_ersem( ld_logchlpftinc, npfts, mld_choice_bgc, & 
     213      &                             k_maxchlinc, logchl_bkginc, logchl_balinc ) 
     214      LOGICAL :: ld_logchlpftinc 
    215215      INTEGER :: npfts 
    216216      INTEGER :: mld_choice_bgc 
     217      REAL    :: k_maxchlinc 
    217218      REAL    :: logchl_bkginc(:,:,:) 
    218219      REAL    :: logchl_balinc(:,:,:,:) 
Note: See TracChangeset for help on using the changeset viewer.