Changeset 9332


Ignore:
Timestamp:
2018-02-19T14:30:30+01:00 (2 years ago)
Author:
jpalmier
Message:

JPALM — MEDUSA Mocsy output check

File:
1 edited

Legend:

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

    r9329 r9332  
    258258               !!          must be associated to air-sea flux or air xCO2... 
    259259               !!          Check MOCSY outputs 
    260               IF ( (f_pco2w(ji,jj) > 1.E4 ).OR.(f_pco2w(ji,jj) < 0.0 ) .OR.     & 
    261                    (f_fco2w(ji,jj) > 1.E4 ).OR.(f_fco2w(ji,jj) < 0.0 ) .OR.     &    
    262                    (f_fco2atm(ji,jj) > 1.E4 ).OR.(f_fco2atm(ji,jj) < 0.0 ) .OR.     & 
    263                    (f_co2flux(ji,jj) > 1.E-1 ).OR.(f_co2flux(ji,jj) < -1.E-1 ) .OR.     & 
    264                    (f_dpco2(ji,jj) > 1.E4 ).OR.(f_dpco2(ji,jj) < -1.E4 ) ) THEN 
    265                  IF(lwp) THEN  
    266                      WRITE(numout,*) ' surface T = ',ztmp(ji,jj) 
    267                      WRITE(numout,*) ' surface S = ',zsal(ji,jj) 
    268                      WRITE(numout,*) ' surface ALK = ',zalk(ji,jj) 
    269                      WRITE(numout,*) ' surface DIC = ',zdic(ji,jj) 
    270                      WRITE(numout,*) ' KW660 = ',f_kw660(ji,jj) 
    271                      WRITE(numout,*) ' atm xCO2 = ',f_xco2a(ji,jj)    
    272                      WRITE(numout,*) ' surface pco2w  = ',f_pco2w(ji,jj) 
    273                      WRITE(numout,*) ' surface fco2w  = ',f_fco2w(ji,jj) 
    274                      WRITE(numout,*) ' surface fco2a  = ',f_fco2atm(ji,jj) 
    275                      WRITE(numout,*) ' surface co2flx = ',f_co2flux(ji,jj) 
    276                      WRITE(numout,*) ' surface dpco2  = ',f_dpco2(ji,jj) 
    277                      WRITE(numout,*) ' MOCSY output: ji =', mig(ji),' jj = ', mjg(jj),  & 
    278                                        ' kt = ', kt      
    279                      WRITE(numout,*) 'MEDUSA - Air-Sea OUTPUT: unrealistic surface Carb. Chemistry' 
    280                  ENDIF      
    281                  CALL ctl_stop( 'MEDUSA - Air-Sea OUTPUT: ',            & 
    282                                 'unrealistic surface Carb. Chemistry -- OUTPUTS' ) 
    283               ENDIF      
     260               !!=================== 
     261               !! Jpalm -- 19-02-2018 -- remove the cap - only check MOCSY inputs 
     262               !!       because of specific area in arabic sea where strangely 
     263               !!       with core 2 forcing, ALK is lower than DIC and result in  
     264               !!       Enormous dpco2 - even if all carb chem caract are OK. 
     265               !!       and this check stops the model. 
     266               !!       ==> let's run the model without it and see how it goes. 
     267               !!=================== 
     268               !!IF ( (f_pco2w(ji,jj) > 1.E4 ).OR.(f_pco2w(ji,jj) < 0.0 ) .OR.     & 
     269               !!    (f_fco2w(ji,jj) > 1.E4 ).OR.(f_fco2w(ji,jj) < 0.0 ) .OR.     &    
     270               !!    (f_fco2atm(ji,jj) > 1.E4 ).OR.(f_fco2atm(ji,jj) < 0.0 ) .OR.     & 
     271               !!    (f_co2flux(ji,jj) > 1.E-1 ).OR.(f_co2flux(ji,jj) < -1.E-1 ) .OR.     & 
     272               !!    (f_dpco2(ji,jj) > 1.E4 ).OR.(f_dpco2(ji,jj) < -1.E4 ) ) THEN 
     273               !!  IF(lwp) THEN  
     274               !!      WRITE(numout,*) ' surface T = ',ztmp(ji,jj) 
     275               !!      WRITE(numout,*) ' surface S = ',zsal(ji,jj) 
     276               !!      WRITE(numout,*) ' surface ALK = ',zalk(ji,jj) 
     277               !!      WRITE(numout,*) ' surface DIC = ',zdic(ji,jj) 
     278               !!      WRITE(numout,*) ' KW660 = ',f_kw660(ji,jj) 
     279               !!      WRITE(numout,*) ' atm xCO2 = ',f_xco2a(ji,jj)    
     280               !!      WRITE(numout,*) ' surface pco2w  = ',f_pco2w(ji,jj) 
     281               !!      WRITE(numout,*) ' surface fco2w  = ',f_fco2w(ji,jj) 
     282               !!      WRITE(numout,*) ' surface fco2a  = ',f_fco2atm(ji,jj) 
     283               !!      WRITE(numout,*) ' surface co2flx = ',f_co2flux(ji,jj) 
     284               !!      WRITE(numout,*) ' surface dpco2  = ',f_dpco2(ji,jj) 
     285               !!      WRITE(numout,*) ' MOCSY output: ji =', mig(ji),' jj = ', mjg(jj),  & 
     286               !!                        ' kt = ', kt      
     287               !!      WRITE(numout,*) 'MEDUSA - Air-Sea OUTPUT: unrealistic surface Carb. Chemistry' 
     288               !!  ENDIF      
     289               !!  CALL ctl_stop( 'MEDUSA - Air-Sea OUTPUT: ',            & 
     290               !!                 'unrealistic surface Carb. Chemistry -- OUTPUTS' ) 
     291               !!ENDIF      
    284292            ENDIF 
    285293         ENDDO 
Note: See TracChangeset for help on using the changeset viewer.