Changeset 10253 for branches/UKMO/dev_r5518_AMM15_package/NEMOGCM/NEMO/TOP_SRC/trcini.F90

Timestamp:

2018-10-29T15:55:40+01:00 (5 years ago)

Author:

kingr

Message:

Merged AMM15_v3_6_STABLE_package_collate@10237

File:

• branches/UKMO/dev_r5518_AMM15_package/NEMOGCM/NEMO/TOP_SRC/trcini.F90 (modified) (6 diffs)

branches/UKMO/dev_r5518_AMM15_package/NEMOGCM/NEMO/TOP_SRC/trcini.F90

@@ -24 +24 @@
    USE trcini_c14b     ! C14 bomb initialisation
    USE trcini_my_trc   ! MY_TRC   initialisation
+   ! +++>>> FABM
+   USE trcsms_fabm     ! FABM initialisation
+   USE trcini_fabm     ! FABM initialisation
+   ! FABM <<<FABM
    USE trcdta          ! initialisation from files
    USE daymod          ! calendar manager
@@ -32 +36 @@
    USE sbc_oce
    USE trcice          ! tracers in sea ice
+   USE trcbc,   only : trc_bc_init ! generalized Boundary Conditions
  
    IMPLICIT NONE
@@ -69 +74 @@
       IF(lwp) WRITE(numout,*) 'trc_init : initial set up of the passive tracers'
       IF(lwp) WRITE(numout,*) '~~~~~~~'
+      ! +++>>> FABM
+      ! Allow FABM to update numbers of biogeochemical tracers, diagnostics (jptra etc.)
+      IF( lk_fabm ) CALL nemo_fabm_init
+      ! FABM <<<+++
 
       CALL top_alloc()              ! allocate TOP arrays
@@ -101 +110 @@
       IF( lk_c14b    )       CALL trc_ini_c14b         ! C14 bomb  tracer
       IF( lk_my_trc  )       CALL trc_ini_my_trc       ! MY_TRC  tracers
+      ! +++>>> FABM
+      IF( lk_fabm    )       CALL trc_ini_fabm         ! FABM    tracers
+      ! FABM <<<+++
 
       CALL trc_ice_ini                                 ! Tracers in sea ice
@@ -110 +122 @@
       ENDIF
 
+      ! Initialisation of tracers Initial Conditions
       IF( ln_trcdta )      CALL trc_dta_init(jptra)
 
@@ -144 +157 @@
         ! 
       ENDIF
+      ! --->>> FABM
+! Initialisation of tracers Boundary Conditions  - here so that you can use initial condition as boundary
+      !IF( lk_my_trc )     CALL trc_bc_init(jptra)
+      ! FABM <<<---
+      ! FABM +++>>>
+! Initialisation of FABM diagnostics and tracer boundary conditions (so that you can use initial condition as boundary)
+      IF( lk_fabm )     THEN
+          wndm=0._wp !uninitiased field at this point
+          qsr=0._wp !uninitiased field at this point
+          CALL compute_fabm ! only needed to set-up diagnostics
+          CALL trc_bc_init(jptra)
+      ENDIF
+      ! FABM <<<+++
  
       tra(:,:,:,:) = 0._wp