Changeset 9245 for branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_9163/NEMOGCM/NEMO/TOP_SRC/MEDUSA/air_sea.F90
- Timestamp:
- 2018-01-16T16:57:50+01:00 (6 years ago)
- 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 92 92 # endif 93 93 USE lib_mpp, ONLY: ctl_stop 94 USE trcstat 94 USE trcstat, ONLY: trc_rst_dia_stat 95 95 96 96 !!* Substitution … … 157 157 !!! 158 158 !!! Jpalm test on atm xCO2 159 IF ( (f_xco2a(ji,jj) > 1500 ).OR.(f_xco2a(ji,jj) < 100 ) ) THEN159 IF ( (f_xco2a(ji,jj) > 1500.0 ).OR.(f_xco2a(ji,jj) < 100.0 ) ) THEN 160 160 IF(lwp) THEN 161 161 WRITE(numout,*) ' atm xCO2 = ',f_xco2a(ji,jj), & 162 162 ' -- ji =', mig(ji),' jj = ', mjg(jj) 163 CALL ctl_stop( 'MEDUSA - Air-Sea :', 'unrealistic atm xCO2 ' )164 ENDIF163 ENDIF 164 CALL ctl_stop( 'MEDUSA - Air-Sea :', 'unrealistic atm xCO2 ' ) 165 165 ENDIF 166 166 ENDIF … … 207 207 !! Jpalm -- 12-09-2017 -- add extra check after reccurent 208 208 !! carbonate failure in the coupled run. 209 !! must be associated to air-sea flux or air xCO2... i209 !! must be associated to air-sea flux or air xCO2... 210 210 !! Check MOCSY inputs 211 211 IF ( (zsal(ji,jj) > 75.0 ).OR.(zsal(ji,jj) < 0.0 ) .OR. & … … 229 229 ' kt = ', kt 230 230 WRITE(numout,*) 'MEDUSA - Air-Sea INPUT: unrealistic surface Carb. Chemistry' 231 CALL ctl_stop( 'MEDUSA - Air-Sea INPUT: ', &232 'unrealistic surface Carb. Chemistry -- INPUTS' )233 231 ENDIF 232 CALL ctl_stop( 'MEDUSA - Air-Sea INPUT: ', & 233 'unrealistic surface Carb. Chemistry -- INPUTS' ) 234 234 ENDIF 235 235 !! … … 259 259 !! Jpalm -- 12-09-2017 -- add extra check after reccurent 260 260 !! carbonate failure in the coupled run. 261 !! must be associated to air-sea flux or air xCO2... i262 !! Check MOCSY inputs261 !! must be associated to air-sea flux or air xCO2... 262 !! Check MOCSY outputs 263 263 IF ( (f_pco2w(ji,jj) > 1.E4 ).OR.(f_pco2w(ji,jj) < 0.0 ) .OR. & 264 264 (f_fco2w(ji,jj) > 1.E4 ).OR.(f_fco2w(ji,jj) < 0.0 ) .OR. & … … 281 281 ' kt = ', kt 282 282 WRITE(numout,*) 'MEDUSA - Air-Sea OUTPUT: unrealistic surface Carb. Chemistry' 283 CALL ctl_stop( 'MEDUSA - Air-Sea OUTPUT: ', &284 'unrealistic surface Carb. Chemistry -- OUTPUTS' )285 283 ENDIF 284 CALL ctl_stop( 'MEDUSA - Air-Sea OUTPUT: ', & 285 'unrealistic surface Carb. Chemistry -- OUTPUTS' ) 286 286 ENDIF 287 287 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.