Changeset 9130


Ignore:
Timestamp:
2017-12-19T14:15:21+01:00 (3 years ago)
Author:
jpalmier
Message:

JPALM — microboil-carb fail Check and report

Location:
branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_8356/NEMOGCM/NEMO
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_8356/NEMOGCM/NEMO/OPA_SRC/stpctl.F90

    r9073 r9130  
    157157      ztmin = 100.e0  ! arbitrary high minimum value 
    158158      DO jj = 2, jpjm1 
    159          DO ji = 1, jpi 
     159         DO ji = 2, jpim1 
    160160            IF( tmask(ji,jj,1) == 1) THEN 
    161161               ztmax = MAX(ztmax,tsn(ji,jj,1,jp_tem))     ! find local maximum 
  • branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_8356/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

    r9073 r9130  
    821821      !!========================== 
    822822      !! DIC Check 
    823       zmax =  -5.e0  ! arbitrary  low maximum value 
     823      zmax =  -5.0  ! arbitrary  low maximum value 
    824824      zmin =  4.0E4  ! arbitrary high minimum value 
    825825      DO jj = 2, jpjm1 
    826          DO ji = 1, jpi 
     826         DO ji = 2,jpim1 
    827827            IF( tmask(ji,jj,1) == 1) THEN 
    828828               zmax = MAX(zmax,zdic(ji,jj))     ! find local maximum 
     
    834834      IF( lk_mpp )   CALL mpp_min( zmin )       ! min over the global domain 
    835835      ! 
    836       IF( zmax > 4.0E3) THEN  ! we've got a problem 
     836      IF( zmax .GT. 4.0E3) THEN  ! we've got a problem 
    837837         IF (lk_mpp) THEN 
    838838            CALL mpp_maxloc ( zdic(:,:),tmask(:,:,1), zmax, ii,ij ) 
    839839         ELSE 
    840             ilocs = MAXLOC( zdic(:,:), mask = tmask(:,:,1) == 1.e0 ) 
     840            ilocs = MAXLOC( zdic(:,:), mask = tmask(:,:,1) == 1. ) 
    841841            ii = ilocs(1) + nimpp - 1 
    842842            ij = ilocs(2) + njmpp - 1 
     
    851851      ENDIF 
    852852      ! 
    853       IF( zmin <= 0.0) THEN  ! we've got a problem 
     853      IF( zmin .LE. 0.0) THEN  ! we've got a problem 
    854854         IF (lk_mpp) THEN 
    855855            CALL mpp_minloc ( zdic(:,:),tmask(:,:,1), zmin, ii,ij ) 
    856856         ELSE 
    857             ilocs = MINLOC( zdic(:,:), mask = tmask(:,:,1) == 1.e0 ) 
     857            ilocs = MINLOC( zdic(:,:), mask = tmask(:,:,1) == 1. ) 
    858858            ii = ilocs(1) + nimpp - 1 
    859859            ij = ilocs(2) + njmpp - 1 
     
    870870      !!========================== 
    871871      !! ALKALINITY Check 
    872       zmax =  -5.e0  ! arbitrary  low maximum value 
     872      zmax =  -5.0  ! arbitrary  low maximum value 
    873873      zmin =  4.0E4  ! arbitrary high minimum value 
    874874      DO jj = 2, jpjm1 
    875          DO ji = 1, jpi 
     875         DO ji = 2,jpim1 
    876876            IF( tmask(ji,jj,1) == 1) THEN 
    877877               zmax = MAX(zmax,zalk(ji,jj))     ! find local maximum 
     
    883883      IF( lk_mpp )   CALL mpp_min( zmin )       ! min over the global domain 
    884884      ! 
    885       IF( zmax > 4.0E3) THEN  ! we've got a problem 
     885      IF( zmax .GT. 4.0E3) THEN  ! we've got a problem 
    886886         IF (lk_mpp) THEN 
    887887            CALL mpp_maxloc ( zalk(:,:),tmask(:,:,1), zmax, ii,ij ) 
    888888         ELSE 
    889             ilocs = MAXLOC( zalk(:,:), mask = tmask(:,:,1) == 1.e0 ) 
     889            ilocs = MAXLOC( zalk(:,:), mask = tmask(:,:,1) == 1. ) 
    890890            ii = ilocs(1) + nimpp - 1 
    891891            ij = ilocs(2) + njmpp - 1 
     
    900900      ENDIF 
    901901      ! 
    902       IF( zmin <= 0.0) THEN  ! we've got a problem 
     902      IF( zmin .LE. 0.0) THEN  ! we've got a problem 
    903903         IF (lk_mpp) THEN 
    904904            CALL mpp_minloc ( zalk(:,:),tmask(:,:,1), zmin, ii,ij ) 
    905905         ELSE 
    906             ilocs = MINLOC( zalk(:,:), mask = tmask(:,:,1) == 1.e0 ) 
     906            ilocs = MINLOC( zalk(:,:), mask = tmask(:,:,1) == 1. ) 
    907907            ii = ilocs(1) + nimpp - 1 
    908908            ij = ilocs(2) + njmpp - 1 
Note: See TracChangeset for help on using the changeset viewer.