Changeset 8632 for branches/NERC

Timestamp:

2017-10-18T16:13:12+02:00 (7 years ago)

Author:

jpalmier

Message:

JPALM -- Coupled Chl -- set negative values to 0.0

Location:

branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC

Files:

• MEDUSA/bio_medusa_diag_slice.F90 (modified) (1 diff)
• trcrst.F90 (modified) (1 diff)

branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_diag_slice.F90

@@ -73 +73 @@
             IF (chl_out.eq.1) THEN
                !! export and scale surface chl
-               zn_chl_srf(:,:) = (trn(:,:,1,jpchd) + trn(:,:,1,jpchn)) * scl_chl * 1.0E-6  
+               zn_chl_srf(:,:) = MAX( 0.0, (trn(:,:,1,jpchd) + trn(:,:,1,jpchn)) * scl_chl * 1.0E-6 )
                                  !! surf Chl in Kg-chl/m3 as needed for cpl
             ELSEIF (chl_out.eq.2) THEN
                !! export and scale mld chl
-               zn_chl_srf(:,:) = fchl_ml * scl_chl * 1.0E-6  
+               zn_chl_srf(:,:) = MAX( 0.0, fchl_ml * scl_chl * 1.0E-6 )
                                  !! mld Chl in Kg-chl/m3 as needed for cpl
             ENDIF

branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/trcrst.F90

@@ -282 +282 @@
       ELSE
          IF(lwp) WRITE(numout,*) 'Chl surf concentration - setting to zero ...'
-         zn_chl_srf(:,:)  = (trn(:,:,1,jpchn) + trn(:,:,1,jpchd)) * scl_chl * 1.E-6
+         zn_chl_srf(:,:)  = MAX( 0.0, (trn(:,:,1,jpchn) + trn(:,:,1,jpchd)) * scl_chl * 1.E-6 )
       ENDIF
       IF (lk_oasis) THEN