Changeset 7829 for branches/UKMO/dev_r5518_nemo_fabm_ukmo/NEMOGCM/NEMO/TOP_SRC/trcini.F90

Timestamp:

2017-03-23T15:36:14+01:00 (7 years ago)

Author:

dford

Message:

Add a version of the NEMO-FABM coupling code. In theory, this should give equivalent results to PML gitlab commit 2e51db55.

File:

• branches/UKMO/dev_r5518_nemo_fabm_ukmo/NEMOGCM/NEMO/TOP_SRC/trcini.F90 (modified) (7 diffs)

branches/UKMO/dev_r5518_nemo_fabm_ukmo/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)
-
 
       IF( ln_rsttr ) THEN
@@ -140 +152 @@
             CALL wrk_dealloc( jpi, jpj, jpk, ztrcdta )
         ENDIF
-        !
-        trb(:,:,:,:) = trn(:,:,:,:)
-        ! 
-      ENDIF
+      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
@@ -158 +180 @@
       trai(:) = 0._wp                                                   ! initial content of all tracers
       DO jn = 1, jptra
-         trai(jn) = trai(jn) + glob_sum( trn(:,:,:,jn) * cvol(:,:,:)   )
+         !trai(jn) = trai(jn) + glob_sum( trn(:,:,:,jn) * cvol(:,:,:)   )
+         trai(jn) = trai(jn)
       END DO