Changeset 15685

Timestamp:

2022-02-02T17:14:34+01:00 (2 years ago)

Author:

dancopsey

Message:

Change chlorophyll ancil to be in the same units

File:

• NEMO/branches/UKMO/NEMO_4.0.4_change_chlorophyll/src/OCE/TRA/traqsr.F90 (modified) (2 diffs)

NEMO/branches/UKMO/NEMO_4.0.4_change_chlorophyll/src/OCE/TRA/traqsr.F90

@@ -169 +169 @@
                DO jj = 2, jpjm1                       ! Separation in R-G-B depending of the surface Chl
                   DO ji = fs_2, fs_jpim1
-                     zchl    = MIN( 10. , MAX( 0.03, sf_chl(1)%fnow(ji,jj,1) ) )
-                     zCtot   = 40.6  * zchl**0.459
-                     zze     = 568.2 * zCtot**(-0.746)
-                     IF( zze > 102. ) zze = 200.0 * zCtot**(-0.293)
-                     zpsi    = gdepw_n(ji,jj,jk) / zze
+
+                     ! Commenting out this bit of code as I want to load in chlorophyll concentrations
+                     ! with the same units as rn_chl_conc
+                     ! zchl    = MIN( 10. , MAX( 0.03, sf_chl(1)%fnow(ji,jj,1) ) )
+                     ! zCtot   = 40.6  * zchl**0.459
+                     ! zze     = 568.2 * zCtot**(-0.746)
+                     ! IF( zze > 102. ) zze = 200.0 * zCtot**(-0.293)
+                     ! zpsi    = gdepw_n(ji,jj,jk) / zze
                      !
-                     zlogc   = LOG( zchl )
-                     zlogc2  = zlogc * zlogc
-                     zlogc3  = zlogc * zlogc * zlogc
-                     zCb     = 0.768 + 0.087 * zlogc - 0.179 * zlogc2 - 0.025 * zlogc3
-                     zCmax   = 0.299 - 0.289 * zlogc + 0.579 * zlogc2
-                     zpsimax = 0.6   - 0.640 * zlogc + 0.021 * zlogc2 + 0.115 * zlogc3
-                     zdelpsi = 0.710 + 0.159 * zlogc + 0.021 * zlogc2
-                     zCze    = 1.12  * (zchl)**0.803 
+                     ! zlogc   = LOG( zchl )
+                     ! zlogc2  = zlogc * zlogc
+                     ! zlogc3  = zlogc * zlogc * zlogc
+                     ! zCb     = 0.768 + 0.087 * zlogc - 0.179 * zlogc2 - 0.025 * zlogc3
+                     ! zCmax   = 0.299 - 0.289 * zlogc + 0.579 * zlogc2
+                     ! zpsimax = 0.6   - 0.640 * zlogc + 0.021 * zlogc2 + 0.115 * zlogc3
+                     ! zdelpsi = 0.710 + 0.159 * zlogc + 0.021 * zlogc2
+                     ! zCze    = 1.12  * (zchl)**0.803 
                      !
-                     zchl3d(ji,jj,jk) = zCze * ( zCb + zCmax * EXP( -( (zpsi - zpsimax) / zdelpsi )**2 ) )
+                     ! zchl3d(ji,jj,jk) = zCze * ( zCb + zCmax * EXP( -( (zpsi - zpsimax) / zdelpsi )**2 ) )
+
+                     zchl3d(ji,jj,jk) = MIN( 50. , MAX( 0.01, sf_chl(1)%fnow(ji,jj,1) ) )
                   END DO
                   !
@@ -404 +409 @@
          ENDIF
          IF( nqsr == np_RGB ) THEN                 ! constant Chl
-            IF(lwp) WRITE(numout,*) '   ==>>>   Constant Chlorophyll concentration = 0.05'
+            IF(lwp) WRITE(numout,*) '   ==>>>   Constant Chlorophyll concentration = ',rn_chl_conc
          ENDIF
          !