Changeset 2798 for branches/2011/dev_r2784_CMCC1_topbfm/NEMOGCM/NEMO/TOP_SRC

Timestamp:

2011-07-11T16:18:18+02:00 (13 years ago)

Author:

vichi

Message:

Added initial coupling with BFM

Location:

branches/2011/dev_r2784_CMCC1_topbfm/NEMOGCM/NEMO/TOP_SRC

Files:

• par_trc.F90 (modified) (1 diff)
• trcini.F90 (modified) (3 diffs)
• trcnam.F90 (modified) (2 diffs)
• trcstp.F90 (modified) (4 diffs)

branches/2011/dev_r2784_CMCC1_topbfm/NEMOGCM/NEMO/TOP_SRC/par_trc.F90

@@ -17 +17 @@
    USE par_cfc       ! CFC 11 and 12 tracers
    USE par_my_trc    ! user defined passive tracers
-
+   USE par_bfm       ! BFM biogeochemical tracers
+ 
    IMPLICIT NONE
 
    ! Passive tracers : Total size
    ! ---------------               ! total number of passive tracers, of 2d and 3d output and trend arrays
-   INTEGER, PUBLIC,  PARAMETER ::   jptra    =  jp_lobster    + jp_pisces     + jp_cfc     + jp_c14b    + jp_my_trc
-   INTEGER, PUBLIC,  PARAMETER ::   jpdia2d  =  jp_lobster_2d + jp_pisces_2d  + jp_cfc_2d  + jp_c14b_2d + jp_my_trc_2d
-   INTEGER, PUBLIC,  PARAMETER ::   jpdia3d  =  jp_lobster_3d + jp_pisces_3d  + jp_cfc_3d  + jp_c14b_3d + jp_my_trc_3d
+   INTEGER, PUBLIC,  PARAMETER ::   jptra    =  jp_lobster    + jp_pisces     + jp_cfc     + jp_c14b    + jp_my_trc + jp_bfm
+   INTEGER, PUBLIC,  PARAMETER ::   jpdia2d  =  jp_lobster_2d + jp_pisces_2d  + jp_cfc_2d  + jp_c14b_2d + jp_my_trc_2d + jp_bfm_2d
+   INTEGER, PUBLIC,  PARAMETER ::   jpdia3d  =  jp_lobster_3d + jp_pisces_3d  + jp_cfc_3d  + jp_c14b_3d + jp_my_trc_3d + jp_bfm_3d
    !                     ! total number of sms diagnostic arrays
-   INTEGER, PUBLIC,  PARAMETER ::   jpdiabio = jp_lobster_trd + jp_pisces_trd + jp_cfc_trd + jp_c14b_trd + jp_my_trc_trd
+   INTEGER, PUBLIC,  PARAMETER ::   jpdiabio = jp_lobster_trd + jp_pisces_trd + jp_cfc_trd + jp_c14b_trd + jp_my_trc_trd +jp_bfm_trd
    
    !  1D configuration ("key_c1d")

branches/2011/dev_r2784_CMCC1_topbfm/NEMOGCM/NEMO/TOP_SRC/trcini.F90

@@ -25 +25 @@
    USE trcini_c14b     ! C14 bomb initialisation
    USE trcini_my_trc   ! MY_TRC   initialisation
+#if defined key_bfm
+   USE par_bfm
+#endif
    USE trcdta   
    USE daymod
@@ -94 +97 @@
       ENDIF
 
-      IF( ln_dm2dc .AND. ( lk_pisces .OR. lk_lobster ) )    &
-         &       CALL ctl_stop( ' The diurnal cycle is not compatible with PISCES or LOBSTER  ' )
+      IF( ln_dm2dc .AND. ( lk_pisces .OR. lk_lobster .OR. lk_bfm) )    &
+         &       CALL ctl_stop( ' The diurnal cycle is not compatible with PISCES or LOBSTER or BFM ' )
 
       IF( nn_cla == 1 )   &
@@ -118 +121 @@
       IF( lk_my_trc  ) THEN   ;   CALL trc_ini_my_trc       ! MY_TRC  tracers
       ELSE                    ;   IF(lwp) WRITE(numout,*) '          MY_TRC not used'
+      ENDIF
+
+      IF( lk_bfm     ) THEN   ;   CALL trc_ini_bfm          ! BFM  tracers
+      ELSE                    ;   IF(lwp) WRITE(numout,*) '          BFM not used'
       ENDIF
 

branches/2011/dev_r2784_CMCC1_topbfm/NEMOGCM/NEMO/TOP_SRC/trcnam.F90

@@ -26 +26 @@
    USE trcnam_c14b       ! C14 SMS namelist
    USE trcnam_my_trc     ! MY_TRC SMS namelist
+   USE trcnam_bfm        ! BFM namelist 
    USE trdmod_trc_oce
 
@@ -189 +190 @@
       ELSE                    ;   IF(lwp) WRITE(numout,*) '          MY_TRC not used'
       ENDIF
+      
+      IF( lk_bfm  ) THEN      ;   CALL trc_nam_bfm         ! BFM  tracers
+      ELSE                    ;   IF(lwp) WRITE(numout,*) '          BFM not used'
+      ENDIF
       !
    END SUBROUTINE trc_nam

branches/2011/dev_r2784_CMCC1_topbfm/NEMOGCM/NEMO/TOP_SRC/trcstp.F90

@@ -11 +11 @@
    !!----------------------------------------------------------------------
    USE oce_trc          ! ocean dynamics and active tracers variables
+#if ! defined key_bfm
    USE trc
    USE trctrp           ! passive tracers transport
@@ -20 +21 @@
    USE trdmod_trc_oce
    USE trdmld_trc
+#else
+   USE api_bfm, ONLY: bio_calc
+   USE par_bfm
+#ifdef key_obc
+   USE trcobcdta_bfm
+#endif
+#endif
    USE iom
    USE in_out_manager
@@ -49 +57 @@
       !!-------------------------------------------------------------------
 
+#if defined key_bfm
+   !---------------------------------------------
+   ! Check the main BFM flag
+   !---------------------------------------------
+      IF (bio_calc) THEN
+
+                             CALL trc_bfm( kt )           ! main call to BFM
+
+                             CALL trc_trp_bfm( kt )       ! transport of BFM tracers
+
+                             CALL trc_dia_bfm( kt )       ! diagnostic output for BFM
+      END IF
+#else
+
       IF( MOD( kt - 1 , nn_dttrc ) == 0 ) THEN      ! only every nn_dttrc time step
          !
@@ -71 +93 @@
          !
       ENDIF
+#endif
 
    END SUBROUTINE trc_stp