Changeset 8632


Ignore:
Timestamp:
2017-10-18T16:13:12+02:00 (3 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:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_diag_slice.F90

    r8607 r8632  
    7373            IF (chl_out.eq.1) THEN 
    7474               !! export and scale surface chl 
    75                zn_chl_srf(:,:) = (trn(:,:,1,jpchd) + trn(:,:,1,jpchn)) * scl_chl * 1.0E-6   
     75               zn_chl_srf(:,:) = MAX( 0.0, (trn(:,:,1,jpchd) + trn(:,:,1,jpchn)) * scl_chl * 1.0E-6 ) 
    7676                                 !! surf Chl in Kg-chl/m3 as needed for cpl 
    7777            ELSEIF (chl_out.eq.2) THEN 
    7878               !! export and scale mld chl 
    79                zn_chl_srf(:,:) = fchl_ml * scl_chl * 1.0E-6   
     79               zn_chl_srf(:,:) = MAX( 0.0, fchl_ml * scl_chl * 1.0E-6 ) 
    8080                                 !! mld Chl in Kg-chl/m3 as needed for cpl 
    8181            ENDIF 
  • branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/trcrst.F90

    r8607 r8632  
    282282      ELSE 
    283283         IF(lwp) WRITE(numout,*) 'Chl surf concentration - setting to zero ...' 
    284          zn_chl_srf(:,:)  = (trn(:,:,1,jpchn) + trn(:,:,1,jpchd)) * scl_chl * 1.E-6 
     284         zn_chl_srf(:,:)  = MAX( 0.0, (trn(:,:,1,jpchn) + trn(:,:,1,jpchd)) * scl_chl * 1.E-6 ) 
    285285      ENDIF 
    286286      IF (lk_oasis) THEN 
Note: See TracChangeset for help on using the changeset viewer.