Changeset 7880


Ignore:
Timestamp:
2017-04-06T12:59:22+02:00 (3 years ago)
Author:
marc
Message:

Marc 6/4/17. Fix bug where flatx is not set in call to mocsy_interface

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

    r7766 r7880  
    354354      REAL(wp) ::    f_riv_loc_n, f_riv_loc_si, f_riv_loc_c, f_riv_loc_alk 
    355355      !! 
    356       !! horizontal grid location 
    357       REAL(wp) ::    flatx, flonx 
    358       !! 
    359356      !! Jpalm -- 11-10-2015 -- adapt diag to iom_use 
    360357      !! 2D var for diagnostics. 
     
    14011398                     !!                 chemistry package 
    14021399                     CALL mocsy_interface( ztmp, zsal, zalk, zdic, zsil, zpho,         &    ! inputs 
    1403                      f_pp0, fdep2, flatx, f_kw660, f_xco2a, 1,                         &    ! inputs 
     1400                     f_pp0, fdep2, gphit(ji,jj), f_kw660, f_xco2a, 1,                  &    ! inputs 
    14041401                     f_ph, f_pco2w, f_fco2w, f_h2co3, f_hco3, f_co3, f_omarg(ji,jj),   &    ! outputs 
    14051402                     f_omcal(ji,jj), f_BetaD, f_rhosw, f_opres, f_insitut,             &    ! outputs 
     
    15231520               !! AXY (01/03/10): set up level depth (bottom of level) 
    15241521               fdep1 = fdep + fthk 
    1525                !! AXY (17/05/13): where am I? 
    1526                flatx = gphit(ji,jj) 
    1527                flonx = glamt(ji,jj) 
    15281522               !! AXY (28/11/16): local seafloor depth 
    15291523               !!                 previously mbathy(ji,jj) - 1, now mbathy(ji,jj) 
     
    17621756                  !!                 zero in this call 
    17631757                  CALL mocsy_interface( ztmp, zsal, zalk, zdic, zsil, zpho,        &  ! inputs 
    1764                   f_pp0, 0.0, flatx, f_kw660, f_xco2a, 1,                          &  ! inputs 
     1758                  f_pp0, 0.0, gphit(ji,jj), f_kw660, f_xco2a, 1,                   &  ! inputs 
    17651759                  f_ph, f_pco2w, f_fco2w, f_h2co3, f_hco3, f_co3, f_omarg(ji,jj),  &  ! outputs 
    17661760                  f_omcal(ji,jj), f_BetaD, f_rhosw, f_opres, f_insitut,            &  ! outputs 
     
    38953889               !!                 away with if I switch off BGC fluxes in the Gulf of  
    38963890               !!                 Finland, which currently appears the source of trouble 
    3897                if ( flonx.gt.24.7 .and. flonx.lt.27.8 .and. & 
    3898                   &   flatx.gt.59.2 .and. flatx.lt.60.2 ) then 
     3891               if ( glamt(ji,jj).gt.24.7 .and. glamt(ji,jj).lt.27.8 .and. & 
     3892                  &   gphit(ji,jj).gt.59.2 .and. gphit(ji,jj).lt.60.2 ) then 
    38993893                  ibio_switch = 0 
    39003894               endif 
Note: See TracChangeset for help on using the changeset viewer.