Changeset 13524


Ignore:
Timestamp:
2020-09-25T14:30:18+02:00 (5 months ago)
Author:
rlod
Message:

First step of nutrient sources from Antarctic Ice Sheet

Location:
NEMO/branches/2020/dev_r12740_TOP-01_rlod_Antarctic_ice_Sheet_Fe_Source/src
Files:
1 added
6 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12740_TOP-01_rlod_Antarctic_ice_Sheet_Fe_Source/src/OCE/SBC/sbcrnf.F90

    r12489 r13524  
    4242   REAL(wp)                   ::      rn_dep_max        !: depth over which runoffs is spread       (ln_rnf_depth_ini =T) 
    4343   INTEGER                    ::      nn_rnf_depth_file !: create (=1) a runoff depth file or not (=0) 
    44    LOGICAL                   ::   ln_rnf_icb        !: iceberg flux is specified in a file 
     44   LOGICAL           , PUBLIC ::   ln_rnf_icb        !: iceberg flux is specified in a file 
    4545   LOGICAL                    ::   ln_rnf_tem        !: temperature river runoffs attribute specified in a file 
    4646   LOGICAL           , PUBLIC ::   ln_rnf_sal        !: salinity    river runoffs attribute specified in a file 
  • NEMO/branches/2020/dev_r12740_TOP-01_rlod_Antarctic_ice_Sheet_Fe_Source/src/OCE/nemogcm.F90

    r12641 r13524  
    481481                           CALL dyn_spg_init         ! surface pressure gradient 
    482482 
     483      !                                      ! Icebergs 
     484                           CALL icb_init( rn_Dt, nit000)   ! initialise icebergs instance 
     485 
     486                                                ! ice shelf 
     487                           CALL isf_init( Nbb, Nnn, Naa ) 
     488 
    483489#if defined key_top 
    484490      !                                      ! Passive tracers 
    485491                           CALL     trc_init( Nbb, Nnn, Naa ) 
    486492#endif 
     493 
    487494      IF( l_ldfslp     )   CALL ldf_slp_init    ! slope of lateral mixing 
    488  
    489       !                                      ! Icebergs 
    490                            CALL icb_init( rn_Dt, nit000)   ! initialise icebergs instance 
    491  
    492                                                 ! ice shelf 
    493                            CALL isf_init( Nbb, Nnn, Naa ) 
    494495 
    495496      !                                      ! Misc. options 
  • NEMO/branches/2020/dev_r12740_TOP-01_rlod_Antarctic_ice_Sheet_Fe_Source/src/TOP/TRP/trctrp.F90

    r12377 r13524  
    2424   USE trcsbc          ! surface boundary condition          (trc_sbc routine) 
    2525   USE trcbc           ! Tracers boundary condtions          ( trc_bc routine) 
     26   USE trcais          ! Antarctic Ice Sheet tracers         (trc_ais routine) 
    2627   USE zpshde          ! partial step: hor. derivative       (zps_hde routine) 
    2728   USE bdy_oce   , ONLY: ln_bdy 
     
    6566         IF( ln_trcbc .AND. lltrcbc .AND. kt /= nit000 )  & 
    6667                                CALL trc_bc     ( kt,      Kmm, tr, Krhs )      ! tracers: surface and lateral Boundary Conditions  
     68         IF( ln_trcais )        CALL trc_ais    ( kt,      Kmm, tr, Krhs )      ! tracers from Antarctic Ice Sheet (icb, isf) 
    6769         IF( ln_trabbl )        CALL trc_bbl    ( kt, Kbb, Kmm, tr, Krhs )      ! advective (and/or diffusive) bottom boundary layer scheme 
    6870         IF( ln_trcdmp )        CALL trc_dmp    ( kt, Kbb, Kmm, tr, Krhs )      ! internal damping trends 
  • NEMO/branches/2020/dev_r12740_TOP-01_rlod_Antarctic_ice_Sheet_Fe_Source/src/TOP/trc.F90

    r12489 r13524  
    4040   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:  ) ::  trc_o          !: prescribed tracer concentration in ocean for SBC 
    4141   INTEGER             , PUBLIC                            ::  nn_ice_tr      !: handling of sea ice tracers 
     42   INTEGER             , PUBLIC                            ::  nn_ais_tr      !: handling of Antarctic Ice Sheet tracers 
    4243 
    4344   !! interpolated gradient 
     
    6566   LOGICAL             , PUBLIC ::   ln_trcdta          !: Read inputs data from files 
    6667   LOGICAL             , PUBLIC ::   ln_trcbc           !: Enable surface, lateral or open boundaries conditions 
     68   LOGICAL             , PUBLIC ::   ln_trcais          !: Enable Antarctic Ice Sheet nutrient supply 
    6769   LOGICAL             , PUBLIC ::   ln_trcdmp          !: internal damping flag 
    6870   LOGICAL             , PUBLIC ::   ln_trcdmp_clo      !: internal damping flag on closed seas 
     
    9395      LOGICAL           ::   llcbc     ! read in a file or not 
    9496      LOGICAL           ::   llobc     ! read in a file or not 
     97      LOGICAL           ::   llais     ! read in a file or not 
    9598   END TYPE PTRACER 
    9699   ! 
     
    114117   LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   ln_trc_sbc    !: Use surface boundary condition data 
    115118   LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   ln_trc_cbc    !: Use coastal boundary condition data 
     119   LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   ln_trc_ais    !: Use Antarctic Ice Sheet boundary condition data 
    116120   LOGICAL , PUBLIC                                  ::   ln_rnf_ctl    !: remove runoff dilution on tracers 
    117121   REAL(wp), PUBLIC                                  ::   rn_sbc_time   !: Time scaling factor for SBC data (seconds in a day) 
     
    159163         &      ln_trc_ini(jptra)     ,                                                       & 
    160164         &      ln_trc_sbc(jptra)     , ln_trc_cbc(jptra)     , ln_trc_obc(jptra)     ,       & 
     165         &      ln_trc_ais(jptra)     ,                                                       & 
    161166         &      STAT = ierr(1)  ) 
    162167      ! 
  • NEMO/branches/2020/dev_r12740_TOP-01_rlod_Antarctic_ice_Sheet_Fe_Source/src/TOP/trcini.F90

    r12377 r13524  
    2525   USE trcice          ! tracers in sea ice 
    2626   USE trcbc           ! generalized Boundary Conditions 
     27   USE trcais          ! tracers from Antartic Ice Sheet  
    2728  
    2829   IMPLICIT NONE 
     
    158159         ln_trc_cbc(jn) =       sn_tracer(jn)%llcbc 
    159160         ln_trc_obc(jn) =       sn_tracer(jn)%llobc 
     161         ln_trc_ais(jn) =       sn_tracer(jn)%llais 
    160162      END DO 
    161163      ! 
     
    180182         WRITE(numout,*) 'trc_init_sms : Summary for selected passive tracers' 
    181183         WRITE(numout,*) '~~~~~~~~~~~~' 
    182          WRITE(numout,*) '    ID     NAME     INI  SBC  CBC  OBC' 
     184         WRITE(numout,*) '    ID     NAME     INI  SBC  CBC  OBC  AIS' 
    183185         DO jn = 1, jptra 
    184             WRITE(numout,9001) jn, TRIM(ctrcnm(jn)), ln_trc_ini(jn), ln_trc_sbc(jn),ln_trc_cbc(jn),ln_trc_obc(jn) 
     186            WRITE(numout,9001) jn, TRIM(ctrcnm(jn)), ln_trc_ini(jn),ln_trc_sbc(jn),ln_trc_cbc(jn),ln_trc_obc(jn),ln_trc_ais(jn) 
    185187         END DO 
    186188      ENDIF 
     
    189191         WRITE(numout,*) ' Applying tracer boundary conditions ' 
    190192      ENDIF 
     193      ! 
     194      IF( lwp .AND. ln_trcais ) THEN 
     195         WRITE(numout,*) 
     196         WRITE(numout,*) ' Applying tracer from Antarctic Ice Sheet ' 
     197      ENDIF 
    191198      
    192 9001  FORMAT(3x,i3,1x,a10,3x,l2,3x,l2,3x,l2,3x,l2) 
     1999001  FORMAT(3x,i3,1x,a10,3x,l2,3x,l2,3x,l2,3x,l2,3x,l2) 
    193200      ! 
    194201   END SUBROUTINE trc_ini_sms 
     
    240247      ENDIF 
    241248      ! 
     249      IF( ln_trcais ) CALL trc_ais_ini   ! set tracers from Antarctic Ice Sheet 
    242250      ! 
    243251      IF( ln_rsttr ) THEN              ! restart from a file 
  • NEMO/branches/2020/dev_r12740_TOP-01_rlod_Antarctic_ice_Sheet_Fe_Source/src/TOP/trcnam.F90

    r12489 r13524  
    136136      !! 
    137137      NAMELIST/namtrc/jp_bgc, ln_pisces, ln_my_trc, ln_age, ln_cfc11, ln_cfc12, ln_sf6, ln_c14, & 
    138          &            sn_tracer, ln_trcdta, ln_trcbc, ln_trcdmp, ln_trcdmp_clo, jp_dia3d, jp_dia2d 
     138         &            sn_tracer, ln_trcdta, ln_trcbc, ln_trcdmp, ln_trcdmp_clo, jp_dia3d, jp_dia2d, & 
     139         &            ln_trcais 
    139140      !!--------------------------------------------------------------------- 
    140141      ! Dummy settings to fill tracers data structure 
    141       !                  !   name   !   title   !   unit   !   init  !   sbc   !   cbc   !   obc  ! 
    142       sn_tracer = PTRACER( 'NONAME' , 'NOTITLE' , 'NOUNIT' , .false. , .false. , .false. , .false.) 
     142      !                  !   name   !   title   !   unit   !   init  !   sbc   !   cbc   !   obc  !  ais  ! 
     143      sn_tracer = PTRACER( 'NONAME' , 'NOTITLE' , 'NOUNIT' , .false. , .false. , .false. , .false., .false.) 
    143144      ! 
    144145      IF(lwp) WRITE(numout,*) 
     
    209210         WRITE(numout,*) '      Read inputs data from file (y/n)             ln_trcdta     = ', ln_trcdta 
    210211         WRITE(numout,*) '      Enable surface, lateral or open boundaries conditions (y/n)  ln_trcbc  = ', ln_trcbc 
     212         WRITE(numout,*) '      Enable Antarctic Ice Sheet nutrient supply   ln_trcais     = ', ln_trcais 
    211213         WRITE(numout,*) '      Damping of passive tracer (y/n)              ln_trcdmp     = ', ln_trcdmp 
    212214         WRITE(numout,*) '      Restoring of tracer on closed seas           ln_trcdmp_clo = ', ln_trcdmp_clo 
Note: See TracChangeset for help on using the changeset viewer.