Changeset 9309


Ignore:
Timestamp:
2018-02-06T12:22:05+01:00 (3 years ago)
Author:
jpalmier
Message:

JPALM — make clearer in the code and output which atm xCO2 source is used

Location:
branches/NERC/dev_r5518_GO6_CO2_cmip/NEMOGCM/NEMO/TOP_SRC/MEDUSA
Files:
5 edited

Legend:

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

    r9257 r9309  
    153153            IF (tmask(ji,jj,1) == 1) then 
    154154               IF (lk_oasis) THEN 
    155                   !! use 2D atm xCO2 from atm coupling 
    156                   f_xco2a(ji,jj) = PCO2a_in_cpl(ji,jj) 
    157                   !!!  
    158155                  !!! Jpalm test on atm xCO2 
    159156                  IF ( (f_xco2a(ji,jj) > 1500.0 ).OR.(f_xco2a(ji,jj) < 100.0 ) ) THEN 
     
    162159                                       ' -- ji =', mig(ji),' jj = ', mjg(jj) 
    163160                    ENDIF 
    164                     CALL ctl_stop( 'MEDUSA - Air-Sea :', 'unrealistic atm xCO2 ' ) 
     161                    CALL ctl_stop( 'MEDUSA - Air-Sea :', 'unrealistic coupled atm xCO2 ' ) 
    165162                 ENDIF  
    166163               ENDIF 
  • branches/NERC/dev_r5518_GO6_CO2_cmip/NEMOGCM/NEMO/TOP_SRC/MEDUSA/carb_chem.F90

    r8441 r9309  
    103103               !! OPEN wet point IF..THEN loop 
    104104               IF (tmask(ji,jj,jk).eq.1) THEN 
    105                   IF (lk_oasis) THEN 
    106                      !! use 2D atm xCO2 from atm coupling 
    107                      f_xco2a(ji,jj) = PCO2a_in_cpl(ji,jj) 
    108                   ENDIF 
    109105                  !! Do carbonate chemistry 
    110106                  !! 
  • branches/NERC/dev_r5518_GO6_CO2_cmip/NEMOGCM/NEMO/TOP_SRC/MEDUSA/sms_medusa.F90

    r9298 r9309  
    289289   INTEGER  :: co2_rec 
    290290   REAL(wp) :: co2_yinit, co2_yend 
     291#endif 
     292 
     293!!---------------------------------------------------------------------- 
     294!! JPALM -- PI CO2 key 
     295!!---------------------------------------------------------------------- 
     296!! 
     297#if defined key_axy_pi_co2 
     298   LOGICAL , PUBLIC ::   lk_pi_co2 = .TRUE.  !: OASIS used 
     299#else 
     300   LOGICAL , PUBLIC ::   lk_pi_co2 = .FALSE. !: OASIS unused 
    291301#endif 
    292302 
  • branches/NERC/dev_r5518_GO6_CO2_cmip/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

    r9304 r9309  
    9393                                            mpp_min, mpp_minloc,            & 
    9494                                            ctl_stop, ctl_warn, lk_mpp   
    95       USE oce,                        ONLY: tsb, tsn 
     95      USE oce,                        ONLY: tsb, tsn, PCO2a_in_cpl 
    9696      USE par_kind,                   ONLY: wp 
    9797      USE par_medusa,                 ONLY: jpalk, jpchd, jpchn, jpdet,     & 
     
    103103      !! JPALM (27-06-2016): add lk_oasis for CO2 and DMS coupling with atm 
    104104      USE sbc_oce,                    ONLY: lk_oasis 
    105       USE sms_medusa,                 ONLY: hist_pco2, co2_yinit, co2_yend 
     105      USE sms_medusa,                 ONLY: hist_pco2, co2_yinit, co2_yend, & 
     106                                            lk_pi_co2 
    106107      USE trc,                        ONLY: ln_rsttr, nittrc000, trn 
    107108      USE bio_medusa_init_mod,        ONLY: bio_medusa_init 
     
    295296      !!------------------------------------------------------------------ 
    296297      !! 
    297       !! what's atmospheric pCO2 doing? (data start in 1859) 
    298       !! AXY - JPALM new interpolation scheme usinf nyear_len 
    299       this_y = real(nyear) 
    300       this_d = real(nday_year) 
    301       this_s = real(nsec_day) 
    302       !! 
    303       fyear = this_y + ((this_d - 1) + (this_s / (60. * 60. * 24.))) / real(nyear_len(1)) 
    304       !! 
    305       IF ( ( kt == nittrc000 ) .AND. lwp ) THEN  
    306          WRITE(numout,*) ' MEDUSA year      =', this_y 
    307          WRITE(numout,*) ' Year length      =', real(nyear_len(1)) 
    308          WRITE(numout,*) ' MEDUSA nday_year =', this_d 
    309          WRITE(numout,*) ' MEDUSA nsec_day  =', this_s 
    310       ENDIF 
    311       !!  
    312       !! different case test  
    313       IF (fyear .LE. co2_yinit) THEN 
    314          !! before first record -- pre-industrial value 
    315          f_xco2a(:,:) = hist_pco2(1)   
    316       ELSEIF (fyear .GE. co2_yend) THEN 
    317          !! after last record - continue to use the last value 
    318          f_xco2a(:,:) = hist_pco2(int(co2_yend - co2_yinit + 1.) ) 
    319       ELSE 
    320          !! just right 
    321          iyr1 = int(fyear - co2_yinit) + 1 
    322          iyr2 = iyr1 + 1  
    323          fq3 = fyear - real(iyr1) - co2_yinit + 1. 
    324          fq4 = ((1 - fq3) * hist_pco2(iyr1)) + (fq3 * hist_pco2(iyr2)) 
    325          f_xco2a(:,:) = fq4 
     298      IF (lk_oasis) THEN  
     299         !! xCO2 from coupled 
     300         IF ( ( kt == nittrc000 ) .AND. lwp )     & 
     301             WRITE(numout,*) '** MEDUSA Atm xCO2 given by the UM **' 
     302         f_xco2a(:,:) = PCO2a_in_cpl(:,:) 
     303      ELSEIF (lk_pi_co2) THEN 
     304         !! OCMIP pre-industrial xCO2 
     305         IF ( ( kt == nittrc000 ) .AND. lwp )     & 
     306             WRITE(numout,*) '** MEDUSA Atm xCO2 fixed to pre-industrial value **' 
     307         !! f_xco2a(:,:) = 284.725  !! CMIP5 pre-industrial pCO2 
     308         f_xco2a(:,:) = 284.317          !! CMIP6 pre-industrial pCO2  
     309      ELSE       
     310         !! xCO2 from file 
     311         !! AXY - JPALM new interpolation scheme usinf nyear_len 
     312         this_y = real(nyear) 
     313         this_d = real(nday_year) 
     314         this_s = real(nsec_day) 
     315         !! 
     316         fyear = this_y + ((this_d - 1) + (this_s / (60. * 60. * 24.))) / real(nyear_len(1)) 
    326317         !! 
    327318         IF ( ( kt == nittrc000 ) .AND. lwp ) THEN  
    328             WRITE(numout,*) ' MEDUSA year1     =', iyr1 
    329             WRITE(numout,*) ' MEDUSA year2     =', iyr2 
    330             WRITE(numout,*) ' xCO2 year1       =', hist_pco2(iyr1) 
    331             WRITE(numout,*) ' xCO2 year2       =', hist_pco2(iyr2) 
    332             WRITE(numout,*) ' Year2 weight     =', fq3 
     319            WRITE(numout,*) '** MEDUSA Atm xCO2 from file **' 
     320            WRITE(numout,*) ' MEDUSA year      =', this_y 
     321            WRITE(numout,*) ' Year length      =', real(nyear_len(1)) 
     322            WRITE(numout,*) ' MEDUSA nday_year =', this_d 
     323            WRITE(numout,*) ' MEDUSA nsec_day  =', this_s 
     324         ENDIF 
     325         !!  
     326         !! different case test  
     327         IF (fyear .LE. co2_yinit) THEN 
     328            !! before first record -- pre-industrial value 
     329            f_xco2a(:,:) = hist_pco2(1)   
     330         ELSEIF (fyear .GE. co2_yend) THEN 
     331            !! after last record - continue to use the last value 
     332            f_xco2a(:,:) = hist_pco2(int(co2_yend - co2_yinit + 1.) ) 
     333         ELSE 
     334            !! just right 
     335            iyr1 = int(fyear - co2_yinit) + 1 
     336            iyr2 = iyr1 + 1  
     337            fq3 = fyear - real(iyr1) - co2_yinit + 1. 
     338            fq4 = ((1 - fq3) * hist_pco2(iyr1)) + (fq3 * hist_pco2(iyr2)) 
     339            f_xco2a(:,:) = fq4 
     340            !! 
     341            IF ( ( kt == nittrc000 ) .AND. lwp ) THEN  
     342               WRITE(numout,*) ' MEDUSA year1     =', iyr1 
     343               WRITE(numout,*) ' MEDUSA year2     =', iyr2 
     344               WRITE(numout,*) ' xCO2 year1       =', hist_pco2(iyr1) 
     345               WRITE(numout,*) ' xCO2 year2       =', hist_pco2(iyr2) 
     346               WRITE(numout,*) ' Year2 weight     =', fq3 
     347            ENDIF 
    333348         ENDIF 
    334349      ENDIF 
    335 #  if defined key_axy_pi_co2 
    336       !! OCMIP pre-industrial pCO2 
    337       !! f_xco2a(:,:) = 284.725  !! CMIP5 pre-industrial pCO2 
    338       f_xco2a(:,:) = 284.317          !! CMIP6 pre-industrial pCO2  
    339 #  endif 
    340       !! AXY (29/01/14): remove surplus diagnostics 
    341        IF ( ( kt == nittrc000 ) .AND. lwp )   &  
     350      IF ( ( kt == nittrc000 ) .AND. lwp )   &  
    342351            WRITE(numout,*) ' final atm xCO2   =', f_xco2a(1,1) 
    343352# endif 
  • branches/NERC/dev_r5518_GO6_CO2_cmip/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcini_medusa.F90

    r9298 r9309  
    329329      !!                     initially done in trcsms_medusa 
    330330      !!---------------------------------------------------------------------- 
    331       !!      
    332       IF(lwp) WRITE(numout,*) ' trc_ini_medusa: initialisating atm CO2 record' 
    333       CALL trc_ini_medusa_co2atm 
    334  
     331      !!  
     332      IF( ( .NOT.lk_oasis ) .AND. ( .NOT.lk_pi_co2 ) ) THEN 
     333         IF(lwp) WRITE(numout,*) ' trc_ini_medusa: initialisating atm CO2 record' 
     334         CALL trc_ini_medusa_co2atm 
     335      ENDIF 
    335336 
    336337   END SUBROUTINE trc_ini_medusa 
Note: See TracChangeset for help on using the changeset viewer.