Ignore:
Timestamp:
2019-10-18T16:49:33+02:00 (13 months ago)
Author:
jpalmier
Message:

add Changes to enable Alkalinity geo-ingeneering tests within MEDUSA

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_GO6_GeoAlk_MEDUSA/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_update.F90

    r10020 r11726  
    7171# endif 
    7272                                   zpds, zphd, zphn 
    73       USE dom_oce,           ONLY: e3t_0, gphit, mbathy, tmask 
     73      USE dom_oce,           ONLY: e3t_0, gphit, mbathy, tmask, nyear          ! AXY (16/10/19) 
     74 
    7475# if defined key_vvl 
    7576      USE dom_oce,           ONLY: e3t_n 
     
    232233               !!  
    233234               !! riverine flux 
     235               !! (NOTE: this is "protected" for near-surface addition of riverine 
     236               !! by array friver_dep which is zero below a defined model level) 
    234237               if ( jriver_n .gt. 0 ) then 
    235238                  f_riv_loc_n(ji,jj) = f_riv_n(ji,jj) *                      & 
     
    276279               !!  
    277280               !! riverine flux 
     281               !! (NOTE: this is "protected" for near-surface addition of riverine 
     282               !! by array friver_dep which is zero below a defined model level) 
    278283               if ( jriver_si .gt. 0 ) then 
    279284                  f_riv_loc_si(ji,jj) = f_riv_si(ji,jj) *                    & 
     
    362367               !!  
    363368               !! riverine flux 
     369               !! (NOTE: this is "protected" for near-surface addition of riverine 
     370               !! by array friver_dep which is zero below a defined model level) 
    364371               if ( jriver_c .gt. 0 ) then 
    365372                  f_riv_loc_c(ji,jj) = f_riv_c(ji,jj) *                      & 
     
    387394               !!  
    388395               !! riverine flux 
    389                if ( jk .eq. 1 .and. jriver_c .gt. 0 ) then 
    390                   fc_prod = fc_prod + f_riv_c(ji,jj) 
    391                endif 
     396               !! (NOTE: the lines immediately below are incorrect and riverine 
     397               !! flux is already handled above) 
     398               !! if ( jk .eq. 1 .and. jriver_c .gt. 0 ) then 
     399               !!    fc_prod = fc_prod + f_riv_c(ji,jj) 
     400               !! endif 
    392401               !! 
    393402               btra(ji,jj,jpdic_lc) = b0 * ( fc_prod + fc_cons ) 
     
    414423               !!  
    415424               !! riverine flux 
     425               !! (NOTE: this is "protected" for near-surface addition of riverine 
     426               !! by array friver_dep which is zero below a defined model level) 
    416427               if ( jriver_alk .gt. 0 ) then 
    417428                  f_riv_loc_alk(ji,jj) = f_riv_alk(ji,jj) *                  & 
     
    419430                                         fse3t(ji,jj,jk) 
    420431                  fa_prod = fa_prod + f_riv_loc_alk(ji,jj) 
     432               endif 
     433          !! 
     434               !! AXY (16/10/19) 
     435               !! geoengineering flux 
     436               !! (NOTE: this is only added to the surface ocean grid cell, only 
     437               !! within a particular time window, and only within a particular 
     438               !! sub-polar geographical range; these parameters could be passed up 
     439               !! to ROSE if they need different values for different experiments) 
     440          if ( jk .eq. 1 .and. nyear .gt. 2014 .and. nyear .lt. 2050 .and. & 
     441                   gphit(ji,jj) .ge .-60.0 .and. gphit(ji,jj) .le. 70.0 ) then 
     442                  f_geoeng_alk(ji,jj) = 1.17 / fse3t(ji,jj,jk) 
     443                  fa_prod = fa_prod + f_geoeng_alk(ji,jj) 
    421444               endif 
    422445               !!   
Note: See TracChangeset for help on using the changeset viewer.