Changeset 14586

Timestamp:

2021-03-04T16:29:51+01:00 (3 years 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:

• sms_medusa.F90 (modified) (1 diff)
• trcnam_medusa.F90 (modified) (3 diffs)
• trcopt_medusa.F90 (modified) (1 diff)

branches/UKMO/dev_r5518_GO6_package_FOAMv14_readchl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/sms_medusa.F90

@@ -318 +318 @@
    REAL(wp) ::   rpig       !: chla/chla+phea ratio                     (NAMELIST)
    LOGICAL  ::   ln_rgb     !: use RGB light scheme rather than 2-band  (NAMELIST)
+   LOGICAL  ::   ln_chlbio  !: use input chl for bio rgb
                                                         
    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

@@ -94 +94 @@
       &    xthetarem,xo2min 
 #endif
-      NAMELIST/natopt/xkg0,xkr0,xkgp,xkrp,xlg,xlr,rpig,ln_rgb
+      NAMELIST/natopt/xkg0,xkr0,xkgp,xkrp,xlg,xlr,rpig,ln_rgb,ln_chlbio
       INTEGER :: jl, jn
       INTEGER :: ios                 ! Local integer output status for namelist read
@@ -1098 +1098 @@
       rpig  = 0.
       ln_rgb = .false.
+      ln_chlbio = .true.
 
       !READ(numnatm,natopt)
@@ -1122 +1123 @@
          WRITE(numout,*) ' chla/chla+phea ratio            rpig  = ',rpig
          WRITE(numout,*) ' use RGB scheme                ln_rgb  = ',ln_rgb
+         WRITE(numout,*) '                            ln_chlbio  = ',ln_chlbio
          WRITE(numout,*) ' '
 

branches/UKMO/dev_r5518_GO6_package_FOAMv14_readchl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcopt_medusa.F90

@@ -126 +126 @@
             DO jj = 1, jpj
                DO ji = 1, jpi
-                  !totchl = trn(ji,jj,jk,jpchn) + trn(ji,jj,jk,jpchd)
-                  totchl = chl_for_qsr(ji,jj,jk)
+                  IF ( ln_chlbio ) THEN
+                     totchl = chl_for_qsr(ji,jj,jk)
+                  ELSE
+                     totchl = trn(ji,jj,jk,jpchn) + trn(ji,jj,jk,jpchd)
+                  ENDIF
                   totchl = MIN( 10. , MAX( 0.05, totchl ) )
                   irgb = NINT( 41 + 20.* LOG10( totchl ) + 1.e-15 )