New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 9245 for branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_9163/NEMOGCM/NEMO/TOP_SRC/MEDUSA/air_sea.F90 – NEMO

Ignore:
Timestamp:
2018-01-16T16:57:50+01:00 (6 years ago)
Author:
jpalmier
Message:

JPALM -- GMED-371 -- corrections after Richard's review -

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_9163/NEMOGCM/NEMO/TOP_SRC/MEDUSA/air_sea.F90

    r9182 r9245  
    9292# endif 
    9393      USE lib_mpp,           ONLY: ctl_stop 
    94       USE trcstat  
     94      USE trcstat,           ONLY: trc_rst_dia_stat  
    9595 
    9696   !!* Substitution 
     
    157157                  !!!  
    158158                  !!! Jpalm test on atm xCO2 
    159                   IF ( (f_xco2a(ji,jj) > 1500 ).OR.(f_xco2a(ji,jj) < 100 ) ) THEN 
     159                  IF ( (f_xco2a(ji,jj) > 1500.0 ).OR.(f_xco2a(ji,jj) < 100.0 ) ) THEN 
    160160                    IF(lwp) THEN  
    161161                       WRITE(numout,*) ' atm xCO2 = ',f_xco2a(ji,jj),          & 
    162162                                       ' -- ji =', mig(ji),' jj = ', mjg(jj) 
    163                        CALL ctl_stop( 'MEDUSA - Air-Sea :', 'unrealistic atm xCO2 ' ) 
    164                     ENDIF  
     163                    ENDIF 
     164                    CALL ctl_stop( 'MEDUSA - Air-Sea :', 'unrealistic atm xCO2 ' ) 
    165165                 ENDIF  
    166166               ENDIF 
     
    207207               !! Jpalm -- 12-09-2017 -- add extra check after reccurent 
    208208               !!          carbonate failure in the coupled run. 
    209                !!          must be associated to air-sea flux or air xCO2...i 
     209               !!          must be associated to air-sea flux or air xCO2... 
    210210               !!          Check MOCSY inputs 
    211211               IF ( (zsal(ji,jj) > 75.0 ).OR.(zsal(ji,jj) < 0.0 ) .OR.        & 
     
    229229                                       ' kt = ', kt  
    230230                      WRITE(numout,*) 'MEDUSA - Air-Sea INPUT: unrealistic surface Carb. Chemistry' 
    231                       CALL ctl_stop( 'MEDUSA - Air-Sea INPUT: ',             & 
    232                                      'unrealistic surface Carb. Chemistry -- INPUTS' ) 
    233231                  ENDIF      
     232                  CALL ctl_stop( 'MEDUSA - Air-Sea INPUT: ',             & 
     233                                 'unrealistic surface Carb. Chemistry -- INPUTS' ) 
    234234               ENDIF      
    235235               !! 
     
    259259               !! Jpalm -- 12-09-2017 -- add extra check after reccurent 
    260260               !!          carbonate failure in the coupled run. 
    261                !!          must be associated to air-sea flux or air xCO2...i 
    262                !!          Check MOCSY inputs 
     261               !!          must be associated to air-sea flux or air xCO2... 
     262               !!          Check MOCSY outputs 
    263263              IF ( (f_pco2w(ji,jj) > 1.E4 ).OR.(f_pco2w(ji,jj) < 0.0 ) .OR.     & 
    264264                   (f_fco2w(ji,jj) > 1.E4 ).OR.(f_fco2w(ji,jj) < 0.0 ) .OR.     &    
     
    281281                                       ' kt = ', kt      
    282282                     WRITE(numout,*) 'MEDUSA - Air-Sea OUTPUT: unrealistic surface Carb. Chemistry' 
    283                      CALL ctl_stop( 'MEDUSA - Air-Sea OUTPUT: ',            & 
    284                                     'unrealistic surface Carb. Chemistry -- OUTPUTS' ) 
    285283                 ENDIF      
     284                 CALL ctl_stop( 'MEDUSA - Air-Sea OUTPUT: ',            & 
     285                                'unrealistic surface Carb. Chemistry -- OUTPUTS' ) 
    286286              ENDIF      
    287287            ENDIF 
Note: See TracChangeset for help on using the changeset viewer.