Ignore:
Timestamp:
2018-01-18T11:14:03+01:00 (3 years ago)
Author:
frrh
Message:

Commit JP's Met Office GMED ticket 371 for trapping or notifying of
peculiar values of MEDUSA fields arising from transient temperature
spikes in the ocean.

Committed using:
svn merge: 9177:9249 svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_9163

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/MEDUSA/air_sea.F90

    r8521 r9257  
    6262# endif 
    6363                                   zchd, zchn, zdin, zsil 
    64       USE dom_oce,           ONLY: e3t_0, e3t_n, gphit, tmask 
     64      USE dom_oce,           ONLY: e3t_0, e3t_n, gphit, tmask, mig, mjg 
    6565# if defined key_iomput 
    6666      USE iom,               ONLY: lk_iomput 
     
    9191      USE trcoxy_medusa,     ONLY: trc_oxy_medusa 
    9292# endif 
     93      USE lib_mpp,           ONLY: ctl_stop 
     94      USE trcstat,           ONLY: trc_rst_dia_stat  
    9395 
    9496   !!* Substitution 
     
    122124 
    123125# if defined key_roam 
     126      !! init 
     127      f_fco2w(:,:)       = 0.0 
     128      f_fco2atm(:,:)     = 0.0 
     129      f_schmidtco2(:,:)  = 0.0 
     130      f_kwco2(:,:)       = 0.0 
     131      f_co2starair(:,:)  = 0.0 
     132      f_dpco2(:,:)       = 0.0 
     133      f_rhosw(:,:)       = 0.0 
     134      f_K0(:,:)          = 0.0 
     135      !! air pressure (atm); ultimately this will use air  
     136      !! pressure at the base of the UKESM1 atmosphere  
     137      !!                                      
     138      f_pp0(:,:)   = 1.0 
     139 
     140 
    124141      !!----------------------------------------------------------- 
    125142      !! Air-sea gas exchange 
     
    134151         DO ji = 2,jpim1 
    135152            !! OPEN wet point IF..THEN loop 
    136             if (tmask(ji,jj,1) == 1) then 
     153            IF (tmask(ji,jj,1) == 1) then 
    137154               IF (lk_oasis) THEN 
    138155                  !! use 2D atm xCO2 from atm coupling 
    139156                  f_xco2a(ji,jj) = PCO2a_in_cpl(ji,jj) 
     157                  !!!  
     158                  !!! Jpalm test on atm xCO2 
     159                  IF ( (f_xco2a(ji,jj) > 1500.0 ).OR.(f_xco2a(ji,jj) < 100.0 ) ) THEN 
     160                    IF(lwp) THEN  
     161                       WRITE(numout,*) ' atm xCO2 = ',f_xco2a(ji,jj),          & 
     162                                       ' -- ji =', mig(ji),' jj = ', mjg(jj) 
     163                    ENDIF 
     164                    CALL ctl_stop( 'MEDUSA - Air-Sea :', 'unrealistic atm xCO2 ' ) 
     165                 ENDIF  
    140166               ENDIF 
    141167               !! 
     
    162188               'air-sea: carb-chem kt = ', kt 
    163189               CALL flush(numout) 
     190               !! JPALM add carb print: 
     191               call trc_rst_dia_stat(f_xco2a(:,:), 'f_xco2a') 
     192               call trc_rst_dia_stat(wndm(:,:), 'wndm') 
     193               call trc_rst_dia_stat(f_kw660(:,:), 'f_kw660') 
     194               call trc_rst_dia_stat(ztmp(:,:), 'ztmp') 
     195               call trc_rst_dia_stat(zsal(:,:), 'zsal') 
     196               call trc_rst_dia_stat(zalk(:,:), 'zalk') 
     197               call trc_rst_dia_stat(zdic(:,:), 'zdic') 
     198               call trc_rst_dia_stat(zsil(:,:), 'zsil') 
     199               call trc_rst_dia_stat(zpho(:,:), 'zpho') 
    164200#   endif 
    165201      DO jj = 2,jpjm1 
    166202         DO ji = 2,jpim1 
    167203            if (tmask(ji,jj,1) == 1) then 
    168                !! air pressure (atm); ultimately this will use air  
    169                !! pressure at the base of the UKESM1 atmosphere  
    170                !!                                      
    171                f_pp0(ji,jj)   = 1.0 
    172                !! 
    173                !! IF(lwp) WRITE(numout,*) ' MEDUSA ztmp    =', ztmp(ji,jj) 
    174                !! IF(lwp) WRITE(numout,*) ' MEDUSA wndm    =', wndm(ji,jj) 
    175                !! IF(lwp) WRITE(numout,*) ' MEDUSA fr_i    =', fr_i(ji,jj) 
    176204               !! 
    177205#  if defined key_axy_carbchem 
    178206#   if defined key_mocsy 
     207               !! Jpalm -- 12-09-2017 -- add extra check after reccurent 
     208               !!          carbonate failure in the coupled run. 
     209               !!          must be associated to air-sea flux or air xCO2... 
     210               !!          Check MOCSY inputs 
     211               IF ( (zsal(ji,jj) > 75.0 ).OR.(zsal(ji,jj) < 0.0 ) .OR.        & 
     212                    (ztmp(ji,jj) > 50.0 ).OR.(ztmp(ji,jj) < -20.0 ) .OR.      & 
     213                    (zalk(ji,jj) > 35.0E2 ).OR.(zalk(ji,jj) <= 0.0 ) .OR.     & 
     214                    (zdic(ji,jj) > 35.0E2 ).OR.(zdic(ji,jj) <= 0.0 ) .OR.     & 
     215                    (f_kw660(ji,jj) > 1.0E-2 ).OR.(f_kw660(ji,jj) < 0.0 ) ) THEN 
     216                  IF(lwp) THEN  
     217                      WRITE(numout,*) ' surface T = ',ztmp(ji,jj) 
     218                      WRITE(numout,*) ' surface S = ',zsal(ji,jj) 
     219                      WRITE(numout,*) ' surface ALK = ',zalk(ji,jj) 
     220                      WRITE(numout,*) ' surface DIC = ',zdic(ji,jj) 
     221                      WRITE(numout,*) ' KW660 = ',f_kw660(ji,jj) 
     222                      WRITE(numout,*) ' atm xCO2 = ',f_xco2a(ji,jj)    
     223                      WRITE(numout,*) ' surface pco2w  = ',f_pco2w(ji,jj) 
     224                      WRITE(numout,*) ' surface fco2w  = ',f_fco2w(ji,jj) 
     225                      WRITE(numout,*) ' surface fco2a  = ',f_fco2atm(ji,jj) 
     226                      WRITE(numout,*) ' surface co2flx = ',f_co2flux(ji,jj) 
     227                      WRITE(numout,*) ' surface dpco2  = ',f_dpco2(ji,jj) 
     228                      WRITE(numout,*) ' MOCSY input: ji =', mig(ji),' jj = ', mjg(jj),  & 
     229                                       ' kt = ', kt  
     230                      WRITE(numout,*) 'MEDUSA - Air-Sea INPUT: unrealistic surface Carb. Chemistry' 
     231                  ENDIF      
     232                  CALL ctl_stop( 'MEDUSA - Air-Sea INPUT: ',             & 
     233                                 'unrealistic surface Carb. Chemistry -- INPUTS' ) 
     234               ENDIF      
    179235               !! 
    180236               !! AXY (22/06/15): use Orr & Epitalon (2015) MOCSY-2 carbonate 
     
    201257               f_TALK(ji,jj) = (zalk(ji,jj) / f_rhosw(ji,jj)) * 1000. 
    202258               f_dcf(ji,jj)  = f_rhosw(ji,jj) 
     259               !! Jpalm -- 12-09-2017 -- add extra check after reccurent 
     260               !!          carbonate failure in the coupled run. 
     261               !!          must be associated to air-sea flux or air xCO2... 
     262               !!          Check MOCSY outputs 
     263              IF ( (f_pco2w(ji,jj) > 1.E4 ).OR.(f_pco2w(ji,jj) < 0.0 ) .OR.     & 
     264                   (f_fco2w(ji,jj) > 1.E4 ).OR.(f_fco2w(ji,jj) < 0.0 ) .OR.     &    
     265                   (f_fco2atm(ji,jj) > 1.E4 ).OR.(f_fco2atm(ji,jj) < 0.0 ) .OR.     & 
     266                   (f_co2flux(ji,jj) > 1.E-2 ).OR.(f_co2flux(ji,jj) < -1.E-2 ) .OR.     & 
     267                   (f_dpco2(ji,jj) > 1.E4 ).OR.(f_dpco2(ji,jj) < -1.E4 ) ) THEN 
     268                 IF(lwp) THEN  
     269                     WRITE(numout,*) ' surface T = ',ztmp(ji,jj) 
     270                     WRITE(numout,*) ' surface S = ',zsal(ji,jj) 
     271                     WRITE(numout,*) ' surface ALK = ',zalk(ji,jj) 
     272                     WRITE(numout,*) ' surface DIC = ',zdic(ji,jj) 
     273                     WRITE(numout,*) ' KW660 = ',f_kw660(ji,jj) 
     274                     WRITE(numout,*) ' atm xCO2 = ',f_xco2a(ji,jj)    
     275                     WRITE(numout,*) ' surface pco2w  = ',f_pco2w(ji,jj) 
     276                     WRITE(numout,*) ' surface fco2w  = ',f_fco2w(ji,jj) 
     277                     WRITE(numout,*) ' surface fco2a  = ',f_fco2atm(ji,jj) 
     278                     WRITE(numout,*) ' surface co2flx = ',f_co2flux(ji,jj) 
     279                     WRITE(numout,*) ' surface dpco2  = ',f_dpco2(ji,jj) 
     280                     WRITE(numout,*) ' MOCSY output: ji =', mig(ji),' jj = ', mjg(jj),  & 
     281                                       ' kt = ', kt      
     282                     WRITE(numout,*) 'MEDUSA - Air-Sea OUTPUT: unrealistic surface Carb. Chemistry' 
     283                 ENDIF      
     284                 CALL ctl_stop( 'MEDUSA - Air-Sea OUTPUT: ',            & 
     285                                'unrealistic surface Carb. Chemistry -- OUTPUTS' ) 
     286              ENDIF      
    203287            ENDIF 
    204288         ENDDO 
    205289      ENDDO 
    206290 
     291#   if defined key_debug_medusa 
     292               !! JPALM add carb print: 
     293               call trc_rst_dia_stat(f_pco2w(:,:), 'f_pco2w') 
     294               call trc_rst_dia_stat(f_fco2w(:,:), 'f_fco2w') 
     295               call trc_rst_dia_stat(f_fco2atm(:,:), 'f_fco2atm') 
     296               call trc_rst_dia_stat(f_schmidtco2(:,:), 'f_schmidtco2') 
     297               call trc_rst_dia_stat(f_kwco2(:,:), 'f_kwco2') 
     298               call trc_rst_dia_stat(f_co2starair(:,:), 'f_co2starair') 
     299               call trc_rst_dia_stat(f_co2flux(:,:), 'f_co2flux') 
     300               call trc_rst_dia_stat(f_dpco2(:,:), 'f_dpco2') 
     301#   endif 
    207302#   else    
    208303 
Note: See TracChangeset for help on using the changeset viewer.