Changeset 14586


Ignore:
Timestamp:
2021-03-04T16:29:51+01:00 (4 months ago)
Author:
dford
Message:

Allow just physics to use input chlorophyll.

Location:
branches/UKMO/dev_r5518_GO6_package_FOAMv14_readchl/NEMOGCM/NEMO/TOP_SRC/MEDUSA
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_FOAMv14_readchl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/sms_medusa.F90

    r13725 r14586  
    318318   REAL(wp) ::   rpig       !: chla/chla+phea ratio                     (NAMELIST) 
    319319   LOGICAL  ::   ln_rgb     !: use RGB light scheme rather than 2-band  (NAMELIST) 
     320   LOGICAL  ::   ln_chlbio  !: use input chl for bio rgb 
    320321                                                         
    321322   INTEGER , ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   neln    !: number of levels in the euphotic layer 
  • branches/UKMO/dev_r5518_GO6_package_FOAMv14_readchl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90

    r13725 r14586  
    9494      &    xthetarem,xo2min  
    9595#endif 
    96       NAMELIST/natopt/xkg0,xkr0,xkgp,xkrp,xlg,xlr,rpig,ln_rgb 
     96      NAMELIST/natopt/xkg0,xkr0,xkgp,xkrp,xlg,xlr,rpig,ln_rgb,ln_chlbio 
    9797      INTEGER :: jl, jn 
    9898      INTEGER :: ios                 ! Local integer output status for namelist read 
     
    10981098      rpig  = 0. 
    10991099      ln_rgb = .false. 
     1100      ln_chlbio = .true. 
    11001101 
    11011102      !READ(numnatm,natopt) 
     
    11221123         WRITE(numout,*) ' chla/chla+phea ratio            rpig  = ',rpig 
    11231124         WRITE(numout,*) ' use RGB scheme                ln_rgb  = ',ln_rgb 
     1125         WRITE(numout,*) '                            ln_chlbio  = ',ln_chlbio 
    11241126         WRITE(numout,*) ' ' 
    11251127 
  • branches/UKMO/dev_r5518_GO6_package_FOAMv14_readchl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcopt_medusa.F90

    r14388 r14586  
    126126            DO jj = 1, jpj 
    127127               DO ji = 1, jpi 
    128                   !totchl = trn(ji,jj,jk,jpchn) + trn(ji,jj,jk,jpchd) 
    129                   totchl = chl_for_qsr(ji,jj,jk) 
     128                  IF ( ln_chlbio ) THEN 
     129                     totchl = chl_for_qsr(ji,jj,jk) 
     130                  ELSE 
     131                     totchl = trn(ji,jj,jk,jpchn) + trn(ji,jj,jk,jpchd) 
     132                  ENDIF 
    130133                  totchl = MIN( 10. , MAX( 0.05, totchl ) ) 
    131134                  irgb = NINT( 41 + 20.* LOG10( totchl ) + 1.e-15 ) 
Note: See TracChangeset for help on using the changeset viewer.